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)
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)
calling an anonymous function involves an overhead quite comparable to the overhead of an m file function Passing a handle to a built in function is because the interpreter is not involved in the internal loop For a
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
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)