| /petsc/src/mat/impls/submat/ |
| H A D | submat.c | 5 Vec lwork, rwork; /* work vectors inside the scatters */ member
35 PetscCall(VecZeroEntries(Na->rwork)); in MatDiagonalScale_SubMatrix()
36 PetscCall(VecScatterBegin(Na->rprolong, right, Na->rwork, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_SubMatrix()
37 PetscCall(VecScatterEnd(Na->rprolong, right, Na->rwork, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_SubMatrix()
44 PetscCall(MatDiagonalScale(Na->A, left ? Na->lwork : NULL, right ? Na->rwork : NULL)); in MatDiagonalScale_SubMatrix()
53 PetscCall(MatGetDiagonal(Na->A, Na->rwork)); in MatGetDiagonal_SubMatrix()
54 PetscCall(VecScatterBegin(Na->rprolong, Na->rwork, d, INSERT_VALUES, SCATTER_REVERSE)); in MatGetDiagonal_SubMatrix()
55 PetscCall(VecScatterEnd(Na->rprolong, Na->rwork, d, INSERT_VALUES, SCATTER_REVERSE)); in MatGetDiagonal_SubMatrix()
64 PetscCall(VecZeroEntries(Na->rwork)); in MatMult_SubMatrix()
65 PetscCall(VecScatterBegin(Na->rprolong, x, Na->rwork, INSERT_VALUES, SCATTER_FORWARD)); in MatMult_SubMatrix()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex120.c | 16 PetscReal *evals, *rwork; in main() local
141 PetscCall(PetscMalloc1(3 * n - 2, &rwork)); in main()
142 LAPACKsyev_("V", "U", &bn, arrayA, &bn, evals, work, &lwork, rwork, &lierr); in main()
143 PetscCall(PetscFree(rwork)); in main()
155 PetscCall(PetscMalloc1(7 * n + 1, &rwork)); in main()
163 …vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_array, &nn, work, &lwork, rwork, iwork, ifail, &li… in main()
166 PetscCall(PetscFree(rwork)); in main()
170 PetscCall(PetscMalloc1(3 * n + 1, &rwork)); in main()
172 LAPACKsygv_(&one, "V", "U", &bn, arrayA, &bn, arrayB, &bn, evals, work, &lwork, rwork, &lierr); in main()
178 PetscCall(PetscFree(rwork)); in main()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/gltr/ |
| H A D | gltrimpl.h | 29 PetscReal *rwork; /* Real workspace for solver computations */ member
|
| H A D | gltr.c | 625 PetscCall(PetscFree(cg->rwork)); in KSPCGSolve_GLTR()
635 PetscCall(PetscMalloc2(10 * cg->alloced, &cg->rwork, 5 * cg->alloced, &cg->iwork)); in KSPCGSolve_GLTR()
639 t_soln = cg->rwork + 0 * t_size; /* Solution */ in KSPCGSolve_GLTR()
640 t_diag = cg->rwork + 1 * t_size; /* Diagonal of T */ in KSPCGSolve_GLTR()
641 t_offd = cg->rwork + 2 * t_size; /* Off-diagonal of T */ in KSPCGSolve_GLTR()
642 e_valu = cg->rwork + 3 * t_size; /* Eigenvalues of T */ in KSPCGSolve_GLTR()
643 e_vect = cg->rwork + 4 * t_size; /* Eigenvector of T */ in KSPCGSolve_GLTR()
644 e_rwrk = cg->rwork + 5 * t_size; /* Eigen workspace */ in KSPCGSolve_GLTR()
921 if (cg->alloced) PetscCall(PetscFree2(cg->rwork, cg->iwork)); in KSPCGDestroy_GLTR()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexmetric.c | 1091 PetscReal *rwork; in DMPlexMetricModify_Private() local
1092 PetscCall(PetscMalloc1(3 * dim, &rwork)); in DMPlexMetricModify_Private()
1093 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, Mpos, &nb, eigs, work, &lwork, rwork, &lier… in DMPlexMetricModify_Private()
1094 PetscCall(PetscFree(rwork)); in DMPlexMetricModify_Private()
1499 PetscReal *rwork; in DMPlexMetricIntersection_Private() local
1500 PetscCall(PetscMalloc1(3 * dim, &rwork)); in DMPlexMetricIntersection_Private()
1501 …CallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, rwork, &lierr)); in DMPlexMetricIntersection_Private()
1502 PetscCall(PetscFree(rwork)); in DMPlexMetricIntersection_Private()
1540 PetscReal *rwork; in DMPlexMetricIntersection_Private() local
1541 PetscCall(PetscMalloc1(3 * dim, &rwork)); in DMPlexMetricIntersection_Private()
[all …]
|
| /petsc/src/ksp/ksp/interface/ |
| H A D | eige.c | 191 PetscReal *rwork; in KSPComputeEigenvaluesExplicitly() local
198 PetscCall(PetscMalloc1(2 * n, &rwork)); in KSPComputeEigenvaluesExplicitly()
206 …"N", "N", &nb, array, &nb, eigs, &sdummy, &idummy, &sdummy, &idummy, work, &lwork, rwork, &lierr)); in KSPComputeEigenvaluesExplicitly()
211 PetscCall(PetscFree(rwork)); in KSPComputeEigenvaluesExplicitly()
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | gmreig.c | 182 PetscReal *rwork = NULL; in KSPComputeRitz_GMRES() local
189 PetscCall(PetscMalloc1(2 * n, &rwork)); in KSPComputeRitz_GMRES()
190 …ev", LAPACKgeev_("N", "V", &bn, H, &bN, wr, &sdummy, &idummy, Q, &bn, work, &lwork, rwork, &info)); in KSPComputeRitz_GMRES()
191 PetscCall(PetscFree(rwork)); in KSPComputeRitz_GMRES()
|
| /petsc/src/ksp/pc/impls/svd/ |
| H A D | svd.c | 86 PetscReal *rwork, *dd; in PCSetUp_SVD() local
87 PetscCall(PetscMalloc1(5 * nb, &rwork)); in PCSetUp_SVD()
90 …gesvd", LAPACKgesvd_("A", "A", &nb, &nb, a, &nb, dd, u, &nb, v, &nb, work, &lwork, rwork, &lierr)); in PCSetUp_SVD()
92 PetscCall(PetscFree(rwork)); in PCSetUp_SVD()
|
| /petsc/src/snes/impls/ngmres/ |
| H A D | snesngmres.h | 53 PetscReal *rwork; /* the real work vector used for complex */ member
|
| H A D | snesngmres.c | 28 PetscCall(PetscFree(ngmres->rwork)); in SNESDestroy_NGMRES()
69 PetscCall(PetscMalloc1(ngmres->lwork, &ngmres->rwork)); in SNESSetUp_NGMRES()
|
| H A D | ngmresfunc.c | 85 …es->s, &ngmres->rcond, &ngmres->rank, ngmres->work, &ngmres->lwork, ngmres->rwork, &ngmres->info)); in SNESNGMRESFormCombinedSolution_Private()
|
| /petsc/src/ts/impls/glee/ |
| H A D | glee.c | 58 PetscScalar *rwork; /* Scalar work (size of the number of steps)*/ member
440 PetscScalar *ws = glee->swork, *wr = glee->rwork; in TSEvaluateStep_GLEE()
503 PetscScalar *ws = glee->swork, *wr = glee->rwork; in TSStep_GLEE()
640 PetscCall(PetscFree2(glee->swork, glee->rwork)); in TSReset_GLEE()
756 PetscCall(PetscMalloc2(s, &glee->swork, r, &glee->rwork)); in TSSetUp_GLEE()
945 PetscScalar *wr = glee->rwork; in TSGetAuxSolution_GLEE()
962 PetscScalar *wr = glee->rwork; in TSGetTimeError_GLEE()
|
| /petsc/src/ksp/ksp/guess/impls/pod/ |
| H A D | pod.c | 33 PetscReal *rwork; member
78 PetscCall(PetscMalloc1(7 * pod->maxn, &pod->rwork)); in KSPGuessSetUp_POD()
88 …, &rdummy, &idummy, pod->eigs, pod->eigv, &bN, &sdummy, &pod->lwork, pod->rwork, pod->iwork, pod->… in KSPGuessSetUp_POD()
124 PetscCall(PetscFree(pod->rwork)); in KSPGuessDestroy_POD()
315 …dummy, pod->eigs, pod->eigv, &bN, pod->swork + bN * bN, &pod->lwork, pod->rwork, pod->iwork, pod->… in KSPGuessUpdate_POD()
|
| /petsc/src/ksp/ksp/guess/impls/fischer/ |
| H A D | fischer.c | 204 PetscReal *rwork; in KSPGuessFormGuess_Fischer_3() local
217 …PetscCall(PetscCalloc5(m * m, &corr, m, &s_values, blas_lwork, &work, 3 * m - 2, &rwork, m, &scrat… in KSPGuessFormGuess_Fischer_3()
228 …v", LAPACKheev_("V", "L", &blas_m, corr, &blas_m, s_values, work, &blas_lwork, rwork, &blas_info)); in KSPGuessFormGuess_Fischer_3()
274 PetscCall(PetscFree5(corr, s_values, work, rwork, scratch_vec)); in KSPGuessFormGuess_Fischer_3()
|
| /petsc/src/snes/impls/composite/ |
| H A D | snescomposite.c | 39 PetscReal *rwork; /* the real work vector used for complex */ member
263 …beta, &jac->lda, jac->s, &jac->rcond, &jac->rank, jac->work, &jac->lwork, jac->rwork, &jac->info)); in SNESCompositeApply_AdditiveOptimal()
365 PetscCall(PetscMalloc1(jac->lwork, &jac->rwork)); in SNESSetUp_Composite()
388 PetscCall(PetscFree(jac->rwork)); in SNESReset_Composite()
844 jac->rwork = NULL; in SNESCreate_Composite()
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 16 …RangeOrComplement(Mat A, PetscBool range, PetscInt lw, PetscScalar *work, PetscReal *rwork, Mat *B) in MatDenseOrthogonalRangeOrComplement() argument
39 if (!rwork) { in MatDenseOrthogonalRangeOrComplement()
42 sing = rwork; in MatDenseOrthogonalRangeOrComplement()
64 if (!rwork) PetscCall(PetscFree(sing)); in MatDenseOrthogonalRangeOrComplement()
87 …tcol, IS corners, Mat *Gins, Mat *GKins, PetscScalar cvals[2], PetscScalar *work, PetscReal *rwork) in PCBDDCComputeNedelecChangeEdge() argument
112 PetscCall(MatDenseOrthogonalRangeOrComplement(GEd, PETSC_FALSE, 5 * esize, work, rwork, GKins)); in PCBDDCComputeNedelecChangeEdge()
165 PetscReal *rwork; in PCBDDCNedelecSupport() local
1327 PetscCall(PetscMalloc2((5 + cum + maxsize) * maxsize, &work, maxsize, &rwork)); in PCBDDCNedelecSupport()
1481 …ecChangeEdge(lG, eedges[i], extrows[i], extcols[i], cornersis, &Gins, &GKins, cvals, work, rwork)); in PCBDDCNedelecSupport()
1640 PetscCall(PetscFree2(work, rwork)); in PCBDDCNedelecSupport()
[all …]
|
| /petsc/src/dm/dt/fv/interface/ |
| H A D | fv.c | 2063 PetscReal *rwork, *rtau; in PetscFVLeastSquaresPseudoInverseSVD_Static() local
2087 PetscCall(PetscMalloc1(rworkSize, &rwork)); in PetscFVLeastSquaresPseudoInverseSVD_Static()
2090 …Kgelss_(&M, &N, &nrhs, A, &lda, Brhs, &ldb, rtau, &rcond, &irank, tmpwork, &ldwork, rwork, &info)); in PetscFVLeastSquaresPseudoInverseSVD_Static()
2092 PetscCall(PetscFree(rwork)); in PetscFVLeastSquaresPseudoInverseSVD_Static()
|
| /petsc/src/mat/impls/dense/mpi/ |
| H A D | mpidense.c | 1076 PetscInt M = A->rmap->N, N = A->cmap->N, m, n, *rwork, rstart = A->rmap->rstart; in MatTranspose_MPIDense() local
1093 PetscCall(PetscMalloc1(m, &rwork)); in MatTranspose_MPIDense()
1094 for (i = 0; i < m; i++) rwork[i] = rstart + i; in MatTranspose_MPIDense()
1096 PetscCall(MatSetValues(B, 1, &j, m, rwork, v, INSERT_VALUES)); in MatTranspose_MPIDense()
1100 PetscCall(PetscFree(rwork)); in MatTranspose_MPIDense()
|