PetscFE: Added discontinuous Lagrange dual spaces
- Updated PetscFECreateDefault() to accept a continuity argument
- Added PetscDualSpaceLagrangeGet/SetContinuity()

PetscFE: Added cell dofs for Lagrange

PetscFE: Now PetscDualSpaceGetDimension_SingleCell_Lagrange() takes an order argument

PetscFE: Fix array initialization bug

Merge branch 'knepley/feature-plex-neumann-bc'

* knepley/feature-plex-neumann-bc:
SNES ex12: Added support for boundary Jacobian
DMPlex: Added boundary Jacobian to DMPlexComputeJacobianFEM()
D

Merge branch 'knepley/feature-plex-neumann-bc'

* knepley/feature-plex-neumann-bc:
SNES ex12: Added support for boundary Jacobian
DMPlex: Added boundary Jacobian to DMPlexComputeJacobianFEM()
DMPlex: Removed defunct interpolator function
PetscFE: Added PetscFEIntegrateBdJacobian() - Added pointwise functions to PetscFEM - Added an implementation for Basic
DMPlex: Correct boundary iteration for Neumann conditions in residual assembly

show more ...

Merge branch 'knepley/feature-fe-derivatives'

* knepley/feature-fe-derivatives:
Doc+Plex: Small fix for distribution
PetscFE: Correct type for reference space derivatives
SNES ex62: Updated to

Merge branch 'knepley/feature-fe-derivatives'

* knepley/feature-fe-derivatives:
Doc+Plex: Small fix for distribution
PetscFE: Correct type for reference space derivatives
SNES ex62: Updated to new derivative mapping
SNES ex12: Updated to new derivative mapping
PetscFE: Changed mapping of derivatives - Now they are transformed only pointwise, instead of transforming the whole tabulation
SNES ex62: Works with new derivative mapping
SNES ex12: Works with new derivative mapping
PetscFE: Changed PetscFEIntegrateResidual_Basic() derivative mapping - Instead of mapping the basis derivatives, we only map the pointwise derivatives gradU and f_1 Suggested-by: Jed Brown <jedbrown@mcs.anl.gov>
PetscFE: Removed PetscGenerateFEMQuadrature*.py

Conflicts:
config/builder.py

show more ...

Merge branch 'knepley/feature-plex-periodic'

* knepley/feature-plex-periodic:
TS: Allow viewing of TS solutions
DMPlex: Fix comparison for complex
Draw: Fixed bad check
SNES ex62: Enhance ex

Merge branch 'knepley/feature-plex-periodic'

* knepley/feature-plex-periodic:
TS: Allow viewing of TS solutions
DMPlex: Fix comparison for complex
Draw: Fixed bad check
SNES ex62: Enhance ex62 to use tensor-product cells
DMPlex: Now we are correctly creating the Mobius band - Coordinates for uniform refinement are still broken
DMPlex: If domain is periodic, uniform refinement of coordinates will respect this
DMPlex: Allow geometry routines to recognize periodic quads
DMPlex: Added DMPlexLocalizeCoordinates() - This pulls coordinates from vertices to the cell, which allows discontinuities - We use this for periodic problems
DM: Added preliminary support for periodic coordinates - Added DMGet/SetPeriodicity() which is for clipping coordinates to a periodic box
DMPlex: Needed a custom Fortran wrapper for DMPlexInterpolate() because now we require that the incoming mesh be NULL
DMPlex: Fix handling of Fortran NULLs
DMPlex: Updated test output - The DMView was changed
SNES ex12: Cleanup structure - Use PetscFECreateDefault()
PetscFE: Added a view for Basic
SNES ex12: Starting to cleanup test names
SNES ex62: Marked MF tests as broken

show more ...

PetscFE: Correct type for reference space derivatives

PetscFE: Added PetscFEIntegrateBdJacobian()
- Added pointwise functions to PetscFEM
- Added an implementation for Basic

PetscFE: Changed mapping of derivatives
- Now they are transformed only pointwise, instead of transforming the whole tabulation

Merge branch 'barry/add-names-for-all-ascii-viewers'

PetscFE: Changed PetscFEIntegrateResidual_Basic() derivative mapping
- Instead of mapping the basis derivatives, we only map the pointwise derivatives gradU and f_1
Suggested-by: Jed Brown <jedbrown@

PetscFE: Changed PetscFEIntegrateResidual_Basic() derivative mapping
- Instead of mapping the basis derivatives, we only map the pointwise derivatives gradU and f_1
Suggested-by: Jed Brown <jedbrown@mcs.anl.gov>

show more ...

PetscFE: Added a view for Basic

Merge branch 'prbrune/removeunwrappedmathfunctions'

* prbrune/removeunwrappedmathfunctions:
Math: fix typo
Added src/sys/examples/tests/ex25.c for PetscScalar and PetscReal wrapping.
Wrapped t

Merge branch 'prbrune/removeunwrappedmathfunctions'

* prbrune/removeunwrappedmathfunctions:
Math: fix typo
Added src/sys/examples/tests/ex25.c for PetscScalar and PetscReal wrapping.
Wrapped tan(),[sin,cos,tan]h() for Real/Scalar, a[sin,cos,tan]() for Real.
Wrapped math.h and quadmath.h functions used in petsc source exempting examples

