Merge branch 'haplav/fix-dmplexinterpolatepointsf-cannot-handle-face' into 'master'

fix dmplexinterpolatepointsf cannot handle face

See merge request petsc/petsc!2478

doc: add some more cross-references

Fix DMPlexInterpolatePointSF().

Actually DMPlexAddSharedFace_Private().
Cells can't be shared for non-overlapping meshes.

Add support for users to pass in FFLAGS the option for promoting integers to 8 bytes.

Tested for ifort and gfortran

Tested with MPICH and OpenMPI and MPI-Uni

Since MPI does not provide any mechani

Add support for users to pass in FFLAGS the option for promoting integers to 8 bytes.

Tested for ifort and gfortran

Tested with MPICH and OpenMPI and MPI-Uni

Since MPI does not provide any mechanism for integer promotion the code very crudely adjusts the
values in the mpif.h to mark them as 4 byte integers this is not neccessary portable for all MPI implementations

This is an extremely fragile option that we should probably not be supporting.

Users must insure their Fortran codes passes only 4 byte integers to MPI calls

PETSc Fortran examples were all fixed to explicitly use 4 byte integers for MPI calls.

Commit-type: portability-fix, testing-fix, feature
Time: 20 hours
Reported-by: Dmitry Melnichuk <dmitry.melnichuk@geosteertech.com>

show more ...

add DMPlexGetPointHeight()

ex18: test reproducing InterpolateSF fail

DMPlexInterpolatePointSF: improve debug/error message

Merge branch 'haplav/feature-dmplex-ex18-logging-improvements' into 'master'

dmplex ex18 logging improvements

See merge request petsc/petsc!2475

Merge branch 'knepley/feature-fem-pk' into 'master'

Knepley/feature fem pk

See merge request petsc/petsc!2464

Plex test ex18 improvements.

* -interpolate {none,create,after_create,after_distribute}
* add logstages create,distribute,interpolate
* factor out CreateMeshFromFile from CreateMesh
* improve handli

Plex test ex18 improvements.

* -interpolate {none,create,after_create,after_distribute}
* add logstages create,distribute,interpolate
* factor out CreateMeshFromFile from CreateMesh
* improve handling -rotate_interface_{0,1}

show more ...

DM: Handle serial case for Global-to-Natural mapping

Plex: Hack for RemapGeometry()
- Right now, we cannot have a coordinate field on the coordinate DM, so
borrow the one from the DM

DMPlex: Use coordinate dimension to check for field type

Update output files for new Cholesky ordering default

Reported-by: gitlab-ci

Merge branch 'jolivet/fix-pardiso-complex-symmetric' into 'maint'

Activate PARDISO for complex symmetric matrices

See merge request petsc/petsc!2437

Merge remote-tracking branch 'origin/maint'

Merge branch 'patch-2' into 'master'

Fix plexpartition

See merge request petsc/petsc!2445

DMPlexCreateBoxMesh: document -dm_plex_box_faces

Commit-type: documentation
Reported-by: James Wright <jrwrigh.iii@gmail.com>

Apply suggestion to src/dm/impls/plex/plexpartition.c

Doc: Text fixes

DM+Plex: Updated test output for cell types

DM+Plex: Now classify cells by polytope type
- Defined the polytope type of a cell. added names and special shortcut
pointer like depth - Fixed DMLabelMakeValid() to optimize for contiguous
runs of

DM+Plex: Now classify cells by polytope type
- Defined the polytope type of a cell. added names and special shortcut
pointer like depth - Fixed DMLabelMakeValid() to optimize for contiguous
runs of points - Added cell polytope type label to DM - This needs
special processing, like depth - Added DMPlexGetPointDepth(),
DMPlexGetCellType(), DMPlexGetCellTypeLabel() - DMPlexComputeCellTypes()
infers the cell type for all points using dim/depth/coneSize/ghostCells
- DMPlexGetRawFace(Hybrid)_Internal() now takes the cell type. This
means we can unify the functions later. - Replaced switch(dim/coneSize)
constructs with cell type

show more ...

Merge branch 'knepley/feature-project-different-spaces' into 'master'

Knepley/feature project different spaces

See merge request petsc/petsc!2153

Plex ex23: Added tests for subdomain and multifield projection
- Plex ex22: Fix reference loop
- Must have DM set for projection

FE+DS: Introduce PetscTabulation to hold tabulated function data
- This was created to prevent mismatches in sizes I was getting by just
passing around arrays
- FE+FV+DS now hold tabulations
- GetTab

FE+DS: Introduce PetscTabulation to hold tabulated function data
- This was created to prevent mismatches in sizes I was getting by just
passing around arrays
- FE+FV+DS now hold tabulations
- GetTabulation() --> CreateTabulation() and it creates an object rather
than using DMGetWorkArray(). I put in ComputeTabulation() to reuse
space, but there still may be performance regression.
- Default tabulations are now "cell tabulations"
- The replicas in the tabulation are used to index faces in a face
tabulation

show more ...