37 #ifndef TEMPLATE_LAPACK_LASQ5_HEADER 38 #define TEMPLATE_LAPACK_LASQ5_HEADER 42 integer *pp, Treal *tau, Treal *dmin__, Treal *dmin1,
43 Treal *dmin2, Treal *dn, Treal *dnm1, Treal *dnm2,
131 if (*n0 - *i0 - 1 <= 0) {
135 j4 = (*i0 << 2) + *pp - 3;
137 d__ = z__[j4] - *tau;
146 i__1 = ( *n0 - 3 ) << 2;
147 for (j4 = *i0 << 2; j4 <= i__1; j4 += 4) {
148 z__[j4 - 2] = d__ + z__[j4 - 1];
149 temp = z__[j4 + 1] / z__[j4 - 2];
150 d__ = d__ * temp - *tau;
152 z__[j4] = z__[j4 - 1] * temp;
159 i__1 = ( *n0 - 3 ) << 2;
160 for (j4 = *i0 << 2; j4 <= i__1; j4 += 4) {
161 z__[j4 - 3] = d__ + z__[j4];
162 temp = z__[j4 + 2] / z__[j4 - 3];
163 d__ = d__ * temp - *tau;
165 z__[j4 - 1] = z__[j4] * temp;
177 j4 = ( ( *n0 - 2 ) << 2) - *pp;
178 j4p2 = j4 + (*pp << 1) - 1;
179 z__[j4 - 2] = *dnm2 + z__[j4p2];
180 z__[j4] = z__[j4p2 + 2] * (z__[j4p2] / z__[j4 - 2]);
181 *dnm1 = z__[j4p2 + 2] * (*dnm2 / z__[j4 - 2]) - *tau;
186 j4p2 = j4 + (*pp << 1) - 1;
187 z__[j4 - 2] = *dnm1 + z__[j4p2];
188 z__[j4] = z__[j4p2 + 2] * (z__[j4p2] / z__[j4 - 2]);
189 *dn = z__[j4p2 + 2] * (*dnm1 / z__[j4 - 2]) - *tau;
197 i__1 = ( *n0 - 3 ) << 2;
198 for (j4 = *i0 << 2; j4 <= i__1; j4 += 4) {
199 z__[j4 - 2] = d__ + z__[j4 - 1];
203 z__[j4] = z__[j4 + 1] * (z__[j4 - 1] / z__[j4 - 2]);
204 d__ = z__[j4 + 1] * (d__ / z__[j4 - 2]) - *tau;
208 d__1 = emin, d__2 = z__[j4];
213 i__1 = ( *n0 - 3 ) << 2;
214 for (j4 = *i0 << 2; j4 <= i__1; j4 += 4) {
215 z__[j4 - 3] = d__ + z__[j4];
219 z__[j4 - 1] = z__[j4 + 2] * (z__[j4] / z__[j4 - 3]);
220 d__ = z__[j4 + 2] * (d__ / z__[j4 - 3]) - *tau;
224 d__1 = emin, d__2 = z__[j4 - 1];
234 j4 = ( ( *n0 - 2 ) << 2) - *pp;
235 j4p2 = j4 + (*pp << 1) - 1;
236 z__[j4 - 2] = *dnm2 + z__[j4p2];
240 z__[j4] = z__[j4p2 + 2] * (z__[j4p2] / z__[j4 - 2]);
241 *dnm1 = z__[j4p2 + 2] * (*dnm2 / z__[j4 - 2]) - *tau;
247 j4p2 = j4 + (*pp << 1) - 1;
248 z__[j4 - 2] = *dnm1 + z__[j4p2];
252 z__[j4] = z__[j4p2 + 2] * (z__[j4p2] / z__[j4 - 2]);
253 *dn = z__[j4p2 + 2] * (*dnm1 / z__[j4 - 2]) - *tau;
260 z__[(*n0 << 2) - *pp] = emin;
int integer
Definition: template_blas_common.h:40
#define minMACRO(a, b)
Definition: template_blas_common.h:46
int template_lapack_lasq5(integer *i0, integer *n0, Treal *z__, integer *pp, Treal *tau, Treal *dmin__, Treal *dmin1, Treal *dmin2, Treal *dn, Treal *dnm1, Treal *dnm2, logical *ieee)
Definition: template_lapack_lasq5.h:41
bool logical
Definition: template_blas_common.h:41