PetscDualSpace: added PetscDualSpaceGetHeightSubspace()

We already acknowledge that a functional in the dual space may have its
quadrature points limited to a mesh point in the boundary of a cell.
T

PetscDualSpace: added PetscDualSpaceGetHeightSubspace()

We already acknowledge that a functional in the dual space may have its
quadrature points limited to a mesh point in the boundary of a cell.
This goes one step further and reproduces those functionals in terms of
a reference element for the lower-dimensional mesh points. This was
written to be used for function projection for non-conformal meshes:
a global dof may not align with any local cell dof, but may align with a
facet or edge dof. Something like this would also be necessary for
implementing Scott-Zhang interpolation
(http://dx.doi.org/10.1090/S0025-5718-1990-1011446-7).

show more ...

PetscDS: Fix up some functions to work for both FV and FE

PetscFV: Change PetscFVIntegrateRHSFunction() to get info from PetscDS

PetscFE: Do not automatically set tensor space for tensor cell
- This allows it to be set from the command line, but you have to remember to set it
- Also allows discontinuous polynomial spaces on te

PetscFE: Do not automatically set tensor space for tensor cell
- This allows it to be set from the command line, but you have to remember to set it
- Also allows discontinuous polynomial spaces on tensor cells

show more ...

Quadrature+Doc: Add function docs

Quadrature: Add an independent idea of quadrature order
- This is not completely correct yet, since we have to rationalize the creation routines to take this as input

PetscFE: PetscBool fix

DMPlex: make DMPlexCreateReferenceCell

dfte: fix TensorPoint_Internal()

There was a bad recurrence

dfte: fix dual space set up for hexahedral meshes

Plex+FE: Removed all use of PetscCellGeometry, Changed FE integration to use PetscFECellGeom

PetscDualSpace: Remove use of PetscCellGeometry in PetscDualSpaceApply()
- Fixed up calls to it for projection into spaces

PetscFV+TS: Removed use of PetscCellGeometry from PetscFVIntegrateRHSFunction()
- Moving toward cached geometry with arrays of structs
- Changed DMPlexTSGetGeometry() to DMPlexTSGetGeometryFVM()

PetscFE: Fortran bindings

PetscDS: Added Fortran wrappers for tabulation

Quadrature: Added Fortran bindings

PetscDS: Fix Fortran bindings

PetscDS+Doc: Added function docs and made some static

PetscDS: Added Riemann solvers and contexts for each field

PetscFV+Doc: Fixed comments

PetscFV: Now can hold a quadrature object
- This allows projection onto the piecewise-constant space

PetscDualSpace: Added a Simple space
- Its just a collection of functionals

Merge branch 'knepley/feature-dm-dimension-handling'

* knepley/feature-dm-dimension-handling:
DM: Fixed DMDASetDim() in examples
Plex: Fix for isoparametric geometry
Plex: Complex fix
TS ex1

Merge branch 'knepley/feature-dm-dimension-handling'

* knepley/feature-dm-dimension-handling:
DM: Fixed DMDASetDim() in examples
Plex: Fix for isoparametric geometry
Plex: Complex fix
TS ex11: Fixed calls to DMSetCoordinateSection()
DM: FIxed DMClone() to copy dimension
Plex+DA: Now cell geometry routine can take a PetscFE in order to evaluate at quadrature points - We evaluate isoparametrically at the quadrature points - Changed DMDAComputeCellGeometry() --> DMDAComputeCellGeometryFEM() - Changed DMPlexComputeCellGeometry() --> DMPlexComputeCellGeometryFEM() - Added DMPlexComputeCellGeometryAffineFEM()
DM: Added a embedding dimension, different from the topological dimension - This can now be set in DMSetCoordinateSection(), or it can be automatically determined
Plex: Fix output for ex4
DM: Added DMGetDimPoints() - This replicates DMPlex/DAGetDepthStratum(), but uses the dimension which people have said is more intuitive
DM: Move dimension up into DM - Both DA and Plex had this and I think it makes sense

Conflicts:
src/dm/impls/plex/plexcreate.c

show more ...

Sys: PetscCheckPointer PETSC_CHAR checks char*, not char**

PetscFV: Fully enable other limiters