133 |
double *x; |
double *x; |
134 |
double *xupper; /* upper bounds on each variable */ |
double *xupper; /* upper bounds on each variable */ |
135 |
double *xlower; /* lower bounds on each variable */ |
double *xlower; /* lower bounds on each variable */ |
136 |
double *f; |
double *g; /* gradients of f with respect to each variable */ |
137 |
double *g; |
TronMatrix A; /* hessian matrix, passed to TRON? */ |
138 |
TronMatrix A, B, L; |
TronMatrix B; /* some matrix that we get back from TRON */ |
139 |
|
TronMatrix L; /* some other matrix that we get back from TRON */ |
140 |
|
|
141 |
double *xc, *s, *dsave, *wa; |
double *xc, *s, *dsave, *wa; |
142 |
int *indfree, *isave, *iwa; |
int *indfree, *isave, *iwa; |
665 |
double delta = SLV_PARAM_REAL(&(sys->params),TRON_PARAM_GTOL); |
double delta = SLV_PARAM_REAL(&(sys->params),TRON_PARAM_GTOL); |
666 |
/** @TODO fmin should be taken from the model declaration somehow, not a solar parameter. */ |
/** @TODO fmin should be taken from the model declaration somehow, not a solar parameter. */ |
667 |
|
|
668 |
DTRON(&(sys->n),sys->x,sys->xlower,sys->xupper,sys->f,sys->g |
DTRON(&(sys->n),sys->x,sys->xlower,sys->xupper,&(sys->objective),sys->g |
669 |
,TRON_MATRIX_ARG(sys,A) |
,TRON_MATRIX_ARG(sys,A) |
670 |
,&frtol,&fatol,&fmin,&cgtol,&itermax,&delta,sys->task |
,&frtol,&fatol,&fmin,&cgtol,&itermax,&delta,sys->task |
671 |
,TRON_MATRIX_ARG(sys,B) |
,TRON_MATRIX_ARG(sys,B) |