73 #define subscr_1 (mpsubscr_->_1)
74 #define subscr_2 (mpsubscr_->_2)
75 #define subscr_3 (mpsubscr_->_3)
120 C_INT maxfun, msglvl;
217 CTruncatedNewton::lmqn_(ierror, n, &x[1], f, &g[1], &w[1], lw, sfun, &msglvl, &maxit,
218 &maxfun, &eta, &stepmx, &accrcy, &xtol);
271 C_INT maxfun, msglvl;
383 , &ipivot[1], &msglvl, &maxit, &maxfun, &eta, &stepmx, &accrcy, &
415 for (i__ = 1; i__ <= i__1; ++i__)
457 C_FLOAT64 gnorm, ftest, fstop, pnorm, rteps, xnorm;
460 C_INT icycle, nlincg, nfeval;
523 chkucp_(&
subscr_1.lwtest, maxfun, &nwhy, n, &alpha, &epsmch, eta, &peps, &
524 rteps, &rtol, &rtolsq, stepmx, &ftest, xtol, &xnorm, &x[1], lw, &
525 small, &tiny, accrcy);
549 ftest = one + fabs(fnew);
551 if (gtg < epsmch * 1e-4 * ftest * ftest)
559 toleps = rtol + rteps;
560 rtleps = rtolsq + epsmch;
571 for (i__ = 1; i__ <= i__1; ++i__)
585 1], &w[
subscr_1.lzk], n, &w[1], lw, &niter, maxit, &nfeval, &
586 nmodif, &nlincg, &upd1, &yksk, &gsk, &yrsr, &lreset, sfun, &
587 c_false, &ipivot, accrcy, >pnew, &gnorm, &xnorm);
600 reltol = rteps * (xnorm + one) / pe;
601 fabstol = -epsmch * ftest / (oldgtp - epsmch);
606 tnytol = epsmch * (xnorm + one) / pe;
611 alpha =
step1_(&fnew, &fm, &oldgtp, &spe);
616 zero, &spe, &w[
subscr_1.lpk], &oldgtp, &x[1], &fnew, &alpha, &g[1]
617 , &numf, &nwhy, &w[1], lw);
642 if (nwhy == 0 || nwhy == 2)
658 (*sfun)(n, &x[1], &fnew, &g[1]);
666 if (nftotl > *maxfun)
676 difnew = oldf - fnew;
686 if (difnew > difold * 2.)
691 if (difnew < difold * .5)
698 ftest = one + fabs(fnew);
703 if ((alpha * pnorm < toleps * (one + xnorm) &&
704 fabs(difnew) < rtleps * ftest &&
705 gtg < peps * ftest * ftest) ||
706 gtg < *accrcy * 1e-4 * ftest * ftest)
720 for (i__ = 1; i__ <= i__1; ++i__)
722 w[iyk] = g[i__] - w[ioldg];
723 w[isk] = alpha * w[ipk];
736 if (icycle == nm1 || difnew < epsred * (fkeep - fnew))
761 1], &w[
subscr_1.lzk], n, &w[1], lw, &niter, maxit, &nfeval, &
762 nmodif, &nlincg, &upd1, &yksk, &gsk, &yrsr, &lreset, sfun, &
763 c_false, &ipivot, accrcy, >pnew, &gnorm, &xnorm);
861 C_INT icycle, nlincg, nfeval;
910 crash_(n, &x[1], &ipivot[1], &low[1], &up[1], &ier);
943 chkucp_(&
subscr_1.lwtest, maxfun, &nwhy, n, &alpha, &epsmch, eta, &peps, &
944 rteps, &rtol, &rtolsq, stepmx, &ftest, xtol, &xnorm, &x[1], lw, &
945 small, &tiny, accrcy);
964 for (i__ = 1; i__ <= i__1; ++i__)
966 if (ipivot[i__] == 2)
971 if (-ipivot[i__] * g[i__] >= 0.)
981 ztime_(n, &g[1], &ipivot[1]);
986 monit_(n, &x[1], &fnew, &g[1], &niter, &nftotl, &nfeval, &lreset, &
992 ftest = one + fabs(fnew);
994 if (gtg < epsmch * 1e-4 * ftest * ftest)
1002 toleps = rtol + rteps;
1003 rtleps = rtolsq + epsmch;
1014 for (i__ = 1; i__ <= i__1; ++i__)
1025 modet = *msglvl - 3;
1028 1], &w[
subscr_1.lzk], n, &w[1], lw, &niter, maxit, &nfeval, &
1029 nmodif, &nlincg, &upd1, &yksk, &gsk, &yrsr, &lreset, sfun, &
1030 c_true, &ipivot[1], accrcy, >pnew, &gnorm, &xnorm);
1043 pe = pnorm + epsmch;
1047 reltol = rteps * (xnorm + one) / pe;
1048 fabstol = -epsmch * ftest / (oldgtp - epsmch);
1053 tnytol = epsmch * (xnorm + one) / pe;
1054 stpmax_(stepmx, &pe, &spe, n, &x[1], &w[
subscr_1.lpk], &ipivot[1], &low[1]
1059 alpha =
step1_(&fnew, &fm, &oldgtp, &spe);
1064 zero, &spe, &w[
subscr_1.lpk], &oldgtp, &x[1], &fnew, &alpha, &g[1]
1065 , &numf, &nwhy, &w[1], lw);
1068 if ((d__1 = alpha - spe, fabs(d__1)) > epsmch * 10.)
1075 modz_(n, &x[1], &w[
subscr_1.lpk], &ipivot[1], &epsmch, &low[1], &up[1], &
1098 monit_(n, &x[1], &fnew, &g[1], &niter, &nftotl, &nfeval, &lreset, &
1107 if (nwhy == 0 || nwhy == 2)
1123 (*sfun)(n, &x[1], &fnew, &g[1]);
1131 if (nftotl > *maxfun)
1141 difnew = oldf - fnew;
1151 if (difnew > difold * 2.)
1156 if (difnew < difold * .5)
1166 ftest = one + fabs(fnew);
1171 cnvtst_(&conv, &alpha, &pnorm, &toleps, &xnorm, &difnew, &rtleps, &ftest,
1172 >g, &peps, &epsmch, >pnew, &fnew, &flast, &g[1], &ipivot[1],
1180 ztime_(n, &g[1], &ipivot[1]);
1196 for (i__ = 1; i__ <= i__1; ++i__)
1198 w[iyk] = g[i__] - w[ioldg];
1199 w[isk] = alpha * w[ipk];
1212 if (icycle == nm1 || difnew < epsred * (fkeep - fnew))
1236 if (upd1 && *msglvl >= 3)
1243 if (newcon && *msglvl >= 3)
1250 modet = *msglvl - 3;
1253 1], &w[
subscr_1.lzk], n, &w[1], lw, &niter, maxit, &nfeval, &
1254 nmodif, &nlincg, &upd1, &yksk, &gsk, &yrsr, &lreset, sfun, &
1255 c_true, &ipivot[1], accrcy, >pnew, &gnorm, &xnorm);
1304 monit_(n, &x[1], f, &g[1], &niter, &nftotl, &nfeval, &ireset, &ipivot[
1343 for (i__ = 1; i__ <= i__1; ++i__)
1345 if (ipivot[i__] != 0)
1350 gtg += g[i__] * g[i__];
1383 for (i__ = 1; i__ <= i__1; ++i__)
1385 if (ipivot[i__] != 0)
1409 --up; --low; --ipivot; --p; --x;
1412 *spe = *stepmx / *pe;
1416 for (i__ = 1; i__ <= i__1; ++i__)
1418 if (ipivot[i__] != 0)
goto L10;
1420 if (p[i__] == 0.)
goto L10;
1422 if (p[i__] > 0.)
goto L5;
1424 t = low[i__] - x[i__];
1426 if (t > *spe * p[i__]) *spe = t / p[i__];
1430 t = up[i__] - x[i__];
1432 if (t < *spe * p[i__]) *spe = t / p[i__];
1464 for (i__ = 1; i__ <= i__1; ++i__)
1466 if (ipivot[i__] != 0)
1481 tol = *epsmch * 10. * ((d__1 = low[i__], fabs(d__1)) + 1.);
1483 if (x[i__] - low[i__] > tol)
1493 tol = *epsmch * 10. * ((d__1 = up[i__], fabs(d__1)) + 1.);
1495 if (up[i__] - x[i__] > tol)
1536 ltest = *flast - *fnew <= *gtpnew * -.5;
1539 for (i__ = 1; i__ <= i__1; ++i__)
1541 if (ipivot[i__] == 0 || ipivot[i__] == 2)
1546 t = -ipivot[i__] * g[i__];
1583 if ((*alpha * *pnorm >= *toleps * (one + *xnorm) || fabs(*difnew) >= *
1584 rtleps * *ftest || *gtg >= *peps * *ftest * *ftest) && *gtg >= *
1585 accrcy * 1e-4 * *ftest * *ftest)
1623 for (i__ = 1; i__ <= i__1; ++i__)
1625 if (x[i__] < low[i__])
1630 if (x[i__] > up[i__])
1637 if (x[i__] == low[i__])
1642 if (x[i__] == up[i__])
1647 if (up[i__] == low[i__])
1652 if (low[i__] > up[i__])
1764 CTruncatedNewton::initpc_(&diagb[1], &emat[1], n, &w[1], lw, modet, upd1, yksk, gsk, yrsr,
1768 for (i__ = 1; i__ <= i__1; ++i__)
1782 for (k = 1; k <= i__1; ++k)
1798 ztime_(n, &r__[1], &ipivot[1]);
1801 CTruncatedNewton::msolve_(&r__[1], &zk[1], n, &w[1], lw, upd1, yksk, gsk, yrsr, lreset,
1806 ztime_(n, &zk[1], &ipivot[1]);
1811 if (rz / rhsnrm < tol)
1828 for (i__ = 1; i__ <= i__2; ++i__)
1830 v[i__] = zk[i__] + beta * v[i__];
1836 ztime_(n, &v[1], &ipivot[1]);
1840 &delta, accrcy, xnorm);
1844 ztime_(n, &gv[1], &ipivot[1]);
1850 if (vgv / rhsnrm < tol)
1855 ndia3_(n, &emat[1], &v[1], &gv[1], &r__[1], &vgv, modet);
1879 qnew = (*gtp + pr) * .5;
1880 qtest = k * (1. - qold / qnew);
1941 CTruncatedNewton::msolve_(&g[1], &zsol[1], n, &w[1], lw, upd1, yksk, gsk, yrsr, lreset, &
1947 ztime_(n, &zsol[1], &ipivot[1]);
1982 ztime_(n, &zsol[1], &ipivot[1]);
2027 for (i__ = 1; i__ <= i__1; ++i__)
2029 e[i__] = e[i__] - r__[i__] * r__[i__] / vr + gv[i__] * gv[i__] / *vgv;
2070 for (i__ = 1; i__ <= i__1; ++i__)
2108 ybnd = *scxbd + *xmin;
2159 d__ = (d__1 = *fnew - *fm, fabs(d__1));
2162 if (d__ * 2. <= -(*gtp) && d__ >= epsmch)
2164 alpha = d__ * -2. / *gtp;
2186 ret_val = std::numeric_limits< C_FLOAT64 >::epsilon();
2207 *small = *epsmch * *epsmch;
2210 *rteps = sqrt(*epsmch);
2213 if (fabs(*rtol) < *accrcy)
2215 *rtol = *rteps * 10.;
2220 if (*lw < *lwtest || *n < 1 || *rtol < 0. || *eta >= 1. || *eta < 0. || *
2221 stepmx < *rtol || *maxfun < 1)
2230 *rtolsq = *rtol * *rtol;
2231 *peps = pow(*accrcy,
c_b246);
2255 (*sfun)(n, &x[1], fnew, &g[1]);
2294 *delta = sqrt(*accrcy) * (*xnorm + 1.);
2301 for (i__ = 1; i__ <= i__1; ++i__)
2303 w[ihg] = x[i__] + *delta * v[i__];
2308 (*sfun)(n, &w[
subscr_2.lhg], &f, &gv[1]);
2311 for (i__ = 1; i__ <= i__1; ++i__)
2313 gv[i__] = (gv[i__] - g[i__]) * dinv;
2335 gsk, yrsr, lreset, first);
2351 C_FLOAT64 ykhyk = 0.0, ykhyr = 0.0, yrhyr = 0.0, rdiagb;
2390 for (i__ = 1; i__ <= i__1; ++i__)
2392 rdiagb = 1. / diagb[i__];
2393 hg[i__] = g[i__] * rdiagb;
2397 hyk[i__] = yk[i__] * rdiagb;
2402 hyr[i__] = yr[i__] * rdiagb;
2426 ssbfgs_(n, &one, &sr[1], &yr[1], &hg[1], &hyr[1], yrsr, &yrhyr, &gsr, &
2431 ssbfgs_(n, &one, &sr[1], &yr[1], &hyk[1], &hyr[1], yrsr, &yrhyr, &
2432 yksr, &ykhyr, &hyk[1]);
2437 ssbfgs_(n, &one, &sk[1], &yk[1], &hg[1], &hyk[1], yksk, &ykhyk, gsk, &
2446 for (i__ = 1; i__ <= i__1; ++i__)
2448 rdiagb = 1. / diagb[i__];
2449 hg[i__] = g[i__] * rdiagb;
2453 hyk[i__] = yk[i__] * rdiagb;
2465 ssbfgs_(n, &one, &sk[1], &yk[1], &hg[1], &hyk[1], yksk, &ykhyk, gsk, &
2471 for (i__ = 1; i__ <= i__1; ++i__)
2474 y[i__] = g[i__] / diagb[i__];
2503 delta = (*gamma * *yjhyj / *yjsj + 1.) * *vsj / *yjsj - *gamma * *vhyj / *
2505 beta = -(*gamma) * *vsj / *yjsj;
2508 for (i__ = 1; i__ <= i__1; ++i__)
2510 hjp1v[i__] = *gamma * hjv[i__] + delta * sj[i__] + beta * hjyj[i__];
2529 initp3_(&diagb[1], &emat[1], n, lreset, yksk, yrsr, &w[
subscr_2.lhyk], &w[
2581 for (i__ = 1; i__ <= i__1; ++i__)
2583 bsk[i__] = diagb[i__] * sr[i__];
2592 for (i__ = 1; i__ <= i__1; ++i__)
2595 bsk[i__] = td * sk[i__] - bsk[i__] * srds / sds + yr[i__] * yrsk / *
2597 emat[i__] = td - td * td * sr[i__] * sr[i__] / sds + yr[i__] * yr[i__]
2605 for (i__ = 1; i__ <= i__1; ++i__)
2607 emat[i__] = emat[i__] - bsk[i__] * bsk[i__] / sds + yk[i__] * yk[i__]
2616 for (i__ = 1; i__ <= i__1; ++i__)
2618 bsk[i__] = diagb[i__] * sk[i__];
2625 for (i__ = 1; i__ <= i__1; ++i__)
2628 emat[i__] = td - td * td * sk[i__] * sk[i__] / sds + yk[i__] * yk[i__]
2647 for (i__ = 1; i__ <= i__1; ++i__)
2677 for (i__ = 1; i__ <= 14; ++i__)
2679 subscr_3.lsub[i__ - 1] = (i__ - 1) * *n + 1;
2736 rtsmll = sqrt(*small);
2742 fpresn = *epsmch * 10.;
2767 xbnd, &u, &fu, &gu, &xmin, &fmin, &gmin, &xw, &fw, &gw, &a, &b, &
2768 oldf, &b1, &scxbnd, &e, &step, &factor, &braktd, >est1, >est2,
2769 &tol, &ientry, &itest);
2772 if (lsprnt >= nprnt)
2774 lsout_(&ientry, &itest, &xmin, &fmin, &gmin, &xw, &fw, &gw, &u, &a, &
2775 b, &tol, reltol, &scxbnd, xbnd);
2790 for (i__ = 1; i__ <= i__1; ++i__)
2792 w[l] = x[i__] + ualpha * p[i__];
2797 (*sfun)(n, &w[lx], &fu, &w[lg]);
2805 if (fu <= fmin && fu <= oldf - ualpha * gtest1)
2830 for (i__ = 1; i__ <= i__1; ++i__)
2832 x[i__] += *alpha * p[i__];
2855 C_FLOAT64 half, abgw, fabsr, five, zero, p, q, r__, s, scale,
2856 denom, three, a1, d1, d2, sumsq, point1, abgmin, chordm, eleven,
2898 if (*u <= zero || *xbnd <= *tnytol || *gu > zero)
2905 if (*xbnd < *fabstol)
2911 twotol = *tol + *tol;
2940 *b = *scxbnd + *reltol * fabs(*scxbnd) + *fabstol;
2946 *gtest1 = -(*rmu) * *gu;
2947 *gtest2 = -(*eta) * *gu;
2968 chordu = *oldf - (*xmin + *u) * *gtest1;
2981 chordm = *oldf - *xmin * *gtest1;
2983 denom = chordm - *fmin;
2985 if (fabs(denom) >= 1e-15)
2992 if (chordm - *fmin < 0.)
3001 *gu = *gmin * (chordu - *fu) / denom;
3004 *fu = half * *u * (*gmin + *gu) + *fmin;
3034 *tol = fabs(*xmin) * *reltol + *fabstol;
3057 twotol = *tol + *tol;
3058 xmidpt = half * (*a + *b);
3062 convrg = fabs(xmidpt) <= twotol - half * (*b - *a) ||
3063 (fabs(*gmin) <= *gtest2 &&
3065 ((d__1 = *xmin - *xbnd, fabs(d__1)) > * tol ||
3087 if ((d__1 = *oldf - *fw, fabs(d__1)) <= *fpresn * (one + fabs(*oldf)))
3092 *tol = point1 * *tol;
3099 *reltol = point1 * *reltol;
3100 *fabstol = point1 * *fabstol;
3101 twotol = point1 * twotol;
3110 if (fabs(*e) <= *tol)
3117 r__ = three * (*fmin - *fw) / *xw + *gmin + *gw;
3121 if (*gw == zero || *gmin == zero)
3130 abgmin = fabs(*gmin);
3131 s = sqrt(abgmin) * sqrt(abgw);
3133 if (*gw / abgw * *gmin > zero)
3159 sumsq = one + d__1 * d__1;
3169 if (fabsr > *rtsmll)
3171 p = fabsr * *rtsmll;
3178 sumsq = one + d__1 * d__1;
3183 sumsq = sqrt(sumsq);
3186 if (scale < *big / sumsq)
3196 q = sqrt((d__1 = r__ + s, fabs(d__1))) * sqrt((d__2 = r__ - s, fabs(d__2)));
3198 if (r__ >= s || r__ <= -s)
3216 s = *xw * (*gmin - r__ - q);
3217 q = *gw - *gmin + q + q;
3231 if (*b1 != *step || *braktd)
3248 *step = -(*factor) * *xw;
3250 if (*step > *scxbnd)
3255 if (*step != *scxbnd)
3257 *factor = five * *factor;
3267 if ((*a != zero || *xw >= zero) && (*b != zero || *xw <= zero))
3286 *step = five * d2 * (point1 + one / *u) / eleven;
3290 *step = half * d2 * sqrt(*u);
3311 if (fabs(s) <= (d__1 = half * q * r__, fabs(d__1)) || s <= q * a1 || s >= q
3323 if (*step - *a >= twotol && *b - *step >= twotol)
3346 if (*step < *scxbnd)
3355 *scxbnd -= (*reltol * fabs(*xbnd) + *fabstol) / (one + *reltol);
3359 if (fabs(*step) < *tol && *step < zero)
3364 if (fabs(*step) < *tol && *step >= zero)
3390 C_INT i__, m, ix, iy, mp1;
3406 if (*incx == 1 && *incy == 1)
3419 ix = (-(*n) + 1) * *incx + 1;
3424 iy = (-(*n) + 1) * *incy + 1;
3429 for (i__ = 1; i__ <= i__1; ++i__)
3453 for (i__ = 1; i__ <= i__1; ++i__)
3468 for (i__ = mp1; i__ <= i__1; i__ += 4)
3471 dy[i__ + 1] += dx[i__ + 1];
3472 dy[i__ + 2] += dx[i__ + 2];
3473 dy[i__ + 3] += dx[i__ + 3];
3482 return(pow(*ap, *bp));
C_FLOAT64 step1_(C_FLOAT64 *fnew, C_FLOAT64 *fm, C_FLOAT64 *gtp, C_FLOAT64 *smax)
int lsout_(C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *)
int linder_(C_INT *, FTruncatedNewton *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *, C_FLOAT64 *, C_INT *)
int daxpy_(integer *n, doublereal *da, doublereal *dx, integer *incx, doublereal *dy, integer *incy)
int dcopy_(integer *n, doublereal *dx, integer *incx, doublereal *dy, integer *incy)
int modz_(C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *)
int negvec_(C_INT *, C_FLOAT64 *)
int ndia3_(C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *)
int getptc_(C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *)
int initpc_(C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_INT *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *)
int crash_(C_INT *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_INT *)
int ssbfgs_(C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *)
int stpmax_(C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_FLOAT64 *)
doublereal dnrm2_(integer *n, doublereal *x, integer *incx)
doublereal ddot_(integer *n, doublereal *dx, integer *incx, doublereal *dy, integer *incy)
int dxpy_(C_INT *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_INT *)
int initp3_(C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *)
int mslv_(C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *)
int modlnp_(C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_INT *, C_INT *, C_INT *, C_INT *, C_INT *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, FTruncatedNewton *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *)
int monit_(C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *, C_INT *, C_INT *, C_INT *)
int tn_(C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, FTruncatedNewton *)
int setucr_(C_FLOAT64 *, C_INT *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, FTruncatedNewton *, C_FLOAT64 *, C_FLOAT64 *)
int lmqnbc_(C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, FTruncatedNewton *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *)
int chkucp_(C_INT *, C_INT *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *)
int gtims_(C_FLOAT64 *v, C_FLOAT64 *gv, C_INT *n, C_FLOAT64 *x, C_FLOAT64 *g, C_FLOAT64 *w, C_INT *, FTruncatedNewton *sfun, C_INT *first, C_FLOAT64 *delta, C_FLOAT64 *accrcy, C_FLOAT64 *xnorm)
int ztime_(C_INT *, C_FLOAT64 *, C_INT *)
int lmqn_(C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, FTruncatedNewton *, C_INT *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *)
C_FLOAT64 pow_dd(C_FLOAT64 *ap, C_FLOAT64 *bp)
int msolve_(C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_FLOAT64 *, C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *)
int cnvtst_(C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, C_INT *, C_FLOAT64 *)
int tnbc_(C_INT *, C_INT *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_FLOAT64 *, C_INT *, FTruncatedNewton *, C_FLOAT64 *, C_FLOAT64 *, C_INT *)