Undo change of pushfoward/pullback interface that wasn't backed up by implementation yet

Pushforwards/pullbacks of form degrees in [1,d) have different sizes depending on the embedding dimension
(a

Undo change of pushfoward/pullback interface that wasn't backed up by implementation yet

Pushforwards/pullbacks of form degrees in [1,d) have different sizes depending on the embedding dimension
(a 1-form has d degrees of freedom in d dimensions). I changed the interface to reflect this, but didn't
have to do the implementation. I'm reverting to the previous interface but putting in a note saying
that the pushforwards and pullbacks are only valid if the embedding dimension is the same as the reference
dimension.

show more ...

Rework the PetscDualSpace interface

Mention form degree in viewer output of dual space

Set from degree from options for dual space

Merge branch 'knepley/feature-dm-remove-hybrid' into 'master'

Knepley/feature dm remove hybrid

See merge request petsc/petsc!2614

DM+Plex: Removed the idea of hybrid and ghost cells, in favor of cell type
- Improved Plex tests to work for other cell types
- Plex creation, interpolation, submesh, and refinement routines now set

DM+Plex: Removed the idea of hybrid and ghost cells, in favor of cell type
- Improved Plex tests to work for other cell types
- Plex creation, interpolation, submesh, and refinement routines now set the cell type
- Distribution sorts strata by cell type
- Fixed orientation of prisms from input formats
- Added Plex ex40 to test regular refinement
- CellRefiner is a PetscObject
- dm_plex_cell_refiner chooses the refinement strategy
- Removed SimplexToTensor stuff in favor of ToHex refiner
- Most functions take DMPolytopeType instead of CellRefiner
- CellRefinerInCellTest_Internal()
- CellRefinerGetAffineTransforms_Internal()
- CellRefinerGetAffineFaceTransforms_Internal()
- DMForest now responds to DMPlexGetOverlap() and propagates its overlap to the underlying Plex
- Use DMPlexComputeCellType_Internal() to classify a cell based only on topology
- Added DMPlexGetSimplexOrBoxCells() to replace DMPlexGetInteriorCellStratum()
- Added DM_POLYTOPE_POINT_PRISM_TENSOR
- Added DMPolytopTypeGetDim() and DMPolyTopeTypeGetConeSize()
- Added DMLabelSetStratumBounds()
- Added DMPlexSetCellType()
- Added celltypes to DMView()
- Added DMPolytopeTypeGetNumVertices()
- Added DMPlexCreateReferenceCellByType()
- Added DMPolytopeTypeFromGmsh()
- Made DMPlexRefineUniform() public

show more ...

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

Knepley/feature project different spaces

See merge request petsc/petsc!2153

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 ...

Plex: Now for projection allow 3 DMs: input, output, and auxiliary
- We get the input DM from the input Vec
- Now we can have more general relations between input/output/auxiliary meshes in projectio

Plex: Now for projection allow 3 DMs: input, output, and auxiliary
- We get the input DM from the input Vec
- Now we can have more general relations between input/output/auxiliary meshes in projection and integration
- Added DMEnclosureType
- Added DMGetEnclosureRelation() and DMGetEnclosurePoint()
- Removed DMPlexGetSubpoint() and DMPlexGetAuxiliaryPoint()
- In projection, convert output mesh points to input mesh points
- Coordinates must be reordered after depth label is constructed
- DMPlexLabelAddFaceCells() only adds cells adjacent to faces
- Put in extra debugging for missing point enclosures
- In boundary projection, calculate face data for geometry
- p4est: use known point maps instead of projection to interp from base to coarse

show more ...

Merge branch 'knepley/fix-fv-ghost-cells' into 'master'

Knepley/fix fv ghost cells

See merge request petsc/petsc!2132

Plex: Use new API to handle hybrid and FV ghost cells

Merge branch 'jed/deprecate-dmgetdefaultsection' [PR #1969]

* jed/deprecate-dmgetdefaultsection:
DMGetDefaultSection et al: add PETSC_DEPRECATED_FUNCTION
These functions were deprecated in v3.9, but

Merge branch 'jed/deprecate-dmgetdefaultsection' [PR #1969]

* jed/deprecate-dmgetdefaultsection:
DMGetDefaultSection et al: add PETSC_DEPRECATED_FUNCTION
These functions were deprecated in v3.9, but no PETSC_DEPRECATED_FUNCTION was added at that time so many uses crept back into the source code. Here we mark it deprecated and fix all existing uses in PETSc.

DMGetLocalSection/DMSetLocalSection: add explicit names for DMGetSection/DMSetSection
These names are for symmetry with DMGetLocalVector() versus DMGetGlobalVector().

DMGetSection: replace all internal usage with DMGetLocalSection
Same for DMSetSection -> DMSetLocalSection. We retain the interface for DMGetSection/DMSetSection without deprecation warnings (for now).

show more ...

DMGetSection: replace all internal usage with DMGetLocalSection

Same for DMSetSection -> DMSetLocalSection. We retain the interface
for DMGetSection/DMSetSection without deprecation warnings (for n

DMGetSection: replace all internal usage with DMGetLocalSection

Same for DMSetSection -> DMSetLocalSection. We retain the interface
for DMGetSection/DMSetSection without deprecation warnings (for now).

Suggested-by: Barry Smith <bsmith@mcs.anl.gov>

show more ...

Merge remote-tracking branch 'origin/barry/add-petscarrayxxx-ops' [PR #1771]

Add PetscArraycmp(), PetscArraycpy(), PetscArrayzero(), PetscArraymove(), update code to use new forms

Replace a couple multi-mallocs with a single malloc
Replace backwards loops in MatSetValues_*()

Add PetscArraycmp(), PetscArraycpy(), PetscArrayzero(), PetscArraymove(), update code to use new forms

Replace a couple multi-mallocs with a single malloc
Replace backwards loops in MatSetValues_*() with PetscArraymove()

Commit-type: style-fix, feature

show more ...

Merged in knepley/feature-dualspace-bdm (pull request #1699)

Knepley/feature dualspace bdm

fix -Wmaybe-uninitialized warnings with arch-c-exodus-dbg-builder

Plex+FE: Fixed improper reference for coordinates

PetscFE+Plex: Use cell and face geometry in boundary integration

Merged in knepley/feature-plex-pullback (pull request #1632)

Plex: Implemented pullback/pushforward

Approved-by: BarryFSmith <bsmith@mcs.anl.gov>

Plex: Left in a C++ comment

Plex+FE: Moved internal functions from header to fe.c
Suggested-by: Jed Brown <jed@jedbrown.org>

Plex: Implemented pullback/pushforward - Added k to PetscDualSpace. It indicates what kind of transformation we use. - DualSpace does pullback and FE integration routines do pushforward - All integra

Plex: Implemented pullback/pushforward - Added k to PetscDualSpace. It indicates what kind of transformation we use. - DualSpace does pullback and FE integration routines do pushforward - All integration now takes a PetscFEGeom argument and pushes forward correctly - Added PetscDualSpaceTransformType - Removed TransformF() because we now transform the basis function instead - Chaged workspace management in DS - Updated L2 diff functions - Updated interpolation functions - Added way to output exact solution in DMSNESCheck()

show more ...

Merged in knepley/feature-dm-global-basis (pull request #1582)

Knepley/feature dm global basis

Approved-by: BarryFSmith <bsmith@mcs.anl.gov>

DM: Correct bugs in last commit
- I am changing editors and messed up

DM: Add ability to specify different global and local bases
- The global basis is for solving and applying BC
- The local basis is for assembly
- We transform bases in DMLocalToGlobal() and DMGlobalT

DM: Add ability to specify different global and local bases
- The global basis is for solving and applying BC
- The local basis is for assembly
- We transform bases in DMLocalToGlobal() and DMGlobalToLocal()
- Added a default implementation which does a rotation
- Transform functions used in BC and exact solutions, which we assume operate in the global basis
- Transform functions used in projection, which we assume operate in the global basis
- Transform Jacobian values to the global basis before insertion
- Transform coordinates to global basis when calling out projection, BC, and exact solution functions
- Copy transforms when forming a hierarchy
- Added many tiny LA things (Eigen is no longer a good dependency)

show more ...

Merged in knepley/feature-bdfunc-normal (pull request #1500)

Plex: Feed the facet normal to a type of projection function

Approved-by: BarryFSmith <bsmith@mcs.anl.gov>