| /petsc/src/snes/tutorials/ |
| H A D | ex75.h | 3 …Real *kn, PetscReal *km, PetscReal *B, PetscScalar vx[41][41], PetscScalar vz[41][41], PetscScalar… in SolKxData5()
1774 vz[0][0] = 0.00000000000000000000e+00; in SolKxData5()
1775 vz[0][1] = 1.11779031706136957277e-06; in SolKxData5()
1776 vz[0][2] = 0.00000000000000000000e+00; in SolKxData5()
1777 vz[0][3] = -1.11779031706136957277e-06; in SolKxData5()
1778 vz[0][4] = 0.00000000000000000000e+00; in SolKxData5()
1779 vz[0][5] = 1.11779031706136957277e-06; in SolKxData5()
1780 vz[0][6] = 0.00000000000000000000e+00; in SolKxData5()
1781 vz[0][7] = -1.11779031706136957277e-06; in SolKxData5()
1782 vz[0][8] = 0.00000000000000000000e+00; in SolKxData5()
[all …]
|
| H A D | ex75.c | 54 … z, PetscReal kn, PetscReal km, PetscReal B, PetscScalar *vx, PetscScalar *vz, PetscScalar *p, Pet… in SolKxSolution() argument
370 …*vz = UU * _C1 * t3 + UU * _C2 * t6 - _C1 * t6 * Rm + _C2 * t3 * Rm - VV * _C3 * t3 - VV * _C4 * t… in SolKxSolution()
482 *vz *= sin(km * z); /* Vz */ in SolKxSolution()
509 PetscScalar vx, vz, sxx, sxz, szz; in SolKxWrapperP() local
512 SolKxSolution(x[0], x[1], kn, km, B, &vx, &vz, &v[0], &sxx, &sxz, &szz); in SolKxWrapperP()
531 PetscScalar vx, vz, p, sxx, sxz, szz; in MapleTest() local
534 PetscCall(SolKxSolution(x[i], z[j], kn, km, B, &vx, &vz, &p, &sxx, &sxz, &szz)); in MapleTest()
535 …norm = PetscSqrt(PetscSqr(PetscAbsScalar(vx - vxMaple[i][j])) + PetscSqr(PetscAbsScalar(vz - vzMap… in MapleTest()
536 …e)z[j], (double)PetscAbsScalar(vx - vxMaple[i][j]), (double)PetscAbsScalar(vz - vzMaple[i][j]), (d… in MapleTest()
|
| /petsc/src/ts/utils/dmplexlandau/tutorials/ |
| H A D | ex1.c | 51 PetscReal vz = PetscRealPart(constants[0]); in f0_v2_par_shift() local
52 …if (dim == 2) *f0 = u[0] * (2. * PETSC_PI * x[0]) * (x[1] - vz) * (x[1] - vz); /* n r v^2_par|perp… in f0_v2_par_shift()
53 else *f0 = u[0] * (x[2] - vz) * (x[2] - vz); in f0_v2_par_shift()
68 PetscReal vz = PetscRealPart(constants[0]); in f0_v2_shift() local
69 if (dim == 2) f0[0] = u[0] * (2. * PETSC_PI * x[0]) * (x[0] * x[0] + (x[1] - vz) * (x[1] - vz)); in f0_v2_shift()
70 else f0[0] = u[0] * (x[0] * x[0] + x[1] * x[1] + (x[2] - vz) * (x[2] - vz)); in f0_v2_shift()
|
| H A D | ex2.c | 98 PetscReal vz = numConstants > 0 ? PetscRealPart(constants[0]) : 0; in f0_ve_shift() local
99 …= u[0] * 2. * PETSC_PI * x[0] * PetscSqrtReal(x[0] * x[0] + (x[1] - vz) * (x[1] - vz)); /* n r v */ in f0_ve_shift()
101 *f0 = u[0] * PetscSqrtReal(x[0] * x[0] + x[1] * x[1] + (x[2] - vz) * (x[2] - vz)); /* n v */ in f0_ve_shift()
153 …tscScalar user[2] = {0., ctx->charges[0]}, q[LANDAU_MAX_SPECIES], tt[LANDAU_MAX_SPECIES], vz; in testSpitzer() local
196 vz = ctx->n_0 * PetscRealPart(tt[0]) / n_e; /* non-dimensional */ in testSpitzer()
197 } else vz = 0; in testSpitzer()
200 PetscCall(PetscDSSetConstants(prob, 1, &vz)); in testSpitzer()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexcreate.c | 718 PetscInt v, vx, vy, vz; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() local
781 for (vz = 0; vz < faces[2]; vz++) { in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
783 voffset = numFaces + vz * vertices[0] * vertices[1] + vx; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
799 for (vz = 0; vz < faces[2]; vz++) { in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
801 voffset = numFaces + vz * vertices[0] * vertices[1] + vertices[0] * (vertices[1] - 1) + vx; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
817 for (vz = 0; vz < faces[2]; vz++) { in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
819 voffset = numFaces + vz * vertices[0] * vertices[1] + vy * vertices[0]; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
835 for (vz = 0; vz < faces[2]; vz++) { in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
837 voffset = numFaces + vz * vertices[0] * vertices[1] + vy * vertices[0] + faces[0]; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
872 for (vz = 0; vz <= faces[2]; ++vz) { in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/pipecg2/ |
| H A D | pipecg2.c | 82 static PetscErrorCode VecMergedOpsShort_Private(Vec vx, Vec vr, Vec vz, Vec vw, Vec vp, Vec vq, Vec… in VecMergedOpsShort_Private() argument
92 PetscCall(VecGetArray(vz, (PetscScalar **)&pz)); in VecMergedOpsShort_Private()
282 PetscCall(VecRestoreArray(vz, (PetscScalar **)&pz)); in VecMergedOpsShort_Private()
304 static PetscErrorCode VecMergedOps_Private(Vec vx, Vec vr, Vec vz, Vec vw, Vec vp, Vec vq, Vec vc, … in VecMergedOps_Private() argument
314 PetscCall(VecGetArray(vz, (PetscScalar **)&pz)); in VecMergedOps_Private()
511 PetscCall(VecRestoreArray(vz, (PetscScalar **)&pz)); in VecMergedOps_Private()
|
| /petsc/src/tao/pde_constrained/tutorials/ |
| H A D | parabolic.c | 665 PetscScalar v, vx, vy, vz; in ParabolicInitialize() local
726 vz = h * (k + 0.5); in ParabolicInitialize()
729 PetscCall(VecSetValues(ZZ, 1, &linear_index, &vz, INSERT_VALUES)); in ParabolicInitialize()
730 …if ((vx < 0.6) && (vx > 0.4) && (vy < 0.6) && (vy > 0.4) && (vy < 0.6) && (vz < 0.6) && (vz > 0.4)… in ParabolicInitialize()
|
| H A D | elliptic.c | 552 PetscScalar v, vx, vy, vz; in EllipticInitialize() local
699 vz = h * (k + 0.5); in EllipticInitialize()
702 PetscCall(VecSetValues(ZZ, 1, &linear_index, &vz, INSERT_VALUES)); in EllipticInitialize()
719 … PetscSinScalar(2 * PETSC_PI * (vy + 0.25 * js)) * PetscSinScalar(2 * PETSC_PI * (vz + 0.25 * ks)); in EllipticInitialize()
|