DM: Replace 0 -> NULL for pointers

DTDS: Fixed bug in SelectEquations() and added SelectDiscretizations()

Plex: Need to add BC for the solution time derivative in TS problems
- Added DMPlexInsertTimeDerivativeBoundaryValues()
- Added time derivative of the solution argument to DMAddBoundary(),
DMGetBou

Plex: Need to add BC for the solution time derivative in TS problems
- Added DMPlexInsertTimeDerivativeBoundaryValues()
- Added time derivative of the solution argument to DMAddBoundary(),
DMGetBoundary(), PetscDSAddBoundary(), PetscDSUpdateBoundary(),
PetscDSGetBoundary()

show more ...

TS+Plex: Convergence estimation now working of time dependent PDEs
- Add exact solution time derivative to DS
- Put monitor in ConvEst for TS
- Fixed error determination in ConvEst for TS

Co-Authore

TS+Plex: Convergence estimation now working of time dependent PDEs
- Add exact solution time derivative to DS
- Put monitor in ConvEst for TS
- Fixed error determination in ConvEst for TS

Co-Authored-By: Brad Aagaard <baagaard@usgs.gov>

show more ...

Merge remote-tracking branch 'origin/maint'

return derivative tabulation, not function tabulation twice

typo

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

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 remote-tracking branch 'origin/maint'

Single petscdir.mk

Fix bad access in PetscDTPKDEvalJet()

It has been passing tests because as long as the misaccessed bits are a PetscReal,
the added terms were zero. This slipped through the CI tests because VALGRIN

Fix bad access in PetscDTPKDEvalJet()

It has been passing tests because as long as the misaccessed bits are a PetscReal,
the added terms were zero. This slipped through the CI tests because VALGRIND=1
isn't running valgrind on all tests in a test set

show more ...

Merge remote-tracking branch 'origin/maint'

gcc/gfortran-10: fix errors/warnings

The primary issue now is - gfortran defaults to flagging argument mismatch as error.
We can work-around this with '-fallow-argument-mismatch' option - however th

gcc/gfortran-10: fix errors/warnings

The primary issue now is - gfortran defaults to flagging argument mismatch as error.
We can work-around this with '-fallow-argument-mismatch' option - however this switches from errors to warnings.

And the only way to disable this warning is '-w' i.e disable all warnings.

Fixes involve adding interface definitions for all prototypes that get used. [and work-around mpich by using a(1) for scalars]

show more ...

PetscDefined: remove many uses of #ifdef PETSC_USE_DEBUG

Merge branch 'knepley/feature-dm-addboundary-collective' into 'master'

Knepley/feature dm addboundary collective

See merge request petsc/petsc!2659

DM: Label completion moved from DMCreateLocalSection() to DMAddBoundary() - Labels must be complete to support BC for higher order spaces - DMAddBoundary() is now collective - PetscDSAddBoundary() is

DM: Label completion moved from DMCreateLocalSection() to DMAddBoundary() - Labels must be complete to support BC for higher order spaces - DMAddBoundary() is now collective - PetscDSAddBoundary() is also collective
- Labels are also completed in DMCopyDS()

show more ...