Changeset 3987

Timestamp:

May 19, 2005, 12:50:05 PM (21 years ago)

Author:

gusciora

Message:

Added comment.

Location:

trunk/psLib/test/dataManip

Files:

• tst_psFunc08.c (modified) (1 diff)
• tst_psMinimize06.c (modified) (1 diff)

trunk/psLib/test/dataManip/tst_psFunc08.c

@@ -4 +4 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.2 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2005-04-07 20:27:42 $
+*  @version  $Revision: 1.3 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2005-05-19 22:50:05 $
+*
+*  XXX: Probably should test single- and multi-dimensional polynomials in
+*  which one diminsion is constant (n == 1).
 *
 * Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii

trunk/psLib/test/dataManip/tst_psMinimize06.c

@@ -51 +51 @@
     return(sum);
 }
+
+/*****************************************************************************
+myFuncEAM(): This routine subtracts the associate value in expectedParm[] from
+each parameter and then squares it, then sums that for all parameters, then
+adds MIN_VALUE to it.  The minimum for this function will be MIN_VALUE, and
+will occur when each parameter equals the associated value in expectedParm[].
+ 
+This procedure ignores the coordinates, other than to ensure that they were
+passed correctly from psMinimizePowell().
+ 
+ *****************************************************************************/
+psVector *myFuncEAM(psImage *myDeriv,
+                    psVector *myParams,
+                    psArray *myCoords)
+{
+    psVector *sum = psVectorAlloc(myCoords->n, PS_TYPE_F32);
+    psS32 i;
+    psS32 j;
+
+
+    if (myDeriv == NULL) {
+        myDeriv = psImageAlloc(myParams->n, myCoords->n, PS_TYPE_F32);
+        psError(PS_ERR_UNKNOWN, true, "myDeriv is NULL.\n");
+    }
+
+    for (i=0;i<N;i++) {
+        sum->data.F32[i] = MIN_VALUE;
+        for (j=0;j<NUM_PARAMS;j++) {
+            sum->data.F32[i]+= (myParams->data.F32[j] - expectedParm[j]) *
+                               (myParams->data.F32[j] - expectedParm[j]);
+
+            myDeriv->data.F32[i][j] = (2.0 * myParams->data.F32[j]) -
+                                      (2.0 * expectedParm[j]);
+        }
+    }
+
+    return(sum);
+}
+
 
 psS32 t01()