Silent overzelous warnings "may be used uninitialized in this function [-Wmaybe-uninitialized]"

Merge branch 'knepley/fix-plex-l2-hybrid' into 'master'

Plex: L2FieldDiff should only integrate the hybrid field on those cells

See merge request petsc/petsc!2842

Plex ex37: Nozzle test works
- Fixed reader for quads and bi-segments
- Added nozzle mesh

Plex+EGADS: Update to geometry interpolation

Plex: All refinement now looks at the CAD model

Plex: L2FieldDiff should only integrate the hybrid field on those cells

Merge remote-tracking branch 'origin/maint'

removed extra getpointsf

Plex ex5: Turn on tensor discretization

Plex: Labels stored in DM fields and DSes need to be refined during uniform refinement

Plex+SNES: Large update to support multiple DSes
- Now allow simple divisions of the domain with DMSetField()
- DMPlexComputeResidual_Internal() supports multiple DSes
- DMPlexComputeJacobian_Interna

Plex+SNES: Large update to support multiple DSes
- Now allow simple divisions of the domain with DMSetField()
- DMPlexComputeResidual_Internal() supports multiple DSes
- DMPlexComputeJacobian_Internal() supports multiple DSes
- DMSNESCheckDisretization() supports multiple DSes
- DMSNESCheckResidual() supports multiple DSes
- DMSNESCheckJacoiban() supports multiple DSes
- DMComputeL2FieldDiff() supports multiple DSes
- PetscConveEst supports multiple DSes
- Fixed bug refering to cStart from a cellIS
- Added PetscDSGetQuadrature()
- Added DMSetField_Internal()
- Added DMSetRegionNumDS()
- Protect PetscSectionView() from empty bcIndices
- ISDifference now allows a NULL second argument

show more ...

Merge remote-tracking branch 'origin/maint'

Plex ex5: Hybrid residuals are now working
- Make functions static
- Put in a realistic initial guess
- Added hybrid Jacobian test
- Remove use of hybrid bounds and update outupt

Plex: More hybrid fixes
- Converted hybrid assembly to use PetscTabulation
- Protect for NULL discretizations, which happen with no dofs for a field
- Added PetscDiscType enum
- Changed PetscDSIsFE_I

Plex: More hybrid fixes
- Converted hybrid assembly to use PetscTabulation
- Protect for NULL discretizations, which happen with no dofs for a field
- Added PetscDiscType enum
- Changed PetscDSIsFE_Internal() to PetscDSGetDiscType_Internal()
- Must inject coordinates from hybrid cells into the ambient space

show more ...

DM: Handle multiple DSes
- Fixed CreateSubDM()
- Made DMCreateDS() idempotent
- When creating section, must read BC from all DSes

Plex+FE: Updates for hybrid integration
- Added support for hybrid Jacobians
- Added PetscDSHasBdJacobian()
- Added boundary Jacobian preconditioner
- Needed a path to integrate the hybrid Jacobian p

Plex+FE: Updates for hybrid integration
- Added support for hybrid Jacobians
- Added PetscDSHasBdJacobian()
- Added boundary Jacobian preconditioner
- Needed a path to integrate the hybrid Jacobian preconditioner
- Fixed indexing in hybrid element mat update
- Be careful about the difference between intrinsic dimension and embedding dimension for hybrid integration
- Need to use rectangular transformation of gradient in case they are different
- Fix use of embedding dimension
- Make sure FEGeom struct has the dimensions set
- Pointwise function derivatives are in the embedding space
- Use PetscArrayzero()
- Change TransitiveClosure() for hybrid cells
- We get the closures of the endcaps, and then the hybrid faces, this allows degenerate cells that we need in PyLith

show more ...

Propagate errors from HPDDM into KSPHPDDM

DMPlexComputeProjection3Dto2D: simplify implementation; make safe with -ffast-math

I was getting NaN with -ffast-math due to speculative divide by
zero/early exit (seems like a miscompile), but this

DMPlexComputeProjection3Dto2D: simplify implementation; make safe with -ffast-math

I was getting NaN with -ffast-math due to speculative divide by
zero/early exit (seems like a miscompile), but this implementation is
much simpler and avoids that issue.

show more ...

Merge branch 'knepley/feature-plex-interior-ghost' into 'master'

Knepley/feature plex interior ghost

See merge request petsc/petsc!2789

Plex+Doc: Added additions to Changes and more documentation

Plex: Fix cell FE geometry

Plex: Allow projection between bases on the same mesh

Plex: Abstract getting dof indices on closure
- Added DMPlexMatSetClosureGeneral()
- Changed interface for existing DMPlexGetClosureIndices()

Plex: Draw viewer for higher order geometry in 2D

Plex: This fixes a pernicious bug with projection - When setting points along a label, we must filter out points in the closure that are not also in the label.