Merge branch 'master' into knepley/feature-pc-patch

* master: (562 commits)
fix clang-analyzer warning
Developers manual never use system randomizers in PETSc and example source code.
missing

Merge branch 'master' into knepley/feature-pc-patch

* master: (562 commits)
fix clang-analyzer warning
Developers manual never use system randomizers in PETSc and example source code.
missing CHKERRQ()
DT: Reverting the tgamma change broke MPFR - I guess we need a nightly build with it in
Add a test for ex21
Fix error in ex21
Add missing fortran interface that handles null objects
DMComposite: add missing guards
MatPreallocate: fix overpreallocation of local rows
style fixes
MatView_MPI_DA: do not print twice the class prefix
MatSetOperation: stash the original MatView_XXX in the viewnative function pointer
DMCreateInterpolation_DA gracefully error instead of catching FPE signals
MatIncreaseOverlap: preserve communicator of index sets
MatMPIAIJ: fix indentation and missing CHKERRQ
ISBuildTwoSided: remove tabs and fix style
Fix documentation
PetscViewerDraw: fix unininitialized warning from valgrind
MatConvert_SeqAIJ_SeqBAIJ: remove error for non square matrices and add MAT_REUSE_MATRIX case
Mat: MAT_INPLACE_MATRIX can be used for transposition with rectangular matrices
...

Conflicts:
src/dm/dt/interface/dtfe.c
src/dm/impls/plex/plex.c
src/dm/impls/plex/plexfem.c
src/snes/examples/tutorials/ex13.c
src/snes/examples/tutorials/ex62.c
src/snes/examples/tutorials/makefile
src/snes/examples/tutorials/output/ex62_18.out
src/snes/examples/tutorials/output/ex62_19.out
src/snes/examples/tutorials/output/ex62_20.out
src/snes/examples/tutorials/output/ex62_20_parmetis.out
src/snes/examples/tutorials/output/ex62_21.out
src/snes/examples/tutorials/output/ex62_22.out
src/snes/examples/tutorials/output/ex62_23.out
src/snes/examples/tutorials/output/ex62_24.out
src/snes/examples/tutorials/output/ex62_25.out
src/snes/examples/tutorials/output/ex62_26.out
src/snes/examples/tutorials/output/ex62_27.out
src/snes/examples/tutorials/output/ex62_28.out
src/snes/examples/tutorials/output/ex62_29.out
src/snes/examples/tutorials/output/ex62_2d_p3_1.out
src/snes/examples/tutorials/output/ex62_30.out
src/snes/examples/tutorials/output/ex62_31.out
src/snes/examples/tutorials/output/ex62_32.out
src/snes/examples/tutorials/output/ex62_33.out
src/snes/examples/tutorials/output/ex62_34.out
src/snes/examples/tutorials/output/ex62_35.out
src/snes/examples/tutorials/output/ex62_36.out
src/snes/examples/tutorials/output/ex62_pc_simple.out
src/snes/examples/tutorials/output/ex62_pc_simplec.out
src/snes/examples/tutorials/output/ex62_quad_q2p1_full.out
src/snes/examples/tutorials/output/ex62_quad_q2q1_full.out
src/snes/utils/dmplexsnes.c

show more ...

Plex: Fixed evaluation of boundary functionals in projection

Plex: Use correct subspace for aux evaluation in projection

Merge branch 'knepley/feature-error-logging' of bitbucket.org:petsc/petsc into knepley/feature-error-logging

Merge branch 'master' into knepley/feature-error-logging

* master: (324 commits)
DT: Reverting the tgamma change broke MPFR - I guess we need a nightly build with it in
Revert "tgamma: icc on wi

Merge branch 'master' into knepley/feature-error-logging

