Lines Matching refs:volume
110 …w, PetscReal *phi, PetscReal *dphidx, PetscReal *jacobian, PetscReal *ijacobian, PetscReal *volume) in Compute_Lagrange_Basis_1D_Internal() argument
117 PetscAssertPointer(volume, 11); in Compute_Lagrange_Basis_1D_Internal()
133 jacobian[0] = ijacobian[0] = volume[0] = 0.0; in Compute_Lagrange_Basis_1D_Internal()
141 *volume = jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
143 jxw[j] *= *volume; in Compute_Lagrange_Basis_1D_Internal()
162 jacobian[0] = ijacobian[0] = volume[0] = 0.0; in Compute_Lagrange_Basis_1D_Internal()
170 *volume = jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
172 if (jxw) jxw[j] *= *volume; in Compute_Lagrange_Basis_1D_Internal()
220 …PetscReal *dphidx, PetscReal *dphidy, PetscReal *jacobian, PetscReal *ijacobian, PetscReal *volume) in Compute_Lagrange_Basis_2D_Internal() argument
227 PetscAssertPointer(volume, 12); in Compute_Lagrange_Basis_2D_Internal()
265 PetscCall(DMatrix_Invert_2x2_Internal(jacobian, ijacobian, volume)); in Compute_Lagrange_Basis_2D_Internal()
266 …volume >= 1e-12, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Quadrangular element has zero volume… in Compute_Lagrange_Basis_2D_Internal()
268 if (jxw) jxw[j] *= *volume; in Compute_Lagrange_Basis_2D_Internal()
293 PetscCall(DMatrix_Invert_2x2_Internal(jacobian, ijacobian, volume)); in Compute_Lagrange_Basis_2D_Internal()
294 …volume >= 1e-12, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Triangular element has zero volume: … in Compute_Lagrange_Basis_2D_Internal()
376 …PetscReal *dphidy, PetscReal *dphidz, PetscReal *jacobian, PetscReal *ijacobian, PetscReal *volume) in Compute_Lagrange_Basis_3D_Internal() argument
383 PetscAssertPointer(volume, 13); in Compute_Lagrange_Basis_3D_Internal()
437 PetscCall(DMatrix_Invert_3x3_Internal(jacobian, ijacobian, volume)); in Compute_Lagrange_Basis_3D_Internal()
438 …volume >= 1e-8, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Hexahedral element has zero volume: %… in Compute_Lagrange_Basis_3D_Internal()
440 if (jxw) jxw[j] *= (*volume); in Compute_Lagrange_Basis_3D_Internal()
470 PetscCall(DMatrix_Invert_3x3_Internal(jacobian, ijacobian, volume)); in Compute_Lagrange_Basis_3D_Internal()
471 …volume >= 1e-8, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Tetrahedral element has zero volume: … in Compute_Lagrange_Basis_3D_Internal()
474 …* 3] - coords[2 + 3 * 3]) - coords[1 + 3 * 3] * (coords[2 + 1 * 3] - coords[2 + 2 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
475 …* 3] - coords[2 + 3 * 3]) - coords[1 + 3 * 3] * (coords[2 + 0 * 3] - coords[2 + 2 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
476 …* 3] - coords[2 + 3 * 3]) - coords[1 + 3 * 3] * (coords[2 + 0 * 3] - coords[2 + 1 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
478 …* 3] - coords[2 + 3 * 3]) + coords[0 + 3 * 3] * (coords[2 + 1 * 3] - coords[2 + 2 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
479 …* 3] - coords[2 + 3 * 3]) + coords[0 + 3 * 3] * (coords[2 + 0 * 3] - coords[2 + 2 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
480 …* 3] - coords[2 + 3 * 3]) + coords[0 + 3 * 3] * (coords[2 + 0 * 3] - coords[2 + 1 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
482 …* 3] - coords[1 + 1 * 3]) + coords[0 + 3 * 3] * (coords[1 + 2 * 3] - coords[1 + 1 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
483 …* 3] - coords[1 + 3 * 3]) - coords[0 + 3 * 3] * (coords[1 + 0 * 3] - coords[1 + 2 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
484 …* 3] - coords[1 + 3 * 3]) - coords[0 + 3 * 3] * (coords[1 + 0 * 3] - coords[1 + 1 * 3])) / *volume; in Compute_Lagrange_Basis_3D_Internal()
494 if (jxw) jxw[j] *= *volume; in Compute_Lagrange_Basis_3D_Internal()
560 PetscReal jacobian[9], ijacobian[9], volume; in DMMoabFEMComputeBasis() local
575 …ht_product, fe_basis, compute_der ? fe_basis_derivatives[0] : NULL, jacobian, ijacobian, &volume)); in DMMoabFEMComputeBasis()
578 …erivatives[0] : NULL, compute_der ? fe_basis_derivatives[1] : NULL, jacobian, ijacobian, &volume)); in DMMoabFEMComputeBasis()
581 …erivatives[1] : NULL, compute_der ? fe_basis_derivatives[2] : NULL, jacobian, ijacobian, &volume)); in DMMoabFEMComputeBasis()
678 …tscReal *phypts, PetscReal *phibasis, PetscReal *jacobian, PetscReal *ijacobian, PetscReal *volume) in FEMComputeBasis_JandF() argument
683 …al(nverts, coordinates, 1, quadrature, phypts, NULL, phibasis, NULL, jacobian, ijacobian, volume)); in FEMComputeBasis_JandF()
686 …rts, coordinates, 1, quadrature, phypts, NULL, phibasis, NULL, NULL, jacobian, ijacobian, volume)); in FEMComputeBasis_JandF()
689 …oordinates, 1, quadrature, phypts, NULL, phibasis, NULL, NULL, NULL, jacobian, ijacobian, volume)); in FEMComputeBasis_JandF()
725 PetscReal phibasis[8], jacobian[9], ijacobian[9], volume; in DMMoabPToRMapping() local
744 …teBasis_JandF(dim, nverts, coordinates, natparam, phypts, phibasis, jacobian, ijacobian, &volume)); in DMMoabPToRMapping()
780 …teBasis_JandF(dim, nverts, coordinates, natparam, phypts, phibasis, jacobian, ijacobian, &volume)); in DMMoabPToRMapping()