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;
static const real gu
Definition: fun-pz81.c:71
#define absMACRO(x)
Definition: template_blas_common.h:45
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:39
int integer
Definition: template_blas_common.h:38
#define maxMACRO(a, b)
Definition: template_blas_common.h:43
Treal template_blas_log(Treal x)
Treal template_lapack_lamch(const char *cmach, Treal dummyReal)
Definition: template_lapack_lamch.h:199