impls/pounders/gqt.c

115 c       subroutine gqt(n,a,lda,b,delta,rtol,atol,itmax,
124 c       a is a double precision array of dimension (lda,n).
129 c       lda is an integer variable.
130 c         On entry lda is the leading dimension of the array a.
131 c         On exit lda is unchanged.
214 PetscErrorCode gqt(PetscInt n, PetscReal *a, PetscInt lda, PetscReal *b, PetscReal delta, PetscReal rtol, PetscReal atol, PetscInt itmax, PetscReal *retpar, PetscReal *retf, PetscReal *x, PetscInt *retinfo, PetscInt *retits, PetscReal *z, PetscReal *wa1, PetscReal *wa2)
224   PetscCall(PetscBLASIntCast(lda, &blaslda));
225   PetscCall(PetscBLASIntCast(lda + 1, &blasldap1));
239     PetscCallBLAS("BLAScopy", BLAScopy_(&iblas, &a[j + lda * (j + 1)], &blaslda, &a[j + 1 + lda * j], &blas1));
246     PetscCallBLAS("BLASasum", wa2[j] = BLASasum_(&blasn, &a[0 + lda * j], &blas1));
286       PetscCallBLAS("BLAScopy", BLAScopy_(&iblas, &a[j + 1 + j * lda], &blas1, &a[j + (j + 1) * lda], &blaslda));
288     for (j = 0; j < n; j++) a[j + j * lda] = wa1[j] + par;
314       PetscCall(estsv(n, a, lda, &rznorm, z));
367         PetscCallBLAS("BLAScopy", BLAScopy_(&iblas, &a[indef - 1 + 0 * lda], &blaslda, &a[0 + (indef - 1) * lda], &blas1));
368         a[indef - 1 + (indef - 1) * lda] = wa1[indef - 1] + par;
371         PetscCallBLAS("BLAScopy", BLAScopy_(&iblas, &a[0 + (indef - 1) * lda], &blas1, wa2, &blas1));
374         PetscCallBLAS("BLAScopy", BLAScopy_(&iblas, wa2, &blas1, &a[0 + (indef - 1) * lda], &blas1));
375         PetscCallBLAS("BLASnrm2", temp = BLASnrm2_(&iblas, &a[0 + (indef - 1) * lda], &blas1));
377         a[indef - 1 + (indef - 1) * lda] -= temp * temp;
385       parc = -a[indef - 1 + (indef - 1) * lda] / (temp * temp);
419         PetscCallBLAS("BLAScopy", BLAScopy_(&iblas, &a[j + 1 + j * lda], &blas1, &a[j + (j + 1) * lda], &blaslda));
421       PetscCall(PetscBLASIntCast(lda + 1, &iblas));