| /petsc/src/ksp/pc/impls/sor/ |
| H A D | sor.c | 7 PetscInt its; /* inner iterations, number of sweeps */ member
33 PetscCall(MatSOR(pc->pmat, x, jac->omega, (MatSORType)flag, jac->fshift, jac->its, jac->lits, y)); in PCApply_SOR()
47 PetscCall(MatSOR(pc->pmat, x, jac->omega, (MatSORType)flag, jac->fshift, jac->its, jac->lits, y)); in PCApplyTranspose_SOR()
52 …, Vec y, Vec w, PetscReal rtol, PetscReal abstol, PetscReal dtol, PetscInt its, PetscBool guesszer… in PCApplyRichardson_SOR() argument
59 PetscCall(MatSOR(pc->pmat, b, jac->omega, stype, jac->fshift, its * jac->its, jac->lits, y)); in PCApplyRichardson_SOR()
61 *outits = its; in PCApplyRichardson_SOR()
77 …"-pc_sor_its", "number of inner SOR iterations", "PCSORSetIterations", jac->its, &jac->its, NULL)); in PCSetFromOptions_SOR()
116 …_FMT ", local iterations = %" PetscInt_FMT ", omega = %g\n", sortype, jac->its, jac->lits, (double… in PCView_SOR()
140 static PetscErrorCode PCSORSetIterations_SOR(PC pc, PetscInt its, PetscInt lits) in PCSORSetIterations_SOR() argument
145 jac->its = its; in PCSORSetIterations_SOR()
[all …]
|
| /petsc/src/ksp/pc/impls/ksp/ |
| H A D | pcksp.c | 7 PetscInt its; /* total number of iterations KSP uses */ member
34 PetscInt its; in PCApply_KSP() local
45 PetscCall(KSPGetIterationNumber(jac->ksp, &its)); in PCApply_KSP()
46 jac->its += its; in PCApply_KSP()
52 PetscInt its; in PCMatApply_KSP() local
63 PetscCall(KSPGetIterationNumber(jac->ksp, &its)); in PCMatApply_KSP()
64 jac->its += its; in PCMatApply_KSP()
70 PetscInt its; in PCApplyTranspose_KSP() local
81 PetscCall(KSPGetIterationNumber(jac->ksp, &its)); in PCApplyTranspose_KSP()
82 jac->its += its; in PCApplyTranspose_KSP()
[all …]
|
| /petsc/src/ts/tutorials/output/ |
| H A D | ex30_0_dirk.out | 23 step 0 stage 0 nonlinear its 0, linear its 0
26 step 0 stage 1 nonlinear its 1, linear its 1
29 step 0 stage 2 nonlinear its 1, linear its 1
|
| H A D | ex30_0_dirk_fieldsplit.out | 23 step 0 stage 0 nonlinear its 0, linear its 0
26 step 0 stage 1 nonlinear its 1, linear its 1
29 step 0 stage 2 nonlinear its 1, linear its 1
|
| H A D | ex30_hdf5_diagnostic.out | 17 step 0 stage 0 nonlinear its 1, linear its 1
19 step 1 stage 0 nonlinear its 1, linear its 1
|
| H A D | ex30_annulus.out | 17 step 0 stage 0 nonlinear its 1, linear its 1
19 step 1 stage 0 nonlinear its 1, linear its 1
|
| H A D | ex30_vtk_diagnostic.out | 17 step 0 stage 0 nonlinear its 1, linear its 1
19 step 1 stage 0 nonlinear its 1, linear its 1
|
| H A D | ex30_0_dirk_mg.out | 25 step 0 stage 0 nonlinear its 0, linear its 0
30 step 0 stage 1 nonlinear its 1, linear its 1
35 step 0 stage 2 nonlinear its 1, linear its 1
|
| /petsc/src/snes/interface/ |
| H A D | snesut.c | 29 PetscErrorCode SNESMonitorSolution(SNES snes, PetscInt its, PetscReal fgnorm, PetscViewerAndFormat … in SNESMonitorSolution() argument
66 PetscErrorCode SNESMonitorResidual(SNES snes, PetscInt its, PetscReal fgnorm, PetscViewerAndFormat … in SNESMonitorResidual() argument
102 PetscErrorCode SNESMonitorSolutionUpdate(SNES snes, PetscInt its, PetscReal fgnorm, PetscViewerAndF… in SNESMonitorSolutionUpdate() argument
298 PetscErrorCode SNESMonitorDefault(SNES snes, PetscInt its, PetscReal fgnorm, PetscViewerAndFormat *… in SNESMonitorDefault() argument
325 …T " SNES Function norm %14.12e, Update norm %14.12e, Objective %14.12e\n", its, (double)fgnorm, (d… in SNESMonitorDefault()
327 …r, "%3" PetscInt_FMT " SNES Function norm %14.12e, Update norm %14.12e\n", its, (double)fgnorm, (d… in SNESMonitorDefault()
330 …ViewerASCIIPrintf(viewer, "%3" PetscInt_FMT " SNES Function norm %14.12e\n", its, (double)fgnorm)); in SNESMonitorDefault()
339 if (!its) PetscCall(PetscDrawLGReset(lg)); in SNESMonitorDefault()
340 x = (PetscReal)its; in SNESMonitorDefault()
344 if (its <= 20 || !(its % 5) || snes->reason) { in SNESMonitorDefault()
[all …]
|
| /petsc/src/ksp/ksp/impls/tsirm/ |
| H A D | tsirm.c | 40 PetscInt i, *ind_row, first_iteration = 1, its = 0, total = 0, col = 0; in KSPSolve_TSIRM() local
64 ksp->its = 0; in KSPSolve_TSIRM()
65 PetscCall(KSPConvergedDefault(ksp, ksp->its, norm, &ksp->reason, ksp->cnvP)); in KSPSolve_TSIRM()
66 PetscCall(KSPMonitor(ksp, ksp->its, norm)); in KSPSolve_TSIRM()
72 PetscCall(KSPGetIterationNumber(sub_ksp, &its)); in KSPSolve_TSIRM()
73 total += its; in KSPSolve_TSIRM()
82 ksp->its++; in KSPSolve_TSIRM()
83 PetscCall(KSPConvergedDefault(ksp, ksp->its, norm, &ksp->reason, ksp->cnvP)); in KSPSolve_TSIRM()
84 PetscCall(KSPMonitor(ksp, ksp->its, norm)); in KSPSolve_TSIRM()
114 } while (ksp->its < ksp->max_it && !ksp->reason); in KSPSolve_TSIRM()
[all …]
|
| /petsc/src/ksp/ksp/impls/tcqmr/ |
| H A D | tcqmr.c | 18 ksp->its = 0; in KSPSolve_TCQMR()
56 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_TCQMR()
58 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPSolve_TCQMR()
59 ksp->its++; in KSPSolve_TCQMR()
102 if (ksp->its > 2) { in KSPSolve_TCQMR()
106 if (ksp->its > 1) { in KSPSolve_TCQMR()
138 rnorm = rnorm0 * PetscSqrtReal((PetscReal)ksp->its + 2.0) * PetscRealPart(sprod); in KSPSolve_TCQMR()
141 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_TCQMR()
142 if (ksp->its >= ksp->max_it) { in KSPSolve_TCQMR()
147 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPSolve_TCQMR()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex53.c | 15 PetscInt i, n = 10, col[3], its; in main() local
90 PetscCall(KSPGetIterationNumber(ksp, &its)); in main()
91 …C_COMM_WORLD, "1. Norm of error for Ax=b: %g, Iterations %" PetscInt_FMT "\n", (double)norm, its)); in main()
100 PetscCall(KSPGetIterationNumber(ksp, &its)); in main()
101 …OMM_WORLD, "2. Norm of error for A^T x=b: %g, Iterations %" PetscInt_FMT "\n", (double)norm, its)); in main()
119 PetscCall(KSPGetIterationNumber(ksp, &its)); in main()
120 …RLD, "3. Norm of error for (A+Delta) x=b: %g, Iterations %" PetscInt_FMT "\n", (double)norm, its)); in main()
|
| /petsc/src/ksp/ksp/impls/cg/ |
| H A D | cgls.c | 48 ksp->its = 0; in KSPSolve_CGLS()
57 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_CGLS()
73 ksp->its++; in KSPSolve_CGLS()
75 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPSolve_CGLS()
76 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_CGLS()
81 } while (ksp->its < ksp->max_it); in KSPSolve_CGLS()
83 if (ksp->its >= ksp->max_it) ksp->reason = KSP_DIVERGED_ITS; in KSPSolve_CGLS()
|
| /petsc/src/ksp/ksp/impls/ibcgs/ |
| H A D | ibcgs.c | 154 for (ksp->its = 1; ksp->its < ksp->max_it + 1; ksp->its++) { in KSPSolve_IBCGS()
156 if (ksp->its == 1) deltan = rhon; in KSPSolve_IBCGS()
161 …T_CONVERGED, "KSPSolve has not converged due to taun is zero, iteration %" PetscInt_FMT, ksp->its); in KSPSolve_IBCGS()
227 if (ksp->lagnorm && ksp->its > 1) { in KSPSolve_IBCGS()
233 if (ksp->lagnorm && ksp->its > 1 && ksp->normtype != KSP_NORM_NONE) { in KSPSolve_IBCGS()
246 …if (ksp->lagnorm && ksp->its > 1 && ksp->normtype != KSP_NORM_NONE) rnorm = PetscSqrtReal(PetscRea… in KSPSolve_IBCGS()
249 …CONVERGED, "KSPSolve has not converged due to kappan is zero, iteration %" PetscInt_FMT, ksp->its); in KSPSolve_IBCGS()
254 …CONVERGED, "KSPSolve has not converged due to thetan is zero, iteration %" PetscInt_FMT, ksp->its); in KSPSolve_IBCGS()
276 if (!ksp->lagnorm && ksp->chknorm < ksp->its && ksp->normtype != KSP_NORM_NONE) { in KSPSolve_IBCGS()
284 PetscCall(KSPMonitor(ksp, ksp->its, rnorm)); in KSPSolve_IBCGS()
[all …]
|
| /petsc/src/ksp/ksp/impls/gcr/ |
| H A D | gcr.c | 41 if (ctx->modifypc) PetscCall((*ctx->modifypc)(ksp, ksp->its, k, ksp->rnorm, ctx->modifypc_ctx)); in KSPSolve_GCR_cycle()
58 if (ksp->its > ksp->chknorm && ksp->normtype != KSP_NORM_NONE) { in KSPSolve_GCR_cycle()
63 ksp->its++; in KSPSolve_GCR_cycle()
67 PetscCall(KSPMonitor(ksp, ksp->its, norm_r)); in KSPSolve_GCR_cycle()
69 if (ksp->its - 1 > ksp->chknorm) { in KSPSolve_GCR_cycle()
70 PetscCall((*ksp->converged)(ksp, ksp->its, norm_r, &ksp->reason, ksp->cnvP)); in KSPSolve_GCR_cycle()
74 if (ksp->its >= ksp->max_it) { in KSPSolve_GCR_cycle()
103 ksp->its = 0; in KSPSolve_GCR()
107 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm0)); in KSPSolve_GCR()
108 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm0, &ksp->reason, ksp->cnvP)); in KSPSolve_GCR()
[all …]
|
| /petsc/src/ksp/ksp/impls/minres/ |
| H A D | minres.c | 140 ksp->its = 0; in KSPSolve_MINRES()
179 ksp->its++; in KSPSolve_MINRES()
184 if (ksp->its > 1) PetscCall(VecAXPY(R3, -beta / betal, R1)); in KSPSolve_MINRES()
190 d[ksp->its - 1] = alpha; in KSPSolve_MINRES()
191 e[ksp->its - 1] = beta; in KSPSolve_MINRES()
238 if (ksp->its > 2 && minres->nutol > 0.0) { in KSPSolve_MINRES()
243 …PetscCall(PetscInfo(ksp, "it = %" PetscInt_FMT ": inexact check %g (%g / %g)\n", ksp->its - 2, (do… in KSPSolve_MINRES()
245 ksp->its--; in KSPSolve_MINRES()
259 if (ksp->its > 2) { in KSPSolve_MINRES()
271 if (ksp->its > 1) { in KSPSolve_MINRES()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/nash/ |
| H A D | nash.c | 51 ksp->its = 0; in KSPCGSolve_NASH()
108 ++ksp->its; in KSPCGSolve_NASH()
143 ++ksp->its; in KSPCGSolve_NASH()
174 PetscCall(KSPMonitor(ksp, ksp->its, norm_r)); in KSPCGSolve_NASH()
177 PetscCall((*ksp->converged)(ksp, ksp->its, norm_r, &ksp->reason, ksp->cnvP)); in KSPCGSolve_NASH()
185 ++ksp->its; in KSPCGSolve_NASH()
221 ++ksp->its; in KSPCGSolve_NASH()
296 ++ksp->its; in KSPCGSolve_NASH()
418 PetscCall(KSPMonitor(ksp, ksp->its, norm_r)); in KSPCGSolve_NASH()
421 PetscCall((*ksp->converged)(ksp, ksp->its, norm_r, &ksp->reason, ksp->cnvP)); in KSPCGSolve_NASH()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/stcg/ |
| H A D | stcg.c | 48 ksp->its = 0; in KSPCGSolve_STCG()
105 ++ksp->its; in KSPCGSolve_STCG()
140 ++ksp->its; in KSPCGSolve_STCG()
171 PetscCall(KSPMonitor(ksp, ksp->its, norm_r)); in KSPCGSolve_STCG()
174 PetscCall((*ksp->converged)(ksp, ksp->its, norm_r, &ksp->reason, ksp->cnvP)); in KSPCGSolve_STCG()
182 ++ksp->its; in KSPCGSolve_STCG()
218 ++ksp->its; in KSPCGSolve_STCG()
294 ++ksp->its; in KSPCGSolve_STCG()
417 PetscCall(KSPMonitor(ksp, ksp->its, norm_r)); in KSPCGSolve_STCG()
420 PetscCall((*ksp->converged)(ksp, ksp->its, norm_r, &ksp->reason, ksp->cnvP)); in KSPCGSolve_STCG()
[all …]
|
| /petsc/src/tao/snes/ |
| H A D | taosnes.c | 11 PetscInt its; in TaoSolve_SNES() local
18 PetscCall(SNESGetIterationNumber(taosnes->snes, &its)); in TaoSolve_SNES()
19 PetscCall(TaoSetIterationNumber(tao, its)); in TaoSolve_SNES()
60 static PetscErrorCode TAOSNESMonitor(SNES snes, PetscInt its, PetscReal fnorm, PetscCtx ctx) in TAOSNESMonitor() argument
69 PetscCall(TaoSetIterationNumber(tao, its)); in TAOSNESMonitor()
70 PetscCall(TaoMonitor(tao, its, obj, fnorm, 0.0, 0.0)); in TAOSNESMonitor()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex16f.F90 | 21 PetscInt its, ione, ifive, izero
67 PetscCallA(KSPGetIterationNumber(ksp1, its, ierr))
69 write (6, 100) norm, its
88 PetscCallA(KSPGetIterationNumber(ksp2, its, ierr))
89 write (6, 100) norm, its
|
| /petsc/src/ksp/ksp/impls/gmres/pgmres/ |
| H A D | pgmres.c | 38 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPPGMRESCycle()
45 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPPGMRESCycle()
54 …if (it < pgmres->max_k + 1 && ksp->its + 1 < PetscMax(2, ksp->max_it)) { /* We don't know whether … in KSPPGMRESCycle()
72 ksp->its++; in KSPPGMRESCycle()
76 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPPGMRESCycle()
80 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPPGMRESCycle()
89 if (!(it < pgmres->max_k + 1 && ksp->its < ksp->max_it)) break; in KSPPGMRESCycle()
156 …if (ksp->reason == KSP_CONVERGED_ITERATING && ksp->its == ksp->max_it) ksp->reason = KSP_DIVERGED_… in KSPPGMRESCycle()
159 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPPGMRESCycle()
166 PetscInt its, itcount; in KSPSolve_PGMRES() local
[all …]
|
| /petsc/src/snes/interface/ftn-custom/ |
| H A D | zsnesf90.c | 15 PetscInt *its, N; in snesgetconvergencehistory_() local
18 *ierr = SNESGetConvergenceHistory(*snes, &hist, &its, &N); in snesgetconvergencehistory_()
22 *ierr = F90Array1dCreate(its, MPIU_INT, 1, N, fits PETSC_F90_2PTR_PARAM(ptrd2)); in snesgetconvergencehistory_()
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_ksp.py | 56 ksp.its = 1
57 self.assertEqual(ksp.its, 1)
58 ksp.its = 0
59 self.assertEqual(ksp.its, 0)
144 def monitor(ksp, its, rnorm): argument
146 reshist[its] = {'r': rnorm, 'o': ksp.getCGObjectiveValue()}
148 reshist[its] = rnorm
165 def converged(ksp, its, rnorm): argument
166 if its > 10:
177 def converged(ksp, its, rnorm): argument
|
| /petsc/src/snes/tests/ |
| H A D | ex12f.F90 | 18 PetscInt :: its, lag
100 subroutine FormMonitor(snes, its, norm, snesm, ierr) argument
103 PetscInt :: its, one, mone
112 if (mod(snesm%its, snesm%lag) == 0) then
119 snesm%its = snesm%its + 1
198 PetscInt ii, its, i1, i0, i3
262 snesm%its = 0
270 PetscCallA(SNESGetIterationNumber(snes, its, ierr))
|
| /petsc/src/dm/dt/fe/tests/ |
| H A D | ex1.c | 10 PetscInt its; /* Number of replications for timing */ member
19 options->its = 1; in ProcessOptions()
25 …sInt("-its", "The number of replications for timing", "ex1.c", options->its, &options->its, NULL)); in ProcessOptions()
180 static PetscErrorCode TestIntegration(DM dm, PetscInt cbs, PetscInt its) in TestIntegration() argument
210 for (i = 0; i < its; ++i) { in TestIntegration()
241 PetscInt N = (cEnd - cStart) * Nf * its; in TestIntegration()
247 …nt_FMT " reps\n Cell rate: %.2f/s flop rate: %.2f MF/s\n", title, N, Nch, its, (double)cellRate, … in TestIntegration()
252 static PetscErrorCode TestIntegration2(DM dm, PetscInt cbs, PetscInt its) in TestIntegration2() argument
263 for (i = 0; i < its; ++i) PetscCall(DMPlexSNESComputeResidualFEM(dm, X, F, NULL)); in TestIntegration2()
304 PetscCall(TestIntegration(dm, ctx.cbs, ctx.its)); in main()
[all …]
|