5 SUBROUTINE dmatd(NEQ,X,Y,YPRIME,DELTA,CJ,H,IER,EWT,E,
6 * wm,iwm,res,ires,uround,jacd,rpar,ipar)
64 IMPLICIT DOUBLE PRECISION(a-
h,o-z)
69 parameter(lml=1, lmu=2, lmtype=4, lnre=12, lnpd=22, llciwp=30)
74 go
to(100,200,300,400,500),mtype
82 CALL jacd(x,y,yprime,wm,cj,rpar,ipar)
99 yprime(i)=yprime(i)+cj*
del
100 iwm(lnre)=iwm(lnre)+1
101 CALL res(x,y,yprime,cj,
e,ires,rpar,ipar)
102 IF (ires .LT. 0)
RETURN
105 220 wm(nrow+l)=(
e(l)-delta(l))*delinv
114 230 CALL dgetrf( neq, neq, wm, neq, iwm(lipvt), ier)
128 CALL jacd(x,y,yprime,wm,cj,rpar,ipar)
129 meband=2*iwm(lml)+iwm(lmu)+1
135 500 mband=iwm(lml)+iwm(lmu)+1
137 meband=mband+iwm(lml)
148 wm(ipsave+k)=yprime(n)
154 510 yprime(n)=yprime(n)+cj*
del
155 iwm(lnre)=iwm(lnre)+1
156 CALL res(x,y,yprime,cj,
e,ires,rpar,ipar)
157 IF (ires .LT. 0)
RETURN
161 yprime(n)=wm(ipsave+k)
167 i1=max0(1,(n-iwm(lmu)))
168 i2=min0(neq,(n+iwm(lml)))
171 520 wm(ii+i)=(
e(i)-delta(i))*delinv
178 550 CALL dgbtrf(neq, neq, iwm(lml), iwm(lmu), wm, meband,
void del(const std::string &file)
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)
charNDArray max(char d, const charNDArray &m)
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup abs(local error in x(i))<
subroutine dmatd(NEQ, X, Y, YPRIME, DELTA, CJ, H, IER, EWT, E, WM, IWM, RES, IRES, UROUND, JACD, RPAR, IPAR)
octave_value sqrt(void) const