DMPlex: added constraints and constraint matrix

DMPlex: fix recent constraint additions

DMPlex: doc correction for DMPlexSetConstraintMatrix()

DMPlex: use PETSC_COMM_SELF for local error messages

DMPlex: added constraints and constraint matrix

DMPlex: fix recent constraint additions

DMPlex: doc correction for DMPlexSetConstraintMatrix()

DMPlex: use PETSC_COMM_SELF for local error messages

DMPlex: replace ISGetSize with ISGetLocalSize where appropriate for anchors

DMPlex: fix DMPlexGlobalToLocalHook_Plex_constraints

show more ...

DMPlex: make DMPlexCreateReferenceCell

Plex: Added DMPlexComputeGeometryFEM() and PetscFECellGeom structure

Plex: Missing declarations from bad merge

DMTS: Reorganized DMTS for Plex
- Now FVM computation functions look just like FEM
- Riemann solvers are specified to DS, just like f0/f1
- Removed explicit DMPlexTS structure, now just attach to DMT

DMTS: Reorganized DMTS for Plex
- Now FVM computation functions look just like FEM
- Riemann solvers are specified to DS, just like f0/f1
- Removed explicit DMPlexTS structure, now just attach to DMTS
- Changed DMPlexTSSetRHSFunctionLocal()
- Added DMPlexTSComputeRHSFunctionFVM()
- Added minradius to DMTS

show more ...

Plex: Fix DMCreateDefaultSection_Plex() to handle FV discretizations
- I will check CLASSID until a superclass for FE/FV becomes clear

Merge branch 'knepley/feature-dm-dimension-handling'

* knepley/feature-dm-dimension-handling:
DM: Fixed DMDASetDim() in examples
Plex: Fix for isoparametric geometry
Plex: Complex fix
TS ex1

Merge branch 'knepley/feature-dm-dimension-handling'

* knepley/feature-dm-dimension-handling:
DM: Fixed DMDASetDim() in examples
Plex: Fix for isoparametric geometry
Plex: Complex fix
TS ex11: Fixed calls to DMSetCoordinateSection()
DM: FIxed DMClone() to copy dimension
Plex+DA: Now cell geometry routine can take a PetscFE in order to evaluate at quadrature points - We evaluate isoparametrically at the quadrature points - Changed DMDAComputeCellGeometry() --> DMDAComputeCellGeometryFEM() - Changed DMPlexComputeCellGeometry() --> DMPlexComputeCellGeometryFEM() - Added DMPlexComputeCellGeometryAffineFEM()
DM: Added a embedding dimension, different from the topological dimension - This can now be set in DMSetCoordinateSection(), or it can be automatically determined
Plex: Fix output for ex4
DM: Added DMGetDimPoints() - This replicates DMPlex/DAGetDepthStratum(), but uses the dimension which people have said is more intuitive
DM: Move dimension up into DM - Both DA and Plex had this and I think it makes sense

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

show more ...

Plex+DA: Now cell geometry routine can take a PetscFE in order to evaluate at quadrature points
- We evaluate isoparametrically at the quadrature points
- Changed DMDAComputeCellGeometry() --> DMDACo

Plex+DA: Now cell geometry routine can take a PetscFE in order to evaluate at quadrature points
- We evaluate isoparametrically at the quadrature points
- Changed DMDAComputeCellGeometry() --> DMDAComputeCellGeometryFEM()
- Changed DMPlexComputeCellGeometry() --> DMPlexComputeCellGeometryFEM()
- Added DMPlexComputeCellGeometryAffineFEM()

show more ...

DM: Move dimension up into DM
- Both DA and Plex had this and I think it makes sense

Merge branch 'knepley/feature-quad-fas'

* knepley/feature-quad-fas:
DMPlex: Proper preallocation for interpolators
DMPlex: Added events for interpolator and injector construction - Fixed event n

Merge branch 'knepley/feature-quad-fas'

