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;
int integer
Definition: template_blas_common.h:38
int template_lapack_larrc(const char *jobt, const integer *n, const Treal *vl, const Treal *vu, Treal *d__, Treal *e, Treal *pivmin, integer *eigcnt, integer *lcnt, integer *rcnt, integer *info)
Definition: template_lapack_larrc.h:39
bool logical
Definition: template_blas_common.h:39
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:44