Home
last modified time | relevance | path

Searched refs:PetscSqrtReal (Results 1 – 25 of 235) sorted by relevance

12345678910

/petsc/src/ts/tests/
H A Dex2.c214 …PetscSqrtReal(2.0)) * t) / 2.0 - PetscExpReal((2.0 - PetscSqrtReal(2.0)) * t) / (2.0 * PetscSqrtRe… in solx()
219 …PetscSqrtReal(2.0)) * t) / 2.0 - PetscExpReal((2.0 - PetscSqrtReal(2.0)) * t) / PetscSqrtReal(2.0)… in soly()
224 …PetscSqrtReal(2.0)) * t) / 2.0 - PetscExpReal((2.0 - PetscSqrtReal(2.0)) * t) / (2.0 * PetscSqrtRe… in solz()
H A Dex35.c36 r = ck[0] * ((PetscSqrtReal(PETSC_PI) / 2.) * x); in erfinv()
43 r += (ck[k] / temp) * PetscPowReal((PetscSqrtReal(PETSC_PI) / 2.) * x, 2. * k + 1.); in erfinv()
130 mag = PetscSqrtReal(-2.0 * PetscLogReal(u1)); in main()
156 speed = PetscSqrtReal(PetscSqr(v[p * dim]) + PetscSqr(v[p * dim + 1])); in main()
/petsc/src/vec/vec/tutorials/
H A Dex1.c78 v = norm - 2.0 * PetscSqrtReal((PetscReal)n); in main()
84 v = norm - 2.0 * PetscSqrtReal((PetscReal)n); in main()
90 v = norm - 8.0 * PetscSqrtReal((PetscReal)n); in main()
96 v = norm - 18.0 * PetscSqrtReal((PetscReal)n); in main()
102 v = norm - 2.0 * PetscSqrtReal((PetscReal)n); in main()
106 v = norm - 18.0 * PetscSqrtReal((PetscReal)n); in main()
112 v = norm - 38.0 * PetscSqrtReal((PetscReal)n); in main()
118 v = norm - 36.0 * PetscSqrtReal((PetscReal)n); in main()
124 v = norm - 9.0 * PetscSqrtReal((PetscReal)n); in main()
133 v = norm - 9.0 * PetscSqrtReal((PetscReal)n); in main()
[all …]
/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()
101 sqn = PetscSqrtReal(nsi); in KSPSolve_PIPECGRR()
162 if (ksp->normtype == KSP_NORM_NATURAL) dp = PetscSqrtReal(PetscAbsScalar(gamma)); in KSPSolve_PIPECGRR()
193 …errncr = PetscSqrtReal(Anorm * xnp + 2.0 * Anorm * PetscAbsScalar(alphap) * dpp + rnp + 2.0 * Pets… in KSPSolve_PIPECGRR()
194 …errncw = PetscSqrtReal(Anorm * unp + 2.0 * Anorm * PetscAbsScalar(alphap) * dq + wnp + 2.0 * Petsc… in KSPSolve_PIPECGRR()
197 …errncs = PetscSqrtReal(Anorm * unp + 2.0 * Anorm * PetscAbsScalar(betap) * pnp + wnp + 2.0 * Petsc… in KSPSolve_PIPECGRR()
198 …errncz = PetscSqrtReal((mnz * sqn + 2) * Anorm * dm + 2.0 * Anorm * PetscAbsScalar(betap) * qnp + … in KSPSolve_PIPECGRR()
203 …errr = PetscSqrtReal((mnz * sqn + 1) * Anorm * xnp + db) * eps + PetscSqrtReal(PetscAbsScalar(alph… in KSPSolve_PIPECGRR()
204 errs = PetscSqrtReal(mnz * sqn * Anorm * dpp) * eps; in KSPSolve_PIPECGRR()
205 …errw = PetscSqrtReal(mnz * sqn * Anorm * unp) * eps + PetscSqrtReal(PetscAbsScalar(alphap) * mnz *… in KSPSolve_PIPECGRR()
[all …]
/petsc/src/ksp/ksp/impls/cg/stcg/
H A Dstcg.c89 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
91 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
124 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
126 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
158 norm_r = PetscSqrtReal(rr); /* norm_r = |r| */ in KSPCGSolve_STCG()
162 norm_r = PetscSqrtReal(rz); /* norm_r = |r|_M */ in KSPCGSolve_STCG()
202 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
204 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
259 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_STCG()
277 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
[all …]
/petsc/src/ksp/ksp/impls/cg/nash/
H A Dnash.c92 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
94 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
127 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
129 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
161 norm_r = PetscSqrtReal(rr); /* norm_r = |r| */ in KSPCGSolve_NASH()
165 norm_r = PetscSqrtReal(rz); /* norm_r = |r|_M */ in KSPCGSolve_NASH()
205 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
207 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
262 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_NASH()
280 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
[all …]
/petsc/src/dm/dt/tests/
H A Dex3.c26 const PetscReal u = PetscSqrtReal(2.0 + x * x); in func4()
34 else *val = PetscSqrtReal(x) * PetscLogReal(x); in func5()
40 *val = PetscSqrtReal(1 - x * x); in func6()
47 else *val = PetscSqrtReal(x) / PetscSqrtReal(1 - x * x); in func7()
67 *val = PetscSqrtReal(PetscTanReal(x)); in func10()
81 else *val = PetscExpReal(1 - 1 / x) / PetscSqrtReal(x * x * x - x * x * x * x); in func12()
H A Dex7.c57 diffMat = PetscSqrtReal(diffMat); in CheckPullback()
58 normMat = PetscSqrtReal(normMat); in CheckPullback()
195 diff = PetscSqrtReal(PetscSqr(wvcheck - wv)); in main()
221 diff = PetscSqrtReal(PetscSqr(wvcheck - wv)); in main()
345 diff = PetscSqrtReal(diff); in main()
346 norm = PetscSqrtReal(norm); in main()
416 diffMat = PetscSqrtReal(diffMat); in main()
417 normMat = PetscSqrtReal(normMat); in main()
430 diffMat = PetscSqrtReal(diffMat); in main()
431 normMat = PetscSqrtReal(normMat); in main()
[all …]
/petsc/src/tao/unconstrained/tutorials/
H A Dminsurf2.c226 f2 = PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction()
227 f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormFunction()
238 ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormFunction()
245 ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormFunction()
252 ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction()
259 ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction()
265 ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2); in FormFunction()
270 ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2); in FormFunction()
385 f1 = PetscSqrtReal(1.0 + d1 * d1 + d7 * d7); in FormFunctionGradient()
386 f2 = PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient()
[all …]
/petsc/src/ksp/ksp/impls/cg/gltr/
H A Dgltr.c143 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
145 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_GLTR()
174 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
176 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_GLTR()
198 cg->norm_r[0] = PetscSqrtReal(rz); /* norm_r = |r|_M */ in KSPCGSolve_GLTR()
206 norm_r = PetscSqrtReal(rr); /* norm_r = |r| */ in KSPCGSolve_GLTR()
243 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
245 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_GLTR()
298 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_GLTR()
310 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
[all …]
/petsc/src/ts/utils/dmplexlandau/
H A Dland_tensors.h3 #define LANDAU_INVSQRT(q) (1. / PetscSqrtReal(q))
89 sqrt_1s = PetscSqrtReal(1. + s); in LandauTensor2D()
120 inorm = PetscSqrtReal(inorm2); in LandauTensor3D()
128 …#define GAMMA3(_x, _c02) PetscSqrtReal(1.0 + ((_x[0] * _x[0]) + (_x[1] * _x[1]) + (_x[2] * _x[2]))…
151 …fact = -rsq / (g1 * g2 * PetscSqrtReal(wsq)); /* flip sign. papers use du/dt = C, PETSc uses form … in LandauTensor3DRelativistic()
153 …for (j = 0; j < 3; ++j) U[i][j] = fact * (-diff[i] * diff[j] / wsq + (PetscSqrtReal(rsq) - 1.) * (… in LandauTensor3DRelativistic()
/petsc/src/snes/impls/gs/
H A Dgssecant.c101 if (k == 0) ft1 = PetscSqrtReal(ft); in SNESComputeNGSDefaultSecant()
104 if (stol * PetscSqrtReal(xt) > PetscSqrtReal(dxt)) isdone = PETSC_TRUE; in SNESComputeNGSDefaultSecant()
105 if (PetscSqrtReal(ft) < atol) isdone = PETSC_TRUE; in SNESComputeNGSDefaultSecant()
106 if (rtol * ft1 > PetscSqrtReal(ft)) isdone = PETSC_TRUE; in SNESComputeNGSDefaultSecant()
/petsc/src/dm/dt/interface/
H A Ddtprob.c27 p[0] = PetscSqrtReal(2. / PETSC_PI) * PetscExpReal(-0.5 * PetscSqr(x[0])); in PetscPDFMaxwellBoltzmann1D()
115 p[0] = PetscSqrtReal(2. / PETSC_PI) * PetscSqr(x[0]) * PetscExpReal(-0.5 * PetscSqr(x[0])); in PetscPDFMaxwellBoltzmann3D()
137 …p[0] = PetscErfReal(x[0] / PETSC_SQRT2) - PetscSqrtReal(2. / PETSC_PI) * x[0] * PetscExpReal(-0.5 … in PetscCDFMaxwellBoltzmann3D()
160 …p[0] = PetscSqrtReal(1. / (2. * PETSC_PI)) * PetscExpReal(-0.5 * PetscSqr(x[0] / … in PetscPDFGaussian1D()
202 r = ck[0] * (PetscSqrtReal(PETSC_PI) / 2.) * q; in PetscPDFSampleGaussian1D()
211 r += (ck[k] / temp) * PetscPowReal((PetscSqrtReal(PETSC_PI) / 2.) * q, 2. * k + 1.); in PetscPDFSampleGaussian1D()
214 x[0] = PetscSqrtReal(PETSC_PI * 0.5) * r; in PetscPDFSampleGaussian1D()
259 const PetscReal mag = PetscSqrtReal(-2.0 * PetscLogReal(p[0])); in PetscPDFSampleGaussian2D()
283 …p[0] = (1. / PETSC_PI * PetscSqrtReal(PETSC_PI)) * PetscExpReal(-0.5 * (PetscSqr(x[0]) + PetscSqr(… in PetscPDFGaussian3D()
943 speed[p] = PetscSqrtReal(mag); in PetscProbComputeKSStatisticMagnitude()
/petsc/src/ts/tutorials/eimex/
H A Dallen_cahn.c187 x[i] = PetscTanhReal((x_map - 0.8) / (2. * PetscSqrtReal(user->param))); in FormInitialSolution()
189 x[i] = PetscTanhReal((0.613 - x_map) / (2. * PetscSqrtReal(user->param))); in FormInitialSolution()
191 x[i] = PetscTanhReal((x_map - 0.36) / (2. * PetscSqrtReal(user->param))); in FormInitialSolution()
193 x[i] = PetscTanhReal((0.2 - x_map) / (2. * PetscSqrtReal(user->param))); in FormInitialSolution()
195 x[i] = PetscTanhReal((x_map + 0.9) / (2. * PetscSqrtReal(user->param))); in FormInitialSolution()
/petsc/src/ksp/ksp/impls/cg/
H A Dcg.c182 dp = PetscSqrtReal(PetscAbsScalar(beta)); /* dp <- r'*z = r'*B*r = e'*A'*B*A*e */ in KSPSolve_CG()
247 e[i] = PetscSqrtReal(PetscAbsScalar(b)) / a; in KSPSolve_CG()
267 a = PetscSqrtReal(r2 / norm_p); in KSPSolve_CG()
273 a = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPSolve_CG()
290 if (eigs) d[i] = PetscSqrtReal(PetscAbsScalar(b)) * e[i] + 1.0 / a; in KSPSolve_CG()
297 dp = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPSolve_CG()
318 dp = PetscSqrtReal(PetscAbsScalar(beta)); in KSPSolve_CG()
420 dp = PetscSqrtReal(PetscAbsScalar(beta)); /* dp <- r'*z = r'*B*r = e'*A'*B*A*e */ in KSPSolve_CG_SingleReduction()
465 e[i] = PetscSqrtReal(PetscAbsScalar(b)) / a; in KSPSolve_CG_SingleReduction()
487 if (eigs) d[i] = PetscSqrtReal(PetscAbsScalar(b)) * e[i] + 1.0 / a; in KSPSolve_CG_SingleReduction()
[all …]
/petsc/src/ksp/ksp/impls/gmres/agmres/
H A Dagmresdeflation.c39 abs_pivot = PetscSqrtReal(pivot_r * pivot_r + pivot_i * pivot_i); in KSPAGMRESQuickSort()
43 abs_val = PetscSqrtReal(val_r[R] * val_r[R] + val_i[R] * val_i[R]); in KSPAGMRESQuickSort()
46 abs_val = PetscSqrtReal(val_r[R] * val_r[R] + val_i[R] * val_i[R]); in KSPAGMRESQuickSort()
54 abs_val = PetscSqrtReal(val_r[L] * val_r[L] + val_i[L] * val_i[L]); in KSPAGMRESQuickSort()
57 abs_val = PetscSqrtReal(val_r[L] * val_r[L] + val_i[L] * val_i[L]); in KSPAGMRESQuickSort()
H A Dagmresleja.c36 pd = pd * PetscSqrtReal(rd * rd + id * id); in KSPAGMRESLejaCfpdMax()
69 for (i = 0; i < m; i++) n_cmpl[i] = PetscSqrtReal(re[i] * re[i] + im[i] * im[i]); in KSPAGMRESLejaOrdering()
/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/vec/vec/impls/mpi/
H A Dpvecimpl.h175 … == NORM_FROBENIUS || type == NORM_1_AND_2) z[type == NORM_1_AND_2] = PetscSqrtReal(z[type == NORM… in VecNorm_MPI_Default()
204 *norm = PetscSqrtReal(loc[0]); in VecErrorWeightedNorms_MPI_Default()
205 *norma = PetscSqrtReal(loc[1]); in VecErrorWeightedNorms_MPI_Default()
206 *normr = PetscSqrtReal(loc[2]); in VecErrorWeightedNorms_MPI_Default()
/petsc/src/ksp/ksp/impls/tfqmr/
H A Dtfqmr.c78 if (!m) w = PetscSqrtReal(dp * dpold); in KSPSolve_TFQMR()
81 cm = 1.0 / PetscSqrtReal(1.0 + psi * psi); in KSPSolve_TFQMR()
92 dpest = PetscSqrtReal(2 * i + m + 2.0) * tau; in KSPSolve_TFQMR()
/petsc/src/snes/tutorials/
H A Dex34.c58 const PetscReal r = PetscSqrtReal(PetscSqr(x[0]) + PetscSqr(x[1])); in obstacle_ball()
59 const PetscReal psi_0 = PetscSqrtReal(1. - PetscSqr(r_0)); in obstacle_ball()
64 if (r < r_0) u[0] = PetscSqrtReal(1.0 - PetscSqr(r)); in obstacle_ball()
75 const PetscReal r = PetscSqrtReal(PetscSqr(x[0]) + PetscSqr(x[1])); in exactSol_ball()
H A Dex8.c21 obj[0] = PetscSqrtReal(PetscAbsScalar(1 + ux2 + uy2)); in objective_2d()
29 const PetscScalar v = 1. / PetscSqrtReal(PetscAbsScalar(1 + ux2 + uy2)); in gradient_1_2d()
41 const PetscScalar v1 = 1. / PetscSqrtReal(PetscAbsScalar(1 + ux2 + uy2)); in hessian_11_2d()
/petsc/src/ts/utils/dmplexlandau/tutorials/
H A Dex1.c112 shifts[0] = 0.5 * PetscSqrtReal(ctx->masses[0] / ctx->masses[1]); in SetMaxwellians()
114 …data[0].shift = ctx->electronShift * shifts[grid] * PetscSqrtReal(data[0].kT_m) / ctx->v_0; // shi… in SetMaxwellians()
142 …v_bar_ab = 1.8e-19 * PetscSqrtReal(m_cgs[0] * m_cgs[1]) * n_cm3[0] * ctx->lambdas[0][1] * PetscPow… in FormFunction()
148 if (AA < 0) ff = PetscAtanhReal(PetscSqrtReal(-AA)) / PetscSqrtReal(-AA); in FormFunction()
149 else ff = PetscAtanReal(PetscSqrtReal(AA)) / PetscSqrtReal(AA); in FormFunction()
152 …vTe = 2 * PetscSqrtReal(PETSC_PI / m_cgs[ii]) * PetscSqr(PetscSqr(e_cgs)) * n_cm3[0] * ctx->lambda… in FormFunction()
/petsc/src/dm/interface/
H A Ddmgeommodel.c20 norm = PetscSqrtReal(norm); in DMSnapToGeomModelBall()
40 norm = PetscSqrtReal(norm); in DMSnapToGeomModelCylinder()
/petsc/src/snes/linesearch/impls/secant/
H A Dlinesearchsecant.c113 …ambda_mid, (double)lambda_old, (double)PetscSqrtReal(fnrm), (double)PetscSqrtReal(fnrm_mid), (doub… in SNESLineSearchApply_Secant()

12345678910