Home
last modified time | relevance | path

Searched refs:PetscAbs (Results 1 – 25 of 48) sorted by relevance

12

/petsc/src/tao/leastsquares/impls/pounders/
H A Dgqt.c14 e = PetscAbs(r[0]); in estsv()
22 if (z[i] >= 0.0) e = -PetscAbs(e); in estsv()
23 else e = PetscAbs(e); in estsv()
25 if (PetscAbs(e - z[i]) > PetscAbs(r[i + ldr * i])) { in estsv()
26 temp = PetscMin(0.01, PetscAbs(r[i + ldr * i])) / PetscAbs(e - z[i]); in estsv()
40 s = PetscAbs(e - z[i]); in estsv()
41 sm = PetscAbs(e + z[i]); in estsv()
42 for (j = i + 1; j < n; j++) sm += PetscAbs(z[j] + wm * r[i + ldr * j]); in estsv()
61 if (PetscAbs(z[j]) > PetscAbs(r[j + ldr * j])) { in estsv()
62 temp = PetscMin(0.01, PetscAbs(r[j + ldr * j] / z[j])); in estsv()
[all …]
/petsc/src/ts/tutorials/
H A Dex9.c54 return (PetscAbs(a) > PetscAbs(b)) ? a : b; in MaxAbs()
58 return (PetscAbs(a) < PetscAbs(b)) ? a : b; in MinAbs()
62 return (a * b < 0) ? 0 : Sgn(a) * PetscMin(PetscAbs(a), PetscAbs(b)); in MinMod2()
66 return (a * b < 0) ? 0 : Sgn(a) * PetscMax(PetscAbs(a), PetscAbs(b)); in MaxMod2()
70 …eturn (a * b < 0 || a * c < 0) ? 0 : Sgn(a) * PetscMin(PetscAbs(a), PetscMin(PetscAbs(b), PetscAbs… in MinMod3()
343 u[0] = PetscAbs(x0); in PhysicsSample_Advect()
435 *maxspeed = (PetscAbs(uL[0]) > PetscAbs(uR[0])) ? uL[0] : uR[0]; in PhysicsRiemann_Burgers_Exact()
445 flux[0] = 0.25 * (PetscSqr(uL[0]) + PetscSqr(uR[0])) - 0.5 * PetscAbs(speed) * (uR[0] - uL[0]); in PhysicsRiemann_Burgers_Roe()
471 c = PetscMax(PetscAbs(uL[0]), PetscAbs(uR[0])); in PhysicsRiemann_Burgers_Rusanov()
575 …flux[0] = 0.5 * (TrafficFlux(a, uL[0]) + TrafficFlux(a, uR[0])) - 0.5 * PetscAbs(speed) * (uR[0]… in PhysicsRiemann_Traffic_Roe()
[all …]
/petsc/src/dm/dt/tests/
H A Dex2.c43 if (PetscAbs(R[j]) < 10 * PETSC_MACHINE_EPSILON) R[j] = 0; in main()
H A Dex1.c197 if (PetscAbs(zeroth) < 1e-10) zeroth = 0.; in main()
198 if (PetscAbs(first) < 1e-10) first = 0.; in main()
199 if (PetscAbs(second) < 1e-10) second = 0.; in main()
/petsc/src/vec/vec/tests/
H A Dex62.c41 …PetscCheck(PetscAbs(nrm - ans) < PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Diff is too big, … in main()
/petsc/src/mat/tests/
H A Dex30.c152 …PetscCheck(PetscAbs(norm2 - norm2_inplace) <= tol, PETSC_COMM_SELF, PETSC_ERR_PLIB, "ILU(0) %g and… in main()
207 …PetscCheck(PetscAbs(norm2 - norm2_inplace) <= tol, PETSC_COMM_SELF, PETSC_ERR_PLIB, "ICC(0) %g and… in main()
H A Dex84.c90 …for (i = 0; i < NNORMS; i++) PetscCheck(PetscAbs(norms0[i] - norms1[i]) <= PETSC_SMALL, PETSC_COMM… in main()
/petsc/src/sys/classes/draw/impls/image/
H A Ddrawimage.h56 int dx = PetscAbs(x_2 - x_1), sx = (x_2 - x_1) >= 0 ? +1 : -1; in PetscImageDrawLine()
57 int dy = PetscAbs(y_2 - y_1), sy = (y_2 - y_1) >= 0 ? +1 : -1; in PetscImageDrawLine()
H A Ddrawimage.c140 if (x_1 == x_2 && PetscAbs(y_1 - y_2) > 7) { in PetscDrawArrow_Image()
149 if (y_1 == y_2 && PetscAbs(x_1 - x_2) > 7) { in PetscDrawArrow_Image()
/petsc/src/ts/tutorials/multirate/
H A Dfinitevolume1d.c24 return (PetscAbs(a) < PetscAbs(b)) ? a : b; in MinAbs()
28 return (a * b < 0) ? 0 : Sgn(a) * PetscMin(PetscAbs(a), PetscAbs(b)); in MinMod2()
32 return (a * b < 0) ? 0 : Sgn(a) * PetscMax(PetscAbs(a), PetscAbs(b)); in MaxMod2()
36 …eturn (a * b < 0 || a * c < 0) ? 0 : Sgn(a) * PetscMin(PetscAbs(a), PetscMin(PetscAbs(b), PetscAbs… in MinMod3()
H A Dex4.c50 return (PetscAbs(a) > PetscAbs(b)) ? a : b; in MaxAbs()
110 u[0] = PetscAbs(x0); in PhysicsSample_Advect()
254 …s = PetscMax(PetscAbs(L.u) + PetscSqrtScalar(g * L.h), PetscAbs(R.u) + PetscSqrtScalar(g *… in PhysicsRiemann_Shallow_Rusanov()
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()
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()
H A Dex7.c129 u[0] = PetscAbs(x0); in PhysicsSample_Advect()
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()
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()
H A Dex6.c91 u[0] = PetscAbs(x0); in PhysicsSample_Advect()
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()
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()
H A Dex8.c84 u[0] = PetscAbs(x0); in PhysicsSample_Advect()
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()
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()
1305 … for (k = 0; k < dof; k++) nrm1 = nrm1 + hf * PetscAbs(ptr_X[i * dof + k] - ptr_XR[i * dof + k]); in main()
/petsc/src/sys/tests/
H A Dex15.cxx35 …if (check) PetscCheck(PetscAbs(ans - res) < 1e-5, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Expected %g, b… in TestComplexOperators()
/petsc/src/ts/tutorials/phasefield/
H A Dheat.c348 max = PetscMax(max, PetscAbs(yy[0])); in MyMonitor()
351 max = PetscMax(max, PetscAbs(yy[1])); in MyMonitor()
H A Dbiharmonic.c571 max = PetscMax(max, PetscAbs(yy[0])); in MyMonitor()
612 max = PetscMax(max, PetscAbs(yy[1])); in MyMonitor()
616 max = PetscMax(max, PetscAbs(yy[1] + yy[2])); in MyMonitor()
/petsc/src/ksp/ksp/tutorials/
H A Dex76.c291 …PetscCheck(reason[0] == reason[1] && PetscAbs(iterations[0] - iterations[1]) <= 3, PetscObjectComm… in main()
301 …PetscCheck(reason[0] == reason[1] && PetscAbs(iterations[0] - iterations[2]) <= 3, PetscObjectComm… in main()
/petsc/src/snes/impls/qn/
H A Dqn.c228 …if (PetscAbs(PetscRealPart(DolddotD)) > qn->powell_gamma * PetscAbs(PetscRealPart(DolddotDold))) p… in SNESSolve_QN()
/petsc/src/sys/classes/draw/impls/x/
H A Dxops.c110 if (x_1 == x_2 && PetscAbs(y_1 - y_2) > 7) { in PetscDrawArrow_X()
119 if (y_1 == y_2 && PetscAbs(x_1 - x_2) > 7) { in PetscDrawArrow_X()
157 w = PetscAbs(w); in PetscDrawEllipse_X()
160 h = PetscAbs(h); in PetscDrawEllipse_X()
/petsc/src/ts/utils/dmplexlandau/kokkos/
H A Dlandau.kokkos.cxx616 …PetscReal Ud[2][2], Uk[2][2], mask = (PetscAbs(vj[0] - x) < 100 * PETSC_SQRT_MACHINE_EPSILON && Pe… in LandauKokkosJacobian()
619 …z[ipidx], mask = (PetscAbs(vj[0]-x) < 100*PETSC_SQRT_MACHINE_EPSILON && PetscAbs(vj[1]-y) < 100*PE… in LandauKokkosJacobian()
/petsc/src/dm/impls/plex/tutorials/
H A Dex15.c148 …if (PetscAbs(array_serial[i] - array_load_bcast[i]) > tol) PetscCall(PetscPrintf(comm, "DoF %" Pet… in VerifyLoadedSolution()
/petsc/src/dm/impls/plex/tests/
H A Dex8.c144 #define RelativeError(a, b) PetscAbs(a - b) / (1.0 + PetscMax(PetscAbs(a), PetscAbs(b)))
/petsc/src/snes/tutorials/
H A Dex56.c274 … if (PetscAbs(faceCoord - b) < PETSC_SMALL) { /* domain have not been set yet, still [0,1]^3 */ in main()
451 …disp[iter] - mdisp[iter - 1]), (double)err[iter], (double)(PetscLogReal(PetscAbs(err[iter - 1] / e… in main()
/petsc/src/ksp/pc/impls/gamg/
H A Dgeo.c355 if ((v = PetscAbs(PetscRealPart(alpha[tt]) - 0.5)) > worst) worst = v; in triangulateAndFormProl()
384 if (PetscAbs(PetscRealPart(shp)) > 1.e-6) { in triangulateAndFormProl()

12