Home
last modified time | relevance | path

Searched refs:det (Results 1 – 22 of 22) sorted by relevance

/petsc/src/dm/impls/moab/
H A Ddmmbfem.cxx13 PetscReal det = DMatrix_Determinant_2x2_Internal(inmat); in DMatrix_Invert_2x2_Internal() local
15 outmat[0] = inmat[3] / det; in DMatrix_Invert_2x2_Internal()
16 outmat[1] = -inmat[1] / det; in DMatrix_Invert_2x2_Internal()
17 outmat[2] = -inmat[2] / det; in DMatrix_Invert_2x2_Internal()
18 outmat[3] = inmat[0] / det; in DMatrix_Invert_2x2_Internal()
20 if (determinant) *determinant = det; in DMatrix_Invert_2x2_Internal()
31 PetscReal det = DMatrix_Determinant_3x3_Internal(inmat); in DMatrix_Invert_3x3_Internal() local
33 outmat[0] = (inmat[8] * inmat[4] - inmat[7] * inmat[5]) / det; in DMatrix_Invert_3x3_Internal()
34 outmat[1] = -(inmat[8] * inmat[1] - inmat[7] * inmat[2]) / det; in DMatrix_Invert_3x3_Internal()
35 outmat[2] = (inmat[5] * inmat[1] - inmat[4] * inmat[2]) / det; in DMatrix_Invert_3x3_Internal()
[all …]
/petsc/src/dm/impls/plex/
H A Dplexmetric.c1169 PetscScalar *met, *det; in DMPlexMetricEnforceSPD() local
1200 PetscCall(VecGetArray(determinant, &det)); in DMPlexMetricEnforceSPD()
1201 PetscCall(DMPlexMetricModify_Private(1, h_min, h_max, a_max, met, det)); in DMPlexMetricEnforceSPD()
1202 PetscCall(VecRestoreArray(determinant, &det)); in DMPlexMetricEnforceSPD()
1211 PetscCall(VecGetArray(determinant, &det)); in DMPlexMetricEnforceSPD()
1215 PetscCall(DMPlexPointLocalRef(dmDet, v, det, &vdet)); in DMPlexMetricEnforceSPD()
1218 PetscCall(VecRestoreArray(determinant, &det)); in DMPlexMetricEnforceSPD()
1267 PetscScalar *met, *det, integral, constants[1]; in DMPlexMetricNormalize() local
1296 PetscCall(VecGetArray(determinant, &det)); in DMPlexMetricNormalize()
1297 PetscCall(VecSet(tmp, det[0])); in DMPlexMetricNormalize()
[all …]
H A Dplexgeometry.c3499 PetscScalar det, idet; in DMPlexCoordinatesToReference_NewtonUpdate() local
3506 det = J[0] * J[3] - J[1] * J[2]; in DMPlexCoordinatesToReference_NewtonUpdate()
3507 idet = 1. / det; in DMPlexCoordinatesToReference_NewtonUpdate()
3517 det = invJ[0] * J[0] + invJ[1] * J[3] + invJ[2] * J[6]; in DMPlexCoordinatesToReference_NewtonUpdate()
3518 idet = 1. / det; in DMPlexCoordinatesToReference_NewtonUpdate()
/petsc/src/binding/petsc4py/demo/legacy/dmplex/
H A Danisotropic_adaptation.py91 det = dmHess.metricDeterminantCreate() variable
92 dmHess.metricNormalize(H, metric, det)
/petsc/src/dm/interface/
H A Ddmceed.c102 const CeedScalar det = J[0][0] * J[1][1] - J[0][1] * J[1][0]; in Geometry2D() local
104 qdata[i + Q * 0] = det * w[i]; in Geometry2D()
107 qdata[i + Q * 3] = J[1][1] / det; in Geometry2D()
108 qdata[i + Q * 4] = -J[1][0] / det; in Geometry2D()
109 qdata[i + Q * 5] = -J[0][1] / det; in Geometry2D()
110 qdata[i + Q * 6] = J[0][0] / det; in Geometry2D()
127 …const CeedScalar det = J[0][0] * (J[1][1] * J[2][2] - J[1][2] * J[2][1]) + J[0][1] * (J[1][2] * J[… in Geometry3D() local
129 qdata[i + Q * 0] = det * w[i]; /* det J * weight */ in Geometry3D()
/petsc/src/ts/tutorials/power_grid/stability_9bus/
H A Dex9busdmnetwork.c604 PetscScalar IGr, IGi, Zdq_inv[4], det; in FormIFunction() local
642 det = Rs * Rs + Xdp * Xqp; in FormIFunction()
644 Zdq_inv[0] = Rs / det; in FormIFunction()
645 Zdq_inv[1] = Xqp / det; in FormIFunction()
646 Zdq_inv[2] = -Xdp / det; in FormIFunction()
647 Zdq_inv[3] = Rs / det; in FormIFunction()
848 PetscScalar Vd, Vq, IGr, IGi, Zdq_inv[4], det; in AlgFunction() local
876 det = Rs * Rs + Xdp * Xqp; in AlgFunction()
878 Zdq_inv[0] = Rs / det; in AlgFunction()
879 Zdq_inv[1] = Xqp / det; in AlgFunction()
[all …]
H A Dex9busadj.c214 PetscScalar Zdq_inv[4], det; in ResidualFunction() local
264 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualFunction()
266 Zdq_inv[0] = Rs[i] / det; in ResidualFunction()
267 Zdq_inv[1] = Xqp[i] / det; in ResidualFunction()
268 Zdq_inv[2] = -Xdp[i] / det; in ResidualFunction()
269 Zdq_inv[3] = Rs[i] / det; in ResidualFunction()
442 PetscScalar Zdq_inv[4], det; in ResidualJacobian() local
517 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualJacobian()
519 Zdq_inv[0] = Rs[i] / det; in ResidualJacobian()
520 Zdq_inv[1] = Xqp[i] / det; in ResidualJacobian()
[all …]
H A Dex9busoptfd.c216 PetscScalar Zdq_inv[4], det; in ResidualFunction() local
265 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualFunction()
267 Zdq_inv[0] = Rs[i] / det; in ResidualFunction()
268 Zdq_inv[1] = Xqp[i] / det; in ResidualFunction()
269 Zdq_inv[2] = -Xdp[i] / det; in ResidualFunction()
270 Zdq_inv[3] = Rs[i] / det; in ResidualFunction()
443 PetscScalar Zdq_inv[4], det; in ResidualJacobian() local
518 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualJacobian()
520 Zdq_inv[0] = Rs[i] / det; in ResidualJacobian()
521 Zdq_inv[1] = Xqp[i] / det; in ResidualJacobian()
[all …]
H A Dex9bus.c442 PetscScalar Zdq_inv[4], det; in ResidualFunction() local
491 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualFunction()
493 Zdq_inv[0] = Rs[i] / det; in ResidualFunction()
494 Zdq_inv[1] = Xqp[i] / det; in ResidualFunction()
495 Zdq_inv[2] = -Xdp[i] / det; in ResidualFunction()
496 Zdq_inv[3] = Rs[i] / det; in ResidualFunction()
700 PetscScalar Zdq_inv[4], det; in ResidualJacobian() local
776 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualJacobian()
778 Zdq_inv[0] = Rs[i] / det; in ResidualJacobian()
779 Zdq_inv[1] = Xqp[i] / det; in ResidualJacobian()
[all …]
H A Dex9busopt.c361 PetscScalar Zdq_inv[4], det; in ResidualFunction() local
410 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualFunction()
412 Zdq_inv[0] = Rs[i] / det; in ResidualFunction()
413 Zdq_inv[1] = Xqp[i] / det; in ResidualFunction()
414 Zdq_inv[2] = -Xdp[i] / det; in ResidualFunction()
415 Zdq_inv[3] = Rs[i] / det; in ResidualFunction()
590 PetscScalar Zdq_inv[4], det; in ResidualJacobian() local
661 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualJacobian()
663 Zdq_inv[0] = Rs[i] / det; in ResidualJacobian()
664 Zdq_inv[1] = Xqp[i] / det; in ResidualJacobian()
[all …]
/petsc/src/dm/dt/tests/
H A Dex7.c199 PetscReal det, wvcheck, diff; in main() local
203 det = 1.; in main()
206 det = v[0]; in main()
209 det = v[0] * v[3] - v[1] * v[2]; in main()
212 …det = v[0] * (v[4] * v[8] - v[5] * v[7]) + v[1] * (v[5] * v[6] - v[3] * v[8]) + v[2] * (v[3] * v[7… in main()
215 …det = v[0] * (v[5] * (v[10] * v[15] - v[11] * v[14]) + v[6] * (v[11] * v[13] - v[9] * v[15]) + v[7… in main()
220 wvcheck = det * w[0]; in main()
/petsc/doc/tutorials/
H A Dguide_to_examples_by_physics.md65 J = \mathrm{det}(F) = \mathrm{det}\left(\nabla u\right)
/petsc/src/tao/complementarity/tutorials/
H A Dminsurf1.c429 PetscReal fnorm, det, hx, hy, xt = 0, yt = 0; in MSA_BoundaryConditions() local
483 det = njac11 * njac22 - njac21 * njac12; in MSA_BoundaryConditions()
484 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in MSA_BoundaryConditions()
485 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in MSA_BoundaryConditions()
/petsc/src/tao/unconstrained/tutorials/
H A Dminsurf1.c477 PetscReal fnorm, det, hx, hy, xt = 0, yt = 0; in MSA_BoundaryConditions() local
531 det = njac11 * njac22 - njac21 * njac12; in MSA_BoundaryConditions()
532 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in MSA_BoundaryConditions()
533 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in MSA_BoundaryConditions()
H A Dminsurf2.c694 PetscReal fnorm, det, hx, hy, xt = 0, yt = 0; in MSA_BoundaryConditions() local
753 det = njac11 * njac22 - njac21 * njac12; in MSA_BoundaryConditions()
754 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in MSA_BoundaryConditions()
755 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in MSA_BoundaryConditions()
/petsc/src/snes/tutorials/
H A Dex58.c466 PetscScalar det, hx, hy, xt = 0, yt = 0; in FormBoundaryConditions() local
533 det = njac11 * njac22 - njac21 * njac12; in FormBoundaryConditions()
534 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in FormBoundaryConditions()
535 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in FormBoundaryConditions()
H A Dex4.c462 PetscScalar det, hx, hy, xt = 0, yt = 0; in FormBoundaryConditions_Enneper() local
523 det = njac11 * njac22 - njac21 * njac12; in FormBoundaryConditions_Enneper()
524 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in FormBoundaryConditions_Enneper()
525 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in FormBoundaryConditions_Enneper()
H A Dex16.c228 …const PetscReal det = PetscRealPart(a * (e * i - f * h) - b * (i * d - f * g) + c * (d * h - e *… in InvertTensor() local
229 const PetscReal di = 1. / det; in InvertTensor()
250 if (dett) *dett = det; in InvertTensor()
/petsc/src/binding/petsc4py/test/
H A Dtest_dmplex.py211 det = self.plex.metricDeterminantCreate()
217 self.plex.metricEnforceSPD(metric, metric1, det[0])
225 metric, metric1, det[0], restrictSizes=False, restrictAnisotropy=False
/petsc/src/tao/bound/tutorials/
H A Dplate2f.F90519 PetscReal scl, fnorm, det, xt
600 det = njac11*njac22 - njac21*njac12
601 u1 = u1 - (njac22*nf1 - njac12*nf2)/det
602 u2 = u2 - (njac11*nf2 - njac21*nf1)/det
H A Dplate2.c610 PetscReal fnorm, det, hx, hy, xt = 0, yt = 0; in MSA_BoundaryConditions() local
675 det = njac11 * njac22 - njac21 * njac12; in MSA_BoundaryConditions()
676 u1 = u1 - (njac22 * nf1 - njac12 * nf2) / det; in MSA_BoundaryConditions()
677 u2 = u2 - (njac11 * nf2 - njac21 * nf1) / det; in MSA_BoundaryConditions()
/petsc/share/petsc/datafiles/meshes/
H A Dtestcase3D.cas6283 (mesh/bad/repair-grad/hinv/det-threshold 0.001)