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)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
charNDArray min(char d, const charNDArray &m)