| /petsc/src/vec/vec/tests/ |
| H A D | ex34.c | 11 PetscReal nrm1, nrm2, nrm3, nrm4; in main() local
30 PetscCall(VecNorm(V, NORM_1, &nrm1)); in main()
32 PetscCall(PetscPrintf(comm, "Original: norm1=%e,norm2=%e\n", (double)nrm1, (double)nrm2)); in main()
35 PetscCall(VecNorm(V, NORM_1, &nrm1)); in main()
37 PetscCall(PetscPrintf(comm, "cached: norm1=%e,norm2=%e\n", (double)nrm1, (double)nrm2)); in main()
45 PetscCall(VecNorm(V, NORM_1, &nrm1)); in main()
47 PetscCall(PetscPrintf(comm, "Altered: norm1=%e,norm2=%e\n", (double)nrm1, (double)nrm2)); in main()
50 PetscCall(VecNorm(V, NORM_1, &nrm1)); in main()
52 PetscCall(PetscPrintf(comm, "recomputed: norm1=%e,norm2=%e\n", (double)nrm1, (double)nrm2)); in main()
60 PetscCall(VecNorm(V, NORM_1, &nrm1)); in main()
[all …]
|
| /petsc/src/ts/tutorials/multirate/ |
| H A D | ex7.c | 632 static PetscErrorCode SolutionErrorNorms(FVCtx *ctx, DM da, PetscReal t, Vec X, PetscReal *nrm1) in SolutionErrorNorms() argument
649 …if (i < count_slow / 2 || i > count_slow / 2 + count_fast - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - p… in SolutionErrorNorms()
650 else *nrm1 += hf * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms()
798 PetscReal nrm1 = 0; in main() local
799 PetscCall(SolutionErrorNorms(&ctx, da, ptime, X, &nrm1)); in main()
800 PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %g\n", (double)nrm1)); in main()
803 PetscReal nrm1 = 0; in main() local
818 …if (i < count_slow / 2 || i > count_slow / 2 + count_fast - 1) nrm1 = nrm1 + hs * PetscAbs(ptr_X[i… in main()
819 else nrm1 = nrm1 + hf * PetscAbs(ptr_X[i] - ptr_XR[i]); in main()
823 PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %g\n", (double)nrm1)); in main()
|
| H A D | ex6.c | 180 …PetscErrorCode SolutionErrorNorms_2WaySplit(FVCtx *ctx, DM da, PetscReal t, Vec X, PetscReal *nrm1) in SolutionErrorNorms_2WaySplit() argument
196 if (i < ctx->sf || i > ctx->fs - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_2WaySplit()
197 else *nrm1 += hf * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_2WaySplit()
912 PetscReal nrm1 = 0; in main() local
913 PetscCall(SolutionErrorNorms_2WaySplit(&ctx, da, ptime, X, &nrm1)); in main()
914 PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %g\n", (double)nrm1)); in main()
917 PetscReal nrm1 = 0; in main() local
933 … for (k = 0; k < dof; k++) nrm1 = nrm1 + hs * PetscAbs(ptr_X[i * dof + k] - ptr_XR[i * dof + k]); in main()
935 … for (k = 0; k < dof; k++) nrm1 = nrm1 + hf * PetscAbs(ptr_X[i * dof + k] - ptr_XR[i * dof + k]); in main()
939 PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %g\n", (double)nrm1)); in main()
|
| H A D | ex8.c | 193 …PetscErrorCode SolutionErrorNorms_3WaySplit(FVCtx *ctx, DM da, PetscReal t, Vec X, PetscReal *nrm1) in SolutionErrorNorms_3WaySplit() argument
210 if (i < ctx->sm || i > ctx->ms - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_3WaySplit()
211 else if (i < ctx->mf || i > ctx->fm - 1) *nrm1 += hm * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_3WaySplit()
212 else *nrm1 += hf * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_3WaySplit()
1280 PetscReal nrm1 = 0; in main() local
1281 PetscCall(SolutionErrorNorms_3WaySplit(&ctx, da, ptime, X, &nrm1)); in main()
1282 PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %g\n", (double)nrm1)); in main()
1285 PetscReal nrm1 = 0; in main() local
1301 … for (k = 0; k < dof; k++) nrm1 = nrm1 + hs * PetscAbs(ptr_X[i * dof + k] - ptr_XR[i * dof + k]); in main()
1303 … for (k = 0; k < dof; k++) nrm1 = nrm1 + hm * PetscAbs(ptr_X[i * dof + k] - ptr_XR[i * dof + k]); in main()
[all …]
|
| H A D | ex4.c | 524 …PetscErrorCode SolutionErrorNorms_2WaySplit(FVCtx *ctx, DM da, PetscReal t, Vec X, PetscReal *nrm1) in SolutionErrorNorms_2WaySplit() argument
540 if (i < ctx->sf || i > ctx->fs - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_2WaySplit()
541 else *nrm1 += hf * PetscAbs(ptr_X[i] - ptr_Y[i]); in SolutionErrorNorms_2WaySplit()
1355 PetscReal nrm1 = 0; in main() local
1356 PetscCall(SolutionErrorNorms_2WaySplit(&ctx, da, ptime, X, &nrm1)); in main()
1357 PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %g\n", (double)nrm1)); in main()
1360 PetscReal nrm1 = 0; in main() local
1376 … for (k = 0; k < dof; k++) nrm1 = nrm1 + hs * PetscAbs(ptr_X[i * dof + k] - ptr_XR[i * dof + k]); in main()
1378 … for (k = 0; k < dof; k++) nrm1 = nrm1 + hf * PetscAbs(ptr_X[i * dof + k] - ptr_XR[i * dof + k]); in main()
1382 PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %g\n", (double)nrm1)); in main()
|
| H A D | ex5.c | 777 PetscReal nrm1, nrmsup; in main() local
787 PetscCall(VecNorm(XR, NORM_1, &nrm1)); in main()
789 nrm1 /= Mx; in main()
790 …PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %g ||x-x_e||_sup %g\n", (double)nrm1, (double)nrms… in main()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex9.c | 1370 static PetscErrorCode SolutionErrorNorms(FVCtx *ctx, DM da, PetscReal t, Vec X, PetscReal *nrm1, Pe… in SolutionErrorNorms() argument
1380 PetscCall(VecNorm(Y, NORM_1, nrm1)); in SolutionErrorNorms()
1382 *nrm1 /= Mx; in SolutionErrorNorms()
1506 PetscReal nrm1, nrmsup; in main() local
1515 PetscCall(SolutionErrorNorms(&ctx, da, ptime, X, &nrm1, &nrmsup)); in main()
1516 …PetscCall(PetscPrintf(comm, "Error ||x-x_e||_1 %8.4e ||x-x_e||_sup %8.4e\n", (double)nrm1, (doubl… in main()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex43.c | 1566 PetscReal nrm1[3], nrm2[3], nrmI[3]; in solve_stokes_2d_coupled() local
1586 PetscCall(VecStrideNorm(X_analytic, 0, NORM_1, &nrm1[0])); in solve_stokes_2d_coupled()
1590 PetscCall(VecStrideNorm(X_analytic, 1, NORM_1, &nrm1[1])); in solve_stokes_2d_coupled()
1594 PetscCall(VecStrideNorm(X_analytic, 2, NORM_1, &nrm1[2])); in solve_stokes_2d_coupled()
|
| /petsc/src/vec/vec/impls/seq/kokkos/ |
| H A D | veckok.kokkos.cxx | 1068 PetscCallCXX(*z = KokkosBlas::nrm1(exec, xv)); in VecNorm_SeqKokkos()
|