DMIsBoundary: don't get link until list is populated

Merge branch 'tisaac/dm-boundary-in-ds'

Boundary conditions are better conceptually matched to the PetscDS,
which holds the equations, than to the DM. This makes it easier to
share the same problem

Merge branch 'tisaac/dm-boundary-in-ds'

Boundary conditions are better conceptually matched to the PetscDS,
which holds the equations, than to the DM. This makes it easier to
share the same problem between different meshes, such as in a hierarchy.

* tisaac/dm-boundary-in-ds:
PetscDS: fix examples affected by moving boundaries
builder.py: skip empty targets
DMPlex: don't copy coordinates to self
SNES ex69: fix recent PetscDS changes
DMSetDS(): reference, then destroy
PetscDS: fix _Basic creation and destruction
PETSC_DEPRECATED in source files doesn't play well with fortran interfaces
DM: switch to using PetscDS boundary in examples
PetscDS: add missing PetscDSCopyBoundary, fix copy loops
DM: set PetscDS from options in DMSetFromOptions()
DM: protect against changing prob when populating boundary
DM, PetscDS: remove redundant boundary representations
DM: transfer boundary methods to PetscDS
PetscDS: start creating parallel for DMBoundary

show more ...

DM+Doc: Added some for DMAddBoundary()

DMLabel: added DMLabelSetStratumIS()

replace uses of PetscMalloc() with PetscNew() or PetscMalloc1() where appropriate

DM: If no coordinates exist on vertices, set coordinate dimension to PETSC_DEFAULT
- If you leave it -1, then you get weird problems with inconsistent dims on different procs

DMSetDS(): reference, then destroy

PETSC_DEPRECATED in source files doesn't play well with fortran interfaces

PetscDS: add missing PetscDSCopyBoundary, fix copy loops

DM: set PetscDS from options in DMSetFromOptions()

DM: protect against changing prob when populating boundary

DM, PetscDS: remove redundant boundary representations

The DMBoundary is now just a container for the DSBoundary

DM: transfer boundary methods to PetscDS

PetscDS: start creating parallel for DMBoundary

The way we interface with a DMBoundary, it is essentially a mesh-free description of the boundary condition,
so when trying to define the same set of

PetscDS: start creating parallel for DMBoundary

The way we interface with a DMBoundary, it is essentially a mesh-free description of the boundary condition,
so when trying to define the same set of equations on multiple meshes, it makes sense for the boundary conditions
to be a part of the same object that contains the rest of the equations, the PetscDS

show more ...

Merge branch 'tisaac/ts-ex11-amr'

We now have an example (ts/examples/tutorials/ex11.c) where we can
use either DMPlex or DMP4est to handle adaptivity when time stepping.
We have tried to minimize t

Merge branch 'tisaac/ts-ex11-amr'

We now have an example (ts/examples/tutorials/ex11.c) where we can
use either DMPlex or DMP4est to handle adaptivity when time stepping.
We have tried to minimize the amount of code that is specialized to the
different DMs, and in the future more work can be done to improve this
interface.

* tisaac/ts-ex11-amr: (62 commits)
DMP4est: type fix
Plex: Fixes for analyzer
DMP4est: update snes ex12 regressions
DMPlex: fix parent/child interpolation with numFields > 0
DMP4est: also compute anyChange for uniformly refined/coarsened meshes
TS ex11: fix command line for AMR without computed gradients
TS ex11: remove gradient reconstruction from adaptivity requirements
DMPlexTree: some compilers want this var initialized
Plex: Changed Get/RestoreFaceFields() and Get/RestoreFaceGeometry() to return the number of active faces - This fixes a long-standing bad access in FV, although the values were never actually used.
DMPlexTS: removed unused variables
DMPlex/DMP4est: small fixes for next warnings and errors
SNES ex36 and ex52: These do not work anymore
DMPforestGetTransferSF_Point(): try to match parents
PetscViewerHDF5: destroy old file data if setting new filename
DMPlex + FV: remove leaky references in projection functions
Plex: Fixed incredibly stupid bug in Jacobian computation
PetscViewer HDF5: added get name method
DMPlexTree: do not skip out of child matrix loops early
PetscFV: correct field offsets when limiting gradients
DMP4estTree: ignore negative values when getting/setting vector values for transfer
...

show more ...

DM: Fixed error check in DMLocatePoints()

In order to be set up properly, MATIS needs to know the l2gmap before any preallocation

The communicator for the l2gmap provided by the DM should be the same of the DM

Merge branch 'dmay/feature-dm-swarm' into knepley/feature-plex-particle-example

* dmay/feature-dm-swarm: (69 commits)
DMSwarm: * Check for size compatabaility between vectors created for DMSwarm f

Merge branch 'dmay/feature-dm-swarm' into knepley/feature-plex-particle-example

