Changeset 2206 for trunk/psLib/src/math/psMinimize.h

Timestamp:

Oct 26, 2004, 3:15:47 PM (22 years ago)

Author:

gusciora

Message:

...

File:

• trunk/psLib/src/math/psMinimize.h (modified) (2 diffs)

trunk/psLib/src/math/psMinimize.h

@@ -8 +8 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.32 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-10-27 00:57:31 $
+ *  @version $Revision: 1.33 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-10-27 01:15:47 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -88 +88 @@
                                   const psArray *x);
 
+/*
+psVector *psMinimizeLMChi2Gauss1D(psImage *deriv,
+                                  const psVector *params,
+                                  const psArray *coords)
+{
+    psVector *x = (psVector *) coords->data[0];
+    psVector *y = (psVector *) coords->data[1];
+    psVector *yErr = NULL;
+    int i;
+    int j;
+    float normalization = params->data.F32[0];
+    float center = params->data.F32[1];
+    float stdev = params->data.F32[2];
+    psVector *out = psVectorAlloc(x->n, PS_TYPE_F32);
+ 
+ 
+    if (deriv == NULL) {
+        deriv = psImageAlloc(params->n, x->n, PS_TYPE_F32);
+    }
+ 
+    if (coords->n == 2) {
+        yErr = (psVector *) coords->data[2];
+    }
+ 
+    for (i=0;i<x->n;i++) {
+        out->data.F32[i] = psGaussian(x->data.F32[i], center, stdev, 1);
+    }
+ 
+    for (i=0;i<x->n;i++) {
+        for (j=0;j<params->n;j++) {
+            deriv->data.F32[i][j] = the derivative of the gaussian w.r.t. the j-th parameter at the coord x[i]
+        }
+    }
+    return(out);
+}
+*/
+
 psBool psMinimizeLMChi2(psMinimization *min,
                         psImage *covar,