* master: (324 commits)
DT: Reverting the tgamma change broke MPFR - I guess we need a nightly build with it in
Revert "tgamma: icc on windows requires mathimf.h for tgamma."
setup routine no longer calls PetscCalloc5 if the arrays have already been allocated in a previous setup call
clarify SNESGetFunction()
Move MOAB docs to submansec so as not to clutter DM mansec
doc: remove matlab from doc generation dependency - i.e remove 'alldoc3' and use 'html' to generate corresponding '*.m.html' files.
Remove matlab pounders related cruft that is not required. Document the matlab code better and include an appropriate message in the help text.
Fixed typo.
Remove the inclusion of matlab in testing.
Finalizing comments and the work for the test problems.
Finished the testing and generation of plots. Some simple bug fixes addressed.
Set -Wl,-rpath for Mac OS Matlab mkl directory; don't pass -nojvm to matlab since it causes crashes on Barry's machine and is not needed
Added list of known problems with the code.
Added a comment on ignoring the return from Tao.
Finish with initial implementation of matlab interface to test TAO pounders.
Working through the matlab interface for pounders testing.
Fix typo in manual Reported-by: Jaroslaw Piwonski (CAU) <jpi@informatik.uni-kiel.de>
Fixed mkl_sparse_sypr() calls that were using an unsupported value in the matrix descriptor.
test harness: simplify code that checks for sequential-only tests
Fix the doc for call-backs drdyf drdpf
...

show more ...

Merge branch 'knepley/fix-plex-residual-loop'

* knepley/fix-plex-residual-loop:
Plex+SNES: Protect Destroy
Plex+SNES/TS: Moved zeroing of Jacobian out of DMPlexComputeJacobian_Internal() - Shoul

Merge branch 'knepley/fix-plex-residual-loop'

* knepley/fix-plex-residual-loop:
Plex+SNES: Protect Destroy
Plex+SNES/TS: Moved zeroing of Jacobian out of DMPlexComputeJacobian_Internal() - Should probably go in DMLocal versions, but now its in the top-level Plex FEM versions
SNES+Plex: Needed cell offset, not cell number
Sys: Added a check for gnome-terminal, which never uses -display - Added some docs
Plex: Made a public interface for ComputeBdResidualSingle()
Plex: The first argument to a pointfunc should be the coordinate dimension

show more ...

Merge branch 'master' into knepley/feature-error-logging

* master: (221 commits)
dashboard: count and list example compile failures like:
fix fix-pr879-jczhang-optimize-vecscatter-multi-copy
a

Merge branch 'master' into knepley/feature-error-logging

* master: (221 commits)
dashboard: count and list example compile failures like:
fix fix-pr879-jczhang-optimize-vecscatter-multi-copy
add manual page cross links for SNESLineSearchSetPostCheck() and related routines
fix help manual of MatCreateSubmatrices() Reported-by: Tim Steinhoff <kandanovian@gmail.com>
Add manual pages for PetscRealPart() and PetscImaginaryPart()
fix compile error with missing MPI_C_DOBULE_COMPLEX
fix typo Reported-by: Febrian Setianto <febrian.setianto@gmail.com>
plex tutorial ex5: updated test requirements
DMPlexLoad_HDF5_Xdmf_Internal: fix for complex PetscScalar.
DMPlexLoad_HDF5_Xdmf_Internal: fix for 64-bit PetscInt.
DMPlexBuildFromCellList_Internal: use DMPlexInvertCell_Internal (with PetscInt* arg)
petscviewer.h: removed notice which is no longer relevant
DMPlexCreateFromFile: document option
check only if PETSC_USE_DEBUG
DMPlexBuildFromCellList_Internal optionally inverts cells.
ex5: new test 4 with -compare without repartitioning
ex5: new test 3 using -compare
ex5: HDF5 write+read moved to separate function
ex5: option -format processed in ProcessOptions()
ex5: added missing hdf5 requirement
...

Conflicts:
include/petsclog.h

show more ...

Merge branch 'master' into knepley/feature-snes-ex13-spectral

* master: (2350 commits)
migrate warning fix from next-tmp
Revert "Remove dead code."
clarifying comment in the code
Remove dead

