Lines Matching refs:nverts
110 static PetscErrorCode Compute_Lagrange_Basis_1D_Internal(const PetscInt nverts, const PetscReal *co… in Compute_Lagrange_Basis_1D_Internal() argument
120 PetscCall(PetscArrayzero(dphidx, npts * nverts)); in Compute_Lagrange_Basis_1D_Internal()
122 if (nverts == 2) { /* Linear Edge */ in Compute_Lagrange_Basis_1D_Internal()
125 const PetscInt offset = j * nverts; in Compute_Lagrange_Basis_1D_Internal()
134 for (i = 0; i < nverts; ++i) { in Compute_Lagrange_Basis_1D_Internal()
146 for (i = 0; i < nverts; i++) { in Compute_Lagrange_Basis_1D_Internal()
150 } else if (nverts == 3) { /* Quadratic Edge */ in Compute_Lagrange_Basis_1D_Internal()
153 const PetscInt offset = j * nverts; in Compute_Lagrange_Basis_1D_Internal()
163 for (i = 0; i < nverts; ++i) { in Compute_Lagrange_Basis_1D_Internal()
175 for (i = 0; i < nverts; i++) { in Compute_Lagrange_Basis_1D_Internal()
179 …nvalid. Currently only support EDGE2 and EDGE3 basis evaluations in 1-D : %" PetscInt_FMT, nverts); in Compute_Lagrange_Basis_1D_Internal()
220 static PetscErrorCode Compute_Lagrange_Basis_2D_Internal(const PetscInt nverts, const PetscReal *co… in Compute_Lagrange_Basis_2D_Internal() argument
231 PetscCall(PetscArrayzero(dphidx, npts * nverts)); in Compute_Lagrange_Basis_2D_Internal()
232 PetscCall(PetscArrayzero(dphidy, npts * nverts)); in Compute_Lagrange_Basis_2D_Internal()
234 if (nverts == 4) { /* Linear Quadrangle */ in Compute_Lagrange_Basis_2D_Internal()
237 const PetscInt offset = j * nverts; in Compute_Lagrange_Basis_2D_Internal()
251 for (i = 0; i < nverts; ++i) { in Compute_Lagrange_Basis_2D_Internal()
272 for (i = 0; i < nverts; i++) { in Compute_Lagrange_Basis_2D_Internal()
280 } else if (nverts == 3) { /* Linear triangle */ in Compute_Lagrange_Basis_2D_Internal()
300 const PetscInt offset = j * nverts; in Compute_Lagrange_Basis_2D_Internal()
332 …invalid. Currently only support QUAD4 and TRI3 basis evaluations in 2-D : %" PetscInt_FMT, nverts); in Compute_Lagrange_Basis_2D_Internal()
376 static PetscErrorCode Compute_Lagrange_Basis_3D_Internal(const PetscInt nverts, const PetscReal *co… in Compute_Lagrange_Basis_3D_Internal() argument
388 PetscCall(PetscArrayzero(dphidx, npts * nverts)); in Compute_Lagrange_Basis_3D_Internal()
389 PetscCall(PetscArrayzero(dphidy, npts * nverts)); in Compute_Lagrange_Basis_3D_Internal()
390 PetscCall(PetscArrayzero(dphidz, npts * nverts)); in Compute_Lagrange_Basis_3D_Internal()
393 if (nverts == 8) { /* Linear Hexahedra */ in Compute_Lagrange_Basis_3D_Internal()
396 const PetscInt offset = j * nverts; in Compute_Lagrange_Basis_3D_Internal()
418 for (i = 0; i < nverts; ++i) { in Compute_Lagrange_Basis_3D_Internal()
443 for (i = 0; i < nverts; ++i) { in Compute_Lagrange_Basis_3D_Internal()
451 } else if (nverts == 4) { /* Linear Tetrahedra */ in Compute_Lagrange_Basis_3D_Internal()
489 const PetscInt offset = j * nverts; in Compute_Lagrange_Basis_3D_Internal()
502 for (i = 0; i < nverts; ++i) { in Compute_Lagrange_Basis_3D_Internal()
529 … invalid. Currently only support HEX8 and TET4 basis evaluations in 3-D : %" PetscInt_FMT, nverts); in Compute_Lagrange_Basis_3D_Internal()
555 PetscErrorCode DMMoabFEMComputeBasis(const PetscInt dim, const PetscInt nverts, const PetscReal *co… in DMMoabFEMComputeBasis() argument
575 …PetscCall(Compute_Lagrange_Basis_1D_Internal(nverts, coordinates, npoints, quadpts, phypts, jacobi… in DMMoabFEMComputeBasis()
578 …PetscCall(Compute_Lagrange_Basis_2D_Internal(nverts, coordinates, npoints, quadpts, phypts, jacobi… in DMMoabFEMComputeBasis()
581 …PetscCall(Compute_Lagrange_Basis_3D_Internal(nverts, coordinates, npoints, quadpts, phypts, jacobi… in DMMoabFEMComputeBasis()
604 PetscErrorCode DMMoabFEMCreateQuadratureDefault(const PetscInt dim, const PetscInt nverts, PetscQua… in DMMoabFEMCreateQuadratureDefault() argument
614 PetscCall(PetscDTStroudConicalQuadrature(1, nc, nverts, 0, 1.0, quadrature)); in DMMoabFEMCreateQuadratureDefault()
618 if (nverts == 3) { in DMMoabFEMCreateQuadratureDefault()
640 } else PetscCall(PetscDTGaussTensorQuadrature(dim, nc, nverts, 0.0, 1.0, quadrature)); in DMMoabFEMCreateQuadratureDefault()
644 if (nverts == 4) { in DMMoabFEMCreateQuadratureDefault()
669 } else PetscCall(PetscDTGaussTensorQuadrature(dim, nc, nverts, 0.0, 1.0, quadrature)); in DMMoabFEMCreateQuadratureDefault()
678 static PetscErrorCode FEMComputeBasis_JandF(const PetscInt dim, const PetscInt nverts, const PetscR… in FEMComputeBasis_JandF() argument
683 …PetscCall(Compute_Lagrange_Basis_1D_Internal(nverts, coordinates, 1, quadrature, phypts, NULL, phi… in FEMComputeBasis_JandF()
686 …PetscCall(Compute_Lagrange_Basis_2D_Internal(nverts, coordinates, 1, quadrature, phypts, NULL, phi… in FEMComputeBasis_JandF()
689 …PetscCall(Compute_Lagrange_Basis_3D_Internal(nverts, coordinates, 1, quadrature, phypts, NULL, phi… in FEMComputeBasis_JandF()
719 PetscErrorCode DMMoabPToRMapping(const PetscInt dim, const PetscInt nverts, const PetscReal *coordi… in DMMoabPToRMapping() argument
738 PetscCall(PetscArrayzero(phibasis, nverts)); in DMMoabPToRMapping()
744 …PetscCall(FEMComputeBasis_JandF(dim, nverts, coordinates, natparam, phypts, phibasis, jacobian, ij… in DMMoabPToRMapping()
780 …PetscCall(FEMComputeBasis_JandF(dim, nverts, coordinates, natparam, phypts, phibasis, jacobian, ij… in DMMoabPToRMapping()
796 if (phi) PetscCall(PetscArraycpy(phi, phibasis, nverts)); in DMMoabPToRMapping()