Home
last modified time | relevance | path

Searched refs:du (Results 1 – 24 of 24) sorted by relevance

/petsc/src/ts/tutorials/power_grid/
H A Dex3.c43 PetscReal du[2] = {0.0, 0.0}; in main() local
101 PetscCall(PetscOptionsRealArray("-du", "Perturbation in initial solution", "", du, &n, &flg2)); in main()
102 u[0] += du[0]; in main()
103 u[1] += du[1]; in main()
145 for (du[1] = -2.5; du[1] <= .01; du[1] += .1) { in main()
149 u[0] += du[0]; in main()
150 u[1] += du[1]; in main()
H A Dex9.c101 PetscReal du[2] = {0.0, 0.0}; in main() local
156 PetscCall(PetscOptionsRealArray("-du", "Perturbation in initial solution", "", du, &n, &flg2)); in main()
157 u[0] += du[0]; in main()
158 u[1] += du[1]; in main()
193 for (du[1] = -2.5; du[1] <= .01; du[1] += .1) { in main()
197 u[0] += du[0]; in main()
198 u[1] += du[1]; in main()
H A Dex2.c122 PetscReal du[2] = {0.0, 0.0}; in main() local
176 PetscCall(PetscOptionsRealArray("-du", "Perturbation in initial solution", "", du, &n, &flg2)); in main()
177 u[0] += du[0]; in main()
178 u[1] += du[1]; in main()
214 for (du[1] = -2.5; du[1] <= .01; du[1] += .1) { in main()
218 u[0] += du[0]; in main()
219 u[1] += du[1]; in main()
H A Dex3sa.c32 PetscReal du[2] = {0.0, 0.0}; in main() local
107 PetscCall(PetscOptionsRealArray("-du", "Perturbation in initial solution", "", du, &n, &flg2)); in main()
108 u[0] += du[0]; in main()
109 u[1] += du[1]; in main()
183 for (du[1] = -2.5; du[1] <= .01; du[1] += .1) { in main()
187 u[0] += du[0]; in main()
188 u[1] += du[1]; in main()
H A Dex9adj.c188 PetscReal du[2] = {0.0, 0.0}; in main() local
260 PetscCall(PetscOptionsRealArray("-du", "Perturbation in initial solution", "", du, &n, &flg2)); in main()
261 u[0] += du[0]; in main()
262 u[1] += du[1]; in main()
321 for (du[1] = -2.5; du[1] <= .01; du[1] += .1) { in main()
325 u[0] += du[0]; in main()
326 u[1] += du[1]; in main()
/petsc/src/snes/tests/
H A Dex20.c139 PetscScalar tleft, tright, beta, td, ad, dd, fd = 0.0, tu, au, du = 0.0, fu = 0.0; in FormFunction() local
201 du = PetscPowScalar(au, beta); in FormFunction()
202 fu = du * (tu - t0); in FormFunction()
246 du = PetscPowScalar(au, beta); in FormFunction()
247 fu = du * (tu - t0); in FormFunction()
294 du = PetscPowScalar(au, beta); in FormFunction()
295 fu = du * (tu - t0); in FormFunction()
331 du = PetscPowScalar(au, beta); in FormFunction()
332 fu = du * (tu - t0); in FormFunction()
368 du = PetscPowScalar(au, beta); in FormFunction()
[all …]
/petsc/src/snes/tutorials/
H A Dex48.c710 … PetscScalar *PETSC_RESTRICT u, PetscScalar *PETSC_RESTRICT v, PetscScalar du[PETSC_RESTRICT], Pet… in PointwiseNonlinearity()
715 du[0] = du[1] = du[2] = 0; in PointwiseNonlinearity()
723 du[ll] += dphi[l][ll] * n[l].u; in PointwiseNonlinearity()
727 …gam = PetscSqr(du[0]) + PetscSqr(dv[1]) + du[0] * dv[1] + 0.25 * PetscSqr(du[1] + dv[0]) + 0.25 * … in PointwiseNonlinearity()
731 …l phi[], PetscReal dphi[4][2], PetscScalar *u, PetscScalar *v, PetscScalar du[], PetscScalar dv[],… in PointwiseNonlinearity2D() argument
736 du[0] = du[1] = 0; in PointwiseNonlinearity2D()
744 du[ll] += dphi[l][ll] * n[l].u; in PointwiseNonlinearity2D()
748 gam = PetscSqr(du[0]) + PetscSqr(dv[1]) + du[0] * dv[1] + 0.25 * PetscSqr(du[1] + dv[0]); in PointwiseNonlinearity2D()
792 PetscScalar du[3], dv[3], u, v; in THIFunctionLocal() local
795 PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta); in THIFunctionLocal()
[all …]
/petsc/src/ts/utils/
H A Ddmplexts.c407 Vec du, uhat, uhat_t, r, rhat, df; in DMTSCheckJacobian() local
415 PetscCall(VecDuplicate(u, &du)); in DMTSCheckJacobian()
416 PetscCall(VecSetRandom(du, rand)); in DMTSCheckJacobian()
419 PetscCall(MatMult(J, du, df)); in DMTSCheckJacobian()
430 PetscCall(VecWAXPY(uhat, h, du, u)); in DMTSCheckJacobian()
431 PetscCall(VecWAXPY(uhat_t, h * shift, du, u_t)); in DMTSCheckJacobian()
445 PetscCall(VecDestroy(&du)); in DMTSCheckJacobian()
/petsc/include/
H A Dpetscfe.h178 const CeedScalar *u = in[0], *du = in[1], *qdata = in[2]; \
192 …edScalar u_x[2] = {invJ[0][0] * du[i + Q * 0] + invJ[1][0] * du[i + Q * 1], invJ[0][1] * du[i + Q …
/petsc/src/ksp/ksp/impls/cg/pipecgrr/
H A Dpipecgrr.c25 …PetscReal ds = 0.0, dz = 0.0, dx = 0.0, dpp = 0.0, dq = 0.0, dm = 0.0, du = 0.0, dw = 0.0, db = … in KSPSolve_PIPECGRR() local
112 unp = du; in KSPSolve_PIPECGRR()
135 PetscCall(VecNormBegin(U, NORM_2, &du)); in KSPSolve_PIPECGRR()
158 PetscCall(VecNormEnd(U, NORM_2, &du)); in KSPSolve_PIPECGRR()
211 errw = PetscSqrtReal(mnz * sqn * Anorm * du) * eps; in KSPSolve_PIPECGRR()
/petsc/src/binding/petsc4py/test/
H A Dtest_ts_py.py18 def function(self, ts, t, u, du, F): argument
21 f = du + u * u
24 def jacobian(self, ts, t, u, du, a, J, P): argument
H A Dtest_ts.py40 def ifunction(self, ts, t, u, du, F): argument
43 f = du + u * u
46 def ijacobian(self, ts, t, u, du, a, J, P): argument
/petsc/src/ts/tutorials/
H A Dex14.c729 …trict 8][3], PetscScalar *restrict u, PetscScalar *restrict v, PetscScalar du[restrict 3], PetscSc… in PointwiseNonlinearity()
734 du[0] = du[1] = du[2] = 0; in PointwiseNonlinearity()
742 du[ll] += dphi[l][ll] * n[l].u; in PointwiseNonlinearity()
746 …gam = Sqr(du[0]) + Sqr(dv[1]) + du[0] * dv[1] + 0.25 * Sqr(du[1] + dv[0]) + 0.25 * Sqr(du[2]) + 0.… in PointwiseNonlinearity()
790 PetscScalar du[3], dv[3], u, v, udot = 0, vdot = 0; in THIFunctionLocal_3D() local
797 PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta); in THIFunctionLocal_3D()
804 …n[l]->u += dp[0] * jw * eta * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp[… in THIFunctionLocal_3D()
805 …fn[l]->v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp… in THIFunctionLocal_3D()
1129 PetscScalar du[3], dv[3], u, v; in THIJacobianLocal_Momentum() local
1132 PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta); in THIJacobianLocal_Momentum()
[all …]
H A Dex9.c825 PetscScalar ubar, du[2], a[2], fL[2], fR[2], lam[2], ustar[2], R[2][2]; in PhysicsRiemann_IsoGas_Roe() local
831 du[0] = uR[0] - uL[0]; in PhysicsRiemann_IsoGas_Roe()
832 du[1] = uR[1] - uL[1]; in PhysicsRiemann_IsoGas_Roe()
833 a[0] = (1 / (2 * c)) * ((ubar + c) * du[0] - du[1]); in PhysicsRiemann_IsoGas_Roe()
834 a[1] = (1 / (2 * c)) * ((-ubar + c) * du[0] + du[1]); in PhysicsRiemann_IsoGas_Roe()
/petsc/src/benchmarks/
H A Dssget177 SIZE=($(du -b ${PATH_INFO[5]}))
185 SIZE=($(du -b ${PATH_INFO[4]}))
/petsc/systems/Apple/OSX/bin/
H A Dmakedmg26 let size=`du -c tempsource | tail -1 | cut -f1`
/petsc/systems/Apple/iOS/bin/
H A Dmakedmg31 let size=`du -c tempsource | tail -1 | cut -f1`
/petsc/src/snes/utils/
H A Ddmplexsnes.c943 Vec du, uhat, r, rhat, df; in DMSNESCheckJacobian() local
951 PetscCall(VecDuplicate(u, &du)); in DMSNESCheckJacobian()
952 PetscCall(VecSetRandom(du, rand)); in DMSNESCheckJacobian()
955 PetscCall(MatMult(J, du, df)); in DMSNESCheckJacobian()
965 PetscCall(VecWAXPY(uhat, h, du, u)); in DMSNESCheckJacobian()
978 PetscCall(VecDestroy(&du)); in DMSNESCheckJacobian()
/petsc/src/binding/petsc4py/demo/legacy/ode/
H A Dvanderpol.py276 def compute_derp(du, dp): argument
278 du[1] * (-10.0 / (81.0 * mu_ * mu_) + 2.0 * 292.0 / (2187.0 * mu_ * mu_ * mu_))
/petsc/src/dm/dt/fv/interface/
H A Dfv.c2125 PetscScalar du;
2129 du = cg1->centroid[0] + 3*cg1->centroid[1] - (cg->centroid[0] + 3*cg->centroid[1]);
2130 grad[0] += fg->grad[!i][0] * du;
2131 grad[1] += fg->grad[!i][1] * du;
/petsc/doc/manual/
H A Dtao.md2198 A_k du = -g_k
2201 to obtain a direction $du$. We need a direction that provides
2208 That is, we require $g_k^T A_k du < 0$.
2211 parameter $\rho_k$ so that $du$ is also a sufficient descent
2217 \displaystyle \min_{\alpha \geq 0} \; \tilde{f}_k(u_k + \alpha du, v_k).
2225 u_{k+\frac{1}{2}} & = & u_k + \alpha_k du \\
2238 $du = -A_k^T g_k$ during the search procedure. However, the
2258 \displaystyle \min_{du,dv} & \tilde{f}_k(u_k+du, v_k+dv) \\
2259 \text{subject to} & A_k du + B_k dv + \alpha_k g_k = 0
2266 du = -A_k^{-1}(B_k dv + \alpha_k g_k).
[all …]
H A Dts.md330 &= \frac{dh}{du} \dot{u} + \frac{\partial h}{\partial t} \\
331 &= \frac{dh}{du} f(t, u, z) + \frac{\partial h}{\partial t}
335 If the Jacobian $\frac{dh}{du} \frac{df}{dz}$ is non-singular, then we have precisely a semi-explic…
/petsc/doc/faq/
H A Dindex.md329 e.g. `find . -name output -type d | xargs du -sh | sort -hr` on a Unix-based system.
/petsc/doc/
H A Dpetsc.bib3325 journal = {Oil and Gas Science and Technology-Revue de l'IFP-Institut Francais du Petrole},
27832 @Article{ du.yu.ea:multileaf,
31427 school = {Universit\'{e} de Montpellier II -- Sciences et Techniques du Languedoc},
35245 title = {{R}\'{e}solution Num\'{e}rique Approach\'{e}e du Probl\`{e}me de Programmation
36664 journal = {Revue d' Economie du developpement},