Home
last modified time | relevance | path

Searched refs:rnorm (Results 1 – 25 of 101) sorted by relevance

12345

/petsc/src/ksp/ksp/impls/rich/
H A Drich.c22 PetscReal rnorm = 0.0, abr; in KSPSolve_Richardson() local
79 PetscCall(VecNorm(r, NORM_2, &rnorm)); /* rnorm <- r'*r */ in KSPSolve_Richardson()
81 PetscCall(VecNorm(z, NORM_2, &rnorm)); /* rnorm <- z'*z */ in KSPSolve_Richardson()
82 } else rnorm = 0.0; in KSPSolve_Richardson()
84 KSPCheckNorm(ksp, rnorm); in KSPSolve_Richardson()
85 ksp->rnorm = rnorm; in KSPSolve_Richardson()
86 PetscCall(KSPMonitor(ksp, i, rnorm)); in KSPSolve_Richardson()
87 PetscCall(KSPLogResidualHistory(ksp, rnorm)); in KSPSolve_Richardson()
88 PetscCall((*ksp->converged)(ksp, i, rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_Richardson()
102 PetscCall(VecNorm(r, NORM_2, &rnorm)); /* rnorm <- r'*r */ in KSPSolve_Richardson()
[all …]
/petsc/src/mat/tests/
H A Dex53.c12 PetscReal s1norm, s2norm, rnorm, tol = 1.e-10; in main() local
50 rnorm = PetscAbsScalar(s2norm - s1norm) / s2norm; in main()
51 …if (rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Error: MatNorm_FROBENIUS()- NormA=%1… in main()
54 rnorm = PetscAbsScalar(s2norm - s1norm) / s2norm; in main()
55 …if (rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Error: MatNorm_INFINITY()- NormA=%16… in main()
58 rnorm = PetscAbsScalar(s2norm - s1norm) / s2norm; in main()
59 …if (rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Error: MatNorm_NORM_1()- NormA=%16.1… in main()
67 PetscCall(VecNorm(s2, NORM_2, &rnorm)); in main()
68 …if (rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Error: MatMult - Norm2=%16.14e bs = … in main()
78 PetscCall(VecNorm(s2, NORM_2, &rnorm)); in main()
[all …]
H A Dex48.c13 PetscReal s1norm, s2norm, rnorm, tol = 1.e-4; in main() local
70 rnorm = PetscAbsReal(s2norm - s1norm) / s2norm; in main()
71 …if (rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "Error: MatNorm_FROBENIUS()- NormA=%16.14e… in main()
74 rnorm = PetscAbsReal(s2norm - s1norm) / s2norm; in main()
75 …if (rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "Error: MatNorm_INFINITY()- NormA=%16.14e … in main()
78 rnorm = PetscAbsReal(s2norm - s1norm) / s2norm; in main()
79 …if (rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "Error: MatNorm_NORM_1()- NormA=%16.14e No… in main()
114 rnorm = s2norm - s1norm; in main()
115 …if (rnorm < -tol || rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "MatMult not equal to MatM… in main()
172 PetscCall(VecNorm(yy, NORM_2, &rnorm)); in main()
[all …]
H A Dex91.c14 PetscReal s1norm, s2norm, rnorm, tol = 10 * PETSC_SMALL; in main() local
94 rnorm = s2norm - s1norm; in main()
95 …if (rnorm < -tol || rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "Error:MatMult - Norm1=%16… in main()
151 rnorm = s2norm - s1norm; in main()
152 …if (rnorm < -tol || rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "Error:MatMult - Norm1=%16… in main()
175 rnorm = s2norm - s1norm; in main()
176 …if (rnorm < -tol || rnorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, "Error:MatMult - Norm1=%16… in main()
/petsc/src/ksp/ksp/impls/tcqmr/
H A Dtcqmr.c12 PetscReal rnorm0, rnorm, dp1, Gamma; in KSPSolve_TCQMR() local
23 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rnorm0; in KSPSolve_TCQMR()
24 else ksp->rnorm = 0; in KSPSolve_TCQMR()
25 PetscCall((*ksp->converged)(ksp, 0, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_TCQMR()
30 rnorm = rnorm0; in KSPSolve_TCQMR()
31 tmp = 1.0 / rnorm; in KSPSolve_TCQMR()
54 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rnorm; in KSPSolve_TCQMR()
55 else ksp->rnorm = 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()
[all …]
/petsc/src/ksp/ksp/impls/cheby/
H A Dcheby.c393 PetscReal rnorm = 0.0, emax, emin; in KSPSolve_Chebyshev_FirstKind() local
443 PetscCall(VecNorm(p[k], NORM_2, &rnorm)); in KSPSolve_Chebyshev_FirstKind()
447 PetscCall(VecNorm(r, NORM_2, &rnorm)); in KSPSolve_Chebyshev_FirstKind()
453 ksp->rnorm = rnorm; in KSPSolve_Chebyshev_FirstKind()
455 PetscCall(KSPLogResidualHistory(ksp, rnorm)); in KSPSolve_Chebyshev_FirstKind()
457 PetscCall(KSPMonitor(ksp, 0, rnorm)); in KSPSolve_Chebyshev_FirstKind()
458 PetscCall((*ksp->converged)(ksp, 0, rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_Chebyshev_FirstKind()
482 PetscCall(VecNorm(p[kp1], NORM_2, &rnorm)); in KSPSolve_Chebyshev_FirstKind()
486 PetscCall(VecNorm(r, NORM_2, &rnorm)); in KSPSolve_Chebyshev_FirstKind()
489 rnorm = 0.0; in KSPSolve_Chebyshev_FirstKind()
[all …]
/petsc/src/ksp/ksp/impls/lsqr/
H A Dlsqr.c72 PetscReal beta, alpha, rnorm; in KSPSolve_LSQR() local
110 PetscCall(VecNorm(U, NORM_2, &rnorm)); in KSPSolve_LSQR()
111 KSPCheckNorm(ksp, rnorm); in KSPSolve_LSQR()
114 ksp->rnorm = rnorm; in KSPSolve_LSQR()
116 PetscCall(KSPLogResidualHistory(ksp, rnorm)); in KSPSolve_LSQR()
117 PetscCall(KSPMonitor(ksp, 0, rnorm)); in KSPSolve_LSQR()
118 PetscCall((*ksp->converged)(ksp, 0, rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_LSQR()
121 beta = rnorm; in KSPSolve_LSQR()
126 KSPCheckNorm(ksp, rnorm); in KSPSolve_LSQR()
215 rnorm = PetscRealPart(phibar); in KSPSolve_LSQR()
[all …]
/petsc/src/ksp/ksp/impls/lcd/
H A Dlcd.c32 PetscReal rnorm = 0.0; in KSPSolve_LCD() local
62 PetscCall(VecNorm(R, NORM_2, &rnorm)); in KSPSolve_LCD()
63 KSPCheckNorm(ksp, rnorm); in KSPSolve_LCD()
65 PetscCall(KSPLogResidualHistory(ksp, rnorm)); in KSPSolve_LCD()
66 PetscCall(KSPMonitor(ksp, 0, rnorm)); in KSPSolve_LCD()
67 ksp->rnorm = rnorm; in KSPSolve_LCD()
70 PetscCall((*ksp->converged)(ksp, 0, rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_LCD()
89 PetscCall(VecNorm(R, NORM_2, &rnorm)); in KSPSolve_LCD()
90 KSPCheckNorm(ksp, rnorm); in KSPSolve_LCD()
93 ksp->rnorm = rnorm; in KSPSolve_LCD()
[all …]
/petsc/src/mat/tutorials/
H A Dex9.c16 PetscReal rnorm; in main() local
65 PetscCall(VecNorm(y, NORM_2, &rnorm)); in main()
66 …if (rnorm > 10000.0 * PETSC_MACHINE_EPSILON) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Error with c… in main()
74 PetscCall(VecNorm(y, NORM_2, &rnorm)); in main()
75 …if (rnorm > 10000.0 * PETSC_MACHINE_EPSILON) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Error with c… in main()
104 PetscCall(VecNorm(y, NORM_2, &rnorm)); in main()
105 …if (rnorm > 10000.0 * PETSC_MACHINE_EPSILON) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Error with c… in main()
127 PetscCall(VecNorm(v2, NORM_2, &rnorm)); in main()
128 …if (rnorm > 10000.0 * PETSC_MACHINE_EPSILON) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Error with c… in main()
/petsc/src/ksp/ksp/interface/
H A Dxmon.c4 PetscErrorCode KSPMonitorLGRange(KSP ksp, PetscInt n, PetscReal rnorm, void *monctx) in KSPMonitorLGRange() argument
16 if (!n) prev = rnorm; in KSPMonitorLGRange()
23 if (rnorm > 0.0) y = PetscLog10Real(rnorm); in KSPMonitorLGRange()
48 y = (prev - rnorm) / prev; in KSPMonitorLGRange()
60 y = (prev - rnorm) / (prev * per); in KSPMonitorLGRange()
66 prev = rnorm; in KSPMonitorLGRange()
H A Diterativ.c34 PetscErrorCode KSPGetResidualNorm(KSP ksp, PetscReal *rnorm) in KSPGetResidualNorm() argument
38 PetscAssertPointer(rnorm, 2); in KSPGetResidualNorm()
39 *rnorm = ksp->rnorm; in KSPGetResidualNorm()
127 PetscErrorCode KSPMonitorResidual(KSP ksp, PetscInt n, PetscReal rnorm, PetscViewerAndFormat *vf) in KSPMonitorResidual() argument
141 …etscViewerASCIIPrintf(viewer, "%3" PetscInt_FMT " KSP Residual norm %14.12e\n", n, (double)rnorm)); in KSPMonitorResidual()
169 PetscErrorCode KSPMonitorResidualView(KSP ksp, PetscInt n, PetscReal rnorm, PetscViewerAndFormat *v… in KSPMonitorResidualView() argument
212 PetscErrorCode KSPMonitorResidualDrawLG(KSP ksp, PetscInt n, PetscReal rnorm, PetscViewerAndFormat … in KSPMonitorResidualDrawLG() argument
226 if (rnorm > 0.0) y = PetscLog10Real(rnorm); in KSPMonitorResidualDrawLG()
338 PetscErrorCode KSPMonitorResidualRange(KSP ksp, PetscInt it, PetscReal rnorm, PetscViewerAndFormat … in KSPMonitorResidualRange() argument
353 if (!it) prev = rnorm; in KSPMonitorResidualRange()
[all …]
/petsc/src/ksp/ksp/impls/tfqmr/
H A Dtfqmr.c39 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = dp; in KSPSolve_TFQMR()
40 else ksp->rnorm = 0.0; in KSPSolve_TFQMR()
43 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPSolve_TFQMR()
44 PetscCall(KSPMonitor(ksp, 0, ksp->rnorm)); in KSPSolve_TFQMR()
45 PetscCall((*ksp->converged)(ksp, 0, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_TFQMR()
94 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = dpest; in KSPSolve_TFQMR()
95 else ksp->rnorm = 0.0; in KSPSolve_TFQMR()
97 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPSolve_TFQMR()
98 PetscCall(KSPMonitor(ksp, i + 1, ksp->rnorm)); in KSPSolve_TFQMR()
99 PetscCall((*ksp->converged)(ksp, i + 1, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_TFQMR()
/petsc/src/binding/petsc4py/src/lib-petsc/
H A Dcustom.h285 PetscErrorCode KSPSetResidualNorm(KSP ksp, PetscReal rnorm) in KSPSetResidualNorm() argument
289 …PetscCheck(rnorm >= 0,PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"residual norm must be nonnegative"… in KSPSetResidualNorm()
290 ksp->rnorm = rnorm; in KSPSetResidualNorm()
295 PetscErrorCode KSPConvergenceTestCall(KSP ksp, PetscInt its, PetscReal rnorm, KSPConvergedReason *r… in KSPConvergenceTestCall() argument
301 …PetscCheck(rnorm >= 0,PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"residual norm must be nonnegative"… in KSPConvergenceTestCall()
302 PetscCall((*ksp->converged)(ksp,its,rnorm,reason,ksp->cnvP)); in KSPConvergenceTestCall()
316 PetscErrorCode KSPConverged(KSP ksp,PetscInt iter,PetscReal rnorm,KSPConvergedReason *reason) in KSPConverged() argument
321 if (!iter) ksp->rnorm0 = rnorm; in KSPConverged()
324 ksp->ttol = PetscMax(rnorm*ksp->rtol,ksp->abstol); in KSPConverged()
327 PetscCall(ksp->converged(ksp,iter,rnorm,&ksp->reason,ksp->cnvP)); in KSPConverged()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/fbcgsr/
H A Dfbcgsr.c86 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rho; in KSPSolve_FBCGSR()
87 else ksp->rnorm = 0; in KSPSolve_FBCGSR()
89 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPSolve_FBCGSR()
90 PetscCall(KSPMonitor(ksp, 0, ksp->rnorm)); in KSPSolve_FBCGSR()
91 PetscCall((*ksp->converged)(ksp, 0, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_FBCGSR()
180 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rho; in KSPSolve_FBCGSR()
181 else ksp->rnorm = 0; in KSPSolve_FBCGSR()
183 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPSolve_FBCGSR()
184 PetscCall(KSPMonitor(ksp, i + 1, ksp->rnorm)); in KSPSolve_FBCGSR()
185 PetscCall((*ksp->converged)(ksp, i + 1, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_FBCGSR()
/petsc/src/ksp/ksp/impls/cg/
H A Dcgls.c49 ksp->rnorm = 0; in KSPSolve_CGLS()
56 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = gamma; in KSPSolve_CGLS()
57 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_CGLS()
74 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = gamma; 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()
/petsc/src/ksp/ksp/impls/minres/
H A Dminres.c105 …PetscReal Acond = 1.0, Acondl = 0.0, rnorml, rnorm, rootl, relAresl, relres, relresl, Arnorml, … in KSPSolve_MINRES() local
141 ksp->rnorm = 0.0; in KSPSolve_MINRES()
163 rnorm = beta1; in KSPSolve_MINRES()
164 if (ksp->normtype == KSP_NORM_PRECONDITIONED) ksp->rnorm = rnorm; in KSPSolve_MINRES()
165 else if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) PetscCall(VecNorm(R2, NORM_2, &ksp->rnorm)); in KSPSolve_MINRES()
166 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPSolve_MINRES()
167 PetscCall(KSPMonitor(ksp, 0, ksp->rnorm)); in KSPSolve_MINRES()
168 …PetscCall((*ksp->converged)(ksp, 0, ksp->rnorm, &ksp->reason, ksp->cnvP)); /* test for convergence… in KSPSolve_MINRES()
171 relres = rnorm / beta1; in KSPSolve_MINRES()
393 rnorml = rnorm; in KSPSolve_MINRES()
[all …]
/petsc/src/binding/petsc4py/test/
H A Dtest_ksp_py.py31 rnorm = r.norm()
32 ksp.setResidualNorm(rnorm)
33 ksp.logConvergenceHistory(rnorm)
34 ksp.monitor(its, rnorm)
35 reason = ksp.callConvergenceTest(its, rnorm)
/petsc/src/ksp/ksp/impls/gmres/pgmres/
H A Dpgmres.c33 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPPGMRESCycle()
34 else ksp->rnorm = 0; in KSPPGMRESCycle()
37 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPPGMRESCycle()
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()
73 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPPGMRESCycle()
74 else ksp->rnorm = 0; in KSPPGMRESCycle()
76 PetscCall((*ksp->converged)(ksp, ksp->its, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPPGMRESCycle()
79 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPPGMRESCycle()
80 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPPGMRESCycle()
[all …]
/petsc/src/ksp/ksp/impls/bcgsl/
H A Dbcgsl.c53 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = zeta0; in KSPSolve_BCGSL()
54 else ksp->rnorm = 0.0; in KSPSolve_BCGSL()
56 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPSolve_BCGSL()
57 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPSolve_BCGSL()
58 PetscCall((*ksp->converged)(ksp, 0, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_BCGSL()
82 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = zeta; in KSPSolve_BCGSL()
83 else ksp->rnorm = 0.0; in KSPSolve_BCGSL()
85 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPSolve_BCGSL()
86 PetscCall(KSPMonitor(ksp, ksp->its, ksp->rnorm)); in KSPSolve_BCGSL()
88 PetscCall((*ksp->converged)(ksp, k, ksp->rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_BCGSL()
[all …]
/petsc/src/ksp/ksp/impls/ibcgs/
H A Dibcgs.c42 PetscReal rnorm = 0.0, rnormin = 0.0; in KSPSolve_IBCGS() local
117 PetscCall(VecNorm(Rn_1, NORM_2, &rnorm)); in KSPSolve_IBCGS()
118 KSPCheckNorm(ksp, rnorm); in KSPSolve_IBCGS()
120 PetscCall(KSPMonitor(ksp, 0, rnorm)); in KSPSolve_IBCGS()
121 PetscCall((*ksp->converged)(ksp, 0, rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_IBCGS()
246 …if (ksp->lagnorm && ksp->its > 1 && ksp->normtype != KSP_NORM_NONE) rnorm = PetscSqrtReal(PetscRea… in KSPSolve_IBCGS()
278 …PetscCallMPI(MPIU_Allreduce(&rnormin, &rnorm, 1, MPIU_REAL, MPIU_SUM, PetscObjectComm((PetscObject… in KSPSolve_IBCGS()
280 rnorm = PetscSqrtReal(rnorm); in KSPSolve_IBCGS()
284 PetscCall(KSPMonitor(ksp, ksp->its, rnorm)); in KSPSolve_IBCGS()
285 PetscCall((*ksp->converged)(ksp, ksp->its, rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_IBCGS()
/petsc/src/ksp/ksp/impls/symmlq/
H A Dsymmlq.c61 ksp->rnorm = 0.0; /* what should we really put here? */ in KSPSolve_SYMMLQ()
89 ksp->rnorm = np; in KSPSolve_SYMMLQ()
165 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = np; in KSPSolve_SYMMLQ()
166 else ksp->rnorm = 0.0; in KSPSolve_SYMMLQ()
167 PetscCall(KSPLogResidualHistory(ksp, ksp->rnorm)); in KSPSolve_SYMMLQ()
168 PetscCall(KSPMonitor(ksp, i + 1, ksp->rnorm)); in KSPSolve_SYMMLQ()
169 …PetscCall((*ksp->converged)(ksp, i + 1, ksp->rnorm, &ksp->reason, ksp->cnvP)); /* test for converg… in KSPSolve_SYMMLQ()
/petsc/src/ksp/ksp/impls/gcr/pipegcr/
H A Dpipegcr.c63 PetscReal rnorm = 0.0, delta, *eta, *etas; in KSPSolve_PIPEGCR_cycle() local
130 …tscCall((*pipegcr->modifypc)(ksp, ksp->its, 0 /* unused argument */, ksp->rnorm, pipegcr->modifypc… in KSPSolve_PIPEGCR_cycle()
189 PetscCall(VecNorm(z, NORM_2, &rnorm)); /* ||r|| <- sqrt(z'*z) */ in KSPSolve_PIPEGCR_cycle()
192 PetscCall(VecNorm(r, NORM_2, &rnorm)); /* ||r|| <- sqrt(r'*r) */ in KSPSolve_PIPEGCR_cycle()
195 rnorm = PetscSqrtReal(PetscAbsScalar(gamma)); /* ||r|| <- sqrt(r,w) */ in KSPSolve_PIPEGCR_cycle()
198 rnorm = 0.0; in KSPSolve_PIPEGCR_cycle()
206 ksp->rnorm = rnorm; in KSPSolve_PIPEGCR_cycle()
208 PetscCall(KSPLogResidualHistory(ksp, rnorm)); in KSPSolve_PIPEGCR_cycle()
209 PetscCall(KSPMonitor(ksp, ksp->its, rnorm)); in KSPSolve_PIPEGCR_cycle()
210 PetscCall((*ksp->converged)(ksp, ksp->its, rnorm, &ksp->reason, ksp->cnvP)); in KSPSolve_PIPEGCR_cycle()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex57f.F9024 subroutine MyKSPMonitor(ksp, n, rnorm, unused, ierr) argument
31 PetscReal rnorm
45 if (rank == 0) write (6, 200) n, rnorm
63 subroutine MyKSPConverged(ksp, n, rnorm, flag, unused, ierr) argument
69 PetscReal rnorm
71 if (rnorm <= .05) then
H A Dexample100.py143 rnorm = r.norm()
144 ksp.setResidualNorm(rnorm)
145 ksp.logConvergenceHistory(rnorm)
146 ksp.monitor(its, rnorm)
147 reason = ksp.callConvergenceTest(its, rnorm)
H A Dex2f.F9024 subroutine MyKSPMonitor(ksp, n, rnorm, vf, ierr) argument
32 PetscReal rnorm
36 PetscCallA(KSPMonitorTrueResidual(ksp, n, rnorm, vf, ierr))
45 if (rank == 0) write (6, 200) n, rnorm
63 subroutine MyKSPConverged(ksp, n, rnorm, flag, unused, ierr) argument
69 PetscReal rnorm
71 if (rnorm <= .05) then

12345