* knepley/feature-quad-fas:
DMPlex: Proper preallocation for interpolators
DMPlex: Added events for interpolator and injector construction - Fixed event names
DM: Added sequence value for each number - Added DMOutputSequenceLoad() - Changed DMGet/SetOutputSequenceNumber() signature
Sys: Made a simple two-int hash - Someone should remove that other IJ crap
PetscDS: Renamed from PetscProblem
DMPlex: Fixed allocation
DMPlex: Support for injection from more complicated discretizations - Added DMPlexMatGetClosureIndicesRefined() - Added DMPlexComputeInjectorFEM() - Fixed check
PetscFE: Use in-cell test specific to cell shape
DMPlex: Added CellRefinerInCellTest_Internal()

show more ...

Merge branch 'master' of bitbucket.org:petsc/petsc

Merge branch 'knepley/feature-problem-class'

* knepley/feature-problem-class: (40 commits)
PetscDS: Renamed from PetscProblem
ftn-types: fix compilers errors with PetscSection [remove duplicate

Merge branch 'knepley/feature-problem-class'

* knepley/feature-problem-class: (40 commits)
PetscDS: Renamed from PetscProblem
ftn-types: fix compilers errors with PetscSection [remove duplicate code from petscvec*.h90]
PetscSection: Fix Fortran datatypes support
DMPlex: Added quads to cell refiner transform support
DMPlex ex3: Updated to Problem interface
DMDA: Fixed FE projection functions
PetscSection: Fixes for Fortran stubs
DMPlex+TS+Doc: Function docs
PetscSection: Fortran bindings
DMPlex+HDF5: Had revered handedness of transformed periodic coordinates
DMPlex+HDF5: Preserve aspect ratio with remapping of periodic coordinates
DMPlex: Added missing declarations
Sys+DMDT: Added PetscValidFunction()
DMPlex: Fix for DMPlexProjectField()
DMPlex+HDF5: Selectively output fields for viz - Fields with NULL name are ignored
PetscSection: Allow NULL names
DMPlex: Fixes for DMPlexProjectField()
PetscProblem: Added scalar integral (objective function) - Added PetscProblemGet/SetObjective() - Used in PetscFEIntegrate() - Reworked PetscProblemEnlarge_Static()
DM+DT: Adding Fortran bindings
PetscViewer: Added HDF_VIZ format to names
...

show more ...

DMPlex: Support for injection from more complicated discretizations
- Added DMPlexMatGetClosureIndicesRefined()
- Added DMPlexComputeInjectorFEM()
- Fixed check

Merge branch 'knepley/feature-problem-class' into knepley/feature-quad-fas

* knepley/feature-problem-class: (38 commits)
PetscSection: Fix Fortran datatypes support
DMPlex: Added quads to cell r

Merge branch 'knepley/feature-problem-class' into knepley/feature-quad-fas

* knepley/feature-problem-class: (38 commits)
PetscSection: Fix Fortran datatypes support
DMPlex: Added quads to cell refiner transform support
DMPlex ex3: Updated to Problem interface
DMDA: Fixed FE projection functions
PetscSection: Fixes for Fortran stubs
DMPlex+TS+Doc: Function docs
PetscSection: Fortran bindings
DMPlex+HDF5: Had revered handedness of transformed periodic coordinates
DMPlex+HDF5: Preserve aspect ratio with remapping of periodic coordinates
DMPlex: Added missing declarations
Sys+DMDT: Added PetscValidFunction()
DMPlex: Fix for DMPlexProjectField()
DMPlex+HDF5: Selectively output fields for viz - Fields with NULL name are ignored
PetscSection: Allow NULL names
DMPlex: Fixes for DMPlexProjectField()
PetscProblem: Added scalar integral (objective function) - Added PetscProblemGet/SetObjective() - Used in PetscFEIntegrate() - Reworked PetscProblemEnlarge_Static()
DM+DT: Adding Fortran bindings
PetscViewer: Added HDF_VIZ format to names
PetscFE+OpenCL: Make more respectable elasticity
PetscFE+OpenCL: Add flag for printing kernel
...

show more ...

Merge branch 'knepley/feature-problem-class' into knepley/fix-fe-gpu

* knepley/feature-problem-class: (34 commits)
PetscSection: Fixes for Fortran stubs
DMPlex+TS+Doc: Function docs
PetscSecti

Merge branch 'knepley/feature-problem-class' into knepley/fix-fe-gpu

* knepley/feature-problem-class: (34 commits)
PetscSection: Fixes for Fortran stubs
DMPlex+TS+Doc: Function docs
PetscSection: Fortran bindings
DMPlex+HDF5: Had revered handedness of transformed periodic coordinates
DMPlex+HDF5: Preserve aspect ratio with remapping of periodic coordinates
DMPlex: Added missing declarations
Sys+DMDT: Added PetscValidFunction()
DMPlex: Fix for DMPlexProjectField()
DMPlex+HDF5: Selectively output fields for viz - Fields with NULL name are ignored
PetscSection: Allow NULL names
DMPlex: Fixes for DMPlexProjectField()
PetscProblem: Added scalar integral (objective function) - Added PetscProblemGet/SetObjective() - Used in PetscFEIntegrate() - Reworked PetscProblemEnlarge_Static()
DM+DT: Adding Fortran bindings
PetscViewer: Added HDF_VIZ format to names
PetscFE+OpenCL: Make more respectable elasticity
PetscFE+OpenCL: Add flag for printing kernel
DMPlex: Fixes to DMPlexProjectField()
PetscFE: Fix work arrays
PetscFE+OpenCL: Fix interface updates
DMDT: Fix complex problems
...

show more ...

DMPlex: Added missing declarations

Merged upstream/master to feature-dmmoab

Merge remote-tracking branch 'upstream/master'

DMPlex: Added DMPlexTSComputeIFunction()
- Fixed docs for the SNES version

Merge branch 'master' into prbrune/snes-fascontinuation

* master: (855 commits)
fix for moab example had wrong blocksize of mapping
docs:remove src/docs/doxygen/manual [per Barry's instructions]

Merge branch 'master' into prbrune/snes-fascontinuation

* master: (855 commits)
fix for moab example had wrong blocksize of mapping
docs:remove src/docs/doxygen/manual [per Barry's instructions]
gfortran: remove 132 char limit from f90 free-form
SNES: Set up matrix-free *after* copying the user Jacobian to snes->pc
DMPlex-Gmsh: Use trueNumCells when building the coordinate section
Sys: Got guard logic reversed
Sys: The #define guard for PetscMalloc*() must match the one for PetscFree*()
PetscFV: I hate fucking complex
removed ltog scatter from DMDA since functionality can be obtained with gtol scatter and less memory
fix test of DMLocalToGlobalBegin/End and add actual test cases
improve object memory logging for DA, ISLocalToGlobal and ISBLOCK
lower high water mark of constructing parallel scatter by allocating only as much room in svalues as needed
owner in VecScatter construction is the rank, therefor can use PetscMPIInt to hold it do not PetscInt
another reuse of allocated array in 2d/3d DMDA create
improved high water mark of DMDASetups by implementing and using PETSC_USE_POINTER for ISCreateBlock()
PetscFV: Another complex fix
DMPlex+TS: Fixing types
ISIdentity() was wrong for block IS! Indices are NOT scaled by bs
PetscFV: Straighten out complex types
TS ex11: Fix casts and remove unused variable
...

show more ...

DMPlex: Assembly routines now use embedded PetscProblem rather than taking an array of PetscFE
- Changed were verified using SNES ex12 and ex62

Merge branch 'master' into barry/saws-options

Merge branch 'knepley/feature-fv-class'

* knepley/feature-fv-class: (44 commits)
PetscFV: Another complex fix
DMPlex+TS: Fixing types
PetscFV: Straighten out complex types
TS ex11: Fix casts

Merge branch 'knepley/feature-fv-class'

* knepley/feature-fv-class: (44 commits)
PetscFV: Another complex fix
DMPlex+TS: Fixing types
PetscFV: Straighten out complex types
TS ex11: Fix casts and remove unused variable
PetscFV: Fix Fortran define
PetscFV: Complex fixes
SNES+DM: Added SNESMonitorFields() - Activated by -snes_monitor_fields
PetscFV: Add Fortran types
SNES ex62: Use new DMPlexAddBoundary() interface - Added declaration of DMPlexComputeL2FieldDiff() - Fixed cast
SNES ex12: Fixed Neumann boundary conditions, and moved to new DMPlexAddBoundary interface() - Now do not automatically mark boundary faces
PetscFV: Fixed Fortran bindings
PetscFV+Plex: Refactored to put reconstruction routines in PetscFV_LeastSquares - Added PetscFVComputeGradient() - Added PetscFVLeastSquaresSetMaxFaces() - Moved last squares fitting routines to PetscFV_LeastSquares, and all the work space
PetscLimiter: Add other implementations
TS ex11: Fixed parallel mesh creation - Distribute after physics - Add output for tests
DMPlex: Fixed bug when using refinement from options - Must copy over the point SF
DMPlex: Distribution should copy over BC info
TS ex11: Added refinement tests
TS ex11: Fix boundary conditions after refinement, set types from command line
DMPlex: Added DMPlexCopyBoundary() - Copy boundary conditions after DMRefine() and DMPlexConstructGhostCells()
TS ex11: Specify boundary label
...

Conflicts:
config/builder.py
include/petscdmplex.h
src/dm/impls/plex/plex.c
src/snes/examples/tutorials/ex62.c

show more ...

Merge branch 'master' into stefano_zampini/pcbddc-multilevel-improvement

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bdd

Merge branch 'master' into stefano_zampini/pcbddc-multilevel-improvement

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bddcprivate.c

show more ...

Merge branch 'knepley/fix-pylith-parallel'

* knepley/fix-pylith-parallel: (22 commits)
DMPlex: In refinement, merge hybrid and non-hybrid SF creation (should be done for everything) - If missing,

Merge branch 'knepley/fix-pylith-parallel'

* knepley/fix-pylith-parallel: (22 commits)
DMPlex: In refinement, merge hybrid and non-hybrid SF creation (should be done for everything) - If missing, define *Max as *End - Fixed bound on hybrid 2D simplicial hybrid cell refinement - Fixes bug in which wrong remote point number was calculated if local process had no hybrid cells, but remote process did so refiners were mismatched
DMPlex: Fix for hybrid refinement - Must correct offsets when neighbor process has no hybrid points
DMPlex: Fix hex hybrid refinement SF - Fixed offset - Fixed number of cell faces - Fixed missing edge offset - Fixed number of cells faces - FIxed missing edge offset
PetscSection: Must initialize remoteOffsets
DMPlex: Fixed cohesive cell creation for clamped boundaries - Explicitly tag cross-edges (edges with two clamped endpoints) with 2 in boundary label - Add cross-edge to vertices when looking for things that need split points replaced - Remove extra check for boundary vertices
DMPlex: Fixes for marking cohesive submesh - Do not mark points as split which are in the boundary label - Mark boundary points from label as unsplit
DMPlex: Changing meaning of boundary label in DMPlexLabelCohesiveComplete() - Now it marks vertices already present in the original label - We mark these as unsplit, and ignore their stars - Remove old code
DMPlex: Fixed stupid typo
DMPlex: Fixed refinement to respect embedded manifolds
DMPlex: Fixed types in DMPlexOrient()
DMPlex: Fixed initialization in DMPlexLabelCohesiveComplete()
DMPlex: Fixed logic in hybrid cell tagging and added missing declaration
DMPlex Fixed huge error with last change set - Default cellHeight is 0
DMPlex: Allow cellHeight argument for DMPlexMarkBoundaryFaces()
DMPlex: Added blabel argument to DMPlexLabelCohesiveComplete() - Lets user designate the boundary of a dividing surface - Added GetSurfaceSize_Static()
DMPlex: Fix for DMPlexOrient() - I was not symmetrizing the orientation graph on processes
DMPlex: Added DMPlexAddConeSize()
DMPlex: Fixed memory overwrite at depth 0 for viewing
PC+Doc: Fixed FieldSplit doc
Sys: Replace Malloc1() in multiple mallocs (fix size)
...

show more ...