5 SUBROUTINE dnsid(X,Y,YPRIME,NEQ,ICOPT,ID,RES,WT,RPAR,IPAR,
6 * delta,r,yic,ypic,wm,iwm,cj,epcon,ratemx,maxit,stptol,
7 * icnflg,icnstr,iernew)
77 IMPLICIT DOUBLE PRECISION(a-h,o-z)
84 parameter(lnni=19, llsoff=35)
96 CALL
dslvd(neq, delta, wm, iwm)
100 delnrm =
ddwnrm(neq,delta,wt,rpar,ipar)
102 IF (fnrm .LE. epcon)
RETURN
107 iwm(lnni) = iwm(lnni) + 1
113 CALL
dlinsd(neq, y, x, yprime, cj, delta, delnrm, wt, lsoff,
114 * stptol, iret, res, ires, wm, iwm, fnrm, icopt, id,
115 * r, yic, ypic, icnflg, icnstr, rlx, rpar, ipar)
120 IF (iret .NE. 0) go to 390
124 IF (fnrm .LE. epcon)
RETURN
130 IF (m .GE. maxit) go to 380
135 CALL dcopy(neq, r, 1, delta, 1)
141 380
IF (rate .LE. ratemx)
THEN
148 390
IF (ires .LE. -2)
THEN
subroutine dnsid(X, Y, YPRIME, NEQ, ICOPT, ID, RES, WT, RPAR, IPAR, DELTA, R, YIC, YPIC, WM, IWM, CJ, EPCON, RATEMX, MAXIT, STPTOL, ICNFLG, ICNSTR, IERNEW)
std::string dimension(void) const
subroutine dlinsd(NEQ, Y, T, YPRIME, CJ, P, PNRM, WT, LSOFF, STPTOL, IRET, RES, IRES, WM, IWM, FNRM, ICOPT, ID, R, YNEW, YPNEW, ICNFLG, ICNSTR, RLX, RPAR, IPAR)
subroutine dslvd(NEQ, DELTA, WM, IWM)
double precision function ddwnrm(NEQ, V, RWT, RPAR, IPAR)