37 #ifndef TEMPLATE_LAPACK_LAEBZ_HEADER
38 #define TEMPLATE_LAPACK_LAEBZ_HEADER
44 const Treal *reltol,
const Treal *pivmin,
const Treal *d__,
const Treal *
45 e,
const Treal *e2,
integer *nval, Treal *ab, Treal *c__,
264 integer nab_dim1, nab_offset, ab_dim1, ab_offset, i__1, i__2, i__3, i__4,
266 Treal d__1, d__2, d__3, d__4;
268 integer itmp1, itmp2, j, kfnew, klnew, kf, ji, kl, jp, jit;
270 #define ab_ref(a_1,a_2) ab[(a_2)*ab_dim1 + a_1]
271 #define nab_ref(a_1,a_2) nab[(a_2)*nab_dim1 + a_1]
274 nab_offset = 1 + nab_dim1 * 1;
277 ab_offset = 1 + ab_dim1 * 1;
289 if (*ijob < 1 || *ijob > 3) {
303 for (ji = 1; ji <= i__1; ++ji) {
304 for (jp = 1; jp <= 2; ++jp) {
305 tmp1 = d__[1] -
ab_ref(ji, jp);
315 for (j = 2; j <= i__2; ++j) {
316 tmp1 = d__[j] - e2[j - 1] / tmp1 -
ab_ref(ji, jp);
347 for (ji = 1; ji <= i__1; ++ji) {
356 for (jit = 1; jit <= i__1; ++jit) {
360 if (kl - kf + 1 >= *nbmin && *nbmin > 0) {
365 for (ji = kf; ji <= i__2; ++ji) {
369 work[ji] = d__[1] - c__[ji];
371 if (work[ji] <= *pivmin) {
374 d__1 = work[ji], d__2 = -(*pivmin);
379 for (j = 2; j <= i__3; ++j) {
380 work[ji] = d__[j] - e2[j - 1] / work[ji] - c__[ji];
381 if (work[ji] <= *pivmin) {
384 d__1 = work[ji], d__2 = -(*pivmin);
398 for (ji = kf; ji <= i__2; ++ji) {
404 i__5 =
nab_ref(ji, 1), i__6 = iwork[ji];
411 if (iwork[ji] ==
nab_ref(ji, 2)) {
418 }
else if (iwork[ji] ==
nab_ref(ji, 1)) {
426 if (klnew <= *mmax) {
433 ab_ref(klnew, 1) = c__[ji];
453 for (ji = kf; ji <= i__2; ++ji) {
454 if (iwork[ji] <= nval[ji]) {
458 if (iwork[ji] >= nval[ji]) {
474 for (ji = kf; ji <= i__2; ++ji) {
479 tmp2 = d__[1] - tmp1;
481 if (tmp2 <= *pivmin) {
484 d__1 = tmp2, d__2 = -(*pivmin);
504 for (j = 2; j <= i__3; ++j) {
505 tmp2 = d__[j] - e2[j - 1] / tmp2 - tmp1;
506 if (tmp2 <= *pivmin) {
509 d__1 = tmp2, d__2 = -(*pivmin);
537 }
else if (itmp1 ==
nab_ref(ji, 1)) {
543 }
else if (klnew < *mmax) {
563 if (itmp1 <= nval[ji]) {
567 if (itmp1 >= nval[ji]) {
584 for (ji = kf; ji <= i__2; ++ji) {
591 d__1 =
maxMACRO(*abstol,*pivmin), d__2 = *reltol * tmp2;
612 nval[ji] = nval[kfnew];
625 for (ji = kf; ji <= i__2; ++ji) {