5 SUBROUTINE ddasid(X,Y,YPRIME,NEQ,ICOPT,ID,RES,JACD,PDUM,H,WT,
6 * jsdum,rpar,ipar,dumsvr,delta,r,yic,ypic,dumpwk,wm,iwm,cj,uround,
7 * dume,dums,dumr,epcon,ratemx,stptol,jfdum,
8 * icnflg,icnstr,iernls)
92 IMPLICIT DOUBLE PRECISION(a-h,o-z)
93 dimension y(*),yprime(*),id(*),wt(*),icnstr(*)
98 parameter(lnre=12, lnje=13, lmxnit=32, lmxnj=33)
111 iwm(lnre) = iwm(lnre) + 1
112 CALL res(x,y,yprime,cj,delta,ires,rpar,ipar)
113 IF (ires .LT. 0) go to 370
129 iwm(lnje)=iwm(lnje)+1
130 CALL
dmatd(neq,x,y,yprime,delta,cj,h,ierj,wt,r,
131 * wm,iwm,res,ires,uround,jacd,rpar,ipar)
132 IF (ires .LT. 0 .OR. ierj .NE. 0) go to 370
136 CALL
dnsid(x,y,yprime,neq,icopt,id,res,wt,rpar,ipar,delta,r,
137 * yic,ypic,wm,iwm,cj,epcon,ratemx,mxnit,stptol,
138 * icnflg,icnstr,iernew)
140 IF (iernew .EQ. 1 .AND. nj .LT. mxnj)
THEN
146 iwm(lnre)=iwm(lnre)+1
147 CALL res(x,y,yprime,cj,delta,ires,rpar,ipar)
148 IF (ires .LT. 0) go to 370
152 IF (iernew .NE. 0) go to 380
161 IF (ires .LE. -2) iernls = -1
164 380 iernls =
min(iernew,2)
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 ddasid(X, Y, YPRIME, NEQ, ICOPT, ID, RES, JACD, PDUM, H, WT, JSDUM, RPAR, IPAR, DUMSVR, DELTA, R, YIC, YPIC, DUMPWK, WM, IWM, CJ, UROUND, DUME, DUMS, DUMR, EPCON, RATEMX, STPTOL, JFDUM, ICNFLG, ICNSTR, IERNLS)
subroutine dmatd(NEQ, X, Y, YPRIME, DELTA, CJ, H, IER, EWT, E, WM, IWM, RES, IRES, UROUND, JACD, RPAR, IPAR)
charNDArray min(char d, const charNDArray &m)