Merge branch 'knepley/feature-plex-high-order-view' into 'main'

Plex: Enable high order viewing

See merge request petsc/petsc!6959

Plex: Enable high order viewing - Refinement recreates affine coordinates - Default for draw viewer - Replace DMProjectCoordinates() with DMSetCoordinateDisc() - Add argument to DMPlexCreateCoordinat

Plex: Enable high order viewing - Refinement recreates affine coordinates - Default for draw viewer - Replace DMProjectCoordinates() with DMSetCoordinateDisc() - Add argument to DMPlexCreateCoordinateSpace()

show more ...

Rename rules.doc and rules.utils because GitLab treats the former as a MS Word document.

Thanks-to: Jed Brown

Plex: Updates for FVM
- Add DMPlexInsertBoundaryValuesFVM()
- Add -dm_plex_print_fvm
- Fix copy of diagnostic flags

DM: Fix destroy

Remove first and last empty lines

Plex: More work adding in CEED backend
- Add DMPlexGetUseCeed() and DMPlexSetUseCeed()
- Need to propagate useCeed flag
- When using CEED, make tensor order the default, but shut it off for matrix in

Plex: More work adding in CEED backend
- Add DMPlexGetUseCeed() and DMPlexSetUseCeed()
- Need to propagate useCeed flag
- When using CEED, make tensor order the default, but shut it off for matrix insertion since this is still done by Plex kernels
- DMPlexMatSetClosureGeneral() now takes flag for closure permutation
- Add DMPlexMatSetClosure_Internal() which takes flag for closure permutation
- Add DMPlexVecGetClosure_Internal() with flag to turn off closure permutation - Turn off closure permutation in L2 diff calculation
- Add DMUseTensorOrder()
- Add DMPlexSNESComputeResidualCEED() - Cleanup SNES ex13 - Do not give CEED negative offsets

show more ...

DM: Helper functions for CEED vectors and error macro

Merge branch 'knepley/fix-swarm-layout' into 'main'

Swarm: Reworked Swarm projection interface and fixed some particle layout

See merge request petsc/petsc!6856

DM: Added DMPrintCellIndices() for debugging

PetscFE: add PETSCFEVECTOR

Reductions in the cost of computing, storing, and applying
a Tabulation can be achieved for a finite element space
that is a vector of "scalar" finite elements (the "scala

PetscFE: add PETSCFEVECTOR

Reductions in the cost of computing, storing, and applying
a Tabulation can be achieved for a finite element space
that is a vector of "scalar" finite elements (the "scalar"
element itself could be vector-valued).

Reductions in the cost of storing and applying
a Tabulation can't be realized until the PetscTabulation structure
itself is changed, but at the very least reduction in the cost of
computing a tabulation is possible with a new PetscFEType,
PETSCFEVECTOR.

To complete the interface, we also add PETSCDUALSPACESUM,
a dual space whose dual basis is the sum of other dual bases,
that parallels PETSCSPACESUM (the space and dual space of
a PETSCFEVECTOR are a PETSCSPACESUM and PETSCDUALSPACESUM,
respectively).

To allow the user to specify whether the
basis vector and value components of the PETSCFEVECTOR
of blockwise or interleaved, we also have to make sure
PETSCSPACESUM and PETSCDUALSPACESUM can have blockwise
or interleaved behavior. We add both capabilities.

PETSCFEVECTOR is now used by detecting compatible
vector-valued dual and primal spaces in PetscFECreateFromSpaces().

show more ...