| /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
148 iu = m; in main()
152 PetscCall(PetscBLASIntCast(0.2 * m, &iu)); in main()
153 …ll(PetscPrintf(PETSC_COMM_WORLD, " LAPACKsyevx: compute %d to %d-th eigensolutions...\n", il, iu)); in main()
163 …LAPACKsyevx_("V", "I", "U", &bn, arrayA, &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_arr… in main()
176 iu = m; in main()
182 PetscCall(PetscBLASIntCast(0.2 * m, &iu)); in main()
183 …all(PetscPrintf(PETSC_COMM_WORLD, " LAPACKsygv: compute %d to %d-th eigensolutions...\n", il, iu)); in main()
192 …LAPACKsygvx_(&one, "V", "I", "U", &bn, arrayA, &bn, arrayB, &bn, &vl, &vu, &il, &iu, &abstol, &nev… in main()
218 PetscCall(CkEigenSolutions(cklvl, A, il - 1, iu - 1, evals, evecs, tols)); in main()
[all …]
|
| H A D | ex118.c | 17 PetscBLASInt n, il = 1, iu = 5, *iblock, *isplit, *iwork, nevs, *ifail, cklvl = 2; in main()
32 nevs = iu - il; in main()
49 …LAPACKstebz_("I", "E", &n, &vl, &vu, &il, &iu, &tol, (PetscReal *)D, (PetscReal *)E, &nevs, &nspli… in main()
90 PetscCall(CkEigenSolutions(cklvl, T, il - 1, iu - 1, evals, evecs, tols)); in main()
125 PetscErrorCode CkEigenSolutions(PetscInt cklvl, Mat A, PetscInt il, PetscInt iu, PetscScalar *eval,… in CkEigenSolutions() argument
134 nev = iu - il; in CkEigenSolutions()
143 for (i = il; i < iu; i++) { in CkEigenSolutions()
145 for (j = il; j < iu; j++) { in CkEigenSolutions()
165 for (i = il; i < iu; i++) { in CkEigenSolutions()
|
| H A D | ex116.c | 26 PetscBLASInt nevs, il, iu, in; in main() local
77 PetscCall(PetscBLASIntCast(m, &iu)); in main()
80 PetscCall(PetscBLASIntCast(0.2 * m, &iu)); in main()
82 …ACKsyevx: compute %" PetscBLASInt_FMT " to %" PetscBLASInt_FMT "-th eigensolutions...\n", il, iu)); in main()
90 …LAPACKsyevx_("V", "I", "U", &bn, arrayA, &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_arr… in main()
113 PetscCall(CkEigenSolutions(cklvl, A, il - 1, iu - 1, evals, evecs, tols)); in main()
196 PetscErrorCode CkEigenSolutions(PetscInt cklvl, Mat A, PetscInt il, PetscInt iu, PetscReal *eval, V… in CkEigenSolutions() argument
203 nev = iu - il; in CkEigenSolutions()
213 for (i = il; i < iu; i++) { in CkEigenSolutions()
215 for (j = il; j < iu; j++) { in CkEigenSolutions()
[all …]
|
| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex5.c | 25 PetscInt iu, ip; in main() local
45 PetscCall(DMStagGetLocationSlot(dmSol, LEFT, 0, &iu)); in main()
49 arrSol[e][iu] = 2 * PRESSURE_CONST; in main()
50 arrRHS[e][iu] = 0.0; in main()
|
| /petsc/src/snes/tests/ |
| H A D | ex8.c | 526 Vec iu, fu, scaling = NULL; in CheckTransfer() local
555 PetscCall(DMGetGlobalVector(idm, &iu)); in CheckTransfer()
561 PetscCall(DMProjectFunction(idm, 0.0, exactFuncs, exactCtxs, INSERT_ALL_VALUES, iu)); in CheckTransfer()
566 PetscCall(MatInterpolate(Interp, iu, fu)); in CheckTransfer()
570 PetscCall(MatRestrict(Interp, iu, fu)); in CheckTransfer()
575 PetscCall(MatRestrict(Interp, iu, fu)); in CheckTransfer()
617 PetscCall(MatInterpolate(InterpAdapt, iu, fu)); in CheckTransfer()
637 PetscCall(DMRestoreGlobalVector(idm, &iu)); in CheckTransfer()
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex1.c | 92 PetscInt iu, ip, ixu, ixp; in main() local
140 PetscCall(DMStagGetLocationSlot(dmSol, LEFT, 0, &iu)); in main()
147 arrSol[e][iu] = a + (b - a - (c / 2.0)) * coordu + (c / 2.0) * coordu * coordu; in main()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 45 PetscInt *jl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR() local
69 PetscCall(PetscShmgetAllocateArray(mbs + 1, sizeof(PetscInt), (void **)&iu)); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
73 iu[0] = mbs + 1; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
112 jmin = iu[prow] + 1; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
113 jmax = iu[prow + 1]; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
137 iu[k + 1] = iu[k] + nzk; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
140 if (iu[k + 1] > umax) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
150 PetscCall(PetscArraycpy(jutmp, ju, iu[k])); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
166 PetscReal af = ((PetscReal)iu[mbs]) / ((PetscReal)ai[mbs]); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
184 PetscCall(PetscShmgetAllocateArray((iu[mbs] + 1) * bs2, sizeof(PetscScalar), (void **)&b->a)); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
[all …]
|
| H A D | sbaijfact2.c | 2285 PetscInt *prowl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd; in MatICCFactorSymbolic_SeqSBAIJ_MSR() local
2305 PetscCall(PetscMalloc1(mbs + 1, &iu)); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2307 iu[0] = mbs + 1; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2348 jmin = iu[prow] + 1; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2349 jmax = iu[prow + 1]; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2376 iu[k + 1] = iu[k] + nzk; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2379 if (iu[k + 1] > umax) { in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2389 PetscCall(PetscArraycpy(jutmp, ju, iu[k])); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2394 PetscCall(PetscArraycpy(jutmp, lev, iu[k] - shift)); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2412 PetscReal af = ((PetscReal)iu[mbs]) / ((PetscReal)ai[mbs]); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
[all …]
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex3.c | 803 Vec iu, fu, scaling; in CheckInterpolation() local
848 PetscCall(DMGetGlobalVector(idm, &iu)); in CheckInterpolation()
854 PetscCall(DMProjectFunction(idm, 0.0, exactFuncs, exactCtxs, INSERT_ALL_VALUES, iu)); in CheckInterpolation()
857 PetscCall(MatRestrict(Interp, iu, fu)); in CheckInterpolation()
859 } else PetscCall(MatInterpolate(Interp, iu, fu)); in CheckInterpolation()
869 PetscCall(DMRestoreGlobalVector(idm, &iu)); in CheckInterpolation()
|
| /petsc/src/ksp/ksp/impls/cg/gltr/ |
| H A D | gltr.c | 77 PetscBLASInt t_size = 0, l_size = 0, il, iu, info; in KSPCGSolve_GLTR()
654 iu = 1; in KSPCGSolve_GLTR()
656 …PetscCallBLAS("LAPACKstebz", LAPACKstebz_("I", "E", &t_size, &vl, &vu, &il, &iu, &cg->eigen_tol, c… in KSPCGSolve_GLTR()
|