Changeset 13124 for trunk/psLib/test/math/tap_psMinimizeLMM.c

Timestamp:

May 1, 2007, 6:20:06 PM (19 years ago)

Author:

gusciora

Message:

Added tabular file of all psLib functions in Chapter 6 of the SDRS (Data
Manipulation). Most of the test file changes in this check-in are fairly
cosmetic.

File:

• trunk/psLib/test/math/tap_psMinimizeLMM.c (modified) (8 diffs)

trunk/psLib/test/math/tap_psMinimizeLMM.c

@@ -2 +2 @@
     This routine must ensure that psMinimizeLM() works correctly.
  
-    XXX: This code needs a lot of additional test case work.
-         The minimization currently (always?) fails and we don't
-         attempt to check the output values.
+    XXX: This code needs a lot of additional test case work.  The minimization
+         currently fails and we don't attempt to check the output values.
+    XXX: Add tests for
+        covar arg set to non-NULL
+        constraint set to non-NULL
+        Set x->vectors to NULL, or use wrong types
+        yWt (errors) vector set to incorrect size, type.
  *****************************************************************************/
 #include <stdio.h>
@@ -19 +23 @@
 float expectedParm[NUM_PARAMS];
 
-
+// y = p2 + p0 * e^( x0^2 + x1^2 / (2 * p1^2) )
 float function(const psVector *params, const psVector *x)
 {
@@ -35 +39 @@
 
 
-/*****************************************************************************
-fitFunc():
-    sum = param[0] * x[0] * x[1] +
-          param[1] * x[0] +
-          param[2] * x[0]^2 +
-          param[3] * x[1] +
-          param[4] * x[1]^2
- 
- *****************************************************************************/
 psF32 fitFunc(psVector *deriv,
               psVector *params,
@@ -56 +51 @@
 }
 
+
 #define NUM_ITER 10
 #define TOL  20.0
@@ -62 +58 @@
     psLogSetFormat("HLNM");
     psLogSetLevel(PS_LOG_INFO);
-    plan_tests(2);
+    plan_tests(34);
+
+
     // Test psMinimizationAlloc()
     {
@@ -79 +77 @@
     }
 
-    /*
-        // Test psMinConstrainAlloc()
-        {
-            psMemId id = psMemGetId();
-            psMinConstrain *tmp = psMinConstrainAlloc(NUM_ITER, TOL);
-            ok(tmp != NULL. "psMinConstrainAlloc() returned non-NULL");
-            skip_start(tmp == NULL, 4, "Skipping tests because psMinConstrainAlloc() failed");
-            ok(tmp->paramMask == NULL, "psMinConstrainAlloc() properly set ->paramMask");
-            ok(tmp->paramMax == NULL, "psMinConstrainAlloc() properly set ->paramMax");
-            ok(tmp->paramMin == NULL, "psMinConstrainAlloc() properly set ->paramMin");
-            ok(tmp->paramDelta == NULL, "psMinConstrainAlloc() properly set ->paramDelta");
-            psFree(tmp);
-            skip_end();
-            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
-        }
-    */
-
-    // Test psMinimizeLMChi2()
+
+    // Test psMinConstraintAlloc()
+    {
+        psMemId id = psMemGetId();
+        psMinConstraint *tmp = psMinConstraintAlloc();
+        ok(tmp != NULL, "psMinConstraintAlloc() returned non-NULL");
+        skip_start(tmp == NULL, 2, "Skipping tests because psMinConstraintAlloc() failed");
+        ok(tmp->paramMask == NULL, "psMinConstraintAlloc() properly set ->paramMask");
+        ok(tmp->checkLimits == NULL, "psMinConstraintAlloc() properly set ->checkLimits");
+        psFree(tmp);
+        skip_end();
+        ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+    }
+
+
+    // Test psMinimizeLMChi2(): unallowed input parameters.
+    {
+        psMemId id = psMemGetId();
+        psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS, 1); // Random number generator; using known seed
+        psMinimization *min = psMinimizationAlloc(NUM_ITERATIONS, ERR_TOL);
+        psArray *ordinates = psArrayAlloc(NUM_DATA_POINTS);
+        psVector *coordinates = psVectorAlloc(NUM_DATA_POINTS, PS_TYPE_F32);
+        psVector *coordinatesF64 = psVectorAlloc(NUM_DATA_POINTS, PS_TYPE_F64);
+        psVector *errors = psVectorAlloc(NUM_DATA_POINTS, PS_TYPE_F32);
+        psVector *params = psVectorAlloc(NUM_PARAMS, PS_TYPE_F32);
+        psVector *paramsF64 = psVectorAlloc(NUM_PARAMS, PS_TYPE_F64);
+        psVector *trueParams = psVectorAlloc(NUM_PARAMS, PS_TYPE_F32);
+
+        // Test with psMinimization set to NULL
+        {
+            psMemId id = psMemGetId();
+            bool tmpBool = psMinimizeLMChi2(NULL, NULL, params, NULL, ordinates,
+                           coordinates, errors, (psMinimizeLMChi2Func)fitFunc);
+            ok(!tmpBool, "psMinimizeLMChi2() returned FALSE with NULL psMinimization");
+            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+        }
+
+
+        // Test with params set to NULL
+        {
+            psMemId id = psMemGetId();
+            bool tmpBool = psMinimizeLMChi2(min, NULL, NULL, NULL, ordinates,
+                           coordinates, errors, (psMinimizeLMChi2Func)fitFunc);
+            ok(!tmpBool, "psMinimizeLMChi2() returned FALSE with NULL params");
+            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+        }
+
+
+        // Test with params wrong type
+        {
+            psMemId id = psMemGetId();
+            bool tmpBool = psMinimizeLMChi2(min, NULL, paramsF64, NULL, ordinates,
+                           coordinates, errors, (psMinimizeLMChi2Func)fitFunc);
+            ok(!tmpBool, "psMinimizeLMChi2() returned FALSE with NULL params");
+            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+        }
+
+
+        // Test with ordinates (x) set to NULL
+        {
+            psMemId id = psMemGetId();
+            bool tmpBool = psMinimizeLMChi2(min, NULL, params, NULL, NULL,
+                           coordinates, errors, (psMinimizeLMChi2Func)fitFunc);
+            ok(!tmpBool, "psMinimizeLMChi2() returned FALSE with NULL ordinates (x)");
+            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+        }
+
+
+        // Test with coordinates (y) set to NULL
+        {
+            psMemId id = psMemGetId();
+            bool tmpBool = psMinimizeLMChi2(min, NULL, params, NULL, ordinates,
+                           NULL, errors, (psMinimizeLMChi2Func)fitFunc);
+            ok(!tmpBool, "psMinimizeLMChi2() returned FALSE with NULL coordinates (y)");
+            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+        }
+
+
+        // Test with coordinates (y) wrong type (F64)
+        {
+            psMemId id = psMemGetId();
+            bool tmpBool = psMinimizeLMChi2(min, NULL, params, NULL, ordinates,
+                           coordinatesF64, errors, (psMinimizeLMChi2Func)fitFunc);
+            ok(!tmpBool, "psMinimizeLMChi2() returned FALSE with coordinates (y) wrong type");
+            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+        }
+
+
+        // Test with ordinates and coordinates wrong size
+        {
+            psMemId id = psMemGetId();
+            coordinates->n--;
+            bool tmpBool = psMinimizeLMChi2(min, NULL, params, NULL, ordinates,
+                           coordinates, errors, (psMinimizeLMChi2Func)fitFunc);
+            ok(!tmpBool, "psMinimizeLMChi2() returned FALSE with NULL psMinimization");
+            coordinates->n++;
+            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+        }
+
+
+        // Test with function set to NULL
+        {
+            psMemId id = psMemGetId();
+            bool tmpBool = psMinimizeLMChi2(min, NULL, params, NULL, ordinates,
+                           coordinates, errors, NULL);
+            ok(!tmpBool, "psMinimizeLMChi2() returned FALSE with NULL fit function");
+            ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+        }
+
+        psFree(min);
+        psFree(params);
+        psFree(paramsF64);
+        psFree(trueParams);
+        psFree(ordinates);
+        psFree(coordinates);
+        psFree(coordinatesF64);
+        psFree(errors);
+        psFree(rng);
+        ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+    }
+
+
+    // Test psMinimizeLMChi2() with legitimate input values
+    // Currently this fails, and we do not attempt to verify output
     {
         psMemId id = psMemGetId();
@@ -158 +262 @@
         }
         printf("Mean relative difference is %f\n", diff/(float)NUM_DATA_POINTS);
-
+        skip_end();
         psFree(min);
         psFree(params);
@@ -166 +270 @@
         psFree(errors);
         psFree(rng);
-        skip_end();
-        ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
-    }
-}
+        ok(!psMemCheckLeaks (id, NULL, NULL, false), "no memory leaks");
+    }
+}