1 SUBROUTINE ddaslv (NEQ, DELTA, WM, IWM)
32 DOUBLE PRECISION DELTA(*), WM(*)
34 EXTERNAL dgbtrs, dgetrs
36 INTEGER LIPVT, LML, LMU, LMTYPE, MEBAND, MTYPE, NPD, INFO
45 go
to(100,100,300,400,400),mtype
48 100 CALL dgetrs(
'N', neq, 1, wm(npd), neq, iwm(lipvt), delta, neq,
57 400 meband=2*iwm(lml)+iwm(lmu)+1
58 CALL dgbtrs(
'N', neq, iwm(lml), iwm(lmu), 1, wm(npd), meband,
59 * iwm(lipvt), delta, neq, inlpck)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
subroutine ddaslv(NEQ, DELTA, WM, IWM)