37 #ifndef TEMPLATE_LAPACK_LARRC_HEADER 38 #define TEMPLATE_LAPACK_LARRC_HEADER 42 const Treal *vu, Treal *d__, Treal *e, Treal *pivmin,
51 Treal sl, su, tmp, tmp2;
139 lpivot = d__[1] - *vl;
140 rpivot = d__[1] - *vu;
148 for (i__ = 1; i__ <= i__1; ++i__) {
152 lpivot = d__[i__ + 1] - *vl - tmp / lpivot;
153 rpivot = d__[i__ + 1] - *vu - tmp / rpivot;
167 for (i__ = 1; i__ <= i__1; ++i__) {
168 lpivot = d__[i__] + sl;
169 rpivot = d__[i__] + su;
176 tmp = e[i__] * d__[i__] * e[i__];
182 sl = sl * tmp2 - *vl;
189 su = su * tmp2 - *vu;
193 lpivot = d__[*n] + sl;
194 rpivot = d__[*n] + su;
202 *eigcnt = *rcnt - *lcnt;
int integer
Definition: template_blas_common.h:40
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:41
bool logical
Definition: template_blas_common.h:41
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46