Searched refs:wa2 (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/tao/leastsquares/impls/pounders/ |
| H A D | gqt.c | 214 …f, PetscReal *x, PetscInt *retinfo, PetscInt *retits, PetscReal *z, PetscReal *wa1, PetscReal *wa2) in gqt() argument
246 PetscCallBLAS("BLASasum", wa2[j] = BLASasum_(&blasn, &a[0 + lda * j], &blas1)); in gqt()
248 anorm = PetscMax(anorm, wa2[j]); in gqt()
250 for (j = 0; j < n; j++) wa2[j] = wa2[j] - PetscAbs(wa1[j]); in gqt()
259 parl = PetscMax(parl, wa1[j] + wa2[j]); in gqt()
260 paru = PetscMax(paru, -wa1[j] + wa2[j]); in gqt()
298 PetscCallBLAS("BLAScopy", BLAScopy_(&blasn, b, &blas1, wa2, &blas1)); in gqt()
299 …PetscCallBLAS("LAPACKtrtrs", LAPACKtrtrs_("U", "T", "N", &blasn, &blas1, a, &blaslda, wa2, &blasn,… in gqt()
301 PetscCallBLAS("BLASnrm2", rxnorm = BLASnrm2_(&blasn, wa2, &blas1)); in gqt()
302 …PetscCallBLAS("LAPACKtrtrs", LAPACKtrtrs_("U", "N", "N", &blasn, &blas1, a, &blaslda, wa2, &blasn,… in gqt()
[all …]
|