| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 25 static PetscErrorCode KSPDGMRESComputeSchurForm(KSP ksp, PetscInt *neig) in KSPDGMRESComputeSchurForm() argument
28 PetscUseMethod(ksp, "KSPDGMRESComputeSchurForm_C", (KSP, PetscInt *), (ksp, neig)); in KSPDGMRESComputeSchurForm()
44 static PetscErrorCode KSPDGMRESImproveEig(KSP ksp, PetscInt neig) in KSPDGMRESImproveEig() argument
47 PetscUseMethod(ksp, "KSPDGMRESImproveEig_C", (KSP, PetscInt), (ksp, neig)); in KSPDGMRESImproveEig()
54 PetscInt neig = dgmres->neig + EIG_OFFSET; in KSPSetUp_DGMRES() local
59 if (!dgmres->neig) PetscFunctionReturn(PETSC_SUCCESS); in KSPSetUp_DGMRES()
62 PetscCall(PetscMalloc1(neig * max_k, &SR)); in KSPSetUp_DGMRES()
264 PetscInt neig1 = dgmres->neig + EIG_OFFSET; in KSPDestroy_DGMRES()
469 …quency of extracted eigenvalues = %" PetscInt_FMT " using Harmonic Ritz values \n", dgmres->neig)); in KSPView_DGMRES()
471 …, " Frequency of extracted eigenvalues = %" PetscInt_FMT " using Ritz values \n", dgmres->neig)); in KSPView_DGMRES()
[all …]
|
| H A D | dgmresimpl.h | 13 PetscInt neig; /* number of eigenvalues to extract at each restart */ \
|
| /petsc/src/ksp/ksp/impls/gmres/agmres/ |
| H A D | agmres.c | 28 PetscInt neig = agmres->neig; in KSPSetUp_AGMRES() local
31 PetscInt lwork = PetscMax(8 * N + 16, 4 * neig * (N - neig)); in KSPSetUp_AGMRES()
56 if (agmres->neig < 1) PetscFunctionReturn(PETSC_SUCCESS); in KSPSetUp_AGMRES()
64 PetscCall(PetscMalloc2(N + 1, &agmres->perm, 2 * neig * N, &agmres->iwork)); in KSPSetUp_AGMRES()
122 …if ((agmres->r == 0) && (agmres->neig > 0)) { /* Compute the eigenvalues for the shifts and the ei… in KSPComputeShifts_DGMRES()
148 …if (agmres->neig > 0) { /* Compute the eigenvalues for the shifts) and the eigenvectors (to augmen… in KSPComputeShifts_DGMRES()
474 if (!ksp->reason && agmres->neig > 0) { in KSPSolve_AGMRES()
478 …gmres->bgv * (ksp->max_it - ksp->its)) && ((agmres->r + 1) < agmres->max_neig))) agmres->neig += 1; in KSPSolve_AGMRES()
479 PetscCall(KSPDGMRESComputeDeflationData_DGMRES(ksp, &agmres->neig)); in KSPSolve_AGMRES()
503 if (agmres->neig) { in KSPDestroy_AGMRES()
[all …]
|
| H A D | agmresdeflation.c | 97 PetscInt neig = agmres->neig; in KSPAGMRESSchurForm() local
121 PetscCall(PetscBLASIntCast(PetscMax(8 * N + 16, 4 * neig * (N - neig)), &lwork)); in KSPAGMRESSchurForm()
122 PetscCall(PetscBLASIntCast(2 * N * neig, &liwork)); in KSPAGMRESSchurForm()
147 while (r < neig) { in KSPAGMRESSchurForm()
200 if (agmres->neig <= 1) PetscFunctionReturn(PETSC_SUCCESS); in KSPAGMRESComputeDeflationData()
|
| /petsc/src/ts/interface/ |
| H A D | tseig.c | 96 PetscInt n, N, nits, neig, i, its = 200; in TSMonitorSPEig() local
143 neig = n; in TSMonitorSPEig()
145 PetscCall(KSPComputeEigenvalues(ksp, N, r, c, &neig)); in TSMonitorSPEig()
148 for (i = 0; i < neig; i++) r[i] = -r[i]; in TSMonitorSPEig()
149 for (i = 0; i < neig; i++) { in TSMonitorSPEig()
|
| /petsc/src/ksp/ksp/impls/cg/ |
| H A D | cgeig.c | 9 …orCode KSPComputeEigenvalues_CG(KSP ksp, PetscInt nmax, PetscReal *r, PetscReal *c, PetscInt *neig) in KSPComputeEigenvalues_CG() argument
19 *neig = n; in KSPComputeEigenvalues_CG()
|
| /petsc/src/ksp/ksp/interface/ |
| H A D | eige.c | 253 PetscErrorCode KSPPlotEigenContours_Private(KSP ksp, PetscInt neig, const PetscReal *r, const Petsc… in KSPPlotEigenContours_Private() argument
271 for (i = 1; i < neig; i++) { in KSPPlotEigenContours_Private()
280 PetscCall(PolyEval(neig, r, c, 0, 0, &px0, &py0)); in KSPPlotEigenContours_Private()
286 PetscCall(PolyEval(neig, r, c, xloc[i], yloc[j], &px, &py)); in KSPPlotEigenContours_Private()
|
| H A D | itfunc.c | 114 …cErrorCode KSPComputeEigenvalues(KSP ksp, PetscInt n, PetscReal r[], PetscReal c[], PetscInt *neig) in KSPComputeEigenvalues() argument
121 PetscAssertPointer(neig, 5); in KSPComputeEigenvalues()
124 if (n && ksp->ops->computeeigenvalues) PetscUseTypeMethod(ksp, computeeigenvalues, n, r, c, neig); in KSPComputeEigenvalues()
125 else *neig = 0; in KSPComputeEigenvalues()
676 PetscInt n, i, neig; in KSPViewEigenvalues_Internal() local
688 neig = n; in KSPViewEigenvalues_Internal()
699 PetscCall(KSPComputeEigenvalues(ksp, n, r, c, &neig)); in KSPViewEigenvalues_Internal()
703 for (i = 0; i < neig; ++i) { in KSPViewEigenvalues_Internal()
712 PetscCall(KSPPlotEigenContours_Private(ksp, neig, r, c)); in KSPViewEigenvalues_Internal()
717 for (i = 0; i < neig; ++i) PetscCall(PetscDrawSPAddPoint(drawsp, r + i, c + i)); in KSPViewEigenvalues_Internal()
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | gmreig.c | 42 …ode KSPComputeEigenvalues_GMRES(KSP ksp, PetscInt nmax, PetscReal *r, PetscReal *c, PetscInt *neig) in KSPComputeEigenvalues_GMRES() argument
57 *neig = n; in KSPComputeEigenvalues_GMRES()
89 *neig = n; in KSPComputeEigenvalues_GMRES()
|
| /petsc/src/ksp/ksp/impls/cheby/ |
| H A D | cheby.c | 20 PetscInt n, neig; in KSPChebyshevComputeExtremeEigenvalues_Private() local
26 PetscCall(KSPComputeEigenvalues(kspest, n, re, im, &neig)); in KSPChebyshevComputeExtremeEigenvalues_Private()
29 for (n = 0; n < neig; n++) { in KSPChebyshevComputeExtremeEigenvalues_Private()
|
| /petsc/src/ksp/ksp/impls/minres/ |
| H A D | minres.c | 872 …de KSPComputeEigenvalues_MINRES(KSP ksp, PetscInt nmax, PetscReal *r, PetscReal *c, PetscInt *neig) in KSPComputeEigenvalues_MINRES() argument
882 *neig = n; in KSPComputeEigenvalues_MINRES()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | KSP.pyx | 2110 cdef PetscInt neig = 0
2116 CHKERR(KSPComputeEigenvalues(self.ksp, its, rdata, idata, &neig))
2117 eigen = empty_c(neig)
2118 eigen.real = r[:neig]
2119 eigen.imag = i[:neig]
|
| /petsc/doc/changes/ |
| H A D | 2024.md | 99 `KSPComputeEigenvalues(KSP,int,double*,double*,int *neig);`
|
| /petsc/src/ksp/ksp/impls/fetidp/ |
| H A D | fetidp.c | 239 …de KSPComputeEigenvalues_FETIDP(KSP ksp, PetscInt nmax, PetscReal *r, PetscReal *c, PetscInt *neig) in KSPComputeEigenvalues_FETIDP() argument
244 PetscCall(KSPComputeEigenvalues(fetidp->innerksp, nmax, r, c, neig)); in KSPComputeEigenvalues_FETIDP()
|
| /petsc/doc/manual/ |
| H A D | ksp.md | 591 KSPComputeEigenvalues(KSP ksp,PetscInt n,PetscReal *realpart,PetscReal *complexpart,PetscInt *neig);
595 inserted into those two arrays. `neig` is the number of eigenvalues
|