29 #define dls001_1 (mdls001_._1)
30 #define dls001_2 (mdls001_._2)
31 #define dls001_3 (mdls001_._3)
33 #define dlsa01_1 (mdlsa01_._1)
34 #define dlsa01_2 (mdlsa01_._2)
35 #define dlsa01_3 (mdlsa01_._3)
48 C_INT *nyh,
double *ewt,
double *ftem,
double *savf,
52 C_INT yh_dim1, yh_offset, i__1, i__2, i__3, i__4;
108 yh_offset = 1 + yh_dim1;
136 for (i__ = 1; i__ <= i__1; ++i__)
146 for (i__ = 1; i__ <= i__1; ++i__)
167 for (j = 1; j <= i__1; ++j)
171 d__1 = srur * fabs(yj), d__2 = r0 / ewt[j];
175 f(&neq[1], &
dls001_1.tn, &y[1], &ftem[1]);
178 for (i__ = 1; i__ <= i__2; ++i__)
181 wm[i__ + j1] = (ftem[i__] - savf[i__]) * fac;
198 for (i__ = 1; i__ <= i__1; ++i__)
227 for (i__ = 1; i__ <= i__1; ++i__)
233 jac(&neq[1], &
dls001_1.tn, &y[1], &ml, &mu, &wm[ml3], &meband);
237 for (i__ = 1; i__ <= i__1; ++i__)
263 for (j = 1; j <= i__1; ++j)
268 for (i__ = j; i__3 < 0 ? i__ >= i__2 : i__ <= i__2; i__ += i__3)
272 d__1 = srur * fabs(yi), d__2 = r0 / ewt[i__];
278 f(&neq[1], &
dls001_1.tn, &y[1], &ftem[1]);
282 for (jj = j; i__2 < 0 ? jj >= i__3 : jj <= i__3; jj += i__2)
284 y[jj] = yh[jj + yh_dim1];
287 d__1 = srur * fabs(yjj), d__2 = r0 / ewt[jj];
296 ii = jj * meb1 - ml + 2;
299 for (i__ = i1; i__ <= i__4; ++i__)
302 wm[ii + i__] = (ftem[i__] - savf[i__]) * fac;
320 for (i__ = 1; i__ <= i__1; ++i__)
void(* evalF)(const C_INT *, const double *, const double *, double *)
double dbnorm_(C_INT *n, double *a, C_INT *nra, C_INT *ml, C_INT *mu, double *w)
void(* evalJ)(const C_INT *, const double *, const double *, const C_INT *, const C_INT *, double *, const C_INT *)
#define dgefa_(__a, __lda, __n, __ipvt, __info)
#define dgbfa_(__ab, __lda, __n, __ml, __mu, __ipvt, __info)
double dfnorm_(C_INT *n, double *a, double *w)
C_INT dprja_(C_INT *neq, double *y, double *yh, C_INT *nyh, double *ewt, double *ftem, double *savf, double *wm, C_INT *iwm, evalF f, evalJ jac)
double dmnorm_(C_INT *n, double *v, double *w)