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) {