Merge branch 'master' into knepley/feature-snes-ex13-spectral

* master: (2350 commits)
migrate warning fix from next-tmp
Revert "Remove dead code."
clarifying comment in the code
Remove dead code.
Add JuMP citation
Copyedits and updates to developers manual
Add more protection to intrinsics
MPI: protect arguments in MPI macros
Add the missing conditionals for 64-bit build
Plex+VTK: Fix field type determination
Fix bug in name cleanup
Bib: Added ref
Plex+VTK: Fix field type determination
Plex: Updated test output - I am now hating this shit
forgot to save file after moving prototypes to petscsys.h
update Apple example project settings per Xcode
Two functions should actually be public, used, for example from Apple MacOS and iOS programs
includes should always be included relative to include directory
Apple removed old locations of docsets
SNES ex12: Updated test output
...

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

show more ...

Merge branch 'knepley/fix-plex-pointfunc-dim' into knepley/fix-plex-residual-loop

* knepley/fix-plex-pointfunc-dim:
Sys: Added a check for gnome-terminal, which never uses -display - Added some do

Merge branch 'knepley/fix-plex-pointfunc-dim' into knepley/fix-plex-residual-loop

* knepley/fix-plex-pointfunc-dim:
Sys: Added a check for gnome-terminal, which never uses -display - Added some docs
Plex: Made a public interface for ComputeBdResidualSingle()
Plex: The first argument to a pointfunc should be the coordinate dimension

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

show more ...

Merge branch 'master' into knepley/fix-plex-residual-loop
Warning: This merge includes real changes

* master: (188 commits)
Add JuMP citation
Copyedits and updates to developers manual
Add mor

Merge branch 'master' into knepley/fix-plex-residual-loop
Warning: This merge includes real changes

* master: (188 commits)
Add JuMP citation
Copyedits and updates to developers manual
Add more protection to intrinsics
MPI: protect arguments in MPI macros
Add the missing conditionals for 64-bit build
Plex+VTK: Fix field type determination
Fix bug in name cleanup
Bib: Added ref
Plex+VTK: Fix field type determination
Plex: Updated test output - I am now hating this shit
forgot to save file after moving prototypes to petscsys.h
update Apple example project settings per Xcode
Two functions should actually be public, used, for example from Apple MacOS and iOS programs
includes should always be included relative to include directory
Apple removed old locations of docsets
SNES ex12: Updated test output
SNES ex62: Set name of solution vector - Update test output
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html
release: set v3.9 strings
move output file to the correct location. i.e fix 305862f9cd0dc399542ae11dd744508e4b6ad2e3
...

Conflicts:
src/snes/utils/dmplexsnes.c
src/ts/utils/dmplexts.c

show more ...

Merge remote-tracking branch 'origin/tisaac/feature-dmfield'

* origin/tisaac/feature-dmfield: (56 commits)
Plex: Updated test output for DMLabel change
DMLabel: make empty ISs strides
DMSwarm:

Merge remote-tracking branch 'origin/tisaac/feature-dmfield'

* origin/tisaac/feature-dmfield: (56 commits)
Plex: Updated test output for DMLabel change
DMLabel: make empty ISs strides
DMSwarm: silence static analyzer warning
PetscFECreateDefault: fix custom Fortran interface to use MPI_Comm
SNES ex69: Updated test output
SNES ex69: Removed all mantle stuff and change fields to constants - Also use dsymutil for executables
SNES ex69: Removed special partition
SNES ex69: Fixed null space creation - Null space locks its vectors
PetscFE: PetscFECreateDefault() now takes MPI_Comm instead of DM
DMField: default initialization of variable
PetscFE: fix integration for openCL
PetscFEGeomComplete(): correct normal vector in 2D.
DMField: c++/complex fixes
DMField: empty processor protections
DMField regression test fixes
DMLabel: try to keep stride IS as stride IS as much as possible
DMPlexRefine(): create depth before hybrid-bounded label
DMField: documentation of interface
DMField: implement shell
DMPlex: rename newly created label
...

