5 SUBROUTINE ddasik(X,Y,YPRIME,NEQ,ICOPT,ID,RES,JACK,PSOL,H,WT,
6 * JSKIP,RPAR,IPAR,SAVR,DELTA,R,YIC,YPIC,PWK,WM,IWM,CJ,UROUND,
7 * EPLI,SQRTN,RSQRTN,EPCON,RATEMX,STPTOL,JFLG,
8 * ICNFLG,ICNSTR,IERNLS)
96 IMPLICIT DOUBLE PRECISION(a-h,o-z)
97 dimension y(*),yprime(*),id(*),wt(*),icnstr(*)
98 dimension savr(*),delta(*),r(*),yic(*),ypic(*),pwk(*)
100 EXTERNAL res, jack, psol
102 parameter(lnre=12, lnje=13, llocwp=29, llciwp=30)
103 parameter(lmxnit=32, lmxnj=33)
119 iwm(lnre) = iwm(lnre) + 1
120 CALL res(x,y,yprime,cj,delta,ires,rpar,ipar)
121 IF (ires .LT. 0)
GO TO 370
135 IF (jflg .EQ. 1 .AND. jskip .EQ. 0)
THEN 137 iwm(lnje)=iwm(lnje)+1
138 CALL jack (res, ires, neq, x, y, yprime, wt, delta, r, h, cj,
139 * wm(lwp), iwm(liwp), ierpj, rpar, ipar)
140 IF (ires .LT. 0 .OR. ierpj .NE. 0)
GO TO 370
146 CALL dnsik(x,y,yprime,neq,icopt,id,res,psol,wt,rpar,ipar,
147 * savr,delta,r,yic,ypic,pwk,wm,iwm,cj,sqrtn,rsqrtn,
148 * eplin,epcon,ratemx,mxnit,stptol,icnflg,icnstr,iernew)
150 IF (iernew .EQ. 1 .AND. nj .LT. mxnj .AND. jflg .EQ. 1)
THEN 157 CALL dcopy (neq, savr, 1, delta, 1)
161 IF (iernew .NE. 0)
GO TO 380
169 IF (ires .LE. -2) iernls = -1
172 380 iernls =
min(iernew,2)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
charNDArray min(char d, const charNDArray &m)