1 SUBROUTINE solsy (WM, IWM, X, TEM)
5 1 icf, ierpj, iersl, jcur, jstart, kflag, l, meth, miter,
6 2 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
7 INTEGER I, MEBAND, ML, MU
8 DOUBLE PRECISION WM, X, TEM
9 DOUBLE PRECISION ROWNS,
10 1 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
11 DOUBLE PRECISION DI, HL0, PHL0, R
13 COMMON /ls0001/ rowns(209),
14 2 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
15 3 iownd(14), iowns(6),
16 4 icf, ierpj, iersl, jcur, jstart, kflag, l, meth, miter,
17 5 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
43 go to(100, 100, 300, 400, 400), miter
44 100 CALL dgetrs(
'N', n, 1, wm(3), n, iwm(21), x, n, inlpck)
50 IF (hl0 .EQ. phl0) go to 330
53 di = 1.0d0 - r*(1.0d0 - 1.0d0/wm(i+2))
54 IF (dabs(di) .EQ. 0.0d0) go to 390
55 320 wm(i+2) = 1.0d0/di
57 340 x(i) = wm(i+2)*x(i)
64 meband = 2*ml + mu + 1
65 CALL dgbtrs(
'N', n, ml, mu, 1, wm(3), meband, iwm(21), x, n,
std::string dimension(void) const
subroutine solsy(WM, IWM, X, TEM)