| /petsc/src/ts/adapt/impls/basic/ |
| H A D | adaptbasic.c | 9 PetscReal enorm = -1; in TSAdaptChoose_Basic() local
20 PetscCall(TSEvaluateWLTE(ts, adapt->wnormtype, &order, &enorm)); in TSAdaptChoose_Basic()
21 …PetscCheck(enorm < 0 || order >= 1, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OUTOFRANGE,… in TSAdaptChoose_Basic()
29 PetscCall(TSErrorWeightedNorm(ts, ts->vec_sol, Y, adapt->wnormtype, &enorm, &enorma, &enormr)); in TSAdaptChoose_Basic()
33 if (enorm < 0) { in TSAdaptChoose_Basic()
41 if (enorm > 1) { in TSAdaptChoose_Basic()
44 …l truncation error %g, accepting because step size %g is at minimum\n", (double)enorm, (double)h)); in TSAdaptChoose_Basic()
47 …on error %g, accepting step of size %g because always_accept is set\n", (double)enorm, (double)h)); in TSAdaptChoose_Basic()
50 …timated scaled local truncation error %g, rejecting step of size %g\n", (double)enorm, (double)h)); in TSAdaptChoose_Basic()
54 …timated scaled local truncation error %g, accepting step of size %g\n", (double)enorm, (double)h)); in TSAdaptChoose_Basic()
[all …]
|
| /petsc/src/ts/adapt/impls/glee/ |
| H A D | adaptglee.c | 12 PetscReal enorm, enorma, enormr, hfac_lte, hfac_ltea, hfac_lter, h_lte, safety; in TSAdaptChoose_GLEE() local
37 PetscCall(TSErrorWeightedENorm(ts, E, X, X, adapt->wnormtype, &enorm, &enorma, &enormr)); in TSAdaptChoose_GLEE()
45 PetscCall(TSErrorWeightedNorm(ts, X, Y, adapt->wnormtype, &enorm, &enorma, &enormr)); in TSAdaptChoose_GLEE()
48 if (enorm < 0) { in TSAdaptChoose_GLEE()
57 if (enorm > 1. || enorma > 1. || enormr > 1.) { in TSAdaptChoose_GLEE()
60 …]] [%g, %g, %g], accepting because step size %g is at minimum\n", (double)enorm, (double)enorma, (… in TSAdaptChoose_GLEE()
63 …, %g], accepting step of size %g because always_accept is set\n", (double)enorm, (double)enorma, (… in TSAdaptChoose_GLEE()
66 … absolute, relative]] [%g, %g, %g], rejecting step of size %g\n", (double)enorm, (double)enorma, (… in TSAdaptChoose_GLEE()
70 …, absolute, relative] [%g, %g, %g], accepting step of size %g\n", (double)enorm, (double)enorma, (… in TSAdaptChoose_GLEE()
82 if (enorm > 0) { in TSAdaptChoose_GLEE()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex112.c | 26 PetscReal enorm, tol = PETSC_SMALL; in main() local
83 PetscCall(VecNorm(z, NORM_1, &enorm)); in main()
84 …if (enorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %g\n", (double)en… in main()
101 PetscCall(VecNorm(z, NORM_1, &enorm)); in main()
102 …if (enorm > tol) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of new |x - z| %g\n", (doubl… in main()
|
| H A D | ex114.c | 15 PetscReal enorm; in main() local
151 PetscCall(VecNorm(e, NORM_INFINITY, &enorm)); in main()
152 …PetscCheck(enorm <= PETSC_MACHINE_EPSILON, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "max+min > PETSC_MACH… in main()
182 PetscCall(VecNorm(e, NORM_INFINITY, &enorm)); in main()
183 …PetscCheck(enorm <= PETSC_MACHINE_EPSILON, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "norm(-max + max_d) %… in main()
190 PetscCall(VecNorm(e, NORM_INFINITY, &enorm)); in main()
191 …PetscCheck(enorm <= PETSC_MACHINE_EPSILON, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "max+min > PETSC_MACH… in main()
197 PetscCall(VecNorm(e, NORM_INFINITY, &enorm)); in main()
198 …tscCheck(enorm <= PETSC_MACHINE_EPSILON, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "norm(-maxabs + maxabs_… in main()
222 PetscCall(VecNorm(e, NORM_INFINITY, &enorm)); in main()
[all …]
|
| H A D | ex143.c | 22 PetscReal enorm; in main() local
80 PetscCall(VecNorm(z, NORM_1, &enorm)); in main()
81 …if (enorm > 1.e-11 && rank == 0) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %… in main()
134 PetscCall(VecNorm(z, NORM_1, &enorm)); in main()
135 …if (enorm > 1.e-9 && rank == 0) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %e… in main()
|
| H A D | ex158.c | 18 PetscReal enorm; in main() local
77 PetscCall(VecNorm(z, NORM_1, &enorm)); in main()
78 …if (enorm > 1.e-11) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %g\n", (double… in main()
144 PetscCall(VecNorm(output, NORM_1, &enorm)); in main()
145 …if (enorm > 1.e-09 && rank == 0) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %… in main()
|
| H A D | ex148.c | 11 PetscReal enorm; in main() local
61 PetscCall(VecNorm(output, NORM_1, &enorm)); in main()
62 …if (enorm > 1.e-11 && rank == 0) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %… in main()
|
| H A D | ex228.c | 20 PetscReal enorm; /* norm for sanity check */ in main() local
73 PetscCall(VecNorm(z1, NORM_1, &enorm)); in main()
74 …if (enorm > 1.e-9) PetscCall(PetscPrintf(PETSC_COMM_WORLD, " Error norm of |x - z1| %g dimension … in main()
|
| H A D | ex155.c | 10 PetscReal enorm; in main() local
66 PetscCall(VecNorm(output, NORM_1, &enorm)); in main()
67 …if (enorm > 1.e-14) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %e\n", enorm)); in main()
|
| H A D | ex149.c | 13 PetscReal enorm; in main() local
66 PetscCall(VecNorm(output, NORM_1, &enorm)); in main()
68 if (rank == 0) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %e\n", enorm)); in main()
|
| H A D | ex153.c | 9 PetscReal enorm; in main() local
66 PetscCall(VecNorm(output, NORM_1, &enorm)); in main()
68 PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %e\n", enorm)); in main()
|
| H A D | ex157.c | 9 PetscReal enorm; in main() local
77 PetscCall(VecNorm(output, NORM_1, &enorm)); in main()
79 PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %e\n", enorm)); in main()
|
| H A D | ex142.c | 25 PetscReal enorm; in main() local
141 PetscCall(VecNorm(z, NORM_1, &enorm)); in main()
142 …if (enorm > 1.e-11) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %g\n", (double… in main()
|
| H A D | ex121.c | 26 PetscReal enorm; in main() local
119 PetscCall(VecNorm(z1, NORM_1, &enorm)); in main()
120 …if (enorm > 1.e-11) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |z1 - z2| %g\n", (doub… in main()
|
| H A D | ex150.c | 12 PetscReal enorm; in main() local
76 PetscCall(VecNorm(output, NORM_1, &enorm)); in main()
78 if (rank == 0) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Error norm of |x - z| %e\n", enorm)); in main()
|
| H A D | ex144.c | 22 PetscReal enorm; in main() local
166 PetscCall(VecNorm(final, NORM_1, &enorm)); in main()
167 …if (enorm > 1.e-10) PetscCall(PetscPrintf(PETSC_COMM_WORLD, " Error norm of |x - z| = %e\n", eno… in main()
|
| H A D | ex146.c | 18 PetscScalar a, *x_arr, *y_arr, *z_arr, enorm; in main() local
158 PetscCall(VecNorm(final, NORM_1, &enorm)); in main()
159 PetscCall(PetscPrintf(PETSC_COMM_WORLD, " Error norm of |x - z| = %e\n", enorm)); in main()
|
| H A D | ex147.c | 18 PetscScalar a, *x_arr, *y_arr, *z_arr, enorm; in main() local
|
| /petsc/src/dm/tests/ |
| H A D | ex27.c | 30 PetscReal norm, enorm; in main() local
154 PetscCall(VecNorm(x, NORM_1, &enorm)); in main()
155 PetscCall(PetscPrintf(PETSC_COMM_SELF, "|x-z| = %g\n", enorm)); in main()
161 PetscCall(VecNorm(xx, NORM_1, &enorm)); in main()
162 PetscCall(PetscPrintf(PETSC_COMM_SELF, "|xx-zz| = %g\n", enorm)); in main()
167 PetscCall(VecNorm(y, NORM_1, &enorm)); in main()
169 PetscCall(PetscPrintf(PETSC_COMM_SELF, "|y-yy| = %g\n", enorm)); in main()
174 PetscCall(VecNorm(z, NORM_1, &enorm)); in main()
176 PetscCall(PetscPrintf(PETSC_COMM_SELF, "|z-zz| = %g\n", enorm)); in main()
|
| H A D | ex28.c | 33 PetscReal norm, enorm; in main() local
178 PetscCall(VecNorm(x, NORM_1, &enorm)); in main()
179 PetscCall(PetscPrintf(PETSC_COMM_SELF, "|x-z| = %g\n", enorm)); in main()
185 PetscCall(VecNorm(xx, NORM_1, &enorm)); in main()
186 PetscCall(PetscPrintf(PETSC_COMM_SELF, "|xx-zz| = %g\n", enorm)); in main()
191 PetscCall(VecNorm(y, NORM_1, &enorm)); in main()
193 PetscCall(PetscPrintf(PETSC_COMM_SELF, "|y-yy| = %g\n", enorm)); in main()
198 PetscCall(VecNorm(z, NORM_1, &enorm)); in main()
200 PetscCall(PetscPrintf(PETSC_COMM_SELF, "|z-zz| = %g\n", enorm)); in main()
|
| /petsc/src/ts/adapt/impls/dsp/ |
| H A D | adaptdsp.c | 72 PetscReal enorm = -1; in TSAdaptChoose_DSP() local
84 PetscCall(TSEvaluateWLTE(ts, adapt->wnormtype, &order, &enorm)); in TSAdaptChoose_DSP()
85 …PetscCheck(enorm < 0 || order >= 1, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OUTOFRANGE,… in TSAdaptChoose_DSP()
96 PetscCall(TSErrorWeightedNorm(ts, ts->vec_sol, Y, adapt->wnormtype, &enorm, &enorma, &enormr)); in TSAdaptChoose_DSP()
99 if (enorm < 0) { in TSAdaptChoose_DSP()
135 enorm = PetscMax(enorm, PETSC_SMALL); in TSAdaptChoose_DSP()
136 ctr0 = PetscPowReal(1 / enorm, 1 / k); in TSAdaptChoose_DSP()
166 *wlte = enorm; in TSAdaptChoose_DSP()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex13.c | 34 PetscReal enorm; in main() local
107 enorm = 0.0; in main()
108 …for (i = 0; i < N; i++) enorm += PetscRealPart(PetscConj(solution[i] - userx[i]) * (solution[i] - … in main()
109 enorm *= PetscRealPart(hx * hy); in main()
110 …PETSC_COMM_WORLD, "m %" PetscInt_FMT " n %" PetscInt_FMT " error norm %g\n", m, n, (double)enorm)); in main()
|
| H A D | ex13f90.F90 | 242 PetscReal enorm
322 enorm = PetscRealPart(cnorm*hx*hy)
323 write (6, 115) m, n, enorm
|
| /petsc/src/snes/interface/noise/ |
| H A D | snesnoise.c | 217 PetscReal yy1n, yy2n, enorm; in JacMatMultCompare() local
272 PetscCall(VecNorm(yy2, NORM_2, &enorm)); in JacMatMultCompare()
273 enorm = enorm / yy1n; in JacMatMultCompare()
274 … PetscCall(PetscFPrintf(comm, stdout, "h = %g: relative error = %g\n", (double)h, (double)enorm)); in JacMatMultCompare()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex30.c | 35 PetscReal rnorm, enorm; in main() local
299 PetscCall(VecNorm(u, NORM_2, &enorm)); in main()
300 PetscCall(PetscPrintf(PETSC_COMM_WORLD, " Error norm %g\n", (double)enorm)); in main()
338 PetscCall(VecNorm(xstar, NORM_2, &enorm)); in main()
339 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Error norm %g\n", (double)enorm)); in main()
|