Index: trunk/psLib/test/math/tap_psSparse.c
===================================================================
--- trunk/psLib/test/math/tap_psSparse.c	(revision 12607)
+++ trunk/psLib/test/math/tap_psSparse.c	(revision 13123)
@@ -10,11 +10,9 @@
     plan_tests(26);
 
-//    diag("psSparse() tests");
-
     // test psSparseSolve for a simple normal example matrix
     {
         psMemId id = psMemGetId();
 
-//        diag ("solve a normalized matrix equation with psSparseSolve");
+        //solve a normalized matrix equation with psSparseSolve
 
         // the basic equation is Ax = b
@@ -28,6 +26,5 @@
         skip_start(matrix == NULL, 5, "Skipping tests because psSparseAlloc() failed");
 
-        for (int i = 0; i < 100; i++)
-        {
+        for (int i = 0; i < 100; i++) {
             psSparseMatrixElement (matrix, i, i, 1.0);
             if (i + 1 < 100) {
@@ -38,5 +35,5 @@
         // incoming matrix elements do not need to be in order; sort before
         // applying sparse matrix
-        psSparseResort (matrix);
+        psSparseResort(matrix);
 
         psVector *xRef = psVectorAlloc (100, PS_TYPE_F32);
@@ -46,9 +43,8 @@
         }
 
-        psVector *bVec = psSparseMatrixTimesVector (NULL, matrix, xRef);
-
-        for (int i = 0; i < 100; i++)
-        {
-            psSparseVectorElement (matrix, i, bVec->data.F32[i]);
+        psVector *bVec = psSparseMatrixTimesVector(NULL, matrix, xRef);
+
+        for (int i = 0; i < 100; i++) {
+            psSparseVectorElement(matrix, i, bVec->data.F32[i]);
         }
 
@@ -59,5 +55,5 @@
 
         // solve for normalization terms (need include local sky?)
-        psVector *xFit = psSparseSolve (NULL, constraint, matrix, 4);
+        psVector *xFit = psSparseSolve(NULL, constraint, matrix, 4);
 
         // measure stdev between xFit and xRes
@@ -89,8 +85,6 @@
         psMemId id = psMemGetId();
 
-//        diag ("solve a non-normalized matrix equation with psSparseSolve");
-
+        //solve a non-normalized matrix equation with psSparseSolve
         // the basic equation is Ax = b
-
         // create a matrix A with diagonals of 1 and a small number of off diagonal elements.
         // construct a vector x, construct the corresponding vector b by multiplication. solve
@@ -101,6 +95,5 @@
         skip_start(matrix == NULL, 5, "Skipping tests because psSparseAlloc() failed");
 
-        for (int i = 0; i < 100; i++)
-        {
+        for (int i = 0; i < 100; i++) {
             psSparseMatrixElement (matrix, i, i, 5.0);
             if (i + 1 < 100) {
@@ -111,5 +104,5 @@
         // incoming matrix elements do not need to be in order; sort before
         // applying sparse matrix
-        psSparseResort (matrix);
+        psSparseResort(matrix);
 
         psVector *xRef = psVectorAlloc (100, PS_TYPE_F32);
@@ -119,9 +112,8 @@
         }
 
-        psVector *bVec = psSparseMatrixTimesVector (NULL, matrix, xRef);
-
-        for (int i = 0; i < 100; i++)
-        {
-            psSparseVectorElement (matrix, i, bVec->data.F32[i]);
+        psVector *bVec = psSparseMatrixTimesVector(NULL, matrix, xRef);
+
+        for (int i = 0; i < 100; i++) {
+            psSparseVectorElement(matrix, i, bVec->data.F32[i]);
         }
 
@@ -132,5 +124,5 @@
 
         // solve for normalization terms (need include local sky?)
-        psVector *xFit = psSparseSolve (NULL, constraint, matrix, 4);
+        psVector *xFit = psSparseSolve(NULL, constraint, matrix, 4);
 
         // measure stdev between xFit and xRes
@@ -173,5 +165,5 @@
         psMemId id = psMemGetId();
 
-//        diag ("solve a simple, small matrix equation ");
+        //solve a simple, small matrix equation
 
         // the basic equation (Ax = b) is:
@@ -192,8 +184,8 @@
         skip_start(matrix == NULL, 5, "Skipping tests because psSparseAlloc() failed");
 
-        psSparseMatrixElement (matrix, 0, 0, 1.0);
-        psSparseMatrixElement (matrix, 1, 1, 1.0);
-        psSparseMatrixElement (matrix, 1, 0, 0.1);
-        psSparseResort (matrix);
+        psSparseMatrixElement(matrix, 0, 0, 1.0);
+        psSparseMatrixElement(matrix, 1, 1, 1.0);
+        psSparseMatrixElement(matrix, 1, 0, 0.1);
+        psSparseResort(matrix);
 
         // border region has a width of 1:
@@ -244,11 +236,9 @@
 
         // supply the fVec and gVec data to the border
-        for (int i = 0; i < Nrows; i++)
-        {
-            psSparseVectorElement (border->sparse, i, fVec->data.F32[i]);
-        }
-        for (int i = 0; i < Nborder; i++)
-        {
-            psSparseBorderElementG (border, i, gVec->data.F32[i]);
+        for (int i = 0; i < Nrows; i++) {
+            psSparseVectorElement(border->sparse, i, fVec->data.F32[i]);
+        }
+        for (int i = 0; i < Nborder; i++) {
+            psSparseBorderElementG(border, i, gVec->data.F32[i]);
         }
 
@@ -262,5 +252,5 @@
         psVector *xFit = NULL;
         psVector *yFit = NULL;
-        psSparseBorderSolve (&xFit, &yFit, constraint, border, 4);
+        psSparseBorderSolve(&xFit, &yFit, constraint, border, 4);
         is_float_tol (xFit->data.F32[0], 1.0, 1e-4, "f(0): %f", xFit->data.F32[0]);
         is_float_tol (xFit->data.F32[1], 1.0, 1e-4, "f(1): %f", xFit->data.F32[1]);
@@ -283,5 +273,5 @@
         psMemId id = psMemGetId();
 
-//        diag ("solve a simple, small matrix equation ");
+        // solve a simple, small matrix equation
 
         // the basic equation (Ax = b) is:
@@ -302,12 +292,11 @@
         skip_start(matrix == NULL, 5, "Skipping tests because psSparseAlloc() failed");
 
-        for (int i = 0; i < Nrows; i++)
-        {
-            psSparseMatrixElement (matrix, i, i, 1.0);
+        for (int i = 0; i < Nrows; i++) {
+            psSparseMatrixElement(matrix, i, i, 1.0);
             if (i + 1 < Nrows) {
-                psSparseMatrixElement (matrix, i + 1, i, 0.1);
-            }
-        }
-        psSparseResort (matrix);
+                psSparseMatrixElement(matrix, i + 1, i, 0.1);
+            }
+        }
+        psSparseResort(matrix);
 
         // border region has a width of 1:
@@ -349,11 +338,9 @@
 
         // supply the fVec and gVec data to the border
-        for (int i = 0; i < Nrows; i++)
-        {
-            psSparseVectorElement (border->sparse, i, fVec->data.F32[i]);
-        }
-        for (int i = 0; i < Nborder; i++)
-        {
-            psSparseBorderElementG (border, i, gVec->data.F32[i]);
+        for (int i = 0; i < Nrows; i++) {
+            psSparseVectorElement(border->sparse, i, fVec->data.F32[i]);
+        }
+        for (int i = 0; i < Nborder; i++) {
+            psSparseBorderElementG(border, i, gVec->data.F32[i]);
         }
 
@@ -367,5 +354,5 @@
         psVector *xFit = NULL;
         psVector *yFit = NULL;
-        psSparseBorderSolve (&xFit, &yFit, constraint, border, 4);
+        psSparseBorderSolve(&xFit, &yFit, constraint, border, 4);
 
         // measure stdev between xFit and xRef