Conflicts:
src/dm/dt/interface/dtfe.c
src/snes/examples/tutorials/ex69.c
src/snes/utils/dmplexsnes.c

show more ...

Merge remote-tracking branch 'origin/knepley/fix-fe-bd-integral'

* origin/knepley/fix-fe-bd-integral:
Plex+SNES: Fix BdResidual - Use DMPlexGetSubpoint() for aux closure, in case it is a boundary

Merge remote-tracking branch 'origin/knepley/fix-fe-bd-integral'

* origin/knepley/fix-fe-bd-integral:
Plex+SNES: Fix BdResidual - Use DMPlexGetSubpoint() for aux closure, in case it is a boundary mesh - Be careful about Plex conversion
FE: Fixed boundary integration - Allow auxiliary dm to be a different dimension, need dimAux and different tabulation
FE: Small confusion about topological vs coordinate dimension
PetscDS: The u_x array must have size dimEmbed because the derivatives are a function of the coordinate dim, not topological dim
Plex: Added DMPlexGetSubpoint() - This returns the point in a submesh corresponding to some point in the original mesh - This replaces custom code in projection

show more ...

Deprecate DMSet/GetDefaultSection() for DMSet/GetSection()

Revert manual page change per Matt's preference

Commit-type: style-fix
Thanks-to: Matthew Knepley

Plex: Added DMPlexGetSubpoint()
- This returns the point in a submesh corresponding to some point in the original mesh
- This replaces custom code in projection

DMPlexSNES: caching of new geometry constructs

The old paradigm cached one geometry object for all cells or all faces,
an approach that doesn't take into account that different quadratures
may requi

DMPlexSNES: caching of new geometry constructs

The old paradigm cached one geometry object for all cells or all faces,
an approach that doesn't take into account that different quadratures
may require different geometry objects. Now the geometry object
is composed with the index set over which the geometry ranges, and the composition
string include the object id of the quadrature. To take advantage of
this caching, we have to be careful not to reconstruct index sets on the
fly when they are known a priori. In particular, the hybrid bounds
have been encapsulated as a label, which makes them accessible as ISs.

show more ...

Use DMField / PetscFEGeom in integration

Added PetscDualSpaceGetAllPointsUnion

Plex: The first argument to a pointfunc should be the coordinate dimension

Merge branch 'barry/tao-to-new-test-harness' into barry/all-new-test-harness-branches

Merge branch 'barry/vec-to-test-harness' into barry/all-new-test-harness-branches

Merge branch 'master' into knepley/feature-plex-adjoint-prelim

* master: (60 commits)
petsc4py: update to use latest 3.8 snapshot
petsc4py: update to snapshot that syncs with 276c550
fix warni

Merge branch 'master' into knepley/feature-plex-adjoint-prelim

* master: (60 commits)
petsc4py: update to use latest 3.8 snapshot
petsc4py: update to snapshot that syncs with 276c550
fix warning
Plex ex26: Protect ExodusII stuff
DM: DMCreateSubDM() field argument should be const
TS ex46: Newlines and some numerical jitter
TS ex45: Newlines and some numerical jitter
TS ex11: Fixed indentation and time step numbering
SNES ex12: Update for DM view changes
Plex tests: Update to new DM output
Plex ex23: Use Nc not dim for vector fields
Plex: Default section now ignores cells hanging off boundary meshes
Plex: Viewer now reports the cell Height
Increase patchlevel to 3.8.3
Plex: Fixed DMPlexInterpolate() bugs - Do not copy associated data if we are not interpolating - Do not interpolate if we already have full depth
fix valgrind errors with using PetscObjectGetName()/SetName() with the same object - resulting in a free() before copy.
ex4: fix compile warning and update makefile
Minor fixes
Improving report_tests
PetscDS: Must initialize the default adjacency
...

