| /petsc/src/snes/tutorials/ |
| H A D | ex30.c | 98 PetscReal dx, dz; member 314 z = (j - grid->jlid - 0.5) * grid->dz; in HorizVelocity() 331 z = (j - grid->jlid) * grid->dz; in VertVelocity() 347 z = (j - grid->jlid - 0.5) * grid->dz; in Pressure() 413 eps22 = (wN - wS) / grid->dz; 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 475 TS = param->potentialT * TInterp(x, i, j - 1) * PetscExpScalar((j - 1.0) * dz * z_scale); in XMomentumResidual() 477 else TN = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale); in XMomentumResidual() 478 TW = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale); in XMomentumResidual() [all …]
|
| H A D | ex48.c | 193 static void HexGrad(const PetscReal dphi[][3], const PetscReal zn[], PetscReal dz[]) in HexGrad() argument 196 dz[0] = dz[1] = dz[2] = 0; in HexGrad() 198 dz[0] += dphi[i][0] * zn[i]; in HexGrad() 199 dz[1] += dphi[i][1] * zn[i]; in HexGrad() 200 dz[2] += dphi[i][2] * zn[i]; in HexGrad() 204 static void HexComputeGeometry(PetscInt q, PetscReal hx, PetscReal hy, const PetscReal dz[PETSC_RES… in HexComputeGeometry() 209 {dz[0], dz[1], dz[2]} in HexComputeGeometry() 791 PetscReal dz[3], phi[8], dphi[8][3], jw, eta, deta; in THIFunctionLocal() local 793 HexGrad(HexQDeriv[q], zn, dz); in THIFunctionLocal() 794 HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw); in THIFunctionLocal() [all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex10.c | 211 PetscReal K[60][60], x, y, z, dx, dy, dz; in Elastic20Stiff() local 221 dz = 2.0; in Elastic20Stiff() 248 xyz[8][2] = z + dz; in Elastic20Stiff() 251 xyz[9][2] = z + dz; in Elastic20Stiff() 254 xyz[10][2] = z + dz; in Elastic20Stiff() 257 xyz[11][2] = z + dz; in Elastic20Stiff() 260 xyz[12][2] = z + 2. * dz; in Elastic20Stiff() 263 xyz[13][2] = z + 2. * dz; in Elastic20Stiff() 266 xyz[14][2] = z + 2. * dz; in Elastic20Stiff() 269 xyz[15][2] = z + 2. * dz; in Elastic20Stiff() [all …]
|
| /petsc/src/ksp/ksp/impls/cg/pipecgrr/ |
| H A D | pipecgrr.c | 25 …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 108 znp = dz; in KSPSolve_PIPECGRR() 129 PetscCall(VecNormBegin(Z, NORM_2, &dz)); in KSPSolve_PIPECGRR() 152 PetscCall(VecNormEnd(Z, NORM_2, &dz)); in KSPSolve_PIPECGRR() 194 …* unp + 2.0 * Anorm * PetscAbsScalar(alphap) * dq + wnp + 2.0 * PetscAbsScalar(alphap) * dz) * eps; in KSPSolve_PIPECGRR()
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex6.c | 510 … const PetscScalar dz = arr_coord_z[ez][slot_coord_next] - arr_coord_z[ez][slot_coord_prev]; in UpdateVelocity_3d() local 513 …y + (arr_stress[ez][ey][ex][slot_tyz_frontdown] - arr_stress[ez][ey][ex][slot_tyz_backdown]) / dz); in UpdateVelocity_3d() 520 … const PetscScalar dz = arr_coord_z[ez][slot_coord_next] - arr_coord_z[ez][slot_coord_prev]; in UpdateVelocity_3d() local 523 …y + (arr_stress[ez][ey][ex][slot_txz_frontleft] - arr_stress[ez][ey][ex][slot_txz_backleft]) / dz); in UpdateVelocity_3d() 530 …const PetscScalar dz = arr_coord_z[ez][slot_coord_element] - arr_coord_z[ez - 1][slot_coord_elemen… in UpdateVelocity_3d() local 533 …_backdown]) / dy + (arr_stress[ez][ey][ex][slot_tzz] - arr_stress[ez - 1][ey][ex][slot_tzz]) / dz); in UpdateVelocity_3d() 705 … const PetscScalar dz = arr_coord_z[ez][slot_coord_next] - arr_coord_z[ez][slot_coord_prev]; in UpdateStress_3d() local 710 … ctx->dt * (arr_velocity[ez][ey][ex][slot_vz_front] - arr_velocity[ez][ey][ex][slot_vz_back]) / dz; in UpdateStress_3d() 712 …ctx->dt * (arr_velocity[ez][ey][ex][slot_vz_front] - arr_velocity[ez][ey][ex][slot_vz_back]) / dz + in UpdateStress_3d() 715 …ity[ez][ey][ex][slot_vz_front] - arr_velocity[ez][ey][ex][slot_vz_back]) / dz + L * ctx->dt * (arr… in UpdateStress_3d() [all …]
|
| /petsc/src/vec/is/sf/impls/basic/ |
| H A D | sfpack.c | 69 for (k = 0; k < opt->dz[r]; k++) \ 112 for (k = 0; k < opt->dz[r]; k++) \ 155 for (k = 0; k < opt->dz[r]; k++) \ 207 for (k = 0; k < srcOpt->dz[0]; k++) \ 1180 PetscInt r, p, start, i, j, k, dx, dy, dz, dydz, m, X, Y; in PetscSFCreatePackOpt() local 1192 opt->dz = opt->array + 4 * n + 2; in PetscSFCreatePackOpt() 1228 dz = m / (dx * dy); in PetscSFCreatePackOpt() 1229 Y = dz > 1 ? (idx[p] - start) / X : dy; in PetscSFCreatePackOpt() 1235 …for (k = 1; k < dz; k++) { /* Go through Z dimension to see if remaining indices follow the patter… in PetscSFCreatePackOpt() 1248 opt->dz[r] = dz; in PetscSFCreatePackOpt() [all …]
|
| H A D | sfpack.h | 48 PetscInt *dx, *dy, *dz; /* [n] Lengths of the submatrix in X, Y, Z dimension. */ member
|
| /petsc/src/snes/tutorials/output/ |
| H A D | ex30_1.out | 8 [ni,nj] = 18, 18 [dx,dz] = 20., 20. km
|
| /petsc/src/ts/tutorials/ |
| H A D | ex14.c | 181 static void HexGrad(const PetscReal dphi[][3], const PetscReal zn[], PetscReal dz[]) in HexGrad() argument 184 dz[0] = dz[1] = dz[2] = 0; in HexGrad() 186 dz[0] += dphi[i][0] * zn[i]; in HexGrad() 187 dz[1] += dphi[i][1] * zn[i]; in HexGrad() 188 dz[2] += dphi[i][2] * zn[i]; in HexGrad() 192 static void HexComputeGeometry(PetscInt q, PetscReal hx, PetscReal hy, const PetscReal dz[restrict]… in HexComputeGeometry() 198 {dz[0], dz[1], dz[2]} in HexComputeGeometry() 789 PetscReal dz[3], phi[8], dphi[8][3], jw, eta, deta; in THIFunctionLocal_3D() local 795 HexGrad(HexQDeriv[q], zn, dz); in THIFunctionLocal_3D() 796 HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw); in THIFunctionLocal_3D() [all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexegads.c | 99 PetscScalar dx, dy, dz, lambda, tolr, obj_old, obj_tmp, target; in DMPlex_EGADS_EDGE_XYZtoUV_Internal() local 114 dz = coords[v * dE + 2] - eval[2]; in DMPlex_EGADS_EDGE_XYZtoUV_Internal() 115 obj_old = dx * dx + dy * dy + dz * dz; in DMPlex_EGADS_EDGE_XYZtoUV_Internal() 127 b = eval[3] * dx + eval[4] * dy + eval[5] * dz; in DMPlex_EGADS_EDGE_XYZtoUV_Internal() 182 PetscScalar dx, dy, dz, lambda, tolr, obj_old, obj_tmp, target; in DMPlex_Geom_EDGE_XYZtoUV_Internal() local 202 dz = coords[v * dE + 2] - eval[2]; in DMPlex_Geom_EDGE_XYZtoUV_Internal() 203 obj_old = dx * dx + dy * dy + dz * dz; in DMPlex_Geom_EDGE_XYZtoUV_Internal() 215 b = eval[3] * dx + eval[4] * dy + eval[5] * dz; in DMPlex_Geom_EDGE_XYZtoUV_Internal() 280 PetscScalar dx, dy, dz, lambda, tolr, denom, obj_old, obj_tmp, target; in DMPlex_EGADS_FACE_XYZtoUV_Internal() local 295 dz = coords[v * dE + 2] - eval[2]; in DMPlex_EGADS_FACE_XYZtoUV_Internal() [all …]
|
| H A D | plexgeometry.c | 2766 const PetscReal dz = e0[0] * e1[1] - e0[1] * e1[0]; in DMPlexComputeGeometryFVM_2D_Internal() local 2767 const PetscReal a = PetscSqrtReal(dx * dx + dy * dy + dz * dz); in DMPlexComputeGeometryFVM_2D_Internal() 2771 n[2] += dz; in DMPlexComputeGeometryFVM_2D_Internal()
|
| /petsc/include/ |
| H A D | petscblaslapack_mangle.h | 93 #define PETSC_BLASLAPACK_RPREFIX_ dz
|
| /petsc/src/dm/impls/moab/ |
| H A D | dmmbutil.cxx | 19 PetscReal dx, dy, dz; member 234 arrays[2][ix] = (z + kk) * genCtx.dz + genCtx.xyzbounds[4]; in DMMoab_GenerateVertices_Private() 431 …genCtx.dz = (genCtx.xyzbounds[5] - genCtx.xyzbounds[4]) / (nelems * genCtx.q); /* distance be… in DMMBUtil_InitializeOptions() 474 PetscCall(PetscInfo(NULL, "Local delxyz:= %g, %g, %g\n", genCtx.dx, genCtx.dy, genCtx.dz)); in DMMBUtil_InitializeOptions()
|
| /petsc/src/ts/tests/ |
| H A D | ex5.c | 442 PetscScalar dz = 1; /* depth of layer between soil surface and deep soil layer */ in calc_gflux() local 447 …*Gflux = (k * (deep_grnd_temp - sfctemp) / dz); … in calc_gflux()
|
| /petsc/src/tao/constrained/impls/ipm/ |
| H A D | pdipm.c | 709 const PetscScalar *dXarr, *dz, *dlambdai; in SNESLineSearch_PDIPM() local 721 dz = dXarr + pdipm->off_z; in SNESLineSearch_PDIPM() 723 if (z[i] - dz[i] < 0.0) alpha_p = PetscMin(alpha_p, 0.9999 * z[i] / dz[i]); in SNESLineSearch_PDIPM()
|
| /petsc/lib/petsc/bin/maint/ |
| H A D | toclapack.sh | 351 i[dz]amax) PR="DOUBLE";; 355 [dz]*) PR="DOUBLE";; 363 ila[dz]l[rc]) PR="DOUBLE";;
|
| /petsc/doc/manual/ |
| H A D | ts.md | 312 where $z$ is a new constraint variable, and the Jacobian $\frac{dh}{dz}$ is non-singular everywhere… 335 If the Jacobian $\frac{dh}{du} \frac{df}{dz}$ is non-singular, then we have precisely a semi-explic…
|