35 #ifndef TEMPLATE_LAPACK_LARRK_HEADER
36 #define TEMPLATE_LAPACK_LARRK_HEADER
40 Treal *
gu, Treal *d__, Treal *e2, Treal *pivmin,
41 Treal *reltol, Treal *w, Treal *werr,
integer *info)
50 Treal mid, eps, tmp1, tmp2,
left, atoli,
right;
149 atoli = *pivmin * 4.;
152 left = *gl - tnorm * 2. * eps * *n - *pivmin * 4.;
153 right = *gu + tnorm * 2. * eps * *n + *pivmin * 4.;
164 d__1 =
maxMACRO(atoli,*pivmin), d__2 = rtoli * tmp2;
176 mid = (left +
right) * .5;
187 for (i__ = 2; i__ <= i__1; ++i__) {
188 tmp1 = d__[i__] - e2[i__ - 1] / tmp1 - mid;
207 *w = (left +
right) * .5;