| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex15.c | 245 PetscReal hx, hy, dv; in CreateSystem2d() local
256 dv = hx * hy; in CreateSystem2d()
315 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem2d()
320 valA[1] = dv * 1.0 / (hy * hy); in CreateSystem2d()
325 valA[2] = dv * 1.0 / (hy * hy); in CreateSystem2d()
331 valA[3] = dv * 1.0 / (hx * hx); in CreateSystem2d()
336 valA[4] = dv * 1.0 / hy; in CreateSystem2d()
341 valA[5] = -dv * 1.0 / hy; in CreateSystem2d()
349 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem2d()
354 valA[1] = dv * 1.0 / (hy * hy); in CreateSystem2d()
[all …]
|
| H A D | ex18.c | 256 PetscReal hx, hy, dv, bogusScale; in CreateSystem() local
283 dv = hx * hy; in CreateSystem()
355 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem()
360 valA[1] = dv * 1.0 / (hy * hy); in CreateSystem()
365 valA[2] = dv * 1.0 / (hy * hy); in CreateSystem()
371 valA[3] = dv * 1.0 / (hx * hx); in CreateSystem()
378 valA[4] = dv * 1.0 / hy; in CreateSystem()
383 valA[5] = -dv * 1.0 / hy; in CreateSystem()
392 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem()
397 valA[1] = dv * 1.0 / (hy * hy); in CreateSystem()
[all …]
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex4.c | 579 PetscReal hx, hy, dv; in CreateSystem2d() local
610 dv = hx * hy; in CreateSystem2d()
674 val_rhs = -ctx->gy * dv * 0.5 * (rho[0] + rho[1]); in CreateSystem2d()
702 val_A[count] = -2.0 * dv * (eta_down + eta_up) / (hy * hy); in CreateSystem2d()
703 if (!left_boundary) val_A[count] += -1.0 * dv * eta_left / (hx * hx); in CreateSystem2d()
704 if (!right_boundary) val_A[count] += -1.0 * dv * eta_right / (hx * hx); in CreateSystem2d()
711 val_A[count] = 2.0 * dv * eta_down / (hy * hy); in CreateSystem2d()
717 val_A[count] = 2.0 * dv * eta_up / (hy * hy); in CreateSystem2d()
724 val_A[count] = dv * eta_left / (hx * hx); in CreateSystem2d()
732 val_A[count] = dv * eta_right / (hx * hx); in CreateSystem2d()
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex48.c | 710 …PetscScalar *PETSC_RESTRICT v, PetscScalar du[PETSC_RESTRICT], PetscScalar dv[PETSC_RESTRICT], Pet… in PointwiseNonlinearity()
716 dv[0] = dv[1] = dv[2] = 0; in PointwiseNonlinearity()
724 dv[ll] += dphi[l][ll] * n[l].v; in PointwiseNonlinearity()
727 …etscSqr(du[0]) + PetscSqr(dv[1]) + du[0] * dv[1] + 0.25 * PetscSqr(du[1] + dv[0]) + 0.25 * PetscSq… in PointwiseNonlinearity()
731 … dphi[4][2], PetscScalar *u, PetscScalar *v, PetscScalar du[], PetscScalar dv[], PetscReal *eta, P… in PointwiseNonlinearity2D() argument
737 dv[0] = dv[1] = 0; in PointwiseNonlinearity2D()
745 dv[ll] += dphi[l][ll] * n[l].v; 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/tutorials/ |
| H A D | ex14.c | 729 …strict u, PetscScalar *restrict v, PetscScalar du[restrict 3], PetscScalar dv[restrict 3], PetscRe… in PointwiseNonlinearity()
735 dv[0] = dv[1] = dv[2] = 0; in PointwiseNonlinearity()
743 dv[ll] += dphi[l][ll] * n[l].v; 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 …fn[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 …v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp[2] * j… 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 …]
|
| /petsc/include/ |
| H A D | petscfe.h | 179 CeedScalar *v = out[0], *dv = out[1]; \
201 dv[i + Q * 0] = qdata[i + Q * 0] * (invJ[0][0] * f1[0] + invJ[0][1] * f1[1]); \
202 dv[i + Q * 1] = qdata[i + Q * 0] * (invJ[1][0] * f1[0] + invJ[1][1] * f1[1]); \
|
| /petsc/doc/manual/ |
| H A D | tao.md | 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).
2273 \displaystyle \min_{dv} & \tilde{f}_k(u_k-A_k^{-1}(B_k dv + \alpha_k g_k), v_k+dv), \\
2281 \displaystyle \min_{dv} & \tilde{f}_k(u_{k+\frac{1}{2}} - A_k^{-1} B_k dv, v_{k+\frac{1}{2}}+dv). \\
2290 \displaystyle \min_{dv} & \frac{1}{2} dv^T \tilde{H}_k dv + \tilde{g}_{k+\frac{1}{2}}^T dv,
2322 dv = -H_k^{-1} \tilde{g}_{k+\frac{1}{2}}
2328 du = -A_k^{-1} B_k dv.
2337 …style \min_{\beta \geq 0} & \tilde{f_k}(u_{k+\frac{1}{2}} + \beta du, v_{k+\frac{1}{2}} + \beta dv)
2347 v_{k+1} & = & v_{k+\frac{1}{2}} + \beta_k dv.
|
| /petsc/include/petsc/private/ |
| H A D | matdensecupmimpl.h | 460 auto dv = VecSeq_CUPM::DeviceArrayWrite(dctx, v); in GetDiagonal() local
471 auto dvp = thrust::device_pointer_cast(dv.data()); in GetDiagonal()
|
| /petsc/src/mat/impls/aij/seq/kokkos/ |
| H A D | aijkok.kokkos.cxx | 1064 ConstPetscScalarKokkosView dv; in MatDiagonalSet_SeqAIJKokkos() local
1069 PetscCall(VecGetKokkosView(D, &dv)); in MatDiagonalSet_SeqAIJKokkos()
1079 if (is == INSERT_VALUES) Aa(Adiag(i)) = dv(i); in MatDiagonalSet_SeqAIJKokkos()
1080 else Aa(Adiag(i)) += dv(i); in MatDiagonalSet_SeqAIJKokkos()
1082 PetscCall(VecRestoreKokkosView(D, &dv)); in MatDiagonalSet_SeqAIJKokkos()
|
| /petsc/src/mat/impls/aij/seq/seqhipsparse/ |
| H A D | aijhipsparse.hip.cxx | 4366 thrust::device_ptr<PetscScalar> dv; in MatSeqAIJCopySubArray_SeqAIJHIPSPARSE() local
4367 if (dmem) dv = thrust::device_pointer_cast(v); in MatSeqAIJCopySubArray_SeqAIJHIPSPARSE()
4370 dv = w->data(); in MatSeqAIJCopySubArray_SeqAIJHIPSPARSE()
4374 …t::make_zip_iterator(thrust::make_tuple(thrust::make_permutation_iterator(dav, widx.begin()), dv)); in MatSeqAIJCopySubArray_SeqAIJHIPSPARSE()
4375 …:make_zip_iterator(thrust::make_tuple(thrust::make_permutation_iterator(dav, widx.end()), dv + n)); in MatSeqAIJCopySubArray_SeqAIJHIPSPARSE()
|
| /petsc/src/mat/impls/aij/seq/seqcusparse/ |
| H A D | aijcusparse.cu | 5047 thrust::device_ptr<PetscScalar> dv; in MatSeqAIJCopySubArray_SeqAIJCUSPARSE() local
5049 dv = thrust::device_pointer_cast(v); in MatSeqAIJCopySubArray_SeqAIJCUSPARSE()
5052 dv = w->data(); in MatSeqAIJCopySubArray_SeqAIJCUSPARSE()
5056 …t::make_zip_iterator(thrust::make_tuple(thrust::make_permutation_iterator(dav, widx.begin()), dv)); in MatSeqAIJCopySubArray_SeqAIJCUSPARSE()
5057 …:make_zip_iterator(thrust::make_tuple(thrust::make_permutation_iterator(dav, widx.end()), dv + n)); in MatSeqAIJCopySubArray_SeqAIJCUSPARSE()
|