137045ce4SJed Brown /* 237045ce4SJed Brown Common tools for constructing discretizations 337045ce4SJed Brown */ 437045ce4SJed Brown #if !defined(__PETSCDT_H) 537045ce4SJed Brown #define __PETSCDT_H 637045ce4SJed Brown 737045ce4SJed Brown #include <petscsys.h> 837045ce4SJed Brown 921454ff5SMatthew G. Knepley /*S 1021454ff5SMatthew G. Knepley PetscQuadrature - Quadrature rule for integration. 1121454ff5SMatthew G. Knepley 1221454ff5SMatthew G. Knepley Level: developer 1321454ff5SMatthew G. Knepley 1421454ff5SMatthew G. Knepley .seealso: PetscQuadratureCreate(), PetscQuadratureDestroy() 1521454ff5SMatthew G. Knepley S*/ 1621454ff5SMatthew G. Knepley typedef struct _p_PetscQuadrature *PetscQuadrature; 1721454ff5SMatthew G. Knepley 1821454ff5SMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscQuadratureCreate(MPI_Comm, PetscQuadrature *); 19*bcede257SMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscQuadratureGetOrder(PetscQuadrature, PetscInt*); 20*bcede257SMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscQuadratureSetOrder(PetscQuadrature, PetscInt); 2121454ff5SMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscQuadratureGetData(PetscQuadrature, PetscInt*, PetscInt*, const PetscReal *[], const PetscReal *[]); 2221454ff5SMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscQuadratureSetData(PetscQuadrature, PetscInt, PetscInt, const PetscReal [], const PetscReal []); 2321454ff5SMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscQuadratureView(PetscQuadrature, PetscViewer); 2421454ff5SMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscQuadratureDestroy(PetscQuadrature *); 25a0845e3aSMatthew G. Knepley 26a0845e3aSMatthew G. Knepley typedef struct { 271828ad4fSMatthew G. Knepley PetscReal *v0, *n, *J, *invJ, *detJ, *vol; 28a0845e3aSMatthew G. Knepley } PetscCellGeometry; 29a0845e3aSMatthew G. Knepley 3037045ce4SJed Brown PETSC_EXTERN PetscErrorCode PetscDTLegendreEval(PetscInt,const PetscReal*,PetscInt,const PetscInt*,PetscReal*,PetscReal*,PetscReal*); 3137045ce4SJed Brown PETSC_EXTERN PetscErrorCode PetscDTGaussQuadrature(PetscInt,PetscReal,PetscReal,PetscReal*,PetscReal*); 32194825f6SJed Brown PETSC_EXTERN PetscErrorCode PetscDTReconstructPoly(PetscInt,PetscInt,const PetscReal*,PetscInt,const PetscReal*,PetscReal*); 33744bafbcSMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscDTGaussTensorQuadrature(PetscInt,PetscInt,PetscReal,PetscReal,PetscQuadrature*); 34a0845e3aSMatthew G. Knepley PETSC_EXTERN PetscErrorCode PetscDTGaussJacobiQuadrature(PetscInt,PetscInt,PetscReal,PetscReal,PetscQuadrature*); 3537045ce4SJed Brown 3637045ce4SJed Brown #endif 37