| /petsc/src/mat/impls/submat/ |
| H A D | submat.c | 5 Vec lwork, rwork; /* work vectors inside the scatters */ member
40 PetscCall(VecZeroEntries(Na->lwork)); in MatDiagonalScale_SubMatrix()
41 PetscCall(VecScatterBegin(Na->lrestrict, left, Na->lwork, INSERT_VALUES, SCATTER_REVERSE)); in MatDiagonalScale_SubMatrix()
42 PetscCall(VecScatterEnd(Na->lrestrict, left, Na->lwork, INSERT_VALUES, SCATTER_REVERSE)); in MatDiagonalScale_SubMatrix()
44 PetscCall(MatDiagonalScale(Na->A, left ? Na->lwork : NULL, right ? Na->rwork : NULL)); in MatDiagonalScale_SubMatrix()
67 PetscCall(MatMult(Na->A, Na->rwork, Na->lwork)); in MatMult_SubMatrix()
68 PetscCall(VecScatterBegin(Na->lrestrict, Na->lwork, y, INSERT_VALUES, SCATTER_FORWARD)); in MatMult_SubMatrix()
69 PetscCall(VecScatterEnd(Na->lrestrict, Na->lwork, y, INSERT_VALUES, SCATTER_FORWARD)); in MatMult_SubMatrix()
82 PetscCall(MatMultAdd(Na->A, Na->rwork, Na->rwork, Na->lwork)); in MatMultAdd_SubMatrix()
84 PetscCall(VecZeroEntries(Na->lwork)); in MatMultAdd_SubMatrix()
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | gmreig.c | 8 PetscBLASInt bn, bN, lwork, idummy, lierr; in KSPComputeExtremeSingularValues_GMRES() local
15 PetscCall(PetscBLASIntCast(5 * N, &lwork)); in KSPComputeExtremeSingularValues_GMRES()
30 …_("N", "N", &bn, &bn, R, &bN, realpart, &sdummy, &idummy, &sdummy, &idummy, work, &lwork, &lierr)); in KSPComputeExtremeSingularValues_GMRES()
32 …", &bn, &bn, R, &bN, realpart, &sdummy, &idummy, &sdummy, &idummy, work, &lwork, realpart + N, &li… in KSPComputeExtremeSingularValues_GMRES()
47 PetscBLASInt bn, bN, lwork, idummy, lierr = -1; in KSPComputeEigenvalues_GMRES() local
54 PetscCall(PetscBLASIntCast(5 * N, &lwork)); in KSPComputeEigenvalues_GMRES()
66 …, "N", &bn, R, &bN, realpart, imagpart, &sdummy, &idummy, &sdummy, &idummy, work, &lwork, &lierr)); in KSPComputeEigenvalues_GMRES()
81 PetscBLASInt bn, bN, lwork, idummy, lierr = -1; in KSPComputeEigenvalues_GMRES()
86 PetscCall(PetscBLASIntCast(5 * N, &lwork)); in KSPComputeEigenvalues_GMRES()
98 …_("N", "N", &bn, R, &bN, eigs, &sdummy, &idummy, &sdummy, &idummy, work, &lwork, gmres->Dsvd, &lie… in KSPComputeEigenvalues_GMRES()
[all …]
|
| /petsc/src/mat/impls/aij/seq/superlu/ |
| H A D | superlu.c | 42 PetscInt lwork; member
79 PetscCall(PetscViewerASCIIPrintf(viewer, " lwork: %" PetscInt_FMT "\n", lu->lwork)); in MatView_Info_SuperLU()
101 if (lu->lwork == -1) PetscFunctionReturn(PETSC_SUCCESS); in MatSolve_SuperLU_Private()
140 …perm_r, lu->etree, lu->equed, lu->R, lu->C, &lu->L, &lu->U, lu->work, lu->lwork, &lu->B, &lu->X, &… in MatSolve_SuperLU_Private()
142 …perm_r, lu->etree, lu->equed, lu->R, lu->C, &lu->L, &lu->U, lu->work, lu->lwork, &lu->B, &lu->X, &… in MatSolve_SuperLU_Private()
146 …perm_r, lu->etree, lu->equed, lu->R, lu->C, &lu->L, &lu->U, lu->work, lu->lwork, &lu->B, &lu->X, &… in MatSolve_SuperLU_Private()
148 …perm_r, lu->etree, lu->equed, lu->R, lu->C, &lu->L, &lu->U, lu->work, lu->lwork, &lu->B, &lu->X, &… in MatSolve_SuperLU_Private()
154 …perm_r, lu->etree, lu->equed, lu->R, lu->C, &lu->L, &lu->U, lu->work, lu->lwork, &lu->B, &lu->X, &… in MatSolve_SuperLU_Private()
156 …perm_r, lu->etree, lu->equed, lu->R, lu->C, &lu->L, &lu->U, lu->work, lu->lwork, &lu->B, &lu->X, &… in MatSolve_SuperLU_Private()
160 …perm_r, lu->etree, lu->equed, lu->R, lu->C, &lu->L, &lu->U, lu->work, lu->lwork, &lu->B, &lu->X, &… in MatSolve_SuperLU_Private()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex116.c | 28 PetscBLASInt *iwork, *ifail, lwork, lierr, bn; in main() local
65 PetscCall(PetscBLASIntCast(8 * n, &lwork)); in main()
68 PetscCall(PetscMalloc1(lwork, &work)); in main()
73 LAPACKsyev_("V", "U", &bn, arrayA, &bn, evals, work, &lwork, &lierr); in main()
90 … &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_array, &in, work, &lwork, iwork, ifail, &li… in main()
137 PetscCall(PetscBLASIntCast(5 * minMN + maxMN, &lwork)); in main()
138 …PetscCall(PetscMalloc4(m * minMN, &arrayU, m * minMN, &arrayVT, m * minMN, &arrayErr, lwork, &work… in main()
153 …vd_("S", "S", &im, &in, arrayA, &im, evals, arrayU, &minMN, arrayVT, &minMN, work, &lwork, &lierr); in main()
|
| H A D | ex120.c | 20 PetscBLASInt nn, nevs, il, iu, *iwork, *ifail, lwork, lierr, bn, one = 1; in main() local
133 PetscCall(PetscBLASIntCast(2 * n, &lwork)); in main()
136 PetscCall(PetscMalloc1(lwork, &work)); in main()
142 LAPACKsyev_("V", "U", &bn, arrayA, &bn, evals, work, &lwork, rwork, &lierr); in main()
163 … &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_array, &nn, work, &lwork, rwork, iwork, ifa… in main()
172 LAPACKsygv_(&one, "V", "U", &bn, arrayA, &bn, arrayB, &bn, evals, work, &lwork, rwork, &lierr); in main()
192 … &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_array, &nn, work, &lwork, rwork, iwork, ifa… in main()
|
| /petsc/src/ksp/ksp/interface/ |
| H A D | eige.c | 158 PetscBLASInt idummy, lwork; in KSPComputeEigenvaluesExplicitly() local
162 PetscCall(PetscBLASIntCast(5 * n, &lwork)); in KSPComputeEigenvaluesExplicitly()
172 …", &bn, array, &bn, realpart, imagpart, &sdummy, &idummy, &sdummy, &idummy, work, &lwork, &lierr)); in KSPComputeEigenvaluesExplicitly()
192 PetscBLASInt idummy, lwork; in KSPComputeEigenvaluesExplicitly() local
196 PetscCall(PetscBLASIntCast(5 * n, &lwork)); in KSPComputeEigenvaluesExplicitly()
206 …"N", "N", &nb, array, &nb, eigs, &sdummy, &idummy, &sdummy, &idummy, work, &lwork, rwork, &lierr)); in KSPComputeEigenvaluesExplicitly()
|
| /petsc/src/ksp/ksp/impls/gmres/agmres/ |
| H A D | agmresdeflation.c | 111 PetscBLASInt lwork, liwork; in KSPAGMRESSchurForm() local
121 PetscCall(PetscBLASIntCast(PetscMax(8 * N + 16, 4 * neig * (N - neig)), &lwork)); in KSPAGMRESSchurForm()
127 …"I", &KspSize, &ilo, &KspSize, A, &ldA, B, &ldB, wr, wi, beta, Q, &N, Z, &N, work, &lwork, &info)); in KSPAGMRESSchurForm()
130 …, NULL, &KspSize, A, &ldA, B, &ldB, &sdim, wr, wi, beta, Q, &N, Z, &N, work, &lwork, NULL, &info)); in KSPAGMRESSchurForm()
158 …ldA, B, &ldB, wr, wi, beta, Q, &N, Z, &N, &r, NULL, NULL, &Dif[0], work, &lwork, iwork, &liwork, &… in KSPAGMRESSchurForm()
|
| H A D | agmres.c | 31 PetscInt lwork = PetscMax(8 * N + 16, 4 * neig * (N - neig)); in KSPSetUp_AGMRES() local
46 PetscCall(PetscMalloc3(N + 1, &agmres->tau, lwork, &agmres->work, N + 1, &agmres->nrs)); in KSPSetUp_AGMRES()
343 PetscBLASInt lwork; in KSPAGMRESBuildSoln() local
349 PetscCall(PetscBLASIntCast(4 * (KspSize + 1), &lwork)); in KSPAGMRESBuildSoln()
358 …", LAPACKgeqrf_(&lC, &KspSize, agmres->hh_origin, &ldH, agmres->tau, agmres->work, &lwork, &info)); in KSPAGMRESBuildSoln()
364 …hs, &KspSize, agmres->hh_origin, &ldH, agmres->tau, agmres->nrs, &N, agmres->work, &lwork, &info)); in KSPAGMRESBuildSoln()
|
| /petsc/src/mat/graphops/order/ |
| H A D | spectral.c | 133 PetscBLASInt bn, bN, lwork = 0, lierr, idummy; in MatGetOrdering_Spectral() local
141 PetscCall(PetscBLASIntCast(5 * n, &lwork)); in MatGetOrdering_Spectral()
143 PetscCall(PetscMalloc4(n, &realpart, n, &imagpart, n * n, &eigvec, lwork, &work)); in MatGetOrdering_Spectral()
145 …_("N", "V", &bn, a, &bN, realpart, imagpart, &sdummy, &idummy, eigvec, &bN, work, &lwork, &lierr)); in MatGetOrdering_Spectral()
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 657 PetscBLASInt lwork; /* size of the working vector */ in KSPDGMRESComputeSchurForm_DGMRES() local
674 PetscCall(PetscBLASIntCast(5 * N, &lwork)); in KSPDGMRESComputeSchurForm_DGMRES()
682 PetscCall(PetscMalloc1(lwork, &work)); in KSPDGMRESComputeSchurForm_DGMRES()
722 …ACKhseqr", LAPACKhseqr_("S", "I", &bn, &ilo, &ihi, A, &ldA, wr, wi, Q, &ldQ, work, &lwork, &info)); in KSPDGMRESComputeSchurForm_DGMRES()
750 lwork = PetscMax(1, 4 * NbrEig * (bn - NbrEig)); in KSPDGMRESComputeSchurForm_DGMRES()
752 PetscCall(PetscMalloc1(lwork, &work)); in KSPDGMRESComputeSchurForm_DGMRES()
758 …elect, &bn, A, &ldA, Q, &ldQ, wr, wi, &NbrEig, &CondEig, &CondSub, work, &lwork, iwork, &liwork, &… in KSPDGMRESComputeSchurForm_DGMRES()
845 PetscBLASInt lwork; /* size of the working vector */ in KSPDGMRESImproveEig_DGMRES() local
882 lwork = 8 * N + 20; /* sizeof the working space */ in KSPDGMRESImproveEig_DGMRES()
890 PetscCall(PetscMalloc1(lwork, &work)); in KSPDGMRESImproveEig_DGMRES()
[all …]
|
| /petsc/src/ksp/pc/impls/svd/ |
| H A D | svd.c | 42 PetscBLASInt nb, lwork; in PCSetUp_SVD() local
69 lwork = 5 * nb; in PCSetUp_SVD()
70 PetscCall(PetscMalloc1(lwork, &work)); in PCSetUp_SVD()
79 …("LAPACKgesvd", LAPACKgesvd_("A", "A", &nb, &nb, a, &nb, d, u, &nb, v, &nb, work, &lwork, &lierr)); in PCSetUp_SVD()
90 …gesvd", LAPACKgesvd_("A", "A", &nb, &nb, a, &nb, dd, u, &nb, v, &nb, work, &lwork, rwork, &lierr)); in PCSetUp_SVD()
|
| /petsc/include/petsc/private/ |
| H A D | cupmsolverinterface.hpp | 101 …lasInt_t /* lda */, cupmScalar_t * /* B */, cupmBlasInt_t /* ldb */, cupmBlasInt_t *lwork) noexcept in cupmSolverXpotrs_bufferSize()
103 *lwork = 0; in cupmSolverXpotrs_bufferSize()
144 …_t * /* devIpiv */, cupmScalar_t * /* B */, cupmBlasInt_t /* ldb */, cupmBlasInt_t *lwork) noexcept in cupmSolverXgetrs_bufferSize()
146 *lwork = 0; in cupmSolverXgetrs_bufferSize()
|
| /petsc/src/ksp/ksp/impls/bcgsl/ |
| H A D | bcgslimpl.h | 28 PetscBLASInt lwork; member
|
| H A D | bcgsl.c | 170 …, &ldMZ, bcgsl->s, bcgsl->u, &bell, bcgsl->v, &bell, bcgsl->work, &bcgsl->lwork, bcgsl->realwork, … in KSPSolve_BCGSL()
172 …1 + ldMZ], &ldMZ, bcgsl->s, bcgsl->u, &bell, bcgsl->v, &bell, bcgsl->work, &bcgsl->lwork, &bierr)); in KSPSolve_BCGSL()
501 PetscCall(PetscBLASIntCast(5 * ell, &bcgsl->lwork)); in KSPSetUp_BCGSL()
502 …PetscCall(PetscMalloc5(bcgsl->lwork, &bcgsl->work, ell, &bcgsl->s, ell * ell, &bcgsl->u, ell * ell… in KSPSetUp_BCGSL()
|
| /petsc/src/ksp/ksp/guess/impls/pod/ |
| H A D | pod.c | 35 PetscBLASInt lwork; member
83 pod->lwork = -1; in KSPGuessSetUp_POD()
86 …mmy, &idummy, &rdummy, &idummy, pod->eigs, pod->eigv, &bN, &sdummy, &pod->lwork, pod->iwork, pod->… in KSPGuessSetUp_POD()
88 …mmy, &idummy, &rdummy, &idummy, pod->eigs, pod->eigv, &bN, &sdummy, &pod->lwork, pod->rwork, pod->… in KSPGuessSetUp_POD()
91 PetscCall(PetscBLASIntCast((PetscInt)PetscRealPart(sdummy), &pod->lwork)); in KSPGuessSetUp_POD()
92 PetscCall(PetscMalloc1(pod->lwork + PetscMax(bN * bN, 6 * bN), &pod->swork)); in KSPGuessSetUp_POD()
313 …y, &reps, &idummy, pod->eigs, pod->eigv, &bN, pod->swork + bN * bN, &pod->lwork, pod->iwork, pod->… in KSPGuessUpdate_POD()
315 …y, &reps, &idummy, pod->eigs, pod->eigv, &bN, pod->swork + bN * bN, &pod->lwork, pod->rwork, pod->… in KSPGuessUpdate_POD()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexmetric.c | 1079 PetscBLASInt lwork; in DMPlexMetricModify_Private() local
1081 PetscCall(PetscBLASIntCast(5 * dim, &lwork)); in DMPlexMetricModify_Private()
1093 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, Mpos, &nb, eigs, work, &lwork, rwork, &lier… in DMPlexMetricModify_Private()
1097 … PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, Mpos, &nb, eigs, work, &lwork, &lierr)); in DMPlexMetricModify_Private()
1486 PetscBLASInt lwork; in DMPlexMetricIntersection_Private() local
1488 PetscCall(PetscBLASIntCast(5 * dim, &lwork)); in DMPlexMetricIntersection_Private()
1501 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, rwork, &li… in DMPlexMetricIntersection_Private()
1505 … PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, &lierr)); in DMPlexMetricIntersection_Private()
1542 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, rwork, &li… in DMPlexMetricIntersection_Private()
1546 … PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, &lierr)); in DMPlexMetricIntersection_Private()
|
| /petsc/src/tao/pde_constrained/tutorials/ |
| H A D | elliptic.c | 58 Vec lwork; member
194 PetscCall(MatMult(user->L, user->uwork, user->lwork)); in FormFunction()
195 PetscCall(VecDot(user->lwork, user->lwork, &d2)); in FormFunction()
215 PetscCall(MatMult(user->L, user->uwork, user->lwork)); in FormGradient()
216 PetscCall(MatMultTranspose(user->L, user->lwork, user->uwork)); in FormGradient()
235 PetscCall(MatMult(user->L, user->uwork, user->lwork)); in FormFunctionGradient()
236 PetscCall(VecDot(user->lwork, user->lwork, &d2)); in FormFunctionGradient()
237 PetscCall(MatMultTranspose(user->L, user->lwork, user->uwork)); in FormFunctionGradient()
946 PetscCall(VecCreate(PETSC_COMM_WORLD, &user->lwork)); in EllipticInitialize()
947 PetscCall(VecSetSizes(user->lwork, PETSC_DECIDE, m + user->mx * user->mx * user->mx)); in EllipticInitialize()
[all …]
|
| H A D | parabolic.c | 54 Vec lwork; member
243 PetscCall(MatMult(user->L, user->uwork, user->lwork)); in FormFunction()
244 PetscCall(VecDot(user->lwork, user->lwork, &d2)); in FormFunction()
279 PetscCall(MatMult(user->L, user->uwork, user->lwork)); in FormGradient()
280 PetscCall(MatMult(user->LT, user->lwork, user->uwork)); in FormGradient()
312 PetscCall(MatMult(user->L, user->uwork, user->lwork)); in FormFunctionGradient()
313 PetscCall(VecDot(user->lwork, user->lwork, &d2)); in FormFunctionGradient()
314 PetscCall(MatMult(user->LT, user->lwork, user->uwork)); in FormFunctionGradient()
896 PetscCall(VecCreate(PETSC_COMM_WORLD, &user->lwork)); in ParabolicInitialize()
897 PetscCall(VecSetSizes(user->lwork, PETSC_DECIDE, m + user->mx * user->mx * user->mx)); in ParabolicInitialize()
[all …]
|
| H A D | hyperbolic.c | 51 Vec lwork; member
230 PetscCall(MatMult(user->L, user->uwork, user->lwork)); in FormFunction()
231 PetscCall(VecDot(user->y, user->lwork, &d2)); in FormFunction()
258 PetscCall(MatMult(user->L, user->vwork, user->lwork)); in FormGradient()
259 PetscCall(VecAXPY(user->ywork, 0.5 * user->alpha, user->lwork)); in FormGradient()
285 PetscCall(MatMult(user->L, user->vwork, user->lwork)); in FormFunctionGradient()
286 PetscCall(VecAXPY(user->ywork, 0.5 * user->alpha, user->lwork)); in FormFunctionGradient()
288 PetscCall(VecDot(user->y, user->lwork, &d2)); in FormFunctionGradient()
998 PetscCall(VecCreate(PETSC_COMM_WORLD, &user->lwork)); in HyperbolicInitialize()
999 PetscCall(VecSetType(user->lwork, VECMPI)); in HyperbolicInitialize()
[all …]
|
| /petsc/src/ksp/ksp/tests/output/ |
| H A D | ex63_1.out | 39 lwork: 0
|
| /petsc/src/snes/impls/ngmres/ |
| H A D | snesngmres.h | 54 PetscBLASInt lwork; /* the size of the work vector */ member
|
| H A D | snesngmres.c | 67 PetscCall(PetscBLASIntCast(12 * msize, &ngmres->lwork)); in SNESSetUp_NGMRES()
69 PetscCall(PetscMalloc1(ngmres->lwork, &ngmres->rwork)); in SNESSetUp_NGMRES()
71 PetscCall(PetscMalloc1(ngmres->lwork, &ngmres->work)); in SNESSetUp_NGMRES()
|
| /petsc/src/snes/impls/composite/ |
| H A D | snescomposite.c | 40 PetscBLASInt lwork; /* the size of the work vector */ member
263 …, jac->beta, &jac->lda, jac->s, &jac->rcond, &jac->rank, jac->work, &jac->lwork, jac->rwork, &jac-… in SNESCompositeApply_AdditiveOptimal()
265 …->lda, jac->beta, &jac->lda, jac->s, &jac->rcond, &jac->rank, jac->work, &jac->lwork, &jac->info)); in SNESCompositeApply_AdditiveOptimal()
363 jac->lwork = 12 * jac->n; in SNESSetUp_Composite()
365 PetscCall(PetscMalloc1(jac->lwork, &jac->rwork)); in SNESSetUp_Composite()
367 PetscCall(PetscMalloc1(jac->lwork, &jac->work)); in SNESSetUp_Composite()
|
| /petsc/src/snes/interface/ |
| H A D | snesut.c | 429 PetscBLASInt nb = 0, lwork; in SNESMonitorJacUpdateSpectrum() local
447 lwork = 3 * nb; in SNESMonitorJacUpdateSpectrum()
450 PetscCall(PetscMalloc1(lwork, &work)); in SNESMonitorJacUpdateSpectrum()
457 …geev", LAPACKgeev_("N", "N", &nb, a, &nb, eigr, eigi, NULL, &nb, NULL, &nb, work, &lwork, &lierr)); in SNESMonitorJacUpdateSpectrum()
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex13.c | 86 PetscBLASInt n, lwork, lierr; in test() local
89 PetscCall(PetscBLASIntCast(5 * Nbpt, &lwork)); in test()
94 …ACKgesvd", LAPACKgesvd_("N", "N", &n, &n, Mcopy, &n, S, NULL, &n, NULL, &n, work, &lwork, &lierr)); in test()
|