42 Treal *
gu, Treal *d__, Treal *e2, Treal *pivmin,
43 Treal *reltol, Treal *w, Treal *werr,
integer *info)
52 Treal mid, eps, tmp1, tmp2, left, atoli, right;
151 atoli = *pivmin * 4.;
154 left = *gl - tnorm * 2. * eps * *n - *pivmin * 4.;
155 right = *
gu + tnorm * 2. * eps * *n + *pivmin * 4.;
161 tmp1 = (d__1 = right - left,
absMACRO(d__1));
166 d__1 =
maxMACRO(atoli,*pivmin), d__2 = rtoli * tmp2;
178 mid = (left + right) * .5;
189 for (i__ = 2; i__ <= i__1; ++i__) {
190 tmp1 = d__[i__] - e2[i__ - 1] / tmp1 - mid;
209 *w = (left + right) * .5;
210 *werr = (d__1 = right - left,
absMACRO(d__1)) * .5;
int template_lapack_larrk(integer *n, integer *iw, Treal *gl, Treal *gu, Treal *d__, Treal *e2, Treal *pivmin, Treal *reltol, Treal *w, Treal *werr, integer *info)
Definition template_lapack_larrk.h:41