37 #ifndef TEMPLATE_LAPACK_LAGTF_HEADER 38 #define TEMPLATE_LAPACK_LAGTF_HEADER 43 Treal *b, Treal *c__,
const Treal *tol, Treal *d__,
143 Treal scale1, scale2;
145 Treal eps, piv1, piv2;
180 for (k = 1; k <= i__1; ++k) {
182 scale2 = (d__1 = c__[k],
absMACRO(d__1)) + (d__2 = a[k + 1],
absMACRO(d__2));
184 scale2 += (d__1 = b[k + 1],
absMACRO(d__1));
189 piv1 = (d__1 = a[k],
absMACRO(d__1)) / scale1;
199 piv2 = (d__1 = c__[k],
absMACRO(d__1)) / scale2;
204 a[k + 1] -= c__[k] * b[k];
210 mult = a[k] / c__[k];
213 a[k + 1] = b[k] - mult * temp;
216 b[k + 1] = -mult * d__[k];
222 if (
maxMACRO(piv1,piv2) <= tl && in[*n] == 0) {
227 if ((d__1 = a[*n],
absMACRO(d__1)) <= scale1 * tl && in[*n] == 0) {
#define absMACRO(x)
Definition: template_blas_common.h:47
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
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:42
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
Treal template_lapack_lamch(const char *cmach, Treal dummyReal)
Definition: template_lapack_lamch.h:202