Changeset 2228 for trunk/psLib/src/dataManip/psMinimize.h
- Timestamp:
- Oct 28, 2004, 12:47:00 PM (22 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/dataManip/psMinimize.h (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/dataManip/psMinimize.h
r2221 r2228 8 8 * @author GLG, MHPCC 9 9 * 10 * @version $Revision: 1.3 4$ $Name: not supported by cvs2svn $11 * @date $Date: 2004-10-2 7 23:31:43$10 * @version $Revision: 1.35 $ $Name: not supported by cvs2svn $ 11 * @date $Date: 2004-10-28 22:46:57 $ 12 12 * 13 13 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 81 81 const psArray *x); 82 82 83 /*84 psVector *psMinimizeLMChi2Gauss1D(psImage *deriv,85 const psVector *params,86 const psArray *coords)87 {88 psVector *x = (psVector *) coords->data[0];89 psVector *y = (psVector *) coords->data[1];90 psVector *yErr = NULL;91 int i;92 int j;93 float normalization = params->data.F32[0];94 float center = params->data.F32[1];95 float stdev = params->data.F32[2];96 psVector *out = psVectorAlloc(x->n, PS_TYPE_F32);97 98 99 if (deriv == NULL) {100 deriv = psImageAlloc(params->n, x->n, PS_TYPE_F32);101 }102 103 if (coords->n == 2) {104 yErr = (psVector *) coords->data[2];105 }106 107 for (i=0;i<x->n;i++) {108 out->data.F32[i] = psGaussian(x->data.F32[i], center, stdev, 1);109 }110 111 for (i=0;i<x->n;i++) {112 for (j=0;j<params->n;j++) {113 deriv->data.F32[i][j] = the derivative of the gaussian w.r.t. the j-th parameter at the coord x[i]114 }115 }116 return(out);117 }118 */119 120 83 psBool psMinimizeLMChi2(psMinimization *min, 121 84 psImage *covar, … … 137 100 psMinimizePowellFunc func); 138 101 102 psVector *psMinimizeLMChi2Gauss1D(psImage *deriv, 103 const psVector *params, 104 const psArray *coords); 105 139 106 typedef 140 107 psVector *(*psMinimizeChi2PowellFunc)(const psVector *params,
Note:
See TracChangeset
for help on using the changeset viewer.