* dmay/feature-dm-swarm: (69 commits)
DMSwarm: * Check for size compatabaility between vectors created for DMSwarm fields * Bug fix for the field register function
DMSwarm DataBucket: * Added guards to ensure that internal arrays are never re-sized if any raw pointers are being currently being accessed
DMSwarm: * Added missing part of macro
DMSwarm: * Added documentation
DMSwarm example: * Added example with DMShell and user point location
DMSwarm Migrate: * Updated to use DMLOCATEPOINT_POINT_NOT_FOUND flag
DM definition: * Added special flag for marking points as being not found within an element during a call to DMLocatePoints
DMSwarm migrate: * Fixed bug with point removal
DMSwarm: * Removed print from data exchanger
DMSwarm: * Use common field names * Re-organized init/finalize fields with DMCreate/DMSetUp * Added methods, flags, checks when using DMSwarmType PIC
SwarmMigrate: * Use common field names
PetscDMSwarm: * Added common names * Corrected extern to PETSC_EXTERN * Added names for migrate
DMSwarm implementation: * Remove unused member
DM implementation: * Added new method locatepointssubdomin
DataBucket: * Print block size
DMSwarm example: * Removed comment
DMSwarm example: * Updated function names
DMSwarm test: * Added new test for variable block sizes
DMSwarm: * Fixed typo in vector name
DMSwarm: * Finishing support for block size != 1
...

Conflicts:
include/petsc/private/dmimpl.h

show more ...

DM: Added mode to DMLocatePoints() to change failure behavior
- Implement to NEAREST mode in Plex for 2D simplices
- Added DMPlexCLosestPoint_Internal()

Merge branch 'master' into tisaac/ts-ex11-amr

* master: (141 commits)
chombo: fix dirname in tarball and update to use snapshot with gcc-6.1 compile fix
petsc4py: fix tarball url [for arbitary c

Merge branch 'master' into tisaac/ts-ex11-amr

* master: (141 commits)
chombo: fix dirname in tarball and update to use snapshot with gcc-6.1 compile fix
petsc4py: fix tarball url [for arbitary commitid] - and dirname in tarball
Plex: Moved ReconstructGradients_Internal() and ApplyLimiter_Internal() to Plex from SNES
CUDA: Added external declaration for cuBLAS handle.
docs: fix reference to KSP_Create in parent commit
Add information on PETSC_EXTERN and PETSC_INTERN to developers manual
Plex: Symbol is shared between DM and SNES, so must be EXTERN
Plex: Remove unused variables
Configure: Fix check for spaces in PETSC_DIR
Plex: Updated examples to new DMPlexCreateBoxMesh()
sick of failure due to incorrect master branch
SNES ex12: Fixed unstructured GMG test
SNES ex12: Added another unstructured GMG test
petsc4py: use the snapshot with the fix for dmplex change
Examples that require Fortran to compile but were listed in tests that did not require Fortran
Plex ex19: Added DMSetFromOptions()
Plex: Store a flag indicating whether the boundary should change during remeshing - Set using -dm_plex_remesh_bd - Changed Remesh_Internal() signature
Always coarsen boundaries.
Plex: Fixed option name
Plex: Protect unused variables
...

Conflicts:
src/dm/impls/plex/makefile
src/snes/utils/dmplexsnes.c

show more ...

DMBoundary: avoid caching the associated label

I would like to be able to add a boundary before the topology is
fixed/determined in DMForest. This means, however, that if a boundary
is declared wit

DMBoundary: avoid caching the associated label

I would like to be able to add a boundary before the topology is
fixed/determined in DMForest. This means, however, that if a boundary
is declared with a name that isn't in the existing labels, that
DMIsBoundaryPoint() can be really bad: searching for a label for every
point. A better solution needs to be implemented.

show more ...

Merge branch 'tisaac/dmp4est-set-periodic-domain' into knepley/fix-plex-mf-prec

* tisaac/dmp4est-set-periodic-domain:
Plex: fix line geometry for localized
Plex: fix section access to last commi

Merge branch 'tisaac/dmp4est-set-periodic-domain' into knepley/fix-plex-mf-prec

* tisaac/dmp4est-set-periodic-domain:
Plex: fix line geometry for localized
Plex: fix section access to last commit
DMPlex: fix rectanlge geometry for localized
DMConvert_p4est_plex: set projection height early, so it passes to coordinate dm
DMLocalizeCoordinates(): make more flexible and targeted
DMP4est: change order of setting projection height and localizing coordinates
DMP4est: prefer [0,L) for periodic coordinates
Plex: more descriptive debugging

show more ...

DMLocalizeCoordinates(): make more flexible and targeted

More flexible because methods that do calculations on facets might need
to localize them as well. This information can be passed to the
coor

DMLocalizeCoordinates(): make more flexible and targeted

More flexible because methods that do calculations on facets might need
to localize them as well. This information can be passed to the
coordinate DM with DMPlexSetMaxProjectionHeight().

More targeted by only localizing cells/points that are needed. If I'm
localizing cells _and_ facets, then that's a lot of indiscriminate
localizing, when I only need it for a lower dimensional manifold.

show more ...

Plex: Must check block size in parallel from DMGetLocalToGlobalMapping()
- Small format cleanup

Plex: Fix support for ADD_BC_VALUES