37 #ifndef TEMPLATE_LAPACK_LARRK_HEADER 38 #define TEMPLATE_LAPACK_LARRK_HEADER 42 Treal *
gu, Treal *d__, Treal *e2, Treal *pivmin,
43 Treal *reltol, Treal *w, Treal *werr,
integer *info)
52 Treal mid, eps, tmp1, tmp2,
left, atoli,
right;
151 atoli = *pivmin * 4.;
154 left = *gl - tnorm * 2. * eps * *n - *pivmin * 4.;
155 right = *
gu + tnorm * 2. * eps * *n + *pivmin * 4.;
166 d__1 =
maxMACRO(atoli,*pivmin), d__2 = rtoli * tmp2;
189 for (i__ = 2; i__ <= i__1; ++i__) {
190 tmp1 = d__[i__] - e2[i__ - 1] / tmp1 - mid;
static const real gu
Definition: fun-pz81.c:68
#define absMACRO(x)
Definition: template_blas_common.h:47
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:41
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
Treal template_blas_log(Treal x)
Treal template_lapack_lamch(const char *cmach, Treal dummyReal)
Definition: template_lapack_lamch.h:202