| /petsc/src/ts/tutorials/network/ |
| H A D | pipeImpls.c | 38 static inline PetscScalar dqdx(PipeField *x, PetscInt i, PetscInt ilast, PetscReal dx) in dqdx() argument
41 return (x[i + 1].q - x[i].q) / dx; in dqdx()
43 return (x[i].q - x[i - 1].q) / dx; in dqdx()
45 return (x[i + 1].q - x[i - 1].q) / (2 * dx); in dqdx()
49 static inline PetscScalar dhdx(PipeField *x, PetscInt i, PetscInt ilast, PetscReal dx) in dhdx() argument
52 return (x[i + 1].h - x[i].h) / dx; in dhdx()
54 return (x[i].h - x[i - 1].h) / dx; in dhdx()
56 return (x[i + 1].h - x[i - 1].h) / (2 * dx); in dhdx()
64 PetscReal dx = pipe->length / (info->mx - 1), dt = pipe->dt; in PipeIFunctionLocal_Lax() local
82 …f[2 * (i - 1) + 2] = (x[i].q - xold_i) + dt * (GRAV * pipe->A * dhdx(xold, i, ilast, dx) + pipe->R… in PipeIFunctionLocal_Lax()
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex9.c | 123 PetscReal **au, dx, dy, x, y; in FormExactSolution() local
126 dx = 4.0 / (PetscReal)(info->mx - 1); in FormExactSolution()
132 x = -2.0 + i * dx; in FormExactSolution()
145 PetscReal **aXl, dx, dy, x, y; in FormBounds() local
150 dx = 4.0 / (PetscReal)(info.mx - 1); in FormBounds()
156 x = -2.0 + i * dx; in FormBounds()
168 PetscReal dx, dy, x, y, ue, un, us, uw; in FormFunctionLocal() local
171 dx = 4.0 / (PetscReal)(info->mx - 1); in FormFunctionLocal()
176 x = -2.0 + i * dx; in FormFunctionLocal()
180 uw = (i - 1 == 0) ? u_exact(x - dx, y) : au[j][i - 1]; in FormFunctionLocal()
[all …]
|
| H A D | ex30.c | 98 PetscReal dx, dz; member
313 x = (i - grid->jlid) * grid->dx; in HorizVelocity()
330 x = (i - grid->jlid - 0.5) * grid->dx; in VertVelocity()
346 x = (i - grid->jlid - 0.5) * grid->dx; in Pressure()
412 eps11 = (uE - uW) / grid->dx; in CalcSecInv()
414 eps12 = 0.5 * ((uN - uS) / grid->dz + (wE - wW) / grid->dx); in CalcSecInv()
466 PetscScalar dx = grid->dx, dz = grid->dz; in XMomentumResidual() local
492 dPdx = (x[j][i + 1].p - x[j][i].p) / dx; in XMomentumResidual()
493 if (j == jlim) dudzN = etaN * (x[j][i].w - x[j][i + 1].w) / dx; in XMomentumResidual()
496 dudxE = etaE * (x[j][i + 1].u - x[j][i].u) / dx; in XMomentumResidual()
[all …]
|
| H A D | ex33.c | 69 PetscReal dx = 1.0 / (PetscReal)(info->mx - 1); in FormFunctionLocal() local
90 PetscScalar K = 2 * dx / (dx / Kappa[i] + dx / Kappa[i - 1]); in FormFunctionLocal()
96 …f[i].s = phi * (u[i].s - uold[i].s) + (dt / dx) * ((lambdaWet / lambda) * u[i].v - (lambdaWetL / l… in FormFunctionLocal()
98 f[i].v = u[i].v + K * lambda * (u[i].p - u[i - 1].p) / dx; in FormFunctionLocal()
|
| /petsc/src/snes/interface/ |
| H A D | snesj.c | 52 PetscScalar dx, *y, wscale; in SNESComputeJacobianDefault() local
112 if (use_wp) dx = PetscSqrtReal(1.0 + unorm); in SNESComputeJacobianDefault()
113 else dx = xx[i - start]; in SNESComputeJacobianDefault()
115 if (PetscAbsScalar(dx) < dx_min) dx = (PetscRealPart(dx) < 0. ? -1. : 1.) * dx_par; in SNESComputeJacobianDefault()
116 dx *= epsilon; in SNESComputeJacobianDefault()
117 wscale = 1.0 / dx; in SNESComputeJacobianDefault()
118 if (x2->ops->setvalues) PetscCall(VecSetValues(x2, 1, &i, &dx, ADD_VALUES)); in SNESComputeJacobianDefault()
121 y[i - start] += dx; in SNESComputeJacobianDefault()
|
| H A D | snesob.c | 144 PetscReal ob, ob1, ob2, ob3, fob, dx, eps = 1e-6; in SNESObjectiveComputeFunctionDefaultFD() local
160 if (fob > 0.) dx = 1e-6 * fob; in SNESObjectiveComputeFunctionDefaultFD()
161 else dx = 1e-6; in SNESObjectiveComputeFunctionDefaultFD()
167 xv = dx; in SNESObjectiveComputeFunctionDefaultFD()
177 xv = 2. * dx; in SNESObjectiveComputeFunctionDefaultFD()
187 xv = -dx; in SNESObjectiveComputeFunctionDefaultFD()
196 fv = (-ob2 + 6. * ob1 - 3. * ob - 2. * ob3) / (6. * dx); in SNESObjectiveComputeFunctionDefaultFD()
|
| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | land_tensors.h | 113 PetscReal dx[3], inorm3, inorm, inorm2, norm2, x2[] = {xp, yp, zp}; in LandauTensor3D() local
116 dx[d] = x2[d] - x1[d]; in LandauTensor3D()
117 norm2 += dx[d] * dx[d]; in LandauTensor3D()
122 for (d = 0; d < 3; ++d) U[d][d] = -(inorm - inorm3 * dx[d] * dx[d]); in LandauTensor3D()
123 U[1][0] = U[0][1] = inorm3 * dx[0] * dx[1]; in LandauTensor3D()
124 U[1][2] = U[2][1] = inorm3 * dx[2] * dx[1]; in LandauTensor3D()
125 U[2][0] = U[0][2] = inorm3 * dx[0] * dx[2]; in LandauTensor3D()
|
| /petsc/src/tao/linesearch/impls/morethuente/ |
| H A D | morethuente.c | 10 static PetscErrorCode Tao_mcstep(TaoLineSearch ls, PetscReal *stx, PetscReal *fx, PetscReal *dx, Pe…
363 static PetscErrorCode Tao_mcstep(TaoLineSearch ls, PetscReal *stx, PetscReal *fx, PetscReal *dx, Pe… in Tao_mcstep() argument
373 …PetscCheck(*dx * (*stp - *stx) < 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "dx * (stp-stx) >… in Tao_mcstep()
377 sgnd = *dp * (*dx / PetscAbsReal(*dx)); in Tao_mcstep()
387 theta = 3 * (*fx - *fp) / (*stp - *stx) + *dx + *dp; in Tao_mcstep()
388 s = PetscMax(PetscAbsReal(theta), PetscAbsReal(*dx)); in Tao_mcstep()
390 gamma1 = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s)); in Tao_mcstep()
393 p = (gamma1 - *dx) + theta; in Tao_mcstep()
394 q = ((gamma1 - *dx) + gamma1) + *dp; in Tao_mcstep()
397 stpq = *stx + ((*dx / ((*fx - *fp) / (*stp - *stx) + *dx)) * 0.5) * (*stp - *stx); in Tao_mcstep()
[all …]
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex11.c | 9 PetscScalar *dx, *grad; in main() local
21 PetscCall(PetscMalloc2(dim * numFaces, &dx, dim * numFaces, &grad)); in main()
27 PetscCall(PetscFVComputeGradient(fvm, numFaces, dx, grad)); in main()
30 PetscCall(PetscFree2(dx, grad)); in main()
|
| /petsc/src/ts/tests/ |
| H A D | ex4.c | 23 PetscReal dx; /* the grid space in x-direction */ member
64 data.dx = 1.0 / (data.m + 1.0); in main()
199 PetscReal x, y, dx, dy; in Initial() local
206 dx = data->dx; in Initial()
219 x = dx * row; in Initial()
277 PetscReal dx, dy, a, epsilon, xc, xl, xr, yl, yr; in RHSJacobian() local
283 dx = data->dx; in RHSJacobian()
288 xc = -2.0 * epsilon * (1.0 / (dx * dx) + 1.0 / (dy * dy)); in RHSJacobian()
289 xl = 0.5 * a / dx + epsilon / (dx * dx); in RHSJacobian()
290 xr = -0.5 * a / dx + epsilon / (dx * dx); in RHSJacobian()
[all …]
|
| /petsc/src/dm/tutorials/ |
| H A D | ex21.c | 182 PetscReal dx = 1.0 / 32.0; in pic_advect() local
186 min[0] = 0.5 * dx; in pic_advect()
187 max[0] = 0.5 * dx + 31.0 * dx; in pic_advect()
188 min[1] = 0.5 * dx; in pic_advect()
189 max[1] = 0.5 * dx; in pic_advect()
193 PetscReal dx = 1.0 / 32.0; in pic_advect() local
197 min[0] = 0.05 * dx; in pic_advect()
198 max[0] = 0.5 * dx; in pic_advect()
199 min[1] = 0.5 * dx; in pic_advect()
200 max[1] = 0.5 * dx + 31.0 * dx; in pic_advect()
|
| H A D | swarm_ex3.c | 17 PetscReal dx, dy; in _DMLocatePoints_DMDARegular_IS() local
31 dx = 2.0 / ((PetscReal)mx); in _DMLocatePoints_DMDARegular_IS()
42 mi = (PetscInt)((coorx - (-1.0)) / dx); in _DMLocatePoints_DMDARegular_IS()
132 PetscReal dx; in ex3_1() local
147 dx = 2.0 / ((PetscReal)mx); in ex3_1()
148 …(DMDASetUniformCoordinates(dmregular, -1.0 + 0.5 * dx, 1.0 - 0.5 * dx, -1.0 + 0.5 * dx, 1.0 - 0.5 … in ex3_1()
189 px = PetscRealPart(LA_coors[2 * p + 0]) + 0.1 * rx * dx; in ex3_1()
190 py = PetscRealPart(LA_coors[2 * p + 1]) + 0.1 * ry * dx; in ex3_1()
|
| H A D | ex10.c | 26 PetscScalar value = 1.0, dx, dy; in main() local
33 dx = (xp - xm) / (Nx - 1); in main()
59 …s; i < ixs + ixm; i++) gauss_ptr[j][i] = PetscExpScalar(-(xm + i * dx) * (xm + i * dx) - (ym + j *… in main()
|
| /petsc/src/ksp/ksp/utils/lmvm/tests/ |
| H A D | solve_performance.c | 11 Vec x, g, dx, df, p; in main() local
26 PetscCall(VecDuplicate(x, &dx)); in main()
47 PetscCall(VecSetRandom(dx, rand)); in main()
49 PetscCall(VecDot(dx, df, &dot)); in main()
59 PetscCall(VecAXPY(x, xscale, dx)); in main()
66 PetscCall(VecAXPY(x, xscale, dx)); in main()
81 PetscCall(VecDestroy(&dx)); in main()
|
| H A D | ex3.c | 51 Vec dx, df; in TestUsability() local
54 PetscCall(VecDuplicate(pattern, &dx)); in TestUsability()
56 PetscCall(VecSetRandom(dx, rand)); in TestUsability()
60 PetscCall(MatMult(B, dx, df)); in TestUsability()
61 PetscCall(MatSolve(B, dx, df)); in TestUsability()
65 PetscCall(VecCopy(dx, df)); in TestUsability()
66 PetscCall(MatLMVMUpdate(B, dx, df)); in TestUsability()
67 PetscCall(MatMult(B, dx, df)); in TestUsability()
68 PetscCall(MatSolve(B, dx, df)); in TestUsability()
72 PetscCall(VecDestroy(&dx)); in TestUsability()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex74.c | 121 PetscReal x, dx, a = ctxt->a, pi = PETSC_PI; in ExactSolution() local
124 dx = (ctxt->xmax - ctxt->xmin) / ((PetscReal)ctxt->imax); in ExactSolution()
128 x = i * dx; in ExactSolution()
148 PetscReal dx, dx2; in RHSJacobian() local
151 dx = (user->xmax - user->xmin) / ((PetscReal)user->imax); in RHSJacobian()
152 dx2 = dx * dx; in RHSJacobian()
164 values[0] = user->a * .5 / dx; in RHSJacobian()
166 values[2] = -user->a * .5 / dx; in RHSJacobian()
|
| /petsc/src/ts/tutorials/power_grid/ |
| H A D | ex6.c | 38 PetscScalar dx; /* x step size */ member
119 …WORLD, "sum(p)*dw*dtheta at t = %3.2f = %3.6f\n", (double)t, (double)(sum * user->dx * user->dy))); in PostStep()
140 user->dx = (user->xmax - user->xmin) / (M - 1); in ini_bou()
183 *p1 = f * (p[j][i + 1] - p[j][i - 1]) / (2 * user->dx); in adv1()
231 …f[j][i] = fwc * (p[j][i + 1] - p[j][i]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j +… in BoundaryConditions()
233 …f[j][i] = fwc * (p[j][i + 1] - p[j][i]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j][… in BoundaryConditions()
235 …f[j][i] = fwc * (p[j][i] - p[j][i - 1]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j +… in BoundaryConditions()
237 …f[j][i] = fwc * (p[j][i] - p[j][i - 1]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j][… in BoundaryConditions()
239 …f[j][i] = fwc * (p[j][i + 1] - p[j][i]) / (user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j… in BoundaryConditions()
241 …f[j][i] = fwc * (p[j][i] - p[j][i - 1]) / (user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j… in BoundaryConditions()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | fdmpiaij.c | 35 PetscScalar dx = 0.0, *w3_array, *dy_i, *dy = coloring->dy; in MatFDColoringApply_BAIJ() local
64 dx = 1.0 / (PetscSqrtReal(1.0 + unorm) * epsilon); in MatFDColoringApply_BAIJ()
69 dx = xx[col]; in MatFDColoringApply_BAIJ()
70 if (PetscAbsScalar(dx) < umin) { in MatFDColoringApply_BAIJ()
71 if (PetscRealPart(dx) >= 0.0) dx = umin; in MatFDColoringApply_BAIJ()
72 else if (PetscRealPart(dx) < 0.0) dx = -umin; in MatFDColoringApply_BAIJ()
74 dx *= epsilon; in MatFDColoringApply_BAIJ()
75 vscale_array[col] = 1.0 / dx; in MatFDColoringApply_BAIJ()
111 w3_array[col] += 1.0 / dx; in MatFDColoringApply_BAIJ()
112 if (i) w3_array[col - 1] -= 1.0 / dx; /* resume original w3[col-1] */ in MatFDColoringApply_BAIJ()
[all …]
|
| /petsc/doc/manual/ |
| H A D | dt.md | 11 … $P(x)$, so that the probability for an event at location in $[x, x+dx]$ is $P(x) dx$. This means …
14 \int_\Omega P(x) dx = 1.
17 …res that the PDF must have units which are the inverse of the volume form $dx$, meaning that it is…
27 \int_\Omega P(x) dx &= \int_\Omega P(\lambda s) \lambda^d ds \\
|
| /petsc/src/dm/impls/stag/tutorials/output/ |
| H A D | ex6_3.out | 2 dx: 0.0833333
5 V_p dt / dx: 0.0207846
|
| H A D | ex6_2.out | 2 dx: 0.0833333
5 V_p dt / dx: 0.0207846
|
| H A D | ex6_4.out | 2 dx: 0.0833333
5 V_p dt / dx: 0.0207846
|
| H A D | ex6_1.out | 2 dx: 0.0833333
5 V_p dt / dx: 0.0207846
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex10.c | 211 PetscReal K[60][60], x, y, z, dx, dy, dz; in Elastic20Stiff() local
219 dx = 2.0; in Elastic20Stiff()
225 xyz[1][0] = x + dx; in Elastic20Stiff()
228 xyz[2][0] = x + 2. * dx; in Elastic20Stiff()
234 xyz[4][0] = x + 2. * dx; in Elastic20Stiff()
240 xyz[6][0] = x + dx; in Elastic20Stiff()
243 xyz[7][0] = x + 2. * dx; in Elastic20Stiff()
249 xyz[9][0] = x + 2. * dx; in Elastic20Stiff()
255 xyz[11][0] = x + 2. * dx; in Elastic20Stiff()
261 xyz[13][0] = x + dx; in Elastic20Stiff()
[all …]
|
| /petsc/src/dm/impls/da/ |
| H A D | dageometry.c | 136 PetscReal gmin_l[2], gmax_l[2], dx[2]; in private_DMDALocatePointsIS_2D_Regular() local
152 dx[0] = (gmax_l[0] - gmin_l[0]) / ((PetscReal)mxlocal); in private_DMDALocatePointsIS_2D_Regular()
153 dx[1] = (gmax_l[1] - gmin_l[1]) / ((PetscReal)mylocal); in private_DMDALocatePointsIS_2D_Regular()
175 for (PetscInt d = 0; d < 2; d++) mi[d] = (PetscInt)((coor_p[d] - gmin[d]) / dx[d]); in private_DMDALocatePointsIS_2D_Regular()
199 PetscReal gmin_l[3], gmax_l[3], dx[3]; in private_DMDALocatePointsIS_3D_Regular() local
218 dx[0] = (gmax_l[0] - gmin_l[0]) / ((PetscReal)mxlocal); in private_DMDALocatePointsIS_3D_Regular()
219 dx[1] = (gmax_l[1] - gmin_l[1]) / ((PetscReal)mylocal); in private_DMDALocatePointsIS_3D_Regular()
220 dx[2] = (gmax_l[2] - gmin_l[2]) / ((PetscReal)mzlocal); in private_DMDALocatePointsIS_3D_Regular()
245 for (PetscInt d = 0; d < 3; d++) mi[d] = (PetscInt)((coor_p[d] - gmin[d]) / dx[d]); in private_DMDALocatePointsIS_3D_Regular()
|