Home
last modified time | relevance | path

Searched refs:PetscAbsScalar (Results 1 – 25 of 206) sorted by relevance

123456789

/petsc/src/vec/vec/tests/
H A Dex43.c45 if (PetscAbsScalar(val_mdot[j] - val_dot[j]) / PetscAbsScalar(val_dot[j]) > 1e-5) { in main()
46 …", val_mdot[j]=%g, val_dot[j]=%g\n", i, j, (double)PetscAbsScalar(val_mdot[j]), (double)PetscAbsSc… in main()
49 if (PetscAbsScalar(tval_mdot[j] - tval_dot[j]) / PetscAbsScalar(tval_dot[j]) > 1e-5) { in main()
50 …tval_mdot[j]=%g, tval_dot[j]=%g\n", i, j, (double)PetscAbsScalar(tval_mdot[j]), (double)PetscAbsSc… in main()
/petsc/src/ksp/ksp/impls/cg/pipecgrr/
H A Dpipecgrr.c83 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- r'*u = r'*B*r = e'*A'*B*A*e */ in KSPSolve_PIPECGRR()
162 if (ksp->normtype == KSP_NORM_NATURAL) dp = PetscSqrtReal(PetscAbsScalar(gamma)); in KSPSolve_PIPECGRR()
193 …ncr = PetscSqrtReal(Anorm * xnp + 2.0 * Anorm * PetscAbsScalar(alphap) * dpp + rnp + 2.0 * PetscAb… in KSPSolve_PIPECGRR()
194 …ncw = PetscSqrtReal(Anorm * unp + 2.0 * Anorm * PetscAbsScalar(alphap) * dq + wnp + 2.0 * PetscAbs… in KSPSolve_PIPECGRR()
197 …ncs = PetscSqrtReal(Anorm * unp + 2.0 * Anorm * PetscAbsScalar(betap) * pnp + wnp + 2.0 * PetscAbs… in KSPSolve_PIPECGRR()
198 …tReal((mnz * sqn + 2) * Anorm * dm + 2.0 * Anorm * PetscAbsScalar(betap) * qnp + 2.0 * PetscAbsSca… in KSPSolve_PIPECGRR()
203 …scSqrtReal((mnz * sqn + 1) * Anorm * xnp + db) * eps + PetscSqrtReal(PetscAbsScalar(alphap) * mnz … in KSPSolve_PIPECGRR()
205 …errw = PetscSqrtReal(mnz * sqn * Anorm * unp) * eps + PetscSqrtReal(PetscAbsScalar(alphap) * mnz *… in KSPSolve_PIPECGRR()
216 …errr = errr + PetscAbsScalar(alphap) * PetscAbsScalar(betap) * errs + PetscAbsScalar(alphap) *… in KSPSolve_PIPECGRR()
217 errs = errw + PetscAbsScalar(betap) * errs + errncs; in KSPSolve_PIPECGRR()
[all …]
/petsc/src/mat/impls/baij/seq/
H A Ddgefa7.c23 …PetscAbsScalar(a[0]) + PetscAbsScalar(a[8]) + PetscAbsScalar(a[16]) + PetscAbsScalar(a[24]) + Pets… in PetscKernel_A_gets_inverse_A_7()
36 max = PetscAbsScalar(aa[0]); in PetscKernel_A_gets_inverse_A_7()
39 tmp = PetscAbsScalar(aa[ll]); in PetscKernel_A_gets_inverse_A_7()
H A Ddgefa6.c23 …(1.e-12 + PetscAbsScalar(a[0]) + PetscAbsScalar(a[7]) + PetscAbsScalar(a[14]) + PetscAbsScalar(a[2… in PetscKernel_A_gets_inverse_A_6()
36 max = PetscAbsScalar(aa[0]); in PetscKernel_A_gets_inverse_A_6()
39 tmp = PetscAbsScalar(aa[ll]); in PetscKernel_A_gets_inverse_A_6()
H A Ddgefa5.c23 …5 * shift * (1.e-12 + PetscAbsScalar(a[0]) + PetscAbsScalar(a[6]) + PetscAbsScalar(a[12]) + PetscA… in PetscKernel_A_gets_inverse_A_5()
36 max = PetscAbsScalar(aa[0]); in PetscKernel_A_gets_inverse_A_5()
39 tmp = PetscAbsScalar(aa[ll]); in PetscKernel_A_gets_inverse_A_5()
H A Ddgefa4.c23 …shift = .25 * shift * (1.e-12 + PetscAbsScalar(a[0]) + PetscAbsScalar(a[5]) + PetscAbsScalar(a[10]… in PetscKernel_A_gets_inverse_A_4()
36 max = PetscAbsScalar(aa[0]); in PetscKernel_A_gets_inverse_A_4()
39 tmp = PetscAbsScalar(aa[ll]); in PetscKernel_A_gets_inverse_A_4()
H A Ddgefa3.c23 …shift = .333 * shift * (1.e-12 + PetscAbsScalar(a[0]) + PetscAbsScalar(a[4]) + PetscAbsScalar(a[8]… in PetscKernel_A_gets_inverse_A_3()
36 max = PetscAbsScalar(aa[0]); in PetscKernel_A_gets_inverse_A_3()
39 tmp = PetscAbsScalar(aa[ll]); in PetscKernel_A_gets_inverse_A_3()
H A Ddgefa2.c23 shift = .25 * shift * (1.e-12 + PetscAbsScalar(a[0]) + PetscAbsScalar(a[3])); in PetscKernel_A_gets_inverse_A_2()
36 max = PetscAbsScalar(aa[0]); in PetscKernel_A_gets_inverse_A_2()
39 tmp = PetscAbsScalar(aa[ll]); in PetscKernel_A_gets_inverse_A_2()
167 max = PetscAbsScalar(aa[0]); in PetscKernel_A_gets_inverse_A_9()
170 tmp = PetscAbsScalar(aa[ll]); in PetscKernel_A_gets_inverse_A_9()
338 max = PetscAbsScalar(aa[0]); in PetscKernel_A_gets_inverse_A_15()
341 tmp = PetscAbsScalar(aa[ll]); in PetscKernel_A_gets_inverse_A_15()
H A Ddgefa.c39 max = PetscAbsScalar(aa[0]); in PetscLINPACKgefa()
42 tmp = PetscAbsScalar(aa[ll]); in PetscLINPACKgefa()
/petsc/src/snes/interface/noise/
H A Dsnesdnest.c134 if (err2 <= PetscAbsScalar(est4) * .1) *fder2 = est4; in SNESNoise_dnest_()
135 else if (err2 < PetscAbsScalar(est4)) *fder2 = est3; in SNESNoise_dnest_()
167 d__3 = PetscAbsScalar(d__1); in SNESNoise_dnest_()
203 d__1 = PetscAbsScalar(f_max); in SNESNoise_dnest_()
204 d__2 = PetscAbsScalar(f_min); in SNESNoise_dnest_()
222 *hopt = PetscSqrtScalar(*fnoise / PetscAbsScalar(*fder2)) * 1.68; in SNESNoise_dnest_()
242 *hopt = PetscSqrtScalar(*fnoise / PetscAbsScalar(*fder2)) * 1.68; in SNESNoise_dnest_()
/petsc/src/mat/utils/
H A Dzerodiag.c90 if (k >= nz || PetscAbsScalar(v[k]) <= abstol) { in MatReorderForNonzeroDiagonal_SeqAIJ()
92 repla = (k >= nz) ? 0.0 : PetscAbsScalar(v[k]); in MatReorderForNonzeroDiagonal_SeqAIJ()
97 if (icol[j[k]] > prow && PetscAbsScalar(v[k]) > repla) { in MatReorderForNonzeroDiagonal_SeqAIJ()
111 if (icol[j[k]] < prow && PetscAbsScalar(v[k]) > repla) { in MatReorderForNonzeroDiagonal_SeqAIJ()
116 if (icol[jj[kk]] == prow && PetscAbsScalar(vv[kk]) > abstol) { in MatReorderForNonzeroDiagonal_SeqAIJ()
133 if (icol[jj[kk]] == prow && PetscAbsScalar(vv[kk]) > abstol) { in MatReorderForNonzeroDiagonal_SeqAIJ()
/petsc/src/snes/linesearch/impls/cp/
H A Dlinesearchcp.c34 if (PetscAbsScalar(fty_old) < atol * ynorm) { in SNESLineSearchApply_CP()
37 …at initial point because dot(F,Y) = %g < atol*||y|| = %g\n", (double)PetscAbsScalar(fty_old), (dou… in SNESLineSearchApply_CP()
70 if (PetscAbsScalar(fty) / PetscAbsScalar(fty_init) < rtol) { in SNESLineSearchApply_CP()
73 …search: abs(fty/fty_init) = %g <= rtol = %g\n", (double)(PetscAbsScalar(fty) / PetscAbsScalar(fty… in SNESLineSearchApply_CP()
80 if (PetscAbsScalar(fty) < atol * ynorm && i > 0) { in SNESLineSearchApply_CP()
83 …tor, " Line search: abs(fty)/||y|| = %g <= atol = %g\n", (double)(PetscAbsScalar(fty) / ynorm),… in SNESLineSearchApply_CP()
/petsc/src/ksp/ksp/impls/symmlq/
H A Dsymmlq.c59 if (PetscAbsScalar(dp) < symmlq->haptol) { in KSPSolve_SYMMLQ()
60 …PetscCall(PetscInfo(ksp, "Detected happy breakdown %g tolerance %g\n", (double)PetscAbsScalar(dp),… in KSPSolve_SYMMLQ()
75 s_prod = PetscAbsScalar(beta1); in KSPSolve_SYMMLQ()
131 if (PetscAbsScalar(dp) < symmlq->haptol) { in KSPSolve_SYMMLQ()
132 …PetscCall(PetscInfo(ksp, "Detected happy breakdown %g tolerance %g\n", (double)PetscAbsScalar(dp),… in KSPSolve_SYMMLQ()
161 s_prod = s_prod * PetscAbsScalar(s); in KSPSolve_SYMMLQ()
163 else np = s_prod / PetscAbsScalar(c); /* residual norm for xc_k (CGNORM) */ in KSPSolve_SYMMLQ()
/petsc/src/ts/tests/
H A Dex13.c72 serr = PetscAbsScalar(func(p, testtime) - aW[0]); in main()
73 derr = PetscAbsScalar(dfunc(p, testtime) - aWdot[0]); in main()
86 serr = PetscAbsScalar(func(p, testtime) - aW[0]); in main()
97 derr = PetscAbsScalar(dfunc(p, testtime) - aWdot[0]); in main()
111 serr = PetscAbsScalar(func(p, testtime) - aW[0]); in main()
112 derr = PetscAbsScalar(dfunc(p, testtime) - aWdot[0]); in main()
150 serr = PetscAbsScalar(func(p, testtime) - aW[0]); in main()
157 derr = PetscAbsScalar(dfunc(p, testtime) - aWdot[0]); in main()
/petsc/src/snes/linesearch/impls/bisection/
H A Dlinesearchbisection.c88 if (PetscAbsScalar(fty) <= atol * ynorm) { in SNESLineSearchApply_Bisection()
91 …r, " Line search: abs(fty)/||y|| = %g <= atol = %g\n", (double)(PetscAbsScalar(fty) / ynorm),… in SNESLineSearchApply_Bisection()
98 if (PetscAbsScalar(fty) / PetscAbsScalar(fty_initial) <= rtol) { in SNESLineSearchApply_Bisection()
101 …rch: abs(fty/fty_initial) = %g <= rtol = %g\n", (double)(PetscAbsScalar(fty) / PetscAbsScalar(fty… in SNESLineSearchApply_Bisection()
/petsc/src/ts/tutorials/network/
H A DpipeImpls.c27 x[i].h = H0 - c * Q0 * PetscAbsScalar(Q0) * coords[i]; in PipeComputeSteadyState()
82 …xold_i) + dt * (GRAV * pipe->A * dhdx(xold, i, ilast, dx) + pipe->R * qavg * PetscAbsScalar(qavg)); in PipeIFunctionLocal_Lax()
90 …a) * (x[start].h - xold[start + 1].h) + dt * R * xold[start + 1].q * PetscAbsScalar(xold[start + 1… in PipeIFunctionLocal_Lax()
91 …a) * (x[ilast].h - xold[ilast - 1].h) + dt * R * xold[ilast - 1].q * PetscAbsScalar(xold[ilast - 1… in PipeIFunctionLocal_Lax()
/petsc/src/mat/tests/
H A Dex178.c116 …PetscAbsScalar(vals[j] - 1) < PETSC_SMALL, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Row %" PetscInt_FMT "… in MatIsDiagonal()
119 …PetscAbsScalar(vals[j]) < PETSC_SMALL, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Row %" PetscInt_FMT " has… in MatIsDiagonal()
H A Dex118.c148 dot = PetscAbsScalar(dot - (PetscScalar)1.0); in CkEigenSolutions()
150 dot = PetscAbsScalar(dot); in CkEigenSolutions()
152 if (PetscAbsScalar(dot) > dot_max) dot_max = PetscAbsScalar(dot); in CkEigenSolutions()
/petsc/src/ksp/ksp/impls/cr/
H A Dcr.c60 dp = PetscSqrtReal(PetscAbsScalar(btop)); /* dp = sqrt(R,AR) */ in KSPSolve_CR()
62 if (PetscAbsScalar(btop) < 0.0) { in KSPSolve_CR()
103 dp = PetscSqrtReal(PetscAbsScalar(btop)); /* dp = sqrt(R,AR) */ in KSPSolve_CR()
114 if (PetscAbsScalar(btop) < 0.0) { in KSPSolve_CR()
/petsc/src/ksp/ksp/impls/gmres/pipefgmres/
H A Dpipefgmres.c184 for (k = 0; k <= loc_it; k++) tt -= PetscAbsScalar(lhh[k]) * PetscAbsScalar(lhh[k]); in KSPPIPEFGMRESCycle()
227 hapbnd = PetscAbsScalar((tt) / *RS(loc_it)); in KSPPIPEFGMRESCycle()
395 hapbnd = PetscMin(PetscAbsScalar(hh[it + 1] / rs[it]), pipefgmres->haptol); in KSPPIPEFGMRESUpdateHessenberg()
396 if (PetscAbsScalar(hh[it + 1]) < hapbnd) { in KSPPIPEFGMRESUpdateHessenberg()
397 …%" PetscInt_FMT ") = %14.12e\n", (double)hapbnd, it + 1, it, (double)PetscAbsScalar(*HH(it + 1, it… in KSPPIPEFGMRESUpdateHessenberg()
424 …PetscReal delta = PetscSqrtReal(PetscSqr(PetscAbsScalar(hh[it])) + PetscSqr(PetscAbsScalar(hh[it +… in KSPPIPEFGMRESUpdateHessenberg()
436 *res = PetscAbsScalar(rs[it + 1]); in KSPPIPEFGMRESUpdateHessenberg()
/petsc/src/sys/objects/device/tests/
H A Dex7.c76 …mory that was not cleared, ptr[%" PetscInt_FMT "] %g != 0", i, (double)PetscAbsScalar(tmp_ptr[i])); in TestAllocate()
93 …ro() did not clear memory, ptr[%" PetscInt_FMT "] %g != 0", i, (double)PetscAbsScalar(tmp_ptr[i])); in TestAllocate()
132 …) was not properly serialized, ptr[%" PetscInt_FMT "] %g != 0", i, (double)PetscAbsScalar(ptr[i])); in TestAsyncCoherence()
143 …s not properly serialized, ptr[%" PetscInt_FMT "] %g != 0", i, (double)PetscAbsScalar(tmp_ptr[i])); in TestAsyncCoherence()
/petsc/src/ksp/pc/impls/tfs/
H A Dtfs.h91 #define MAX_FABS(x, y) (PetscAbsScalar(x) > PetscAbsScalar(y)) ? ((PetscScalar)x) : ((PetscScalar)y)
92 #define MIN_FABS(x, y) (PetscAbsScalar(x) < PetscAbsScalar(y)) ? ((PetscScalar)x) : ((PetscScalar)y)
/petsc/src/dm/dt/fe/tests/
H A Dex4.c47 …RLD, "%s: %g (%g - %g)\n", msg, (double)PetscAbsScalar(a - b), (double)PetscAbsScalar(a), (double)…
/petsc/src/ksp/ksp/impls/cg/cgne/
H A Dcgne.c98 dp = PetscSqrtReal(PetscAbsScalar(beta)); in KSPSolve_CGNE()
129 e[i] = PetscSqrtReal(PetscAbsScalar(b)) / a; in KSPSolve_CGNE()
139 if (eigs) d[i] = PetscSqrtReal(PetscAbsScalar(b)) * e[i] + 1.0 / a; in KSPSolve_CGNE()
153 dp = PetscSqrtReal(PetscAbsScalar(beta)); in KSPSolve_CGNE()
/petsc/src/ksp/ksp/impls/gmres/pgmres/
H A Dpgmres.c253 hapbnd = PetscMin(PetscAbsScalar(hh[it + 1] / rs[it]), pgmres->haptol); in KSPPGMRESUpdateHessenberg()
254 if (PetscAbsScalar(hh[it + 1]) < hapbnd) { in KSPPGMRESUpdateHessenberg()
255 …%" PetscInt_FMT ") = %14.12e\n", (double)hapbnd, it + 1, it, (double)PetscAbsScalar(*HH(it + 1, it… in KSPPGMRESUpdateHessenberg()
282 …PetscReal delta = PetscSqrtReal(PetscSqr(PetscAbsScalar(hh[it])) + PetscSqr(PetscAbsScalar(hh[it +… in KSPPGMRESUpdateHessenberg()
294 *res = PetscAbsScalar(rs[it + 1]); in KSPPGMRESUpdateHessenberg()

123456789