Conflicts:
include/petscmath.h
src/dm/dt/interface/dt.c
src/ksp/ksp/impls/cg/stcg/stcg.c
src/snes/linesearch/interface/linesearch.c
src/ts/impls/explicit/rk/rk.c

show more ...

Merge branch 'knepley/feature-plex-hdf5'

* knepley/feature-plex-hdf5: (63 commits)
DMPlex ex2: Added an example that checks a mesh
DMPlex: #ifdef support functions for HDF5
TS: Fix output sequ

Merge branch 'knepley/feature-plex-hdf5'

* knepley/feature-plex-hdf5: (63 commits)
DMPlex ex2: Added an example that checks a mesh
DMPlex: #ifdef support functions for HDF5
TS: Fix output sequence
DMPlex: Check for HDF5 attributes before writing them
DMPlex: Label with timestep number since we do not have times
Viewer+HDF5: Added PetscViewerHDF5HasAttribute()
TS: Set output sequence for timesteps
DMPlex: Handle output sequences - Added static to some functions
DM: Add idea of an output sequence - Added DMGet/SetOutputSequenceNumber() - Added doc for outputDM
REVERT TS: Timestep gets set into DMTS
DMPlex: Fix up HDF5 segregation
DMPlex: Replace MAXPATHLEN with PETSC_MAX_PATH_LEN
DMPlex: Fix HDF5 timestep handling
DMPlex: Set correct block size on HDF5 field vectors
TS: Timestep gets set into DMTS
DMTS: Added timestep
DMPlex: Fix output of fields - Corrected name - Fixed vector type
DMPlex: Return errors from HDF5 function
DMPlex+Viewer: Now output separate fields for viz - Added the PETSC_VIEWER_HDF5_VIZ format for creating viz fields - Added field type determination
Sys+HDF5: Added PetscDataType string and HDF5 support for it
...

Conflicts:
config/BuildSystem/config/packages/__init__.py
src/ts/interface/ts.c

show more ...

Insure that all ASCII viewers for PETSc objects call PetscObjectPrintClassNamePrefixType()

DMPlex: Now DMPlexInterpolate() allows the new mesh to be inserted into an existing DM
- If the input is non-NULL, no new DM is created
- This allows it to functions within the DMLoad() paradigm

DMPlex+PetscFE: Adding support for TS
- Added DMPlexCompiteIFunctionFEM()
- Added PetscFEIntegrateIFunction() and Bd version
- Added IFunction to FEM struct

Merge branch 'knepley/fix-fem-multifield'

* knepley/fix-fem-multifield:
DMDT: Fixed Fortran PetscQuadrature definition
DMPlex ex3: Added Q2 tests - Made DA tests only 2x2 - Changed output format

Merge branch 'knepley/fix-fem-multifield'

* knepley/fix-fem-multifield:
DMDT: Fixed Fortran PetscQuadrature definition
DMPlex ex3: Added Q2 tests - Made DA tests only 2x2 - Changed output format - Quad tests now actually work
DMDA: Added DMDAComputeL2GradientDiff()
PetscSpace: Fixed polynomial derivative evaluation for quads
DMDA: Fix closure operations in 2D
DMPlex: Interpolation cleanup - Remove unused PetscSpace - Fix uninitialized variable
PetscFE: Fix PetscBool handling
SNES ex62: Updated output
DMPlex ex3: Use PetscFECreateDefault() and correct tensor tests
PetscFE: Fix PetscFECreateDefault() - The space should default to tensor on a non-simplex - The quadrature should be tensor Gauss on a non-simplex
DT: Added PetscDTGaussTensorQuadrature()
DT: Added citation for Gloub+Welsch method
PetscSection: Now global section can also be viewed
DMPlex: Allow NULL function for ProjectFunction()
PetscDualSpace: Added doc for PetscDualSpaceApply() - Cherry-picked from df2d8bb
PetscFE: Fixed PetscFECreateDefault() - Setup FE - Minimum quad order - Destroy quad after setting
PetscFE: PetscFECreateDefault() now takes a number of components and a reference cell type - Cherry-picked from a255314
PetscFE: Added PetscFECreateDefault() - Cherry-picked from 8919f5a
PetscFE: Allow lowest order cell Lagrange dual basis functions
DMPlex: Fix preallocation with cell variables - I had an incorrect loop optimization

Conflicts:
src/dm/impls/plex/plexpreallocate.c

show more ...

PetscSpace: Fixed polynomial derivative evaluation for quads

Merge branch 'maint'

* maint:
DTGaussQuadrature: fix weight averaging for npoints=2 and output

Conflicts:
src/dm/dt/examples/tests/ex1.c

DTGaussQuadrature: fix weight averaging for npoints=2 and output

PetscFE: Fix PetscBool handling

PetscFE: Fix PetscFECreateDefault()
- The space should default to tensor on a non-simplex
- The quadrature should be tensor Gauss on a non-simplex

DT: Added PetscDTGaussTensorQuadrature()