26 #define dls001_1 (mdls001_._1)
27 #define dls001_2 (mdls001_._2)
28 #define dls001_3 (mdls001_._3)
30 #define dlsa01_1 (mdlsa01_._1)
31 #define dlsa01_2 (mdlsa01_._2)
32 #define dlsa01_3 (mdlsa01_._3)
136 for (i__ = 1; i__ <= i__1; ++i__)
138 di = 1. - r__ * (1. - 1. / wm[i__ + 2]);
146 wm[i__ + 2] = 1. / di;
152 for (i__ = 1; i__ <= i__1; ++i__)
155 x[i__] = wm[i__ + 2] * x[i__];
166 meband = (ml << 1) + mu + 1;
C_INT dsolsy_(double *wm, C_INT *iwm, double *x, double *tem)
#define dgesl_(__a, __lda, __n, __ipvt, __b, __job)
#define dgbsl_(__ab, __lda, __n, __ml, __mu, __ipvt, __b, __job)