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)
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)
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)
double precision function ddwnrm(NEQ, V, RWT, RPAR, IPAR)