Changeset 6484 for trunk/psLib/src/math/psMinimizeLMM.c
- Timestamp:
- Feb 24, 2006, 1:43:16 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/math/psMinimizeLMM.c (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/math/psMinimizeLMM.c
r6346 r6484 10 10 * @author EAM, IfA 11 11 * 12 * @version $Revision: 1. 5$ $Name: not supported by cvs2svn $13 * @date $Date: 2006-02- 07 23:14:21$12 * @version $Revision: 1.6 $ $Name: not supported by cvs2svn $ 13 * @date $Date: 2006-02-24 23:43:15 $ 14 14 * 15 15 * Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii … … 239 239 psF64 ymodel; 240 240 psVector *deriv = psVectorAlloc(params->n, PS_TYPE_F32); 241 deriv->n = deriv->nalloc; 241 242 242 243 // zero alpha and beta for summing below … … 370 371 psVector *Beta = psVectorAlloc(params->n, PS_TYPE_F64); 371 372 psVector *Params = psVectorAlloc(params->n, PS_TYPE_F32); 373 beta->n = beta->nalloc; 374 Beta->n = Beta->nalloc; 375 Params->n = Params->nalloc; 372 376 psVector *dy = NULL; 373 377 psF64 Chisq = 0.0; … … 385 389 param_min = psVectorAlloc(params->n, PS_TYPE_F32); 386 390 param_max = psVectorAlloc(params->n, PS_TYPE_F32); 391 beta_lim->n = beta_lim->nalloc; 392 param_min->n = param_min->nalloc; 393 param_max->n = param_max->nalloc; 387 394 for (int i = 0; i < params->n; i++) { 388 395 beta_lim->data.F32[i] = covar->data.F64[0][i]; … … 402 409 } else { 403 410 dy = psVectorAlloc(y->n, PS_TYPE_F32); 411 dy->n = dy->nalloc; 404 412 psVectorInit(dy, 1.0); 405 413 }
Note:
See TracChangeset
for help on using the changeset viewer.