Conflicts:
include/petsc/private/dmimpl.h
include/petscdm.h
include/petscis.h
src/dm/impls/plex/plexcreate.c

show more ...

Merge remote-tracking branch 'origin/knepley/fix-plex-multifield-eval'

Merge branch 'master' into knepley/feature-plex-adjoint-prelim

* master: (251 commits)
yaml: add a backup download URL
Fix documentation of MatSetRandom function
Ignore -malloc-hbw option if m

Merge branch 'master' into knepley/feature-plex-adjoint-prelim

* master: (251 commits)
yaml: add a backup download URL
Fix documentation of MatSetRandom function
Ignore -malloc-hbw option if malloc is already set
Plex: Fixed array initialization
pgi: compile fix
knl: update test to use memkind and remove option that worked arround a bug with intel 17
Check the return value for memkind_posix_memalign
Change the default behavour for -malloc_hbw
Bib: More refs
test: remove runex13f90 as its now converted to new test suite
fix win32 crash. This test was previously run with '-n 4' - so restoring it
Fix petscsfgetgraph_ to work for all Fortran compilers by directly using array of PetscSFNodes
bugfix for MatDuplicate_MPIAIJ_PtAP Reported-by: Samuel Lanthaler <s.lanthaler@gmail.com>
strip out dead sundials usage code; change sundials test to only run 5 timesteps to prevent different number of steps on different systems.
Remove extra lines of output file
Example uses C99 constructs so must mark with requires: c99
fix const missmatch for DMDAVecRestoreArrayRead()
remove conflict between function name in example and new PETSc TS function
builds: attempt to do maint/master/next test builds in separate clones
ex62: it turned out that CHACO partioning is not deterministic. Use a prescribed partition instead
...

Conflicts:
src/dm/interface/dm.c
src/snes/examples/tutorials/ex13.c

show more ...

Merge branch 'master' into knepley/fix-plex-bc-multiple

* master: (347 commits)
tests: add 'PETSC_HAVE_C99' flag - so that 'requires: c99' can be used in the test suite
PGI: fix build error with

Merge branch 'master' into knepley/fix-plex-bc-multiple

* master: (347 commits)
tests: add 'PETSC_HAVE_C99' flag - so that 'requires: c99' can be used in the test suite
PGI: fix build error with updated 17.10 compilers
Fix for test dirs that are not tutorials/tests
fix use of double in PetscAbsReal() which produces compiler warnings with single precision build
elemental: remove info on conflict of -with-cxx-dialect=C++11 on osx - it works now.
Cleanup .bib files
Fix allfortranstubs rule to rm generated files from correct directory, not outdated location
ierr set but not checked
Remove initial refine so test runs faster and does not time out on some machines
Gail's copyedits for bibtex items
Docs: manual: add further description of 3 different current testing options
DM: check for ops->localtolocalbegin/end before calling
PetscSF: fix typo in tutorial ex1
Vec+HDF5: Fix output for complex - Need to check for empty group name
SNES: SNESDestroy() can be called with snes->dm == NULL
SNES: SNESDestroy() can be called with snes->dm == NULL
Vec+HDF5: Fix output for complex - Need to check for empty group name
Docs: Manual: Update tests section
gmakegentest: also exclude .#filenames created by emacs
DMPlex: Fix output of test ex1 after changes to Gmsh reader
...

show more ...

Plex: Fixed projection routines when we iterate over a surface with volume field and surface aux fields
- If we are iterating over a surface with volume input, then the effective DS is the height=1 s

Plex: Fixed projection routines when we iterate over a surface with volume field and surface aux fields
- If we are iterating over a surface with volume input, then the effective DS is the height=1 subspace
- Made DS an argument for projection functions
- Only allow iteration over a surface with volume input when we have a surface auxiliary field with FIELD type function
- Documented DMProjectLocal_Generic_Plex()

show more ...