Home
last modified time | relevance | path

Searched refs:normtype (Results 1 – 25 of 53) sorted by relevance

123

/petsc/src/ksp/ksp/interface/
H A Ditcreate.c152 …if (ksp->converged == KSPConvergedSkip || ksp->normtype == KSP_NORM_NONE) PetscCall(PetscViewerASC… in KSPView()
153 …erASCIIPrintf(viewer, " using %s norm type for convergence test\n", KSPNormTypes[ksp->normtype])); in KSPView()
280 PetscErrorCode KSPSetNormType(KSP ksp, KSPNormType normtype) in KSPSetNormType() argument
284 PetscValidLogicalCollectiveEnum(ksp, normtype, 2); in KSPSetNormType()
285 ksp->normtype = ksp->normtype_set = normtype; in KSPSetNormType()
385 PetscErrorCode KSPSetSupportedNorm(KSP ksp, KSPNormType normtype, PCSide pcside, PetscInt priority) in KSPSetSupportedNorm() argument
389 ksp->normsupporttable[normtype][pcside] = priority; in KSPSetSupportedNorm()
398 ksp->normtype = ksp->normtype_set; in KSPNormSupportTableReset_Private()
402 PetscErrorCode KSPSetUpNorms_Private(KSP ksp, PetscBool errorifnotsupported, KSPNormType *normtype,… in KSPSetUpNorms_Private() argument
410 …if ((ksp->normtype == KSP_NORM_DEFAULT || ksp->normtype == i) && (ksp->pc_side == PC_SIDE_DEFAULT … in KSPSetUpNorms_Private()
[all …]
H A Diterativ.c395 char normtype[256]; in KSPMonitorTrueResidual() local
403 PetscCall(PetscStrncpy(normtype, KSPNormTypes[ksp->normtype], sizeof(normtype))); in KSPMonitorTrueResidual()
404 PetscCall(PetscStrtolower(normtype)); in KSPMonitorTrueResidual()
414 …%s resid norm %14.12e true resid norm %14.12e ||r(i)||/||b|| inf\n", n, normtype, (double)rnorm, (… in KSPMonitorTrueResidual()
416 …esid norm %14.12e true resid norm %14.12e ||r(i)||/||b|| %14.12e\n", n, normtype, (double)rnorm, (… in KSPMonitorTrueResidual()
579 char normtype[256]; in KSPMonitorTrueResidualMax() local
587 PetscCall(PetscStrncpy(normtype, KSPNormTypes[ksp->normtype], sizeof(normtype))); in KSPMonitorTrueResidualMax()
588 PetscCall(PetscStrtolower(normtype)); in KSPMonitorTrueResidualMax()
597 …_FMT " KSP %s true resid max norm %14.12e ||r(i)||/||b|| %14.12e\n", n, normtype, (double)truenorm… in KSPMonitorTrueResidualMax()
1522 KSPNormType normtype; in KSPConvergedDefault() local
[all …]
/petsc/src/ksp/ksp/impls/cr/
H A Dcr.c45 if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_CR()
50 } else if (ksp->normtype == KSP_NORM_NONE) { in KSPSolve_CR()
53 } else if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_CR()
58 } else if (ksp->normtype == KSP_NORM_NATURAL) { in KSPSolve_CR()
61 …bjectComm((PetscObject)ksp), PETSC_ERR_SUP, "KSPNormType of %d not supported", (int)ksp->normtype); in KSPSolve_CR()
96 if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_CR()
101 } else if (ksp->normtype == KSP_NORM_NATURAL) { in KSPSolve_CR()
104 } else if (ksp->normtype == KSP_NORM_NONE) { in KSPSolve_CR()
107 } else if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_CR()
113 …bjectComm((PetscObject)ksp), PETSC_ERR_SUP, "KSPNormType of %d not supported", (int)ksp->normtype); in KSPSolve_CR()
/petsc/src/ksp/ksp/impls/cg/pipecg/
H A Dpipecg.c57 switch (ksp->normtype) { in KSPSolve_PIPECG()
83 SETERRQ(PetscObjectComm((PetscObject)ksp), PETSC_ERR_SUP, "%s", KSPNormTypes[ksp->normtype]); in KSPSolve_PIPECG()
93 if (i > 0 && ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_PIPECG()
95 } else if (i > 0 && ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_PIPECG()
98 if (!(i == 0 && ksp->normtype == KSP_NORM_NATURAL)) PetscCall(VecDotBegin(R, U, &gamma)); in KSPSolve_PIPECG()
105 if (i > 0 && ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_PIPECG()
107 } else if (i > 0 && ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_PIPECG()
110 if (!(i == 0 && ksp->normtype == KSP_NORM_NATURAL)) PetscCall(VecDotEnd(R, U, &gamma)); in KSPSolve_PIPECG()
114 if (ksp->normtype == KSP_NORM_NATURAL) dp = PetscSqrtReal(PetscAbsScalar(gamma)); in KSPSolve_PIPECG()
115 else if (ksp->normtype == KSP_NORM_NONE) dp = 0.0; in KSPSolve_PIPECG()
/petsc/src/ksp/ksp/impls/cg/groppcg/
H A Dgroppcg.c62 switch (ksp->normtype) { in KSPSolve_GROPPCG()
79 SETERRQ(PetscObjectComm((PetscObject)ksp), PETSC_ERR_SUP, "%s", KSPNormTypes[ksp->normtype]); in KSPSolve_GROPPCG()
104 if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_GROPPCG()
106 } else if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_GROPPCG()
114 if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_GROPPCG()
116 } else if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_GROPPCG()
121 if (ksp->normtype == KSP_NORM_NATURAL) { in KSPSolve_GROPPCG()
124 } else if (ksp->normtype == KSP_NORM_NONE) { in KSPSolve_GROPPCG()
/petsc/src/ksp/ksp/impls/cg/
H A Dcg.c168 switch (ksp->normtype) { in KSPSolve_CG()
188 SETERRQ(PetscObjectComm((PetscObject)ksp), PETSC_ERR_SUP, "%s", KSPNormTypes[ksp->normtype]); in KSPSolve_CG()
213 …if (ksp->normtype != KSP_NORM_PRECONDITIONED && (ksp->normtype != KSP_NORM_NATURAL)) PetscCall(KSP… in KSPSolve_CG()
214 if (ksp->normtype != KSP_NORM_NATURAL) { in KSPSolve_CG()
307 if (ksp->normtype == KSP_NORM_PRECONDITIONED && ksp->chknorm < i + 2) { in KSPSolve_CG()
311 } else if (ksp->normtype == KSP_NORM_UNPRECONDITIONED && ksp->chknorm < i + 2) { in KSPSolve_CG()
314 } else if (ksp->normtype == KSP_NORM_NATURAL) { in KSPSolve_CG()
342 …if ((ksp->normtype != KSP_NORM_PRECONDITIONED && ksp->normtype != KSP_NORM_NATURAL) || ksp->chknor… in KSPSolve_CG()
343 if (ksp->normtype != KSP_NORM_NATURAL || ksp->chknorm >= i + 2) { in KSPSolve_CG()
404 switch (ksp->normtype) { in KSPSolve_CG_SingleReduction()
[all …]
H A Dcgls.c56 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = gamma; in KSPSolve_CGLS()
74 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = gamma; in KSPSolve_CGLS()
/petsc/src/ksp/ksp/impls/rich/
H A Drich.c78 if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_Richardson()
80 } else if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_Richardson()
101 if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_Richardson()
103 } else if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_Richardson()
112 …if (ksp->normtype != KSP_NORM_PRECONDITIONED) PetscCall(KSP_PCApply(ksp, r, z)); /* z <- B r … in KSPSolve_Richardson()
117 if (i + 1 < maxit || ksp->normtype != KSP_NORM_NONE) { in KSPSolve_Richardson()
124 if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_Richardson()
126 } else if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_Richardson()
136 if (ksp->normtype != KSP_NORM_NONE) { in KSPSolve_Richardson()
212 if (ksp->normtype == KSP_NORM_NONE) { in KSPBuildResidual_Richardson()
/petsc/src/mat/impls/h2opus/cuda/
H A Dmath2opusutils.cu122 PETSC_INTERN PetscErrorCode MatApproximateNorm_Private(Mat A, NormType normtype, PetscInt normsampl… in MatApproximateNorm_Private() argument
134 switch (normtype) { in MatApproximateNorm_Private()
138 if (normtype == NORM_INFINITY) { in MatApproximateNorm_Private()
161 …PetscCall(PetscInfo(A, "%s norm it %" PetscInt_FMT " -> (%g %g)\n", NormTypes[normtype], i, (doubl… in MatApproximateNorm_Private()
169 …all(PetscInfo(A, "%s norm it %" PetscInt_FMT " -> breakdown (j==jold)\n", NormTypes[normtype], i)); in MatApproximateNorm_Private()
198 …PetscCall(PetscInfo(A, "%s norm it %" PetscInt_FMT " -> %g (%g %g)\n", NormTypes[normtype], i, (do… in MatApproximateNorm_Private()
214 …ERRQ(PetscObjectComm((PetscObject)A), PETSC_ERR_SUP, "%s norm not supported", NormTypes[normtype]); in MatApproximateNorm_Private()
216 …o(A, "%s norm %g computed in %" PetscInt_FMT " iterations\n", NormTypes[normtype], (double)*n, i)); in MatApproximateNorm_Private()
/petsc/src/ksp/ksp/impls/cg/cgne/
H A Dcgne.c91 if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_CGNE()
93 } else if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_CGNE()
95 } else if (ksp->normtype == KSP_NORM_NATURAL) { in KSPSolve_CGNE()
142 if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_CGNE()
150 } else if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_CGNE()
152 } else if (ksp->normtype == KSP_NORM_NATURAL) { in KSPSolve_CGNE()
160 if (ksp->normtype != KSP_NORM_PRECONDITIONED) { in KSPSolve_CGNE()
/petsc/src/ksp/ksp/impls/cheby/
H A Dcheby.c439 if (ksp->normtype) { in KSPSolve_Chebyshev_FirstKind()
440 switch (ksp->normtype) { in KSPSolve_Chebyshev_FirstKind()
450 SETERRQ(PetscObjectComm((PetscObject)ksp), PETSC_ERR_SUP, "%s", KSPNormTypes[ksp->normtype]); in KSPSolve_Chebyshev_FirstKind()
464 …if (ksp->normtype != KSP_NORM_PRECONDITIONED) PetscCall(KSP_PCApply(ksp, r, p[k])); /* p[k] = B^{-… in KSPSolve_Chebyshev_FirstKind()
478 if (ksp->normtype) { in KSPSolve_Chebyshev_FirstKind()
479 switch (ksp->normtype) { in KSPSolve_Chebyshev_FirstKind()
500 …if (ksp->normtype != KSP_NORM_PRECONDITIONED) PetscCall(KSP_PCApply(ksp, r, p[kp1])); /* p[kp1] =… in KSPSolve_Chebyshev_FirstKind()
519 if (ksp->normtype) { in KSPSolve_Chebyshev_FirstKind()
522 switch (ksp->normtype) { in KSPSolve_Chebyshev_FirstKind()
543 if (ksp->normtype != KSP_NORM_NONE) { in KSPSolve_Chebyshev_FirstKind()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecgrr/
H A Dpipecgrr.c58 switch (ksp->normtype) { in KSPSolve_PIPECGRR()
90 SETERRQ(PetscObjectComm((PetscObject)ksp), PETSC_ERR_SUP, "%s", KSPNormTypes[ksp->normtype]); in KSPSolve_PIPECGRR()
119 if (i > 0 && ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_PIPECGRR()
121 } else if (i > 0 && ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_PIPECGRR()
124 if (!(i == 0 && ksp->normtype == KSP_NORM_NATURAL)) PetscCall(VecDotBegin(R, U, &gamma)); in KSPSolve_PIPECGRR()
142 if (i > 0 && ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_PIPECGRR()
144 } else if (i > 0 && ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_PIPECGRR()
147 if (!(i == 0 && ksp->normtype == KSP_NORM_NATURAL)) PetscCall(VecDotEnd(R, U, &gamma)); in KSPSolve_PIPECGRR()
162 if (ksp->normtype == KSP_NORM_NATURAL) dp = PetscSqrtReal(PetscAbsScalar(gamma)); in KSPSolve_PIPECGRR()
163 else if (ksp->normtype == KSP_NORM_NONE) dp = 0.0; in KSPSolve_PIPECGRR()
/petsc/src/ksp/ksp/impls/bicg/
H A Dbicg.c43 if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_BiCG()
45 } else if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_BiCG()
87 if (ksp->normtype == KSP_NORM_PRECONDITIONED) { in KSPSolve_BiCG()
91 } else if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_BiCG()
104 if (ksp->normtype == KSP_NORM_UNPRECONDITIONED) { in KSPSolve_BiCG()
/petsc/src/ksp/ksp/impls/cr/pipecr/
H A Dpipecr.c55 switch (ksp->normtype) { in KSPSolve_PIPECR()
67 SETERRQ(PetscObjectComm((PetscObject)ksp), PETSC_ERR_SUP, "%s", KSPNormTypes[ksp->normtype]); in KSPSolve_PIPECR()
79 if (i > 0 && ksp->normtype == KSP_NORM_PRECONDITIONED) PetscCall(VecNormBegin(U, NORM_2, &dp)); in KSPSolve_PIPECR()
86 if (i > 0 && ksp->normtype == KSP_NORM_PRECONDITIONED) PetscCall(VecNormEnd(U, NORM_2, &dp)); in KSPSolve_PIPECR()
91 if (ksp->normtype == KSP_NORM_NONE) dp = 0.0; in KSPSolve_PIPECR()
/petsc/src/ksp/ksp/impls/cgs/
H A Dcgs.c44 if (ksp->normtype != KSP_NORM_NONE) { in KSPSolve_CGS()
47 if (ksp->normtype == KSP_NORM_NATURAL) dp *= dp; in KSPSolve_CGS()
66 if (ksp->normtype == KSP_NORM_NATURAL) { in KSPSolve_CGS()
98 if (ksp->normtype == KSP_NORM_NATURAL) { in KSPSolve_CGS()
100 } else if (ksp->normtype != KSP_NORM_NONE) { in KSPSolve_CGS()
/petsc/src/ksp/ksp/impls/fcg/
H A Dfcg.c107 switch (ksp->normtype) { in KSPSolve_FCG()
127 SETERRQ(PetscObjectComm((PetscObject)ksp), PETSC_ERR_SUP, "%s", KSPNormTypes[ksp->normtype]); in KSPSolve_FCG()
134 if (ksp->normtype == KSP_NORM_NONE) { in KSPSolve_FCG()
142 …if (ksp->normtype == KSP_NORM_UNPRECONDITIONED || ksp->normtype == KSP_NORM_NONE) PetscCall(KSP_PC… in KSPSolve_FCG()
219 switch (ksp->normtype) { in KSPSolve_FCG()
239 SETERRQ(PetscObjectComm((PetscObject)ksp), PETSC_ERR_SUP, "%s", KSPNormTypes[ksp->normtype]); in KSPSolve_FCG()
260 …if (ksp->normtype == KSP_NORM_UNPRECONDITIONED || ksp->normtype == KSP_NORM_NONE) PetscCall(KSP_PC… in KSPSolve_FCG()
/petsc/src/ksp/ksp/impls/cg/pipecg2/
H A Dpipecg2.c8 static PetscErrorCode VecMergedDot_Private(Vec U, Vec W, Vec R, PetscInt normtype, PetscScalar *ru,… in VecMergedDot_Private() argument
20 if (normtype == KSP_NORM_PRECONDITIONED) { in VecMergedDot_Private()
27 } else if (normtype == KSP_NORM_UNPRECONDITIONED) { in VecMergedDot_Private()
34 } else if (normtype == KSP_NORM_NATURAL) { in VecMergedDot_Private()
82 … vs, Vec va1, Vec vb1, Vec ve, Vec vf, Vec vm, Vec vn, Vec vu, PetscInt normtype, PetscScalar beta… in VecMergedOpsShort_Private() argument
114 if (normtype == KSP_NORM_PRECONDITIONED) { in VecMergedOpsShort_Private()
169 } else if (normtype == KSP_NORM_UNPRECONDITIONED) { in VecMergedOpsShort_Private()
224 } else if (normtype == KSP_NORM_NATURAL) { in VecMergedOpsShort_Private()
304 … vs, Vec va1, Vec vb1, Vec ve, Vec vf, Vec vm, Vec vn, Vec vu, PetscInt normtype, PetscScalar beta… in VecMergedOps_Private() argument
336 if (normtype == KSP_NORM_PRECONDITIONED) { in VecMergedOps_Private()
[all …]
/petsc/src/ksp/ksp/impls/tcqmr/
H A Dtcqmr.c23 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rnorm0; in KSPSolve_TCQMR()
54 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rnorm; in KSPSolve_TCQMR()
139 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rnorm; in KSPSolve_TCQMR()
/petsc/src/ksp/ksp/impls/ibcgs/
H A Dibcgs.c116 if (ksp->normtype != KSP_NORM_NONE) { 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()
276 if (!ksp->lagnorm && ksp->chknorm < ksp->its && ksp->normtype != KSP_NORM_NONE) { in KSPSolve_IBCGS()
/petsc/src/ksp/ksp/impls/bcgs/pipebcgs/
H A Dpipebcgs.c61 if (ksp->normtype != KSP_NORM_NONE) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_PIPEBCGS()
144 …if (ksp->normtype != KSP_NORM_NONE && ksp->chknorm < i + 2) PetscCall(VecNormBegin(R, NORM_2, &dp)… in KSPSolve_PIPEBCGS()
154 … if (ksp->normtype != KSP_NORM_NONE && ksp->chknorm < i + 2) PetscCall(VecNormEnd(R, NORM_2, &dp)); in KSPSolve_PIPEBCGS()
/petsc/src/snes/impls/composite/
H A Dsnescomposite.c693 SNESNormSchedule normtype; in SNESSolve_Composite() local
708 PetscCall(SNESGetNormSchedule(snes, &normtype)); in SNESSolve_Composite()
709 …if (normtype == SNES_NORM_ALWAYS || normtype == SNES_NORM_INITIAL_ONLY || normtype == SNES_NORM_IN… in SNESSolve_Composite()
755 …if ((i == snes->max_its - 1) && (normtype == SNES_NORM_INITIAL_FINAL_ONLY || normtype == SNES_NORM… in SNESSolve_Composite()
774 } else if (normtype == SNES_NORM_ALWAYS) { in SNESSolve_Composite()
/petsc/src/ksp/ksp/impls/lsqr/
H A Dlsqr.c382 char normtype[256]; in KSPLSQRMonitorResidual_LSQR() local
389 PetscCall(PetscStrncpy(normtype, KSPNormTypes[ksp->normtype], sizeof(normtype))); in KSPLSQRMonitorResidual_LSQR()
390 PetscCall(PetscStrtolower(normtype)); in KSPLSQRMonitorResidual_LSQR()
/petsc/src/ksp/ksp/impls/tfqmr/
H A Dtfqmr.c39 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = dp; in KSPSolve_TFQMR()
94 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = dpest; in KSPSolve_TFQMR()
/petsc/src/ksp/ksp/impls/bcgs/fbcgs/
H A Dfbcgs.c56 if (ksp->normtype != KSP_NORM_NONE) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_FBCGS()
121 if (ksp->normtype != KSP_NORM_NONE && ksp->chknorm < i + 2) PetscCall(VecNorm(R, NORM_2, &dp)); in KSPSolve_FBCGS()
/petsc/src/ksp/ksp/impls/lcd/
H A Dlcd.c61 if (ksp->normtype != KSP_NORM_NONE) { in KSPSolve_LCD()
88 if (ksp->normtype != KSP_NORM_NONE) { in KSPSolve_LCD()

123