Lines Matching refs:geom
457 PetscFEGeom *geom; in DMFieldEvaluateFV_DS() local
487 PetscCall(DMFieldCreateFEGeom(coordField, pointIS, quad, PETSC_FEGEOM_BASIC, &geom)); in DMFieldEvaluateFV_DS()
507 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
508 …for (j = 0; j < Nc; j++) sB[i * Nc + j] += geom->detJ[i * Nq + k] * weights[k] * sqB[(i * Nq + k) … in DMFieldEvaluateFV_DS()
521 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
540 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
545 …for (m = 0; m < dim; m++) pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * sqD[((i * Nq… in DMFieldEvaluateFV_DS()
547 …for (l = 0; l < dimC; l++) sD[(i * Nc + j) * dimC + l] += geom->detJ[i * Nq + k] * weights[k] * pD… in DMFieldEvaluateFV_DS()
561 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
566 …for (m = 0; m < dim; m++) pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * rqD[((i * Nq… in DMFieldEvaluateFV_DS()
568 …for (l = 0; l < dimC; l++) rD[(i * Nc + j) * dimC + l] += geom->detJ[i * Nq + k] * weights[k] * pD… in DMFieldEvaluateFV_DS()
587 const PetscReal *invJ = &geom->invJ[(i * Nq + k) * dimC * dimC]; in DMFieldEvaluateFV_DS()
589 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
606 …for (m = 0; m < dimC; m++) sH[(i * Nc + j) * dimC * dimC + l * dimC + m] += geom->detJ[i * Nq + k]… in DMFieldEvaluateFV_DS()
621 const PetscReal *invJ = &geom->invJ[(i * Nq + k) * dimC * dimC]; in DMFieldEvaluateFV_DS()
623 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
640 …for (m = 0; m < dimC; m++) rH[(i * Nc + j) * dimC * dimC + l * dimC + m] += geom->detJ[i * Nq + k]… in DMFieldEvaluateFV_DS()
651 PetscCall(PetscFEGeomDestroy(&geom)); in DMFieldEvaluateFV_DS()
786 …rCode DMFieldComputeFaceData_DS(DMField field, IS pointIS, PetscQuadrature quad, PetscFEGeom *geom) in DMFieldComputeFaceData_DS() argument
796 dim = geom->dim; in DMFieldComputeFaceData_DS()
797 dE = geom->dimEmbed; in DMFieldComputeFaceData_DS()
802 numFaces = geom->numCells; in DMFieldComputeFaceData_DS()
803 Nq = geom->numPoints; in DMFieldComputeFaceData_DS()
823 geom->face[p][s * 2 + 0] = c; in DMFieldComputeFaceData_DS()
824 geom->face[p][s * 2 + 1] = ornt[c]; in DMFieldComputeFaceData_DS()
828 if (geom->face[p][1] < 0) { in DMFieldComputeFaceData_DS()
829 PetscInt Np = geom->numPoints, q, dE = geom->dimEmbed, d; in DMFieldComputeFaceData_DS()
832 for (d = 0; d < dE; ++d) geom->n[(p * Np + q) * dE + d] = -geom->n[(p * Np + q) * dE + d]; in DMFieldComputeFaceData_DS()
893 geom->suppJ[s][p * Nq * dE * dE + q] = cellGeom->J[offset * Nq * dE * dE + q]; in DMFieldComputeFaceData_DS()
894 geom->suppInvJ[s][p * Nq * dE * dE + q] = cellGeom->invJ[offset * Nq * dE * dE + q]; in DMFieldComputeFaceData_DS()
896 for (q = 0; q < Nq; q++) geom->suppDetJ[s][p * Nq + q] = cellGeom->detJ[offset * Nq + q]; in DMFieldComputeFaceData_DS()
1008 for (q = 0; q < Nq; q++) orientPoints[o][q] = -geom->xi[q]; in DMFieldComputeFaceData_DS()
1010 for (q = 0; q < Nq; q++) orientPoints[o][q] = geom->xi[q]; in DMFieldComputeFaceData_DS()
1019 lambda[0] = -(geom->xi[2 * q] + geom->xi[2 * q + 1]) / 2.; in DMFieldComputeFaceData_DS()
1020 lambda[1] = (geom->xi[2 * q] + 1.) / 2.; in DMFieldComputeFaceData_DS()
1021 lambda[2] = (geom->xi[2 * q + 1] + 1.) / 2.; in DMFieldComputeFaceData_DS()
1037 xi[0] = geom->xi[2 * q]; in DMFieldComputeFaceData_DS()
1038 xi[1] = geom->xi[2 * q + 1]; in DMFieldComputeFaceData_DS()
1105 geom->suppJ[s][p * Nq * dE * dE + q] = cellGeom->J[offset * Nq * dE * dE + q]; in DMFieldComputeFaceData_DS()
1106 geom->suppInvJ[s][p * Nq * dE * dE + q] = cellGeom->invJ[offset * Nq * dE * dE + q]; in DMFieldComputeFaceData_DS()
1108 … for (q = 0; q < Nq; q++) geom->suppDetJ[s][p * Nq + q] = cellGeom->detJ[offset * Nq + q]; in DMFieldComputeFaceData_DS()