impls/pounders/gqt.c

8   PetscReal    e, temp, w, wm, ynorm, znorm, s, sm;
26         temp = PetscMin(0.01, PetscAbs(r[i + ldr * i])) / PetscAbs(e - z[i]);
27         PetscCallBLAS("BLASscal", BLASscal_(&blasn, &temp, z, &blas1));
28         e = temp * e;
49         temp = wm - w;
51         if (i < n - 1) PetscCallBLAS("BLASaxpy", BLASaxpy_(&blasnmi, &temp, &r[i + ldr * (i + 1)], &blasldr, &z[i + 1], &blas1));
62         temp = PetscMin(0.01, PetscAbs(r[j + ldr * j] / z[j]));
63         PetscCallBLAS("BLASscal", BLASscal_(&blasn, &temp, z, &blas1));
64         ynorm *= temp;
71       temp = -z[j];
73       PetscCallBLAS("BLASaxpy", BLASaxpy_(&blasj, &temp, &r[0 + ldr * j], &blas1, z, &blas1));
220   PetscReal    alpha, anorm, bnorm, parc, parf, parl, pars, par = *retpar, paru, prod, rxnorm, rznorm = 0.0, temp, xnorm;
324         temp  = (delta - xnorm) * ((delta + xnorm) / delta);
325         alpha = temp / (PetscAbs(prod) + PetscSqrtScalar(prod * prod + temp / delta));
345         temp = 1.0 / xnorm;
346         PetscCallBLAS("BLASscal", BLASscal_(&blasn, &temp, wa2, &blas1));
349         PetscCallBLAS("BLASnrm2", temp = BLASnrm2_(&blasn, wa2, &blas1));
350         parc = (xnorm - delta) / (delta * temp * temp);
375         PetscCallBLAS("BLASnrm2", temp = BLASnrm2_(&iblas, &a[0 + (indef - 1) * lda], &blas1));
377         a[indef - 1 + (indef - 1) * lda] -= temp * temp;
384       PetscCallBLAS("BLASnrm2", temp = BLASnrm2_(&iblas, wa2, &blas1));
385       parc = -a[indef - 1 + (indef - 1) * lda] / (temp * temp);