DMPlex: DMPlexCheckSkeleton() now handles hybrid cells

DMPlex ex4: All refinement tests passing

DMPlex: Added hybrid hex refinement

DMPlex: Bug fix for hybrid tet refinement
- Wrong offset

DMPlex: Fixed errors in hybrid tet refinement
- Use fMax instead of fStart for offset
- Use p instead of fEnd for point value
- 2 faces are created from each hybrid face, not 4

DMPlex: Corrected refinement coordinates for hybrid faces

DMPlex ex4: Hybrid tet refinement now working

DMPlex: Propagate fixes from tet refinement to hybrid tet refinement

DMPlex: Missed malloc conversion

Merge branch 'knepley/feature-plex-hybrid-3d'

* knepley/feature-plex-hybrid-3d: (45 commits)
DMPlex ex5: Check faces on all meshes
DMPlex ex5: Fix consistency of original meshes - Fix orientatio

Merge branch 'knepley/feature-plex-hybrid-3d'

* knepley/feature-plex-hybrid-3d: (45 commits)
DMPlex ex5: Check faces on all meshes
DMPlex ex5: Fix consistency of original meshes - Fix orientations of tet meshes - Redo hex test to automatically interpolate - Replaced output
DMPlex ex5: Check original meshes for consistency
DMPlex: make non-namespaced CellRefiner functions static
DMPlex: Fix hybrid bounds in DMPlexCheckFaces()
DMPlex ex5: Now all test working - Two faults working correctly
DMPlex: Fix for hybrid meshes with multiple surfaces - Move the original hybrid cells to the end of the numbering, so that split+hybrid points are inserted as a block - Use DMPlexShiftPoint_Internal() uniformly
SNES ex4: Fix for DMCheckSkeleton()
DMPlex ex4: All non-hybrid refinement test pass (Fuck Yeah)
DMPlex: In tet refinement, fixed edge orientation in internal face
DMPlex ex4: Correct output after cross-tet edge orientation fix
DMPlex: In tet refinement, correct internal tet face orientations
DMPlex: In tet refinement, correct internal face cones
DMPlex: In tet refinement, correct check for cross-tet edge
DMPlex: Now DMPlexCheckFaces() also checks all subfaces
DMPlex: Added 1D stuff to DMPlexGetRawFaces_Internal()
DMPlex ex4: hex refinement now fixed
DMPlex: Fixed split face supports for tet refinement
DMPlex: Multiple refinement for hex now working - Fixed face orientations for hexes - Fixed split face orientations - Fixed split face supports - Fixed internal face orientations - Fixed internal edge supports
DMPlex: Added GetQuadEdgeInverse_Static() and GetQuadSubfaceInverse_Static() for refinement
...

Conflicts:
src/dm/impls/plex/examples/tests/ex4.c
src/dm/impls/plex/plexsubmesh.c

show more ...

Merge branch 'knepley/fix-plex-submesh-parallel'

* knepley/fix-plex-submesh-parallel:
DMPlex: Fixed leak in SF creation for submesh - Changed name of subpointIS to avoid confusion
DMPlex: Fixed

Merge branch 'knepley/fix-plex-submesh-parallel'

* knepley/fix-plex-submesh-parallel:
DMPlex: Fixed leak in SF creation for submesh - Changed name of subpointIS to avoid confusion
DMPlex: Fixed SF creation in hybrid tet refinement, added better check
DMPlex: Guard for empty submesh in parallel during subSF creation
DMPlex: Make DMPlexCreateCohesiveSubmesh_Interpolated() also create PetscSF
PetscSection: Fix for PetscSectionCreateGlobalSectionCensored() - Offset was wrong
DM: Added check in DMCreateDefaultSF() - Check for overrun of largest offset
PetscSection: Fixed PetscSectionGetOffsetRange() for empty section
DMPlex: Small doc fix for DMPlexConstructCohesiveCells()
DMPlex: Protect another instance of empty IS
DMPlex: Protect submesh selection when some procs have empty labels
DMPlex: Turn off check for cells bordering a cohesive face - In parallel, the negative side cell for a cohesive face can be on another process - We might want to consider communicating these

show more ...

DMPlex ex5: Check faces on all meshes

DMPlex ex5: Fix consistency of original meshes
- Fix orientations of tet meshes
- Redo hex test to automatically interpolate
- Replaced output

DMPlex ex5: Check original meshes for consistency

DMPlex: Slience bogus compiler warning

DMPlex: make non-namespaced CellRefiner functions static

DMPlex: Fix hybrid bounds in DMPlexCheckFaces()

DMPlex ex5: Now all test working
- Two faults working correctly

DMPlex: Fix for hybrid meshes with multiple surfaces
- Move the original hybrid cells to the end of the numbering, so that split+hybrid points are inserted as a block
- Use DMPlexShiftPoint_Internal(

DMPlex: Fix for hybrid meshes with multiple surfaces
- Move the original hybrid cells to the end of the numbering, so that split+hybrid points are inserted as a block
- Use DMPlexShiftPoint_Internal() uniformly

show more ...

SNES ex4: Fix for DMCheckSkeleton()

Merge commit '58723a9' into knepley/feature-plex-hybrid-3d

* commit '58723a9':
DMPlex: Fix DMPlexCheckSkeleton() for submeshes
DMPlex ex5: Added second fault for quad test
DMPlex: Fix setting

Merge commit '58723a9' into knepley/feature-plex-hybrid-3d

* commit '58723a9':
DMPlex: Fix DMPlexCheckSkeleton() for submeshes
DMPlex ex5: Added second fault for quad test
DMPlex: Fix setting of hybrid bounds for multiple faults
DMPlex ex5: Turn on second fault in quad test
DMPlex: Removed compiler warnings
DMPlex ex5: Cleaned up 2D tensor tests to rely on interpolation
DMPlex ex5: Cleaned up 2D simplex tests to rely on interpolation

Conflicts:
include/petscdmplex.h
src/dm/impls/plex/examples/tests/ex4.c
src/dm/impls/plex/plexrefine.c

show more ...

DMPlex: fix stale docs about DMPlexGetCoordinateSection

DMPlex Gmsh: Fixed integer size problem

DMPlex Gmsh: Fixed interface for coordinates

Sys: drop explicit type arguments from PetscNew() and PetscNewLog()

git grep -l 'PetscNew\(Log\)\?(' -- $1 | xargs perl -pi -e '
s@PetscNew\([^,;()]+ *, *@PetscNew(@;
s@PetscNewLog\(([^,;()]+) *,[^,

Sys: drop explicit type arguments from PetscNew() and PetscNewLog()

git grep -l 'PetscNew\(Log\)\?(' -- $1 | xargs perl -pi -e '
s@PetscNew\([^,;()]+ *, *@PetscNew(@;
s@PetscNewLog\(([^,;()]+) *,[^,;()]+, *@PetscNewLog($1,@'

show more ...