35 #ifndef TEMPLATE_LAPACK_LASQ5_HEADER 36 #define TEMPLATE_LAPACK_LASQ5_HEADER 40 integer *pp, Treal *tau, Treal *dmin__, Treal *dmin1,
41 Treal *dmin2, Treal *dn, Treal *dnm1, Treal *dnm2,
129 if (*n0 - *i0 - 1 <= 0) {
133 j4 = (*i0 << 2) + *pp - 3;
135 d__ = z__[j4] - *tau;
144 i__1 = ( *n0 - 3 ) << 2;
145 for (j4 = *i0 << 2; j4 <= i__1; j4 += 4) {
146 z__[j4 - 2] = d__ + z__[j4 - 1];
147 temp = z__[j4 + 1] / z__[j4 - 2];
148 d__ = d__ * temp - *tau;
150 z__[j4] = z__[j4 - 1] * temp;
157 i__1 = ( *n0 - 3 ) << 2;
158 for (j4 = *i0 << 2; j4 <= i__1; j4 += 4) {
159 z__[j4 - 3] = d__ + z__[j4];
160 temp = z__[j4 + 2] / z__[j4 - 3];
161 d__ = d__ * temp - *tau;
163 z__[j4 - 1] = z__[j4] * temp;
175 j4 = ( ( *n0 - 2 ) << 2) - *pp;
176 j4p2 = j4 + (*pp << 1) - 1;
177 z__[j4 - 2] = *dnm2 + z__[j4p2];
178 z__[j4] = z__[j4p2 + 2] * (z__[j4p2] / z__[j4 - 2]);
179 *dnm1 = z__[j4p2 + 2] * (*dnm2 / z__[j4 - 2]) - *tau;
184 j4p2 = j4 + (*pp << 1) - 1;
185 z__[j4 - 2] = *dnm1 + z__[j4p2];
186 z__[j4] = z__[j4p2 + 2] * (z__[j4p2] / z__[j4 - 2]);
187 *dn = z__[j4p2 + 2] * (*dnm1 / z__[j4 - 2]) - *tau;
195 i__1 = ( *n0 - 3 ) << 2;
196 for (j4 = *i0 << 2; j4 <= i__1; j4 += 4) {
197 z__[j4 - 2] = d__ + z__[j4 - 1];
201 z__[j4] = z__[j4 + 1] * (z__[j4 - 1] / z__[j4 - 2]);
202 d__ = z__[j4 + 1] * (d__ / z__[j4 - 2]) - *tau;
206 d__1 = emin, d__2 = z__[j4];
211 i__1 = ( *n0 - 3 ) << 2;
212 for (j4 = *i0 << 2; j4 <= i__1; j4 += 4) {
213 z__[j4 - 3] = d__ + z__[j4];
217 z__[j4 - 1] = z__[j4 + 2] * (z__[j4] / z__[j4 - 3]);
218 d__ = z__[j4 + 2] * (d__ / z__[j4 - 3]) - *tau;
222 d__1 = emin, d__2 = z__[j4 - 1];
232 j4 = ( ( *n0 - 2 ) << 2) - *pp;
233 j4p2 = j4 + (*pp << 1) - 1;
234 z__[j4 - 2] = *dnm2 + z__[j4p2];
238 z__[j4] = z__[j4p2 + 2] * (z__[j4p2] / z__[j4 - 2]);
239 *dnm1 = z__[j4p2 + 2] * (*dnm2 / z__[j4 - 2]) - *tau;
245 j4p2 = j4 + (*pp << 1) - 1;
246 z__[j4 - 2] = *dnm1 + z__[j4p2];
250 z__[j4] = z__[j4p2 + 2] * (z__[j4p2] / z__[j4 - 2]);
251 *dn = z__[j4p2 + 2] * (*dnm1 / z__[j4 - 2]) - *tau;
258 z__[(*n0 << 2) - *pp] = emin;
int integer
Definition: template_blas_common.h:38
#define minMACRO(a, b)
Definition: template_blas_common.h:44
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:39
bool logical
Definition: template_blas_common.h:39