petsc-master 2020-10-26
Report Typos and Errors


Evaluate bases and derivatives at quadrature points for a EDGE2 or EDGE3 element.


#include "petscdt.h" 
#include "petscdmmoab.h"   
PetscErrorCode Compute_Lagrange_Basis_1D_Internal(const PetscInt nverts, const PetscReal *coords, const PetscInt npts, const PetscReal *quad, PetscReal *phypts,
                                                  PetscReal *jxw, PetscReal *phi, PetscReal *dphidx,
                                                  PetscReal *jacobian, PetscReal *ijacobian, PetscReal *volume)
The routine is given the coordinates of the vertices of a linear or quadratic edge element.

The routine evaluates the basis functions associated with each quadrature point provided, and their derivatives with respect to X.


Example Physical Element

    1-------2        1----3----2
      EDGE2             EDGE3

Input Parameter

PetscInt nverts - the number of element vertices
PetscReal coords[3*nverts] - the physical coordinates of the vertices (in canonical numbering)
PetscInt npts - the number of evaluation points (quadrature points)
PetscReal quad[3*npts] - the evaluation points (quadrature points) in the reference space

Output Parameter

PetscReal phypts[3*npts] - the evaluation points (quadrature points) transformed to the physical space
PetscReal jxw[npts] - the jacobian determinant * quadrature weight necessary for assembling discrete contributions
PetscReal phi[npts] - the bases evaluated at the specified quadrature points
PetscReal dphidx[npts] - the derivative of the bases wrt X-direction evaluated at the specified quadrature points




Index of all DMMOAB routines
Table of Contents for all manual pages
Index of all manual pages