| /petsc/src/mat/tests/ |
| H A D | ex120.c | 20 PetscBLASInt nn, nevs, il, iu, *iwork, *ifail, lwork, lierr, bn, one = 1; in main() local
156 PetscCall(PetscMalloc1(5 * n + 1, &iwork)); in main()
163 …&vu, &il, &iu, &abstol, &nevs, evals, evecs_array, &nn, work, &lwork, rwork, iwork, ifail, &lierr); in main()
164 PetscCall(PetscFree(iwork)); in main()
185 PetscCall(PetscMalloc1(6 * n + 1, &iwork)); in main()
186 ifail = iwork + 5 * n; in main()
192 …&vu, &il, &iu, &abstol, &nevs, evals, evecs_array, &nn, work, &lwork, rwork, iwork, ifail, &lierr); in main()
194 PetscCall(PetscFree(iwork)); in main()
|
| H A D | ex118.c | 17 PetscBLASInt n, il = 1, iu = 5, *iblock, *isplit, *iwork, nevs, *ifail, cklvl = 2; in main()
37 PetscCall(PetscMalloc1(3 * n + 1, &iwork)); in main()
49 …cReal *)D, (PetscReal *)E, &nevs, &nsplit, (PetscReal *)evals, iblock, isplit, work, iwork, &info); in main()
55 …scReal *)E, &nevs, (PetscReal *)evals, iblock, isplit, evecs_array, &n, work, iwork, ifail, &info); in main()
102 PetscCall(PetscFree(iwork)); in main()
|
| H A D | ex116.c | 28 PetscBLASInt *iwork, *ifail, lwork, lierr, bn; in main() local
84 PetscCall(PetscMalloc1(6 * n + 1, &iwork)); in main()
85 ifail = iwork + 5 * n; in main()
90 …, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_array, &in, work, &lwork, iwork, ifail, &lierr); in main()
91 PetscCall(PetscFree(iwork)); in main()
|
| /petsc/src/ksp/ksp/guess/impls/pod/ |
| H A D | pod.c | 29 PetscBLASInt *iwork; /* integer work vector */ member
76 …maxn, &pod->eigs, pod->maxn * pod->maxn, &pod->eigv, 6 * pod->maxn, &pod->iwork, pod->maxn * pod->… in KSPGuessSetUp_POD()
86 …mmy, &idummy, pod->eigs, pod->eigv, &bN, &sdummy, &pod->lwork, pod->iwork, pod->iwork + 5 * bN, &l… in KSPGuessSetUp_POD()
88 …, pod->eigs, pod->eigv, &bN, &sdummy, &pod->lwork, pod->rwork, pod->iwork, pod->iwork + 5 * bN, &l… in KSPGuessSetUp_POD()
122 PetscCall(PetscFree6(pod->corr, pod->eigs, pod->eigv, pod->iwork, pod->yhay, pod->low)); in KSPGuessDestroy_POD()
190 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bNen, &bNen, pod->low, &bNen, pod->iwork, &lierr)); in KSPGuessFormGuess_POD()
192 …rs", LAPACKgetrs_(tsolve ? "T" : "N", &bNen, &ione, pod->low, &bNen, pod->iwork, pod->swork, &bNen… in KSPGuessFormGuess_POD()
313 … pod->eigs, pod->eigv, &bN, pod->swork + bN * bN, &pod->lwork, pod->iwork, pod->iwork + 5 * bN, &l… in KSPGuessUpdate_POD()
315 …pod->eigv, &bN, pod->swork + bN * bN, &pod->lwork, pod->rwork, pod->iwork, pod->iwork + 5 * bN, &l… in KSPGuessUpdate_POD()
|
| /petsc/src/mat/impls/aij/seq/bas/ |
| H A D | spbas.c | 737 static PetscErrorCode spbas_mark_row_power(PetscInt *iwork, /* marker-vector */ in spbas_mark_row_power() argument
753 if (minmrk <= j && j < maxmrk && iwork[j] < marker) { in spbas_mark_row_power()
754 …PetscCall(spbas_mark_row_power(iwork, row + in_matrix->icols[row][i], in_matrix, marker / 2, minmr… in spbas_mark_row_power()
755 iwork[j] |= marker; in spbas_mark_row_power()
762 if (minmrk <= j && j < maxmrk) iwork[j] |= 1; in spbas_mark_row_power()
782 PetscInt *iwork; in spbas_power() local
803 PetscCall(PetscCalloc1(PetscMax(ncols, nrows), &iwork)); in spbas_power()
815 PetscCall(spbas_mark_row_power(iwork, i, &in_matrix, marker, i, maxmrk)); in spbas_power()
819 for (j = i; j < maxmrk; j++) nnz += (iwork[j] != 0); in spbas_power()
828 if (iwork[j]) { in spbas_power()
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmresimpl.h | 38 PetscBLASInt *iwork; /* work space for LAPACK functions */ \
88 #define IWORK dgmres->iwork
|
| H A D | dgmres.c | 284 PetscCall(PetscFree(dgmres->iwork)); in KSPDestroy_DGMRES()
663 PetscBLASInt *iwork; in KSPDGMRESComputeSchurForm_DGMRES() local
753 PetscCall(PetscMalloc1(liwork, &iwork)); in KSPDGMRESComputeSchurForm_DGMRES()
758 … &bn, A, &ldA, Q, &ldQ, wr, wi, &NbrEig, &CondEig, &CondSub, work, &lwork, iwork, &liwork, &info)); in KSPDGMRESComputeSchurForm_DGMRES()
770 PetscCall(PetscFree(iwork)); in KSPDGMRESComputeSchurForm_DGMRES()
850 PetscBLASInt liwork, *iwork; in KSPDGMRESImproveEig_DGMRES() local
931 PetscCall(PetscMalloc1(liwork, &iwork)); in KSPDGMRESImproveEig_DGMRES()
937 …A, wr, wi, beta, Q, &N, Z, &N, &NbrEig, NULL, NULL, &Dif[0], work, &lwork, iwork, &liwork, &info)); in KSPDGMRESImproveEig_DGMRES()
975 PetscCall(PetscFree(iwork)); in KSPDGMRESImproveEig_DGMRES()
|
| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | sbaijov.c | 162 …MPIInt proc_end = 0, size, rank, tag1, tag2, *len_s, nrqr, nrqs, *id_r1, *len_r1, flag, *iwork; in MatIncreaseOverlap_MPISBAIJ_Once() local
216 PetscCall(PetscMalloc4(size, &len_s, size, &btable, size, &iwork, size + 1, &Bowners)); in MatIncreaseOverlap_MPISBAIJ_Once()
253 iwork[proc_id] = 0; in MatIncreaseOverlap_MPISBAIJ_Once()
284 *(data1 + proc_id * len + 1 + i) = len_s[proc_id] - iwork[proc_id]; in MatIncreaseOverlap_MPISBAIJ_Once()
285 iwork[proc_id] = len_s[proc_id]; in MatIncreaseOverlap_MPISBAIJ_Once()
326 for (PetscMPIInt proc_id = 0; proc_id < size; proc_id++) len_s[proc_id] = iwork[proc_id] = 0; in MatIncreaseOverlap_MPISBAIJ_Once()
393 PetscCallMPI(MPIU_Allreduce(len_s, iwork, size, MPI_INT, MPI_MAX, comm)); in MatIncreaseOverlap_MPISBAIJ_Once()
394 PetscCall(PetscMalloc1(iwork[rank] + 1, &data2)); in MatIncreaseOverlap_MPISBAIJ_Once()
395 PetscCall(PetscFree4(len_s, btable, iwork, Bowners)); in MatIncreaseOverlap_MPISBAIJ_Once()
|
| /petsc/src/ksp/ksp/impls/cg/gltr/ |
| H A D | gltrimpl.h | 30 PetscBLASInt *iwork; /* Integer workspace for solver computations */ member
|
| H A D | gltr.c | 626 PetscCall(PetscFree(cg->iwork)); in KSPCGSolve_GLTR()
635 PetscCall(PetscMalloc2(10 * cg->alloced, &cg->rwork, 5 * cg->alloced, &cg->iwork)); in KSPCGSolve_GLTR()
646 e_iblk = cg->iwork + 0 * t_size; /* Eigen blocks */ in KSPCGSolve_GLTR()
647 e_splt = cg->iwork + 1 * t_size; /* Eigen splits */ in KSPCGSolve_GLTR()
648 e_iwrk = cg->iwork + 2 * t_size; /* Eigen workspace */ in KSPCGSolve_GLTR()
921 if (cg->alloced) PetscCall(PetscFree2(cg->rwork, cg->iwork)); in KSPCGDestroy_GLTR()
|
| /petsc/include/ |
| H A D | petscdt.h | 410 PetscInt iwork[PETSC_FACTORIAL_MAX]; in PetscDTPermIndex() local
417 for (i = 0; i < n; i++) iwork[i] = i; /* partial permutation inverse */ in PetscDTPermIndex()
421 PetscInt jloc = iwork[j]; in PetscDTPermIndex()
428 iwork[j] = i; in PetscDTPermIndex()
429 iwork[icur] = jloc; in PetscDTPermIndex()
|
| /petsc/src/ksp/ksp/impls/gmres/agmres/ |
| H A D | agmresdeflation.c | 109 PetscBLASInt *iwork = agmres->iwork; in KSPAGMRESSchurForm() local
158 …, &ldB, wr, wi, beta, Q, &N, Z, &N, &r, NULL, NULL, &Dif[0], work, &lwork, iwork, &liwork, &info)); in KSPAGMRESSchurForm()
|
| H A D | agmres.c | 64 PetscCall(PetscMalloc2(N + 1, &agmres->perm, 2 * neig * N, &agmres->iwork)); in KSPSetUp_AGMRES()
|
| /petsc/src/tao/leastsquares/impls/pounders/ |
| H A D | pounders.h | 35 PetscBLASInt *iwork; /* n */ member
|
| H A D | pounders.c | 734 …PetscCallBLAS("LAPACKgesv", LAPACKgesv_(&blasn, &blasm, mfqP->Disp, &blasnpmax, mfqP->iwork, mfqP-… in TaoSolve_POUNDERS()
983 PetscCall(PetscMalloc1(mfqP->n, &mfqP->iwork)); in TaoSetUp_POUNDERS()
1095 PetscCall(PetscFree(mfqP->iwork)); in TaoDestroy_POUNDERS()
|
| /petsc/src/mat/impls/dense/seq/ |
| H A D | dense.c | 1355 PetscInt nnz = m * N, *iwork; in MatView_Dense_Binary() local
1357 PetscCall(PetscMalloc1(nnz, &iwork)); in MatView_Dense_Binary()
1358 for (i = 0; i < m; i++) iwork[i] = N; in MatView_Dense_Binary()
1359 …PetscCall(PetscViewerBinaryWriteAll(viewer, iwork, m, PETSC_DETERMINE, PETSC_DETERMINE, PETSC_INT)… in MatView_Dense_Binary()
1362 for (j = 0; j < N; j++, k++) iwork[k] = j; in MatView_Dense_Binary()
1363 …PetscCall(PetscViewerBinaryWriteAll(viewer, iwork, nnz, PETSC_DETERMINE, PETSC_DETERMINE, PETSC_IN… in MatView_Dense_Binary()
1364 PetscCall(PetscFree(iwork)); in MatView_Dense_Binary()
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/ |
| H A D | dspacelagrange.c | 1499 PetscInt *iwork; in PetscDualSpaceCreateAllDataFromInteriorData() local
1544 …0, dim * dim, &J, dim * dim, &Jinv, Nk * Nk, &L, maxNzforms * Nk, &work, maxNzforms * Nk, &iwork)); in PetscDualSpaceCreateAllDataFromInteriorData()
1616 for (d = 0; d < Nk; d++) iwork[l * Nk + d] = (blockcol + countNodesIn) * Nk + d; in PetscDualSpaceCreateAllDataFromInteriorData()
1625 PetscCall(MatSetValues(allMat, 1, &row, (ncols / pNk) * Nk, iwork, work, INSERT_VALUES)); in PetscDualSpaceCreateAllDataFromInteriorData()
1634 PetscCall(PetscFree7(v0, pv0, J, Jinv, L, work, iwork)); in PetscDualSpaceCreateAllDataFromInteriorData()
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 1432 PetscBLASInt *iwork; in PetscDTSymmetricTridiagonalEigensolve() local
1450 PetscCall(PetscMalloc2(lwork, &work, liwork, &iwork)); in PetscDTSymmetricTridiagonalEigensolve()
1452 …ag, &VL, &VU, &IL, &IU, &abstol, &bm, eigs, V, &ldz, isuppz, work, &lwork, iwork, &liwork, &info)); in PetscDTSymmetricTridiagonalEigensolve()
1455 PetscCall(PetscFree2(work, iwork)); in PetscDTSymmetricTridiagonalEigensolve()
|