DMP4est: use DMCoarsen_Forest

Merge branch 'tisaac/sf-fix-multi-sf-leaves' into tisaac/dmp4est-feature-injection

* tisaac/sf-fix-multi-sf-leaves: (174 commits)
PetscSF: fix PetscSFGetMultiSF() for sparse leaves
PetscSF: make

Merge branch 'tisaac/sf-fix-multi-sf-leaves' into tisaac/dmp4est-feature-injection

* tisaac/sf-fix-multi-sf-leaves: (174 commits)
PetscSF: fix PetscSFGetMultiSF() for sparse leaves
PetscSF: make regression test that fails
few more clang static analyzer fixes
many issues detected by clang static analyzer
Update hypre to its release version
Updates for xSDKTrilinos builds
Added --download-xsdktrilinos
Communicator passed to VecCreateSeq should be PETSC_COMM_SELF
PCMG man page: fix typo in options name
Bib: Added refs
DMLabel: Removed unused variable
Metis tries to use backtrace which requires -lexecinfo on some systems such as freebsd Remove -lexecinfo from freebsd configure files since it is only needed by metis
add alternative output that occurs on some systems due to use of random number generator in partitioner
Plex: When creating cohesive submeshes, fix memory leak
Plex ex11: Updated output for new distribution algorithm
Plex: When explicitly distributing the depth label, we must recreate empty strata
Plex: When partitioning, the cell numbering should include hybrid points
Plex: Allow cell and vertex numberings to include hybrid points
need to show diffs of ex22 if it fails
Support two different output for test example, since due to roundoff this example can produce slightly different convergence history
...

show more ...

Merge branch 'master' into sanderarens/fix-plex-neumann-bc

* master: (28 commits)
few more clang static analyzer fixes
many issues detected by clang static analyzer
Update hypre to its release

Merge branch 'master' into sanderarens/fix-plex-neumann-bc

* master: (28 commits)
few more clang static analyzer fixes
many issues detected by clang static analyzer
Update hypre to its release version
Updates for xSDKTrilinos builds
Added --download-xsdktrilinos
DMLabel: Removed unused variable
Plex: When creating cohesive submeshes, fix memory leak
Plex ex11: Updated output for new distribution algorithm
Plex: When explicitly distributing the depth label, we must recreate empty strata
Plex: When partitioning, the cell numbering should include hybrid points
Plex: Allow cell and vertex numberings to include hybrid points
Plex ex1: Fixed output for empty strata
DMLabel: Added DMLabelHasStratum()
Minor changes and fix to bugs introduced in 82f73ecaa
Modified makefile to fix aijcusparse tests
Add missing code for the case PETSC_HAVE_VECCUDA
Fix bugs introduced in 82f73ecaa
Use VECCUDA with MATAIJCUSPARSE
CUDA: Fixed visibility and disabled Fortran bindings for VecScatters on GPU.
Replace kernels with thrust
...

show more ...

many issues detected by clang static analyzer

Time: .5 hours

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Merge branch 'master' into knepley/feature-plex-handle-func-failure

* master: (378 commits)
SNES ex5: Fix coordinate determination for empty subdomains
fixed petsc.bib entry
p4est: fix syntax

Merge branch 'master' into knepley/feature-plex-handle-func-failure

* master: (378 commits)
SNES ex5: Fix coordinate determination for empty subdomains
fixed petsc.bib entry
p4est: fix syntax error with error message
SNES ex5: Correct Jacobian handling on patches
SNES ex5: Stupid complex fixes
DMP4est: fix transfer of labels to parents
petsc4py: update to latest snapshot. Fixes da example breakage in petsc4py test suite
Plex: Correct setting of options on refined or coarsened meshes
Plex: Add info for grid hashing
Plex: Add logging in interpolator creation
Plex: Correct use of tabs in ASCII viewer
DM: Shorten long event name
Plex ex12: Create boundary label automatically at all levels of mesh hierarchy
Plex ex12: Change PetscViewerSetFormat --> PetscViewerPushFormat - Fix order of SetMode and SetFilename
Plex ex12: Change all functions to static
Since DMDA with periodic boundary conditions does support P0 interpolation remove checks that error out under those conditions
ascem-io.py and amanzi.py need to handle self.hdf5.include being None correctly
fix MPI Uni to recycle communicators so it won't run out like with the petsc4py tests also fix check for location of python shared library on mac to check all cases each time
Switch back to Sherry's version since she accepted by pull request
PetscDraw: Update help messages for -draw_save* options
...

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

show more ...

Merge branch 'tisaac/dmforest'

* tisaac/dmforest: (36 commits)
DMP4est: fix transfer of labels to parents
DMPlexCreateInterpolatorTree: fix handling of global dofs (can be < 0)
DMPlexComputeIn

Merge branch 'tisaac/dmforest'

* tisaac/dmforest: (36 commits)
DMP4est: fix transfer of labels to parents
DMPlexCreateInterpolatorTree: fix handling of global dofs (can be < 0)
DMPlexComputeInterpolatorTree(): fix workspace size
DM: warn when no injection implemented
DMPlexComputeInterpolatorTree(): use row layout to preallocate
DMP4est: all processes agree on root type
DMPlexComputeInterpolatorTree(): valgrind fixes
DMP4est: reset childID if anctecedent is not parent
DMP4est: send labels from children to parents and bcast
DMForest: added DMForestGetAdaptivityPurpose()
DMP4est: enable caching of results of DMPforestGetTransferSF()
Plex: fix empty partition errors
DMPforestGetTransferSF_Internal(): handle different type roots on different processes
DMP4est: broadcast labels when transfering for agreement
PetscSFBasic: fix usage of MPI_Reduce_local
DMPforestLabelsInitialize(): search through parents when finding a point in the closure
DMP4est: try to connect high strata points when transfering
DMP4est: can't use PetscSFCreateSectionSF. It compresses when we don't want to
DMP4est: expand transfer sfs by the CellSF
DMP4est: add DMForestCreateCellSF_p4est
...

show more ...

DMP4est: pass along geometry function to p4est if supplied

DMForest: added interface for mapping base coordinates to refined coordinates

Merged petsc/petsc into master

DMForest: added DMForestGetAdaptivityPurpose()

DMP4est WIP: compute transfer sfs at setup

Merge branch 'mlange/fix-plex-multi-level-overlap'

* mlange/fix-plex-multi-level-overlap: (259 commits)
Plex: Update ex12_5 with correct level-2 overlap
DMLabel: Add missing cleanup for DMLabelG

Merge branch 'mlange/fix-plex-multi-level-overlap'

* mlange/fix-plex-multi-level-overlap: (259 commits)
Plex: Update ex12_5 with correct level-2 overlap
DMLabel: Add missing cleanup for DMLabelGather()
Plex ex12: Fixed variable initialization and some output
Plex ex11: Added error checking
Plex ex14: Now boundary label does not have cells
Plex ex1: Now boundary label does not have cells
Plex ex1: Fix for VTK format
added some debugging code, getting ready to branch
added some print statements, should be an -info...
smome cleaning up
restructured ex11 to allow for problem specific BCs
DMPlex: Update test ex12_5 with correct level-2 overlap
DMPlex: Use global processSF in DMPlexPartitionInvert
DMPlex: Add remote leaf contributions to roots in DMLabelPropagate
DMPlex: Add DMLabelGather to migrate label values from leafs to roots
DMPlex: Refactor DMLabelDistribute to separate point stratum distribution
DMPlex: Add propagation of remote point donations for multi-level overlap
DMPlex: Fix DMLabelDistribute for all-to-all migrations
Plex: make extern exactly match prototype: semantic match is not enough.
DMP4est: typo fix
...

show more ...

DMFOrest: added DMForestGetAdaptivitySF()

It makes most sense to compute transfer SFs when the new forest is being
setup, but we may not need them if we're doing some initiali a prior
refinement wit

DMFOrest: added DMForestGetAdaptivitySF()

It makes most sense to compute transfer SFs when the new forest is being
setup, but we may not need them if we're doing some initiali a prior
refinement with nothing to transfer.

show more ...

DMForest: comments on the interface

DMForest: added DMForestSetAdaptivityPurpose()

DMDestroy() can't do cycle checks with forest->adapt. I could try to
introduce a robust, extensible cycle check into DM, but that's a bigger
project.

DMForest: added DMForestSetAdaptivityPurpose()

DMDestroy() can't do cycle checks with forest->adapt. I could try to
introduce a robust, extensible cycle check into DM, but that's a bigger
project. This lets the user move forest->adapt to dm->coarseDM or
dm->fineDM, where it can be picked up in cycle detection.

show more ...

DMForestTemplate(): copy more attributes

DMForest: less confusing member name

Merge branch 'tisaac/dmforest' into mlange/fix-plex-multi-level-overlap

Conflicts:
src/dm/impls/plex/examples/tests/output/ex12_5.out
src/dm/label/dmlabel.c

DMForest: default refinement and coarsening

DMP4est: simplify DMForestTemplate() interface

Now when you want to create dmAdapted from dm,
just call

DMForestTemplate(dm,MPI_COMM_NULL,&dmAdapted).

If you called DMSetUp(dmAdapted) right af

DMP4est: simplify DMForestTemplate() interface

Now when you want to create dmAdapted from dm,
just call

DMForestTemplate(dm,MPI_COMM_NULL,&dmAdapted).

If you called DMSetUp(dmAdapted) right after this, the topology of
dmAdapted would essentially be a copy of dm, but you can do things like
set adaptivity labels, change the overlap, and change the partitioning
before calling set up, and then all of those changes will be executed in
DMSetUp().

show more ...

WIP: Simplify/correct adaptivity interface for DMForest

Merge branch 'tisaac/dmforest' of bitbucket.org:petsc/petsc into tisaac/dmforest

DMForest: switch adaptivity from an array of markers to label

Now that labels are general, it makes more sense

DMForest: added DMIsForest() and DMForestRegisterType()