| /petsc/src/ksp/pc/impls/tfs/ |
| H A D | bitmask.c | 86 PetscInt PCTFS_div_ceil(PetscInt numer, PetscInt denom) in PCTFS_div_ceil() argument
88 …denom <= 0)) SETERRABORT(PETSC_COMM_SELF, PETSC_ERR_PLIB, "PCTFS_div_ceil() :: numer=%" PetscInt_F… in PCTFS_div_ceil()
89 return PetscCeilInt(numer, denom); in PCTFS_div_ceil()
|
| /petsc/src/ksp/ksp/utils/lmvm/rescale/ |
| H A D | symbrdnrescale.c | 272 PetscReal denom = PetscMax(-b - sqrt_disc, ldb->tol); in SymBroydenRescaleUpdateDiagonal() local
274 sigma = (2 * c) / denom; in SymBroydenRescaleUpdateDiagonal()
276 PetscReal denom = PetscMax(2 * a, ldb->tol); in SymBroydenRescaleUpdateDiagonal() local
278 sigma = (-b - sqrt_disc) / denom; in SymBroydenRescaleUpdateDiagonal()
282 PetscReal denom = PetscMax(2 * a, ldb->tol); in SymBroydenRescaleUpdateDiagonal() local
284 sigma = (-b + sqrt_disc) / denom; in SymBroydenRescaleUpdateDiagonal()
286 PetscReal denom = PetscMax(-b + sqrt_disc, ldb->tol); in SymBroydenRescaleUpdateDiagonal() local
288 sigma = (2 * c) / denom; in SymBroydenRescaleUpdateDiagonal()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfvm.c | 45 PetscReal denom = DMPlex_DotD_Internal(dim, &cgrad[d * dim], v); in DMPlexApplyLimiter_Internal() local
46 PetscReal fact = denom == 0 ? 1.0e+30 : 1 / denom; in DMPlexApplyLimiter_Internal()
|
| H A D | plexgeometry.c | 123 const PetscReal denom = (-s2_x * s1_y + s1_x * s2_y);
128 if (denom != 0.0) {
129 const PetscReal s = (-s1_y * (p0_x - p2_x) + s1_x * (p0_y - p2_y)) / denom;
130 const PetscReal t = (s2_x * (p0_y - p2_y) - s2_y * (p0_x - p2_x)) / denom;
182 const PetscReal denom = -(s0_x * s3_x + s0_y * s3_y + s0_z * s3_z); /* -s0 . (s1 x s2) */
187 if (denom != 0.0) {
188 … const PetscReal t = (s3_x * (p0_x - q0_x) + s3_y * (p0_y - q0_y) + s3_z * (p0_z - q0_z)) / denom;
189 … const PetscReal u = (s4_x * (p0_x - q0_x) + s4_y * (p0_y - q0_y) + s4_z * (p0_z - q0_z)) / denom;
190 … const PetscReal v = (s5_x * (p0_x - q0_x) + s5_y * (p0_y - q0_y) + s5_z * (p0_z - q0_z)) / denom;
|
| H A D | plexegads.c | 280 PetscScalar dx, dy, dz, lambda, tolr, denom, obj_old, obj_tmp, target; in DMPlex_EGADS_FACE_XYZtoUV_Internal() local
312 denom = A[0] * A[3] - A[2] * A[1]; in DMPlex_EGADS_FACE_XYZtoUV_Internal()
313 if (denom == 0.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "denom = 0.0 \n")); in DMPlex_EGADS_FACE_XYZtoUV_Internal()
314 delta[0] = (b[0] * A[3] - b[1] * A[1]) / denom; in DMPlex_EGADS_FACE_XYZtoUV_Internal()
315 delta[1] = (A[0] * b[1] - A[2] * b[0]) / denom; in DMPlex_EGADS_FACE_XYZtoUV_Internal()
377 PetscScalar dx, dy, dz, lambda, tolr, denom, obj_old, obj_tmp, target; in DMPlex_Geom_FACE_XYZtoUV_Internal() local
414 denom = A[0] * A[3] - A[2] * A[1]; in DMPlex_Geom_FACE_XYZtoUV_Internal()
415 if (denom == 0.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "denom = 0.0 \n")); in DMPlex_Geom_FACE_XYZtoUV_Internal()
416 delta[0] = (b[0] * A[3] - b[1] * A[1]) / denom; in DMPlex_Geom_FACE_XYZtoUV_Internal()
417 delta[1] = (A[0] * b[1] - A[2] * b[0]) / denom; in DMPlex_Geom_FACE_XYZtoUV_Internal()
|
| /petsc/src/ts/tests/ |
| H A D | ex35.c | 39 PetscReal denom = (m + 1.) * (2. * m + 1.); in erfinv() local
40 ck[k] += (ck[m] * ck[k - 1 - m]) / denom; in erfinv()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex10.c | 148 PetscScalar sigma_R, denom; in RDDiffusionCoefficient() local
152 denom = 3. * rd->rho * sigma_R + (int)limit * PetscAbsScalar(nx->E) / n->E; in RDDiffusionCoefficient()
157 *D_R = rd->c / denom; in RDDiffusionCoefficient()
159 dD_R->E = -rd->c / PetscSqr(denom) * ddenom.E; in RDDiffusionCoefficient()
160 dD_R->T = -rd->c / PetscSqr(denom) * ddenom.T; in RDDiffusionCoefficient()
163 dxD_R->E = -rd->c / PetscSqr(denom) * dxdenom.E; in RDDiffusionCoefficient()
164 dxD_R->T = -rd->c / PetscSqr(denom) * dxdenom.T; in RDDiffusionCoefficient()
|
| /petsc/src/ksp/ksp/utils/dm/ |
| H A D | dmproject.c | 524 Vec v_field_l, denom_l, coor_l, denom; in DMSwarmProjectField_ApproxQ1_DA_2D() local
539 PetscCall(DMGetGlobalVector(dm, &denom)); in DMSwarmProjectField_ApproxQ1_DA_2D()
542 PetscCall(VecZeroEntries(denom)); in DMSwarmProjectField_ApproxQ1_DA_2D()
602 PetscCall(DMLocalToGlobalBegin(dm, denom_l, ADD_VALUES, denom)); in DMSwarmProjectField_ApproxQ1_DA_2D()
603 PetscCall(DMLocalToGlobalEnd(dm, denom_l, ADD_VALUES, denom)); in DMSwarmProjectField_ApproxQ1_DA_2D()
605 PetscCall(VecPointwiseDivide(v_field, v_field, denom)); in DMSwarmProjectField_ApproxQ1_DA_2D()
609 PetscCall(DMRestoreGlobalVector(dm, &denom)); in DMSwarmProjectField_ApproxQ1_DA_2D()
|
| /petsc/src/dm/dt/interface/ |
| H A D | dtprob.c | 207 PetscReal denom = (m + 1.) * (2. * m + 1.); in PetscPDFSampleGaussian1D() local
209 ck[k] += (ck[m] * ck[k - 1 - m]) / denom; in PetscPDFSampleGaussian1D()
|