1 SUBROUTINE ssolsy (WM, IWM, X, TEM)
49 1 icf, ierpj, iersl, jcur, jstart, kflag, l,
50 2 lyh, lewt, lacor, lsavf, lwm, liwm, meth, miter,
51 3 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
53 1 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
54 COMMON /sls001/ rowns(209),
55 1 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
57 3 icf, ierpj, iersl, jcur, jstart, kflag, l,
58 4 lyh, lewt, lacor, lsavf, lwm, liwm, meth, miter,
59 5 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
60 INTEGER I, MEBAND, ML, MU
65 go to(100, 100, 300, 400, 400), miter
66 100 CALL sgetrs(
'N', n, 1, wm(3), n, iwm(21), x, n, inlpck)
72 IF (hl0 .EQ. phl0) go to 330
75 di = 1.0e0 - r*(1.0e0 - 1.0e0/wm(i+2))
76 IF (
abs(di) .EQ. 0.0e0) go to 390
77 320 wm(i+2) = 1.0e0/di
79 340 x(i) = wm(i+2)*x(i)
86 meband = 2*ml + mu + 1
87 CALL sgbtrs(
'N', n, ml, mu, 1, wm(3), meband, iwm(21), x, n,
subroutine ssolsy(WM, IWM, X, TEM)
std::string dimension(void) const