Home
last modified time | relevance | path

Searched refs:PetscPowScalar (Results 1 – 25 of 39) sorted by relevance

12

/petsc/src/snes/tests/
H A Dex20.c176 dw = PetscPowScalar(aw, beta); in FormFunction()
181 de = PetscPowScalar(ae, beta); in FormFunction()
186 ds = PetscPowScalar(as, beta); in FormFunction()
191 dn = PetscPowScalar(an, beta); in FormFunction()
196 dd = PetscPowScalar(ad, beta); in FormFunction()
201 du = PetscPowScalar(au, beta); in FormFunction()
208 dw = PetscPowScalar(aw, beta); in FormFunction()
213 de = PetscPowScalar(ae, beta); in FormFunction()
219 ds = PetscPowScalar(as, beta); in FormFunction()
228 dn = PetscPowScalar(an, beta); in FormFunction()
[all …]
H A Dex5.c134 v = 6.0 * xp + PetscPowScalar(xp + 1.e-12, 6.0); /* +1.e-12 is to prevent 0^6 */ in main()
H A Dex7.c52 v = 6.0 * xp + PetscPowScalar(xp + 1.e-12, 6.0); /* +1.e-12 is to prevent 0^6 */ in main()
/petsc/src/snes/tutorials/
H A Dex18.c173 dw = PetscPowScalar(aw, beta); in FormFunction()
178 de = PetscPowScalar(ae, beta); in FormFunction()
183 ds = PetscPowScalar(as, beta); in FormFunction()
188 dn = PetscPowScalar(an, beta); in FormFunction()
195 dw = PetscPowScalar(aw, beta); in FormFunction()
200 de = PetscPowScalar(ae, beta); in FormFunction()
206 ds = PetscPowScalar(as, beta); in FormFunction()
213 dn = PetscPowScalar(an, beta); in FormFunction()
221 dw = PetscPowScalar(aw, beta); in FormFunction()
226 de = PetscPowScalar(ae, beta); in FormFunction()
[all …]
H A Dex59.c94 … 1.) * (b - a) * PetscPowScalar(xp, k - 2.) + SQR(a * xp) + SQR(b - a) * PetscPowScalar(xp, 2. * k… in main()
96 v2 = a * xp + (b - a) * PetscPowScalar(xp, k); in main()
H A Dex33.c91 PetscReal lambdaWet = kappaWet * PetscRealPart(PetscPowScalar(u[i].s, alpha)); in FormFunctionLocal()
92 …PetscReal lambda = lambdaWet + kappaNoWet * PetscRealPart(PetscPowScalar(1. - u[i].s, beta)); in FormFunctionLocal()
93 PetscReal lambdaWetL = kappaWet * PetscRealPart(PetscPowScalar(u[i - 1].s, alpha)); in FormFunctionLocal()
94 …PetscReal lambdaL = lambdaWetL + kappaNoWet * PetscRealPart(PetscPowScalar(1. - u[i - 1].s, b… in FormFunctionLocal()
H A Dex2.c133 v = 6.0 * xp + PetscPowScalar(xp + 1.e-12, 6.0); /* +1.e-12 is to prevent 0^6 */ in main()
H A Dex6.c117 v = 6.0 * xp + PetscPowScalar(xp + 1.e-12, 6.0); /* +1.e-12 is to prevent 0^6 */ in main()
/petsc/src/snes/tutorials/network/water/
H A Dwaterreaddata.c19 …f[0] += x[0] - x[1] * PetscPowScalar(flow[i], x[2]) - head[i]; … in PumpHeadCurveResidual()
20 …f[1] += (x[0] - x[1] * PetscPowScalar(flow[i], x[2]) - head[i]) * -1 * PetscPowScalar(flow[i], x[2… in PumpHeadCurveResidual()
21 …f[2] += (x[0] - x[1] * PetscPowScalar(flow[i], x[2]) - head[i]) * -1 * x[1] * x[2] * PetscPowScala… in PumpHeadCurveResidual()
265 …pipe->k = 4.72 * pipe->length / (PetscPowScalar(pipe->roughness, pipe->n) * PetscPowScalar(0.08333… in WaterReadData()
H A Dwaterfunctions.c10 flow_pipe = PetscSign(hf - ht) * PetscPowScalar(PetscAbsScalar(hf - ht) / pipe->k, 1 / pipe->n); in Flow_Pipe()
17 …flow_pump = PetscSign(hf - ht + pump->h0) * PetscPowScalar(PetscAbsScalar(hf - ht + pump->h0) / pu… in Flow_Pump()
/petsc/src/ts/tutorials/power_grid/stability_9bus/
H A Dex9busadj.c300 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap… in ResidualFunction()
301 …for (k = 0; k < ld_nsegsq[i]; k++) QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq… in ResidualFunction()
668 PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap[k]); in ResidualJacobian()
669 …D_dVr += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vr * PetscPo… in ResidualJacobian()
670 …D_dVi += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * PetscPo… in ResidualJacobian()
673 QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq[k]); in ResidualJacobian()
674 …D_dVr += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vr * PetscPo… in ResidualJacobian()
675 …D_dVi += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * PetscPo… in ResidualJacobian()
H A Dex9busoptfd.c301 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap… in ResidualFunction()
302 …for (k = 0; k < ld_nsegsq[i]; k++) QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq… in ResidualFunction()
670 PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap[k]); in ResidualJacobian()
671 …D_dVr += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vr * PetscPo… in ResidualJacobian()
672 …D_dVi += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * PetscPo… in ResidualJacobian()
675 QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq[k]); in ResidualJacobian()
676 …D_dVr += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vr * PetscPo… in ResidualJacobian()
677 …D_dVi += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * PetscPo… in ResidualJacobian()
H A Dex9bus.c529 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap… in ResidualFunction()
530 …for (k = 0; k < ld_nsegsq[i]; k++) QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq… in ResidualFunction()
940 PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap[k]); in ResidualJacobian()
941 …D_dVr += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vr * PetscPo… in ResidualJacobian()
942 …D_dVi += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * PetscPo… in ResidualJacobian()
945 QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq[k]); in ResidualJacobian()
946 …D_dVr += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vr * PetscPo… in ResidualJacobian()
947 …D_dVi += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * PetscPo… in ResidualJacobian()
H A Dex9busopt.c446 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap… in ResidualFunction()
447 …for (k = 0; k < ld_nsegsq[i]; k++) QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq… in ResidualFunction()
812 PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap[k]); in ResidualJacobian()
813 …D_dVr += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vr * PetscPo… in ResidualJacobian()
814 …D_dVi += ld_alphap[k] * ld_betap[k] * PD0[i] * PetscPowScalar(1 / Vm0, ld_betap[k]) * Vi * PetscPo… in ResidualJacobian()
817 QD += ld_alphaq[k] * QD0[i] * PetscPowScalar(Vm / Vm0, ld_betaq[k]); in ResidualJacobian()
818 …D_dVr += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vr * PetscPo… in ResidualJacobian()
819 …D_dVi += ld_alphaq[k] * ld_betaq[k] * QD0[i] * PetscPowScalar(1 / Vm0, ld_betaq[k]) * Vi * PetscPo… in ResidualJacobian()
H A Dex9busdmnetwork.c714 … for (k = 0; k < ld_nsegsp; k++) PD += ld_alphap[k] * PD0 * PetscPowScalar(Vm / Vm0, ld_betap[k]); in FormIFunction()
715 … for (k = 0; k < ld_nsegsq; k++) QD += ld_alphaq[k] * QD0 * PetscPowScalar(Vm / Vm0, ld_betaq[k]); in FormIFunction()
925 … for (k = 0; k < ld_nsegsp; k++) PD += ld_alphap[k] * PD0 * PetscPowScalar(Vm / Vm0, ld_betap[k]); in AlgFunction()
926 … for (k = 0; k < ld_nsegsq; k++) QD += ld_alphaq[k] * QD0 * PetscPowScalar(Vm / Vm0, ld_betaq[k]); in AlgFunction()
/petsc/src/ts/tutorials/
H A Dex11.h35 #define PetscPowScalar(a, b) pow(a, b) macro
445 rarepr2 = gascr4 * 2. * cr * (1. - PetscPowScalar(pst, d__1)); in riem1mdt()
449 rarepr1 = gascl4 * 2. * cl * (1. - PetscPowScalar(pst, d__1)); in riem1mdt()
475 *rstarl = *rl * PetscPowScalar(d__1, d__2); in riem1mdt()
481 *rstarr = *rr * PetscPowScalar(d__1, d__2); in riem1mdt()
504 *rstarr = *rr * PetscPowScalar(d__1, d__2); in riem1mdt()
544 *rstarl = *rl * PetscPowScalar(d__1, d__2); in riem1mdt()
631 rstarr = r0 * PetscPowScalar(d__1, d__2); in riemannsolver()
647 s = p0 / PetscPowScalar(r0, *gam); in riemannsolver()
649 *rx = PetscPowScalar(d__1, gasc4); in riemannsolver()
H A Dex31.c602 fac = PetscPowScalar(y[0] * y[0] + y[1] * y[1], -1.5); in IJacobian_Hull1972B4()
603 fac2 = PetscPowScalar(y[0] * y[0] + y[1] * y[1], -0.5); in IJacobian_Hull1972B4()
690 f[0] = 2. * t * PetscPowScalar(y[1], 1. / 5.) * y[3]; in RHSFunction_Kulikov2013I()
708 m1 = (2. * t * y[3]) / (5. * PetscPowScalar(y[1], 4. / 5.)); in RHSJacobian_Kulikov2013I()
709 m2 = 2. * t * PetscPowScalar(y[1], 1. / 5.); in RHSJacobian_Kulikov2013I()
743 f[0] = 2. * t * PetscPowScalar(y[1], 1. / 5.) * y[3]; in IFunction_Kulikov2013I()
763 m1 = (2. * t * y[3]) / (5. * PetscPowScalar(y[1], 4. / 5.)); in IJacobian_Kulikov2013I()
764 m2 = 2. * t * PetscPowScalar(y[1], 1. / 5.); in IJacobian_Kulikov2013I()
H A Dex30.c161 const PetscScalar fnorm = PetscPowScalar(norm, nexp); in C_0()
184 const PetscScalar fnorm = PetscPowScalar(norm, nexp); in C_0_3d()
202 const PetscScalar fnorm = PetscPowScalar(norm, nexp); in JC_0_c0c0()
203 const PetscScalar dfnorm = nexp * PetscPowScalar(norm, nexp - 1.0); in JC_0_c0c0()
229 const PetscScalar fnorm = PetscPowScalar(norm, nexp); in JC_0_c0c0_3d()
230 const PetscScalar dfnorm = nexp * PetscPowScalar(norm, nexp - 1.0); in JC_0_c0c0_3d()
633 const PetscScalar t2 = alpha / gamma * PetscPowScalar(normC, nexp); in energy()
659 const PetscScalar t2 = alpha / gamma * PetscPowScalar(normC, nexp); in energy()
/petsc/src/tao/linesearch/impls/morethuente/
H A Dmorethuente.c390 gamma1 = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s)); in Tao_mcstep()
413 gamma1 = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s)); in Tao_mcstep()
442 …gamma1 = s * PetscSqrtScalar(PetscMax(0.0, PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s))… in Tao_mcstep()
471 gamma1 = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dy / s) * (*dp / s)); in Tao_mcstep()
/petsc/src/ts/tutorials/power_grid/
H A Dex7.c187 …ho))) * PetscExpScalar(-0.5 / (1 - rho * rho) * (PetscPowScalar((xi - mux) / sigmax, 2) + PetscPow… in ini_bou()
281 …user->disper_coe = PetscPowScalar((user->lambda * user->ws) / (2 * user->H), 2) * user->q * (1.0 -… in IFunction()
331 …c5 = (PetscPowScalar((user->lambda * user->ws) / (2 * user->H), 2) * user->q * (1.0 - Pet… in IJacobian()
H A Dex6.c153 …ho))) * PetscExpScalar(-0.5 / (1 - rho * rho) * (PetscPowScalar((xi - mux) / sigmax, 2) + PetscPow… in ini_bou()
282 …user->disper_coe = PetscPowScalar((user->lambda * user->ws) / (2 * user->H), 2) * user->q * (1.0 -… in IFunction()
436 …c5 = (PetscPowScalar((user->lambda * user->ws) / (2 * user->H), 2) * user->q * (1.0 - Petsc… in IJacobian()
H A Dex5.c126 for (i = 0; i < user->nsamples; i++) x[i] = PetscPowScalar(x[i], 1 / user->kw); in WindSpeeds()
/petsc/src/sys/tests/
H A Dex25.c86 b = PetscPowScalar(a, c); in main()
/petsc/include/
H A Dpetscmath.h519 #define PetscPowScalar(a, b) PetscPowComplex(a, b) macro
543 #define PetscPowScalar(a, b) PetscPowReal(a, b) macro
1106 return PetscPowScalar(base, cpower); in PetscPowScalarReal()
/petsc/src/ts/tests/
H A Dex5.c415 *pottemp = temp * (PetscPowScalar(pressure1 / pavg, kdry)); /* calculates potential temperature */ in potential_temperature()
424 …e = 6.11 * (PetscPowScalar(10, (7.5 * dtemp) / (237.7 + dtemp))); /* converts from dew poin… in calcmixingr()

12