| /petsc/src/ksp/ksp/interface/ |
| H A D | eige.c | 190 PetscScalar *work, *eigs; in KSPComputeEigenvaluesExplicitly() local
199 PetscCall(PetscMalloc1(n, &eigs)); in KSPComputeEigenvaluesExplicitly()
206 …PetscCallBLAS("LAPACKgeev", LAPACKgeev_("N", "N", &nb, array, &nb, eigs, &sdummy, &idummy, &sdummy… in KSPComputeEigenvaluesExplicitly()
214 for (i = 0; i < n; i++) r[i] = PetscRealPart(eigs[i]); in KSPComputeEigenvaluesExplicitly()
217 r[i] = PetscRealPart(eigs[perm[i]]); in KSPComputeEigenvaluesExplicitly()
218 c[i] = PetscImaginaryPart(eigs[perm[i]]); in KSPComputeEigenvaluesExplicitly()
221 PetscCall(PetscFree(eigs)); in KSPComputeEigenvaluesExplicitly()
|
| /petsc/src/ksp/ksp/impls/cg/cgne/ |
| H A D | cgne.c | 43 PetscInt i, stored_max_it, eigs; in KSPSolve_CGNE() local
57 eigs = ksp->calc_sings; in KSPSolve_CGNE()
68 if (eigs) { in KSPSolve_CGNE()
127 if (eigs) { in KSPSolve_CGNE()
139 if (eigs) d[i] = PetscSqrtReal(PetscAbsScalar(b)) * e[i] + 1.0 / a; in KSPSolve_CGNE()
|
| /petsc/src/ksp/ksp/guess/impls/pod/ |
| H A D | pod.c | 25 PetscReal *eigs; /* eigenvalues */ member
76 …PetscCall(PetscCalloc6(pod->maxn * pod->maxn, &pod->corr, pod->maxn, &pod->eigs, pod->maxn * pod->… in KSPGuessSetUp_POD()
86 …od->corr, &bN, &rdummy, &rdummy, &idummy, &idummy, &rdummy, &idummy, pod->eigs, pod->eigv, &bN, &s… in KSPGuessSetUp_POD()
88 …od->corr, &bN, &rdummy, &rdummy, &idummy, &idummy, &rdummy, &idummy, pod->eigs, pod->eigv, &bN, &s… in KSPGuessSetUp_POD()
122 PetscCall(PetscFree6(pod->corr, pod->eigs, pod->eigv, pod->iwork, pod->yhay, pod->low)); in KSPGuessDestroy_POD()
313 …bN, pod->swork, &bN, &reps, &reps, &idummy, &idummy, &reps, &idummy, pod->eigs, pod->eigv, &bN, po… in KSPGuessUpdate_POD()
315 …bN, pod->swork, &bN, &reps, &reps, &idummy, &idummy, &reps, &idummy, pod->eigs, pod->eigv, &bN, po… in KSPGuessUpdate_POD()
323 pod->eigs[i] = PetscMax(pod->eigs[i], 0.0); in KSPGuessUpdate_POD()
324 toten += pod->eigs[i]; in KSPGuessUpdate_POD()
325 if (!pod->eigs[i]) pod->st = i; in KSPGuessUpdate_POD()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/ |
| H A D | cg.c | 121 PetscInt i, stored_max_it, eigs; in KSPSolve_CG() local
135 eigs = ksp->calc_sings; in KSPSolve_CG()
145 if (eigs) { in KSPSolve_CG()
245 if (eigs) { in KSPSolve_CG()
290 if (eigs) d[i] = PetscSqrtReal(PetscAbsScalar(b)) * e[i] + 1.0 / a; in KSPSolve_CG()
366 PetscInt i, stored_max_it, eigs; in KSPSolve_CG_SingleReduction() local
379 eigs = ksp->calc_sings; in KSPSolve_CG_SingleReduction()
389 if (eigs) { in KSPSolve_CG_SingleReduction()
463 if (eigs) { in KSPSolve_CG_SingleReduction()
487 if (eigs) d[i] = PetscSqrtReal(PetscAbsScalar(b)) * e[i] + 1.0 / a; in KSPSolve_CG_SingleReduction()
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | gmreig.c | 80 PetscScalar *R = gmres->Rsvd, *work = R + N * N, *eigs = work + 5 * N, sdummy; in KSPComputeEigenvalues_GMRES()
98 …PetscCallBLAS("LAPACKgeev", LAPACKgeev_("N", "N", &bn, R, &bN, eigs, &sdummy, &idummy, &sdummy, &i… in KSPComputeEigenvalues_GMRES()
103 for (i = 0; i < n; i++) r[i] = PetscRealPart(eigs[i]); in KSPComputeEigenvalues_GMRES()
106 r[i] = PetscRealPart(eigs[perm[i]]); in KSPComputeEigenvalues_GMRES()
107 c[i] = PetscImaginaryPart(eigs[perm[i]]); in KSPComputeEigenvalues_GMRES()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexmetric.c | 1056 …PetscReal *eigs, max_eig, l_min = 1.0 / (h_max * h_max), l_max = 1.0 / (h_min * h_min), la_min =… in DMPlexMetricModify_Private() local
1060 PetscCall(PetscMalloc2(dim * dim, &Mpos, dim, &eigs)); in DMPlexMetricModify_Private()
1074 eigs[0] = PetscRealPart(Mpos[0]); 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()
1118 eigs[i] = PetscMin(l_max, PetscMax(l_min, PetscAbsReal(eigs[i]))); in DMPlexMetricModify_Private()
1119 max_eig = PetscMax(eigs[i], max_eig); in DMPlexMetricModify_Private()
1125 if (a_max >= 1.0) eigs[i] = PetscMax(eigs[i], max_eig * la_min); in DMPlexMetricModify_Private()
1126 *detMp *= eigs[i]; in DMPlexMetricModify_Private()
1133 for (k = 0; k < dim; ++k) Mp[i * dim + j] += Mpos[k * dim + i] * eigs[k] * Mpos[k * dim + j]; in DMPlexMetricModify_Private()
[all …]
|
| /petsc/src/ksp/ksp/impls/minres/ |
| H A D | minres.c | 112 PetscInt stored_max_it, eigs; in KSPSolve_MINRES() local
120 eigs = ksp->calc_sings; in KSPSolve_MINRES()
122 if (eigs) { in KSPSolve_MINRES()
188 if (eigs) { in KSPSolve_MINRES()
488 PetscInt stored_max_it, eigs; in KSPSolve_MINRES_OLD() local
510 eigs = ksp->calc_sings; in KSPSolve_MINRES_OLD()
512 if (eigs) { in KSPSolve_MINRES_OLD()
559 if (eigs) { in KSPSolve_MINRES_OLD()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex30.c | 131 PetscReal eigs[3], s = 0.0; in FIX_C_3d() local
133 Eigenvalues_Sym3x3(C00, C01, C02, C11, C12, C22, eigs); in FIX_C_3d()
134 if (twod) eigs[2] = 1.0; in FIX_C_3d()
135 …if (eigs[0] <= 0 || eigs[1] <= 0 || eigs[2] <= 0) s = -PetscMin(eigs[0], PetscMin(eigs[1], eigs[2]… in FIX_C_3d()
669 PetscReal eigs[3]; in ellipticity_fail_private() local
685 Eigenvalues_Sym3x3(C00, C01, C02, C11, C12, C22, eigs); in ellipticity_fail_private()
686 …if (eigs[0] < 0 || eigs[1] < 0 || eigs[2] < 0) obj[0] = -PetscMin(eigs[0], PetscMin(eigs[1], eigs[… in ellipticity_fail_private()
835 PetscReal eigs[3]; in eigsc() local
849 Eigenvalues_Sym3x3(C00, C01, C02, C11, C12, C22, eigs); in eigsc()
850 PetscCallVoid(PetscSortReal(dim, eigs)); in eigsc()
[all …]
|
| /petsc/src/ksp/ksp/impls/fcg/ |
| H A D | fcg.c | 80 PetscInt eigs = ksp->calc_sings; /* Variables for eigen estimation - START*/ in KSPSolve_FCG() local
94 if (eigs) { in KSPSolve_FCG()
242 if (eigs) { in KSPSolve_FCG()
|
| /petsc/src/ksp/ksp/tutorials/output/ |
| H A D | ex71_bddc_elast_4lev.out | 29 Use estimated eigs for coarse problem: 0
142 Use estimated eigs for coarse problem: 0
261 Use estimated eigs for coarse problem: 0
|
| H A D | ex59_bddc_fetidp_ml_2.out | 29 Use estimated eigs for coarse problem: 0
140 Use estimated eigs for coarse problem: 0
|
| H A D | ex59_bddc_fetidp_ml_1.out | 29 Use estimated eigs for coarse problem: 0
140 Use estimated eigs for coarse problem: 0
|
| H A D | ex71_bddc_elast_3lev.out | 29 Use estimated eigs for coarse problem: 0
142 Use estimated eigs for coarse problem: 0
|
| H A D | ex71_bddc_elast_3lev_alt.out | 29 Use estimated eigs for coarse problem: 0
142 Use estimated eigs for coarse problem: 0
|
| H A D | ex71_bddc_elast_deluxe_layers.out | 29 Use estimated eigs for coarse problem: 0
|
| H A D | ex71_bddc_cusparse.out | 29 Use estimated eigs for coarse problem: 0
|
| H A D | ex71_bddc_elast.out | 29 Use estimated eigs for coarse problem: 0
|
| H A D | ex71_bddc_1.out | 29 Use estimated eigs for coarse problem: 0
|
| H A D | ex71_bddc_2.out | 29 Use estimated eigs for coarse problem: 0
|
| H A D | ex59_bddc_fetidp_ml_3.out | 29 Use estimated eigs for coarse problem: 0
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 3360 PetscReal *eigs, thresh, lthresh, uthresh; in PCBDDCAdaptiveSelection() local
3429 eigs = &rdummy; in PCBDDCAdaptiveSelection()
3439 …B_N, St, &B_N, S, &B_N, &zero, &thresh, B_iwork, B_iwork, &eps, &B_neigs, eigs, eigv, &B_N, &lwork… in PCBDDCAdaptiveSelection()
3441 …B_N, St, &B_N, S, &B_N, &zero, &thresh, B_iwork, B_iwork, &eps, &B_neigs, eigs, eigv, &B_N, &lwork… in PCBDDCAdaptiveSelection()
3453 …PetscCall(PetscMalloc5(mss * mss, &eigv, mss, &eigs, B_lwork, &work, 5 * mss, &B_iwork, mss, &B_if… in PCBDDCAdaptiveSelection()
3597 …", &B_N, St, &B_N, S, &B_N, &lower, &upper, &B_IL, &B_IU, &eps, &B_neigs, eigs, eigv, &B_N, work, … in PCBDDCAdaptiveSelection()
3599 …", &B_N, St, &B_N, S, &B_N, &lower, &upper, &B_IL, &B_IU, &eps, &B_neigs, eigs, eigv, &B_N, work, … in PCBDDCAdaptiveSelection()
3617 …", &B_N, St, &B_N, S, &B_N, &bb[0], &bb[1], &B_IL, &B_IU, &eps, &B_neigs, eigs, eigv, &B_N, work, … in PCBDDCAdaptiveSelection()
3619 …", &B_N, St, &B_N, S, &B_N, &bb[0], &bb[1], &B_IL, &B_IU, &eps, &B_neigs, eigs, eigv, &B_N, work, … in PCBDDCAdaptiveSelection()
3627 …", &B_N, St, &B_N, S, &B_N, &bb[0], &bb[1], &B_IL, &B_IU, &eps, &B_neigs, eigs, eigv, &B_N, work, … in PCBDDCAdaptiveSelection()
[all …]
|
| /petsc/src/snes/tutorials/output/ |
| H A D | ex12_p4est_full_q2_nonconformal_parallel_bddcfas.out | 60 Use estimated eigs for coarse problem: 0
235 Use estimated eigs for coarse problem: 0
411 Use estimated eigs for coarse problem: 0
|
| H A D | ex12_p4est_full_q2_nonconformal_parallel_bddcfas_alt.out | 60 Use estimated eigs for coarse problem: 0
235 Use estimated eigs for coarse problem: 0
411 Use estimated eigs for coarse problem: 0
|
| H A D | ex12_tri_parmetis_bddc.out | 54 Use estimated eigs for coarse problem: 0
|
| H A D | ex12_tri_bddc.out | 57 Use estimated eigs for coarse problem: 0
|