35 #ifndef TEMPLATE_LAPACK_LAGTF_HEADER 36 #define TEMPLATE_LAPACK_LAGTF_HEADER 41 Treal *b, Treal *c__,
const Treal *tol, Treal *d__,
141 Treal scale1, scale2;
143 Treal eps, piv1, piv2;
178 for (k = 1; k <= i__1; ++k) {
180 scale2 = (d__1 = c__[k],
absMACRO(d__1)) + (d__2 = a[k + 1],
absMACRO(d__2));
182 scale2 += (d__1 = b[k + 1],
absMACRO(d__1));
187 piv1 = (d__1 = a[k],
absMACRO(d__1)) / scale1;
197 piv2 = (d__1 = c__[k],
absMACRO(d__1)) / scale2;
202 a[k + 1] -= c__[k] * b[k];
208 mult = a[k] / c__[k];
211 a[k + 1] = b[k] - mult * temp;
214 b[k + 1] = -mult * d__[k];
220 if (
maxMACRO(piv1,piv2) <= tl && in[*n] == 0) {
225 if ((d__1 = a[*n],
absMACRO(d__1)) <= scale1 * tl && in[*n] == 0) {
#define absMACRO(x)
Definition: template_blas_common.h:45
int integer
Definition: template_blas_common.h:38
#define maxMACRO(a, b)
Definition: template_blas_common.h:43
int template_lapack_lagtf(const integer *n, Treal *a, const Treal *lambda, Treal *b, Treal *c__, const Treal *tol, Treal *d__, integer *in, integer *info)
Definition: template_lapack_lagtf.h:40
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:144
Treal template_lapack_lamch(const char *cmach, Treal dummyReal)
Definition: template_lapack_lamch.h:199