Changeset 12607 for trunk/psLib/test/math/tap_psSparse.c
- Timestamp:
- Mar 27, 2007, 12:52:03 PM (19 years ago)
- File:
-
- 1 edited
-
trunk/psLib/test/math/tap_psSparse.c (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/test/math/tap_psSparse.c
r12513 r12607 74 74 dS2 = sqrt (dS2/100.0 - dS*dS); 75 75 76 ok_float_tol (dS2, 0.0, 1e-4, "scatter: %.20f", dS2);76 is_float_tol (dS2, 0.0, 1e-4, "scatter: %.20f", dS2); 77 77 78 78 psFree (matrix); … … 147 147 dS2 = sqrt (dS2/100.0 - dS*dS); 148 148 149 ok_float_tol (dS2, 0.0, 1e-4, "scatter: %.20f", dS2);149 is_float_tol (dS2, 0.0, 1e-4, "scatter: %.20f", dS2); 150 150 151 151 psFree (matrix); … … 219 219 // test the support functions: LowerProduct 220 220 fVec = psSparseBorderLowerProduct (NULL, border, xRef); 221 ok_float (fVec->n, 1.0, "f dimen: %d", fVec->n);222 ok_float (fVec->data.F32[0], 0.3, "f(0): %f", fVec->data.F32[0]);221 is_float (fVec->n, 1.0, "f dimen: %d", fVec->n); 222 is_float (fVec->data.F32[0], 0.3, "f(0): %f", fVec->data.F32[0]); 223 223 psFree (fVec); 224 224 225 225 // test the support functions: Upper Product 226 226 fVec = psSparseBorderUpperProduct (NULL, border, yRef); 227 ok_float (fVec->n, 2.0, "f dimen: %d", fVec->n);228 ok_float (fVec->data.F32[0], 0.05, "f(0): %f", fVec->data.F32[0]);229 ok_float (fVec->data.F32[1], 0.10, "f(0): %f", fVec->data.F32[1]);227 is_float (fVec->n, 2.0, "f dimen: %d", fVec->n); 228 is_float (fVec->data.F32[0], 0.05, "f(0): %f", fVec->data.F32[0]); 229 is_float (fVec->data.F32[1], 0.10, "f(0): %f", fVec->data.F32[1]); 230 230 psFree (fVec); 231 231 232 232 // test the support functions: Square Product 233 233 fVec = psSparseBorderSquareProduct (NULL, border, yRef); 234 ok_float (fVec->n, 1.0, "f dimen: %d", fVec->n);235 ok_float (fVec->data.F32[0], 0.25, "f(0): %f", fVec->data.F32[0]);234 is_float (fVec->n, 1.0, "f dimen: %d", fVec->n); 235 is_float (fVec->data.F32[0], 0.25, "f(0): %f", fVec->data.F32[0]); 236 236 psFree (fVec); 237 237 … … 239 239 psVector *gVec = NULL; 240 240 psSparseBorderMultiply (&fVec, &gVec, border, xRef, yRef); 241 ok_float (fVec->data.F32[0], 1.15, "f(0): %f", fVec->data.F32[0]);242 ok_float (fVec->data.F32[1], 1.20, "f(1): %f", fVec->data.F32[1]);243 ok_float (gVec->data.F32[0], 0.55, "g(0): %f", gVec->data.F32[0]);241 is_float (fVec->data.F32[0], 1.15, "f(0): %f", fVec->data.F32[0]); 242 is_float (fVec->data.F32[1], 1.20, "f(1): %f", fVec->data.F32[1]); 243 is_float (gVec->data.F32[0], 0.55, "g(0): %f", gVec->data.F32[0]); 244 244 245 245 // supply the fVec and gVec data to the border … … 263 263 psVector *yFit = NULL; 264 264 psSparseBorderSolve (&xFit, &yFit, constraint, border, 4); 265 ok_float_tol (xFit->data.F32[0], 1.0, 1e-4, "f(0): %f", xFit->data.F32[0]);266 ok_float_tol (xFit->data.F32[1], 1.0, 1e-4, "f(1): %f", xFit->data.F32[1]);267 ok_float_tol (yFit->data.F32[0], 0.5, 1e-4, "g(0): %f", yFit->data.F32[0]);265 is_float_tol (xFit->data.F32[0], 1.0, 1e-4, "f(0): %f", xFit->data.F32[0]); 266 is_float_tol (xFit->data.F32[1], 1.0, 1e-4, "f(1): %f", xFit->data.F32[1]); 267 is_float_tol (yFit->data.F32[0], 0.5, 1e-4, "g(0): %f", yFit->data.F32[0]); 268 268 269 269 psFree (xFit); … … 321 321 } 322 322 } 323 ok_float(border->Bij->data.F32[0][0], 0.1, "set matrix element %d,%d", 0, 0);323 is_float(border->Bij->data.F32[0][0], 0.1, "set matrix element %d,%d", 0, 0); 324 324 325 325 // construct the Q component: … … 381 381 dS /= Nrows; 382 382 dS2 = sqrt (dS2/Nrows - dS*dS); 383 ok_float_tol (dS2, 0.0, 1e-4, "scatter: %.20f", dS2);383 is_float_tol (dS2, 0.0, 1e-4, "scatter: %.20f", dS2); 384 384 385 385 // measure stdev between yFit and yRef 386 386 float dY = yRef->data.F32[0] - yFit->data.F32[0]; 387 387 // fprintf (stderr, "%5.3f %5.3f %5.3f %5.3f\n", gVec->data.F32[0], yRef->data.F32[0], yFit->data.F32[0], dS); 388 ok_float_tol (dY, 0.0, 2e-4, "scatter: %.20f", dY);388 is_float_tol (dY, 0.0, 2e-4, "scatter: %.20f", dY); 389 389 390 390 psFree (xRef);
Note:
See TracChangeset
for help on using the changeset viewer.
