DM+Plex+FE: Initial CEED integration
- Plex ex42: Surface area CEED example

DM+Plex: Boundary integration functions are now stored in a PetscWeakForm held in the Boundary object
- Boundary residual/Jacobian functions now take a form key
- Now store DMLabel in DSBoundary
- Im

DM+Plex: Boundary integration functions are now stored in a PetscWeakForm held in the Boundary object
- Boundary residual/Jacobian functions now take a form key
- Now store DMLabel in DSBoundary
- Improved Boundary, Weak Form, and DS viewers
- Remove DMCopyBoundary(), DMCopyDisc() is sufficient
- Need DMTransferDS_Internal() to create new DS with label in the correct DM
- Correct all DMAddBoundary() in tests
- Allow DMBoundary to hold just a label name if necessary since p4est is missing labels in a lot of places
- Fix propagating coordinate dimension in DS

show more ...

exodusii: updated test output for ex26 and added output for ex62, ex26f90, and ex62f90

exodusii: added fortran versions of src/dm/impls/plex/tests/ex26.c and ex62.c

exodusii: viewer updates
- finished implementing PETSCVIEWEREXODUSII
- Call the proper exodusii function when VecView/Load and DMView/Load are invoked on a PETSCVIEWEREXODUSII viewer
- Added

exodusii: viewer updates
- finished implementing PETSCVIEWEREXODUSII
- Call the proper exodusii function when VecView/Load and DMView/Load are invoked on a PETSCVIEWEREXODUSII viewer
- Added interface to define second order exodusii meshes
- updated dm/impls/plex/tests/ex26.c

show more ...

DM: Replace object composition with an API for auxiliary data - Now allow data to be tagged with a label+value - Do not separately require a DM and Vec, get the DM from the Vec

Merge branch 'knepley/feature-plex-ds-traversals' into 'main'

DS: Moved management of pointwise functions into PetscWeakForm object - Added...

See merge request petsc/petsc!3726

Plex+DS: Moved management of pointwise functions into PetscWeakForm object
- Added PetscChunkBuffer to manage contiguous buffer space
- WeakForm can index functions by Label/val, field
- WeakForm

Plex+DS: Moved management of pointwise functions into PetscWeakForm object
- Added PetscChunkBuffer to manage contiguous buffer space
- WeakForm can index functions by Label/val, field
- WeakForm allows multiple kernels per key

show more ...

DMPLEX: fix GLVIS regression with mfem > 4.2

DMPlex: Some fixes to the CGNS reader

PetscSFView() never called viewer for the specific type (bug), hence many output files were incorrect.

Added PetscSFView() support for basic MPI implementation and binary viewer saves communication

PetscSFView() never called viewer for the specific type (bug), hence many output files were incorrect.

Added PetscSFView() support for basic MPI implementation and binary viewer saves communication data to binary file.
This new viewer shows exactly how much data is sent from each rank to the other ranks
allowing fiting various performanace models to the data.

This uses MatView() to save the PetscSF data, but since Mat is above PetscSF in the software stack
the code will not work with multiple PETSc libraries. How should this be managed if one wishes to preserve the
traditionl PETSc algorithm stack?

Commit-type: bug-fix, feature
Funded-by: brain project
Project: brain project
spend 1h

Update sfbasic.c save the transpose of the matrix since that indicates the sends for each rank which is easier to understand

Enable PetscSFView_Basic_PatternAndSizes only with single library

show more ...

CHKERRQ() -> CHKERRMPI()

Plex+FE: Now we can turn on second derivatives in an FEM calculation - Now FE interpolation and PetscFEEvaluateFieldJets_Internal() can also do Hessians - Add jet degree to DS - Add PTAP operations f

Plex+FE: Now we can turn on second derivatives in an FEM calculation - Now FE interpolation and PetscFEEvaluateFieldJets_Internal() can also do Hessians - Add jet degree to DS - Add PTAP operations for small matrices

show more ...

Convert MPI error type to PETSc error with string message for all MPI calls

Now PETSc examples will ONLY return PETSc error codes and never MPI error codes directly so we can understand and post-pro

Convert MPI error type to PETSc error with string message for all MPI calls

Now PETSc examples will ONLY return PETSc error codes and never MPI error codes directly so we can understand and post-process their
errors better.

The test harness will now automatically retry tests that fail with MPI, this may help with Intel MPI that produces
seemingly random failures.

Commit-type: error-checking
/spend 30m

show more ...