Searched refs:qqc (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/ksp/pc/impls/gamg/ |
| H A D | agg.c | 559 PetscScalar *qqc, *qqr, *TAU, *WORK; in formProl0() local
573 PetscCall(PetscMalloc5(Mdata * N, &qqc, M * N, &qqr, N, &TAU, LWORK, &WORK, M, &fids)); in formProl0()
594 qqc[jj * Mdata + aggID * bs + ii] = d; in formProl0()
604 for (jj = 0; jj < N; jj++, kk++) qqc[jj * Mdata + ii] = .0; in formProl0()
609 PetscCallBLAS("LAPACKgeqrf", LAPACKgeqrf_(&Mdata, &N, qqc, &LDA, TAU, WORK, &LWORK, &INFO)); in formProl0()
619 if (ii <= jj) data[jj * out_data_stride + ii] = PetscRealPart(qqc[jj * Mdata + ii]); in formProl0()
626 … PetscCallBLAS("LAPACKorgqr", LAPACKorgqr_(&Mdata, &N, &N, qqc, &LDA, TAU, WORK, &LWORK, &INFO)); in formProl0()
630 for (jj = 0; jj < N; jj++) qqr[N * ii + jj] = qqc[jj * Mdata + ii]; in formProl0()
636 PetscCall(PetscFree5(qqc, qqr, TAU, WORK, fids)); in formProl0()
|