Lines Matching refs:fegeom
1464 PetscFEGeom fegeom; in DMPlexComputeL2DiffLocal() local
1475 fegeom.dimEmbed = coordDim; in DMPlexComputeL2DiffLocal()
1504 …coordDim * (Nq + 1), &coords, Nq, &fegeom.detJ, coordDim * coordDim * Nq, &fegeom.J, coordDim * co… in DMPlexComputeL2DiffLocal()
1512 … PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)); in DMPlexComputeL2DiffLocal()
1539 qgeom.dimEmbed = fegeom.dimEmbed; in DMPlexComputeL2DiffLocal()
1540 qgeom.J = &fegeom.J[q * coordDim * coordDim]; in DMPlexComputeL2DiffLocal()
1541 qgeom.invJ = &fegeom.invJ[q * coordDim * coordDim]; in DMPlexComputeL2DiffLocal()
1542 qgeom.detJ = &fegeom.detJ[q]; in DMPlexComputeL2DiffLocal()
1543 …fegeom.detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for eleme… in DMPlexComputeL2DiffLocal()
1555 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2DiffLocal()
1565 …(double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]), (double)Pe… in DMPlexComputeL2DiffLocal()
1566 elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]; in DMPlexComputeL2DiffLocal()
1576 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2DiffLocal()
1591 PetscFEGeom fegeom; in DMComputeL2GradientDiff_Plex() local
1600 fegeom.dimEmbed = coordDim; in DMComputeL2GradientDiff_Plex()
1621 …oordDim * coordDim * Nq, &fegeom.J, coordDim * coordDim * Nq, &fegeom.invJ, numComponents * coordD… in DMComputeL2GradientDiff_Plex()
1628 … PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)); in DMComputeL2GradientDiff_Plex()
1648 qgeom.dimEmbed = fegeom.dimEmbed; in DMComputeL2GradientDiff_Plex()
1649 qgeom.J = &fegeom.J[q * coordDim * coordDim]; in DMComputeL2GradientDiff_Plex()
1650 qgeom.invJ = &fegeom.invJ[q * coordDim * coordDim]; in DMComputeL2GradientDiff_Plex()
1651 qgeom.detJ = &fegeom.detJ[q]; in DMComputeL2GradientDiff_Plex()
1652 …fegeom.detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for eleme… in DMComputeL2GradientDiff_Plex()
1664 PetscCall(PetscFree6(funcVal, coords, fegeom.J, fegeom.invJ, interpolant, fegeom.detJ)); in DMComputeL2GradientDiff_Plex()
1679 …ield, fc, (double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]))); in DMComputeL2GradientDiff_Plex()
1680 elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]; in DMComputeL2GradientDiff_Plex()
1690 PetscCall(PetscFree6(funcVal, coords, fegeom.J, fegeom.invJ, interpolant, fegeom.detJ)); in DMComputeL2GradientDiff_Plex()
1733 PetscFEGeom fegeom; in DMComputeL2FieldDiff_Plex() local
1747 …Nc, &interpolant, dE * (Nq + 1), &coords, Nq, &fegeom.detJ, dE * dE * Nq, &fegeom.J, dE * dE * Nq,… in DMComputeL2FieldDiff_Plex()
1765 …cCall(DMPlexComputeCellGeometryFEM(dm, cone[0], quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)); in DMComputeL2FieldDiff_Plex()
1767 …etscCall(DMPlexComputeCellGeometryFEM(dm, cell, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)); in DMComputeL2FieldDiff_Plex()
1802 qgeom.dimEmbed = fegeom.dimEmbed; in DMComputeL2FieldDiff_Plex()
1803 qgeom.J = &fegeom.J[q * dE * dE]; in DMComputeL2FieldDiff_Plex()
1804 qgeom.invJ = &fegeom.invJ[q * dE * dE]; in DMComputeL2FieldDiff_Plex()
1805 qgeom.detJ = &fegeom.detJ[q]; in DMComputeL2FieldDiff_Plex()
1806 …fegeom.detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for cell … in DMComputeL2FieldDiff_Plex()
1818 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMComputeL2FieldDiff_Plex()
1829 … (double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]))); in DMComputeL2FieldDiff_Plex()
1830 … elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]; in DMComputeL2FieldDiff_Plex()
1845 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMComputeL2FieldDiff_Plex()
1878 PetscFEGeom fegeom; in DMPlexComputeL2DiffVec() local
1916 …lant, coordDim * Nq, &coords, Nq, &fegeom.detJ, coordDim * coordDim * Nq, &fegeom.J, coordDim * co… in DMPlexComputeL2DiffVec()
1923 … PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)); in DMPlexComputeL2DiffVec()
1945 qgeom.dimEmbed = fegeom.dimEmbed; in DMPlexComputeL2DiffVec()
1946 qgeom.J = &fegeom.J[q * coordDim * coordDim]; in DMPlexComputeL2DiffVec()
1947 qgeom.invJ = &fegeom.invJ[q * coordDim * coordDim]; in DMPlexComputeL2DiffVec()
1948 qgeom.detJ = &fegeom.detJ[q]; in DMPlexComputeL2DiffVec()
1949 …fegeom.detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for eleme… in DMPlexComputeL2DiffVec()
1954 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2DiffVec()
1963 … elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]; in DMPlexComputeL2DiffVec()
1973 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2DiffVec()
2006 PetscFEGeom fegeom; in DMPlexComputeL2FluxDiffVecLocal() local
2025 fegeom.dimEmbed = cdim; in DMPlexComputeL2FluxDiffVecLocal()
2037 …ant, cdim * (Nq + 1), &coords, cdim * cdim * Nq, &fegeom.J, cdim * cdim * Nq, &fegeom.invJ, Nq, &f… in DMPlexComputeL2FluxDiffVecLocal()
2044 … PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)); in DMPlexComputeL2FluxDiffVecLocal()
2051 qgeom.dimEmbed = fegeom.dimEmbed; in DMPlexComputeL2FluxDiffVecLocal()
2052 qgeom.J = &fegeom.J[q * cdim * cdim]; in DMPlexComputeL2FluxDiffVecLocal()
2053 qgeom.invJ = &fegeom.invJ[q * cdim * cdim]; in DMPlexComputeL2FluxDiffVecLocal()
2054 qgeom.detJ = &fegeom.detJ[q]; in DMPlexComputeL2FluxDiffVecLocal()
2056 …fegeom.detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for eleme… in DMPlexComputeL2FluxDiffVecLocal()
2064 … fluxElemDiff += PetscSqr(PetscRealPart(interpolant[fc] - minterpolant[fc])) * wt * fegeom.detJ[q]; in DMPlexComputeL2FluxDiffVecLocal()
2072 PetscCall(PetscFree6(interpolant, minterpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2FluxDiffVecLocal()
2152 PetscFEGeom fegeom; in DMPlexComputeClementInterpolant() local
2163 fegeom.dimEmbed = cdim; in DMPlexComputeClementInterpolant()
2188 …Nc, &interpolant, cdim * Nq, &coords, Nq, &fegeom.detJ, cdim * cdim * Nq, &fegeom.J, cdim * cdim *… in DMPlexComputeClementInterpolant()
2206 …PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)… in DMPlexComputeClementInterpolant()
2224 const PetscReal wt = quadWeights[q] * fegeom.detJ[q]; in DMPlexComputeClementInterpolant()
2227 qgeom.dimEmbed = fegeom.dimEmbed; in DMPlexComputeClementInterpolant()
2228 qgeom.J = &fegeom.J[q * cdim * cdim]; in DMPlexComputeClementInterpolant()
2229 qgeom.invJ = &fegeom.invJ[q * cdim * cdim]; in DMPlexComputeClementInterpolant()
2230 qgeom.detJ = &fegeom.detJ[q]; in DMPlexComputeClementInterpolant()
2231 …fegeom.detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for eleme… in DMPlexComputeClementInterpolant()
2255 PetscCall(PetscFree6(valsum, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeClementInterpolant()
2285 PetscFEGeom fegeom; in DMPlexComputeGradientClementInterpolant() local
2296 fegeom.dimEmbed = coordDim; in DMPlexComputeGradientClementInterpolant()
2321 …lant, coordDim * Nq, &coords, Nq, &fegeom.detJ, coordDim * coordDim * Nq, &fegeom.J, coordDim * co… in DMPlexComputeGradientClementInterpolant()
2338 …PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)… in DMPlexComputeGradientClementInterpolant()
2358 qgeom.dimEmbed = fegeom.dimEmbed; in DMPlexComputeGradientClementInterpolant()
2359 qgeom.J = &fegeom.J[q * coordDim * coordDim]; in DMPlexComputeGradientClementInterpolant()
2360 qgeom.invJ = &fegeom.invJ[q * coordDim * coordDim]; in DMPlexComputeGradientClementInterpolant()
2361 qgeom.detJ = &fegeom.detJ[q]; in DMPlexComputeGradientClementInterpolant()
2362 …fegeom.detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for eleme… in DMPlexComputeGradientClementInterpolant()
2368 …= 0; d < coordDim; ++d) grad[fc * coordDim + d] += interpolant[fc * dim + d] * wt * fegeom.detJ[q]; in DMPlexComputeGradientClementInterpolant()
2370 vol += quadWeights[q * qNc] * fegeom.detJ[q]; in DMPlexComputeGradientClementInterpolant()
2397 PetscCall(PetscFree6(gradsum, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeGradientClementInterpolant()