Merge branch 'tisaac/dm-plex-fix-nonconforming'

Changes to DMPlexMigrate() broke the nonconforming examples. The fix was
simple (reintroduce tree/closure traversal when computing partition labels),

Merge branch 'tisaac/dm-plex-fix-nonconforming'

Changes to DMPlexMigrate() broke the nonconforming examples. The fix was
simple (reintroduce tree/closure traversal when computing partition labels),
but the bigger problem was that I didn't have the nonconforming examples in
the test suite. To test in as many environments as possible, I made the
nonconforming tests independent of the partitioning packages available
(several don't support 64bit indices), and I had to make
PetscFEGetTabulation() work with complex scalars.

show more ...

Plex: use DMPlexProjectFunction() in DMPlexRigidBody()

Merge branch 'master' into stefano_zampini/pcbddc-adaptive-selection

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c

Merge branch 'master' into barry/saws-options-merge-master

Totally confusing manual merge, could be wrong

Conflicts:
include/petscoptions.h
src/dm/impls/da/dacreate.c
src/sys/objects/aoptions.c

Merge branch 'master' into barry/saws-options-merge-master

Totally confusing manual merge, could be wrong

Conflicts:
include/petscoptions.h
src/dm/impls/da/dacreate.c
src/sys/objects/aoptions.c
src/sys/threadcomm/interface/threadcomm.c

show more ...

Merge branch 'master' into pr252/dmay/PetscViewerBinarySkipHeaderVecViewLoad
Manual merge of makefile

Conflicts:
src/sys/classes/viewer/examples/tutorials/makefile

Merge branch 'tisaac/feature-gamg-extensible' into barry/experiment

Conflicts:
src/ksp/pc/impls/gamg/gamg.h

Merge branch 'knepley/feature-ts-fem+fvm-example'

* knepley/feature-ts-fem+fvm-example: (32 commits)
Plex: Fix printf
Plex: Use DMGetCoordinateDim() instead of a calculation
Plex: Added routin

Merge branch 'knepley/feature-ts-fem+fvm-example'

* knepley/feature-ts-fem+fvm-example: (32 commits)
Plex: Fix printf
Plex: Use DMGetCoordinateDim() instead of a calculation
Plex: Added routine to view geometry
TS ex18: Added new tests
TS ex18: Updated test output
TS ex18: Add support for functionals
TS ex18: Global field offsets depend on the point since we allow constraints
TS ex18: Move view to after LocalizeCoordinates()
Plex+HDF5: Fixed uninitialized variable
TS ex18: Advection tests working
Plex: Fixed integration with Toby's code that handles sub manifold interpolation - Needed to check for cEndInterior from adding ghost cells
TS ex18: Updated output
TS ex18: Trying out implicit and explicit TS - Added another test with implicit TS
TS ex18: Updated output
SNES ex62: Updated output for new Plex distribution
Plex: Fix leak - No need for extra BC copy
TS ex18: All examples now working - Added ability to use zero velocity - Fixed coordinate localization
TS ex18: Cast boundary values function
TS ex18: Fixes outflow boundary condition - Must index into vector of all field values
TS ex18: Added another advection test
...

show more ...

Merge branch 'knepley/feature-combined-tabulation'

* knepley/feature-combined-tabulation: (38 commits)
Plex: Fix compiler warning
Plex+SNES: Fixes for nonzero cStart
Plex: Preserve periodicity

Merge branch 'knepley/feature-combined-tabulation'

* knepley/feature-combined-tabulation: (38 commits)
Plex: Fix compiler warning
Plex+SNES: Fixes for nonzero cStart
Plex: Preserve periodicity when adding ghost cells
Plex: Fixed DMPlexGetGlobalFieldOffset_Private() in optimized mode
Plex+SNES: Offset should be into field components, not basis functions
Plex+SNES: Add guards for optional arrays
DS: Added PetscDSGetComponentOffset(), Added to petsc.h
Plex+TS: FVM RHS computation uses a global F
Plex+FV: Corrected accumulation of fluxes to cells and added time derivative term
Plex: Added DMPlexPointGlobalFieldRead/Ref()
Plex+SNES: Fixed DMPlexGetFaceFields() - It returns field values, not coefficients - FE evaluation is now done
FV: In PetscFVIntegrateRHSFunction(), the arguments are field values, not coefficients
FE: Added EvaluateFaceFields() to FE fields at face centroids
Plex: Added DMPlexPointLocalFieldRead()
FE: Have a standard tabulation at face centroids
Plex: Copy BC for hierarchical refinement
FE: Allow P0 over 1D cells
Plex+FV: Fixed BC with gradient calculation
Plex ex1: Move overlap computation after SetFromOptions
Plex: DMPlexDistribute() no longer calls DMSetFromOptions() - Fixed examples
...

show more ...

Plex: Fix compiler warning

Plex: Fixed integration with Toby's code that handles sub manifold interpolation
- Needed to check for cEndInterior from adding ghost cells

Merge branch 'knepley/feature-combined-tabulation' into knepley/feature-ts-fem+fvm-example

* knepley/feature-combined-tabulation: (41 commits)
Plex: Fix distribution of uninterpolated meshes - Doe

Merge branch 'knepley/feature-combined-tabulation' into knepley/feature-ts-fem+fvm-example

* knepley/feature-combined-tabulation: (41 commits)
Plex: Fix distribution of uninterpolated meshes - Does SF stratification by depth, not by dimension - The old code numbered cells after vertices
Plex+SNES: Fixes for nonzero cStart
cannot call PetscObjectGetName() on a subset of the MPI processes
MatView() with some matrix formats and ASCII_MATLAB did not print correct matrix name
Plex: Preserve periodicity when adding ghost cells
Plex: Fixed DMPlexGetGlobalFieldOffset_Private() in optimized mode
Plex+SNES: Offset should be into field components, not basis functions
Plex+SNES: Add guards for optional arrays
DS: Added PetscDSGetComponentOffset(), Added to petsc.h
Plex+TS: FVM RHS computation uses a global F
Plex+FV: Corrected accumulation of fluxes to cells and added time derivative term
Plex: Added DMPlexPointGlobalFieldRead/Ref()
Plex+SNES: Fixed DMPlexGetFaceFields() - It returns field values, not coefficients - FE evaluation is now done
FV: In PetscFVIntegrateRHSFunction(), the arguments are field values, not coefficients
FE: Added EvaluateFaceFields() to FE fields at face centroids
Plex: Added DMPlexPointLocalFieldRead()
FE: Have a standard tabulation at face centroids
Plex: Copy BC for hierarchical refinement
FE: Allow P0 over 1D cells
Plex+FV: Fixed BC with gradient calculation
...

Conflicts:
config/builder.py

show more ...

Plex+FV: Fixed BC with gradient calculation

Plex+FV: Hack to cope with BC for output

Plex+FV: When setting FV BC, get the correct per point offset

FV: DMPlexInsertBoundaryValues_FVM_Internal() must put values at the correct field offset

FE+FV: When mixing, FE should ignore hybrid (ghost) cells which FV uses for BC, and FV should evaluate on cells
- When creating default section, FE should not run over hybrid cells
- In DMPlexCompute

FE+FV: When mixing, FE should ignore hybrid (ghost) cells which FV uses for BC, and FV should evaluate on cells
- When creating default section, FE should not run over hybrid cells
- In DMPlexComputeResidual_Internal(), FE should not run over hybrid cells
- In EvaluateFieldJets(), FV does cellwise constant evaluation (should work in linear but demands per cell lookup which should be worked into the interface)
- In DMPlexProjectFunctionLabelLocal(), DMProejctFunctionLocal(), DMProjectFieldLocal(), DMPlexComputeL2Diff(), DMPlexComputeL2GradientDiff(), DMPlexComputeL2FieldDiff(), DMPlexComputeIntegralFEM(), DMPlexCompute InterpolatorFEM(), and DMPlexComputeInjectorFEM(), do not loop over hybrid cells
- Protect use of ghost label in DMPlexComputeGeometryFVM()

show more ...

Merge branch 'master' into mlange/refactor-plex-distribute

Conflicts:
src/dm/impls/plex/examples/tests/output/ex12_1.out
src/dm/impls/plex/examples/tests/output/ex12_2.out

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

Merge remote-tracking branch 'origin/tisaac/plex-flexible-projection'

* origin/tisaac/plex-flexible-projection:
PetscDualSpace: Fixed function signature
Plex: update libraries for ex3
Plex: mo

Merge remote-tracking branch 'origin/tisaac/plex-flexible-projection'

* origin/tisaac/plex-flexible-projection:
PetscDualSpace: Fixed function signature
Plex: update libraries for ex3
Plex: move DMPlexProjectionFunction() back to libdm
Plex: use -dm_plex_max_projection_height in non-conforming examples
Plex: fixed a typo in DMPlexProjectFunctionLocal()
Plex: added DMPlexComputeProjection3Dto1D_Internal()
Plex: fixes to recent merge
PetscFE: added dim, dimEmbed to PetscFECellGeom
Plex: skip empty loops for empty dual spaces
DualSpace: 0-order Lagrange spaces do not have dual spaces for height > 0.
Plex: protect against field projection with height > 0.
Plex: start implementing projection with height > 0.
PetscDualSpace: added PetscDualSpaceGetHeightSubspace()
Plex: added DMPlexGet/SetMaxProjectionHeight()

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

show more ...

Merge branch 'master' into knepley/fix-plex-examples-refinement

* master: (62 commits)
Plex+Label: Fix logic for label distribution - It is restrictive now, but doing it right would mean a ton of

Merge branch 'master' into knepley/fix-plex-examples-refinement

* master: (62 commits)
Plex+Label: Fix logic for label distribution - It is restrictive now, but doing it right would mean a ton of infrastructure I think
Plex ex12: Fix output after label distribution fix
Plex: Only call SetFromOptions() after Distribute, since we cannot avoid it
TS ex11: Since DMPlexDistribute() calls SetFromOptions(), do not call it again, and cleaned up tests
Plex+FV: If using FVM adjacency, do not create ghost cells for faces in the overlap, or set BC on them
Plex: If using FVM adjacency, we must include the closure in the overlap
Plex+Label: All processes with a valid label should participate in label distribution - Formerly it was just proc 0
Plex: Add check that parallel renumbering fully covers the cone points
Plex: For one level adjacency, include the point itself since this makes sense when looking in parallel from the overlap
Plex: Cannot Malloc2() the arrays passed to PetscSFSetGraph()
Plex: Preallocation fixes for indegree > 1 - Again this should be rewritten
SNES ex62: Changes to output for fixes with quads
Plex+Label: Set default output flag to PETSC_TRUE
Plex ex1: Fixed LaTeX output
Plex ex12: Fixed overlap test
Plex: Fixed incorrect usage of PetscSFGatherBegin/End() in preallocation - This should be rewritten using new style
PetscSF: Incorrect check for computed degrees
Builder: Fix bad search and replace
Configure: Fix for Barry's change to organization - I do not know what happened here
Builder: Fix for Barry's change
...

Conflicts:
src/snes/examples/tutorials/ex12.c

show more ...

Merge branch 'mlange/feature-parallel-distribute'

* mlange/feature-parallel-distribute: (33 commits)
Plex+Label: Fix logic for label distribution - It is restrictive now, but doing it right would

Merge branch 'mlange/feature-parallel-distribute'

* mlange/feature-parallel-distribute: (33 commits)
Plex+Label: Fix logic for label distribution - It is restrictive now, but doing it right would mean a ton of infrastructure I think
Plex ex12: Fix output after label distribution fix
Plex: Only call SetFromOptions() after Distribute, since we cannot avoid it
TS ex11: Since DMPlexDistribute() calls SetFromOptions(), do not call it again, and cleaned up tests
Plex+FV: If using FVM adjacency, do not create ghost cells for faces in the overlap, or set BC on them
Plex: If using FVM adjacency, we must include the closure in the overlap
Plex+Label: All processes with a valid label should participate in label distribution - Formerly it was just proc 0
Plex: Add check that parallel renumbering fully covers the cone points
Plex: For one level adjacency, include the point itself since this makes sense when looking in parallel from the overlap
Plex: Cannot Malloc2() the arrays passed to PetscSFSetGraph()
Plex: Preallocation fixes for indegree > 1 - Again this should be rewritten
SNES ex62: Changes to output for fixes with quads
Plex ex1: Fixed LaTeX output
Plex ex12: Fixed overlap test
Plex: Fixed incorrect usage of PetscSFGatherBegin/End() in preallocation - This should be rewritten using new style
PetscSF: Incorrect check for computed degrees
Builder: Fix for Barry's change
PetscSF: Fixes for sparse leaves - When sf->mine exists, the leaves can be numbered non-sequentially. This causes allocation based upon sf->nleaves to underallocate. Replacing sf->nleaves with the maximum leaf number is sufficient in this case.
Plex: Fix leaks
Configure: How did this break?
...

Conflicts:
config/builder.py

show more ...

Plex+FV: If using FVM adjacency, do not create ghost cells for faces in the overlap, or set BC on them

Merge branch 'master' into knepley/fix-configure-petsc-dep

* master: (210 commits)
error out on special cases where DMLocalToGlobal_DA won't work
remove unneed uses of local DM vectors and DMLoc

Merge branch 'master' into knepley/fix-configure-petsc-dep

* master: (210 commits)
error out on special cases where DMLocalToGlobal_DA won't work
remove unneed uses of local DM vectors and DMLocalToGlobalBegin/End in examples removed to prevent cut and paste repeats of their use
openmpi: update to version 1.8.3
added more movies and pubs
more movies
MatGetDiagonalBlock() can be called before matrix is assembled
added a couple more movies/images
PetscLogEventGetId() was missing public prototype
simply output from make streams benchmark test
c2html: fixup so that C2HTML make variable is set [similar to cmake.py]
chaco: current blas test does'nt work with Cray CC [which automatically links with blas] and linux - where -lm needs to be listed explicitly. Switch the test to check if we are using PETSc modifled chaco
builddist: always install c2html/sowing
argonne logo for users manual cover
Remove PETSC_HAVE_CUSP_SMOOTHED_AGGREGATION, use PETSC_HAVE_CUSP instead.
mpich: use older tarball with cygwin/gnu compilers
updated information about packages PETSc can use in users manual
Updated users manual to latest ANL Report format requirements
mpich: update to mpich-3.1.3.tar.gz - remove -download-mpich-mpe [mpich no longer packages mpe] - remove -download-mpich-shared. It was added as a split from --with-mpi-shared to --known-mpi-shared and --download-mpich-shared in 50237a55163f90ac7341b723c6422dfa1dbce792 removing this gives back the control of shared=0/1 to package.py - remove download_solaris as the new tarball works with both gcov and on solaris - remove MPICHInstall() [and use GNUPackage.Install()] as the new tarball fixes 'make clean' requiring sudo
added more images and PETSc pubs
Doc: More refs
...

show more ...

Plex: move DMPlexProjectionFunction() back to libdm

When I first implemented hanging constraints, I moved
DMPlexProjectFunction() to libksp because I used a global solve to
determine the values of d

Plex: move DMPlexProjectionFunction() back to libdm

When I first implemented hanging constraints, I moved
DMPlexProjectFunction() to libksp because I used a global solve to
determine the values of dofs on an anchor point that is not in the
closure of any cell. This was a hack, as those values do not require a
global solve, but can be determined with functionals evaluated on the
anchor point: -dm_plex_max_projection_height > 0 allows functionals to
be evaluated on these anchor points, so the global solve is no longer
necessary.

To see the effect, you can modify the scripts to run ex3 in
src/dm/impls/plex/examples/tests: set -dm_plex_max_projection_height 0
in the runex3_nonconforming_{simplex,tensor}_{2,3} examples, and
projecting a function fails.

I did not move DMPlexProjectField() back into libdm because it evaluates
field jets, not just field values, and doing this on non-cells isn't
straightforward.

show more ...

Plex: fixed a typo in DMPlexProjectFunctionLocal()