Plex: Fix projection for different dimension in input and output - Set height for all three meshes explicitly - Get different DS for input and output DMs - Also check closure over input section - Sho

Plex: Fix projection for different dimension in input and output - Set height for all three meshes explicitly - Get different DS for input and output DMs - Also check closure over input section - Should pass NULL for localU in function projection

show more ...

Merge branch 'knepley/feature-dm-aux-data-keys' into 'main'

Plex: Updates to declarative interface for FEM problems

See merge request petsc/petsc!3972

DM+Plex: Now auxliary data is retrieved by key

Merge branch 'knepley/feature-dm-aux-data' into 'main'

DM: Replace object composition with an API for auxiliary data

See merge request petsc/petsc!3798

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 'jolivet/fix-formatting-bracket' into 'master'

checkbadSource: new rule for wrongfully placed "{"

See merge request petsc/petsc!3655

checkbadSource: new rule for wrongfully placed "{"

Merge branch 'knepley/fix-plex-bd-normal' into 'master'

Plex: Fixed bugs in boundary integration - DMFieldComputeFaceData() now...

See merge request petsc/petsc!3255

Plex: Fixed bugs in boundary integration - DMFieldComputeFaceData() now calculates the supporting cells and orients the normal - Removed code that calculates supporting cells from plexfem.c - PetscFE

Plex: Fixed bugs in boundary integration - DMFieldComputeFaceData() now calculates the supporting cells and orients the normal - Removed code that calculates supporting cells from plexfem.c - PetscFEGeom normals our outward for the first supporting cell - In projection, it looks like I was using h instead of effectiveHeight

show more ...

Merge branch 'dalcinl/fix-null' into 'master'

Fix 0 -> NULL for pointers

See merge request petsc/petsc!3086

DM: Replace 0 -> NULL for pointers

Merge branch 'jed/plex-tensor-closure-depth' into 'master'

Jed/plex tensor closure depth

See merge request petsc/petsc!2802

Plex: During projection, must do field number translation when checking for active fields in a DS

Merge branch 'knepley/feature-plex-project-point' into 'master'

Plex: Have projection report the point it is currently working on

See merge request petsc/petsc!2865

Plex: Have projection report the point it is currently working on

Merge branch 'knepley/feature-plex-check-multiple-ds' into 'master'

Knepley/feature plex check multiple ds

See merge request petsc/petsc!2806

Plex+SNES: Large update to support multiple DSes
- Now allow simple divisions of the domain with DMSetField()
- DMPlexComputeResidual_Internal() supports multiple DSes
- DMPlexComputeJacobian_Interna

Plex+SNES: Large update to support multiple DSes
- Now allow simple divisions of the domain with DMSetField()
- DMPlexComputeResidual_Internal() supports multiple DSes
- DMPlexComputeJacobian_Internal() supports multiple DSes
- DMSNESCheckDisretization() supports multiple DSes
- DMSNESCheckResidual() supports multiple DSes
- DMSNESCheckJacoiban() supports multiple DSes
- DMComputeL2FieldDiff() supports multiple DSes
- PetscConveEst supports multiple DSes
- Fixed bug refering to cStart from a cellIS
- Added PetscDSGetQuadrature()
- Added DMSetField_Internal()
- Added DMSetRegionNumDS()
- Protect PetscSectionView() from empty bcIndices
- ISDifference now allows a NULL second argument

show more ...

Merge branch 'knepley/fix-plex-hybrid-assembly' into 'master'

Knepley/fix plex hybrid assembly

See merge request petsc/petsc!2793

Plex: More hybrid fixes
- Converted hybrid assembly to use PetscTabulation
- Protect for NULL discretizations, which happen with no dofs for a field
- Added PetscDiscType enum
- Changed PetscDSIsFE_I

Plex: More hybrid fixes
- Converted hybrid assembly to use PetscTabulation
- Protect for NULL discretizations, which happen with no dofs for a field
- Added PetscDiscType enum
- Changed PetscDSIsFE_Internal() to PetscDSGetDiscType_Internal()
- Must inject coordinates from hybrid cells into the ambient space

show more ...

Plex+FE: Updates for hybrid integration
- Added support for hybrid Jacobians
- Added PetscDSHasBdJacobian()
- Added boundary Jacobian preconditioner
- Needed a path to integrate the hybrid Jacobian p

Plex+FE: Updates for hybrid integration
- Added support for hybrid Jacobians
- Added PetscDSHasBdJacobian()
- Added boundary Jacobian preconditioner
- Needed a path to integrate the hybrid Jacobian preconditioner
- Fixed indexing in hybrid element mat update
- Be careful about the difference between intrinsic dimension and embedding dimension for hybrid integration
- Need to use rectangular transformation of gradient in case they are different
- Fix use of embedding dimension
- Make sure FEGeom struct has the dimensions set
- Pointwise function derivatives are in the embedding space
- Use PetscArrayzero()
- Change TransitiveClosure() for hybrid cells
- We get the closures of the endcaps, and then the hybrid faces, this allows degenerate cells that we need in PyLith

show more ...

Merge branch 'knepley/fix-plex-projection-label' into 'master'

Plex: This fixes a pernicious bug with projection - When setting points along...

See merge request petsc/petsc!2773

Plex: This fixes a pernicious bug with projection - When setting points along a label, we must filter out points in the closure that are not also in the label.

Merge branch 'knepley/feature-dm-boundary-function-projection' into 'master'

DM: Added DMProjectBdFieldLabelLocal()

See merge request petsc/petsc!2667

DM: Added DMProjectBdFieldLabelLocal()
- Called from DMPlexInsertBoundaryValuesEssentialBdFIeld()
- Added projectbdfieldlabelLocal() to DM
- Added Plex impl
- Added more documentation of the InsertBo

DM: Added DMProjectBdFieldLabelLocal()
- Called from DMPlexInsertBoundaryValuesEssentialBdFIeld()
- Added projectbdfieldlabelLocal() to DM
- Added Plex impl
- Added more documentation of the InsertBoundaryValues*() methods

show more ...

Merge branch 'tisaac/feature-dual-space-rework' into 'master'

Rework PetscDualSpace

See merge request petsc/petsc!2399