35 #ifndef TEMPLATE_LAPACK_LARRC_HEADER
36 #define TEMPLATE_LAPACK_LARRC_HEADER
40 const Treal *vu, Treal *d__, Treal *e, Treal *pivmin,
49 Treal sl, su, tmp, tmp2;
137 lpivot = d__[1] - *vl;
138 rpivot = d__[1] - *vu;
146 for (i__ = 1; i__ <= i__1; ++i__) {
150 lpivot = d__[i__ + 1] - *vl - tmp / lpivot;
151 rpivot = d__[i__ + 1] - *vu - tmp / rpivot;
165 for (i__ = 1; i__ <= i__1; ++i__) {
166 lpivot = d__[i__] + sl;
167 rpivot = d__[i__] + su;
174 tmp = e[i__] * d__[i__] * e[i__];
180 sl = sl * tmp2 - *vl;
187 su = su * tmp2 - *vu;
191 lpivot = d__[*n] + sl;
192 rpivot = d__[*n] + su;
200 *eigcnt = *rcnt - *lcnt;