Plex: After creating GtoN map, set useNatural for parallel DM

Merge remote-tracking branch 'origin/knepley/fix-plex-fvm-adjacency'

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 ...

PetscDS: Can now store default adjacency when it has no fields
- Get it with PETSC_DEFAULT for field

Plex+DS: Moved adjacency information completely to DS
- The Plex versions now just delegate to DS
- Defaults in DS are clearer
- Input to GetAdjacency are optional

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 ...

Merge branch 'maint'

* maint:
Spell check: interger -> integer
petscsys.h : partial man page cleanup
configure exodusii: pass HDF5_DIR to Exodus cmake
Docs: DMPlex man page and source commen

Merge branch 'maint'

* maint:
Spell check: interger -> integer
petscsys.h : partial man page cleanup
configure exodusii: pass HDF5_DIR to Exodus cmake
Docs: DMPlex man page and source comments fixes
Docs: Manual: clarifications for DMPlex chapter

show more ...

Merged in psanan/petsc-docproject/psanan/dmplex-docs-update (pull request #773)

Psanan/dmplex docs update

Approved-by: Matthew Knepley <knepley@gmail.com>

DMPlexDistribute: fix undefined output sf in sequential run

Merged in haplav/petsc/haplav/fix-dmplexdistribute-seq (pull request #775)

fixed bug in DMPlexDistribute for sequential runs

Approved-by: Matthew Knepley <knepley@gmail.com>

DMPlexDistribute: fix undefined output sf in sequential run

Docs: DMPlex man page and source comments fixes

* remove sieve mentions from all DMPlex man pages and source comments
* DMPlexSymmetrize man page: typo fix
* DMPlexCreateExodusFromFile man page: fix

Docs: DMPlex man page and source comments fixes

* remove sieve mentions from all DMPlex man pages and source comments
* DMPlexSymmetrize man page: typo fix
* DMPlexCreateExodusFromFile man page: fix function name
* DMPlexCreateCubeBoundary man page typo fix

show more ...

Merge remote-tracking branch 'origin/pr690/mapdes/feature-plex-user-adjacency/master'

plex: Support user-defined adjacencies via callback

Merge branch 'master' into knepley/fix-fe-default

* master: (297 commits)
Bib: Added refs and fixed one
Doc: Fixed formatting problems for MatCreateAIJ() manual page.
Add to PCSOR manual page

Merge branch 'master' into knepley/fix-fe-default

* master: (297 commits)
Bib: Added refs and fixed one
Doc: Fixed formatting problems for MatCreateAIJ() manual page.
Add to PCSOR manual page information about convergence testing when using KSPRICHARDSON
SNES ex5: Restore old output since this breaks an example I have
Noted that some Cray XC systems will use 'aprun' instead of 'srun'. Thanks to Herr Doktor Hong Zhang for this suggestion.
Added example configure script for Cray XC40 with KNL processors (e.g., NERSC Cori).
Fix possible sprintf -Wformat-overflow
PetscMalloc2-7: fix -Wint-in-bool-context for GCC 7
SNES ex5: make boundary condition handling consistent for symmetry
SNES ex5: use L2 norm and *_monitor_short for output
SNES ex5: the L^2 norm scales with the l^2 norm over sqrt(N)
SNES ex5: refactor MMS capability to avoid extreme code duplication
Documentation fix: typo in PetscLogStageRegister manual page.
nodes -> vertices
rename DMNetworkGetComponentTypeOffset -> DMNetworkGetComponentKeyOffset
rename DMNetworkGetConnectedNodes() -> DMNetworkGetConnectedVertices()
Jacobians for adjoints for explicit methods and constant Jacobians never got initialized
fix complex build error
fix build failure with fortran interface
fix compile error with cdc753b64c658d532d7813e942002e67025667a6
...

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

show more ...

Merge remote-tracking branch 'origin/stefano_zampini/feature-gmsh-periodic'

Merge branch 'master' into stefano_zampini/feature-gmsh-periodic

* master: (1016 commits)
metis: setself.downloaddirnames for change in 6d90219d8fa4f8b2fd242cc6ff9008358d87c26f
fix compile warni

Merge branch 'master' into stefano_zampini/feature-gmsh-periodic

* master: (1016 commits)
metis: setself.downloaddirnames for change in 6d90219d8fa4f8b2fd242cc6ff9008358d87c26f
fix compile warning: /sandbox/petsc/petsc.clone/src/snes/examples/tests/ex2.c:250:8: warning: ‘pointsAllProcs’ may be used uninitialized in this function [-Wmaybe-uninitialized] ierr = DMInterpolationSetUp(interpolator, dm, pointsAllProcs);CHKERRQ(ierr);
migrate 'build requires: !mpiuni' to new format
migrate 'build requires: !complex' to new format
workarround 'division by zero' compiler warning on windows C:\cygwin64\home\petsc\PETSC~3.CLO\src\ts\examples\TUTORI~1\ex11.c(416): warning #39: division by zero if (uL->h < 0 || uR->h < 0) {for (i=0; i<1+dim; i++) flux[i] = 0./0.; return;} /* SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Reconstructed thickness is negative"); */
dos2unix to prevent warnings on opensolaris "/export/home/petsc/petsc.clone/src/ksp/ksp/examples/tests/ex53.c", line 6: warning: invalid white space character in directive
fix output file - corresponding to code changes in efd4aadf157bf1ba2d80c2be092fcf4247860003
disable test on windows
update log from opensolaris so that the diff is clean on linux aswell
update grep in sync with change in efd4aadf157bf1ba2d80c2be092fcf4247860003
Fixed code that did not correctly set current color for block insertion
Update output files that did not get fixed after previous merges
Change MPIUni to not immediately abort on operations it cannot perform, instead return error code
Add filter for example that can produce 1 different iteration count
Example add system calloc, bad formating of %d in print, commented out code
Always use (void) for functions that take no arguments, never ()
always need to caste PetscReal variables to (double) in print statements
updates to manual page for clarity
update help menu for MatRARt() -- MatRARt_mpiaij_mpiaij is supported now, but not efficient.
fix nightlybuild error with c1cbb1def3f8ddca04c3998b9c06618fd1ffef35 on arch-linux-IPL64 src/ksp/pc/examples/tests ex5_1
...

Conflicts:
include/petscdm.h
src/dm/impls/plex/plexcreate.c
src/dm/impls/plex/plexgmsh.c
src/dm/impls/plex/plexhdf5.c

show more ...

Merge branch 'master' into knepley/feature-plex-wedge

* master: (61 commits)
Added ignore
Plex ex7: Fixed mesh name
Plex ex1/ex2: Fixed requirements
PCGAMG: fix spelling in option docs
Pet

Merge branch 'master' into knepley/feature-plex-wedge

* master: (61 commits)
Added ignore
Plex ex7: Fixed mesh name
Plex ex1/ex2: Fixed requirements
PCGAMG: fix spelling in option docs
PetscDS: Fix print
Plex ex19: Output fixes for proper destruction
TSMonitorDrawSolutionPhase: Add missing draw pause
fix warning with gcc-7.1
Docs+PetscFE: Evaluate() operates in reference space
Plex ex7: Added 600 cell test
Plex: Added 600-cell for S^3
Math: Add the golden ratio
Plex: Handle unused variable
Plex ex7: Meshing a sphere
Plex: Changed to DMPlexCreateSphereMesh() - Now can do quads or triangles
Share: Added large 3D mesh file - Its the only one I have gotten to work
Mat: Added example for MatIS
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
PetscDS: Fix update fund for new API
...

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

show more ...

Merge branch 'master' into knepley/fix-fe-fortran

* master: (56 commits)
PCGAMG: fix spelling in option docs
PetscDS: Fix print
Plex ex19: Output fixes for proper destruction
TSMonitorDrawSo

Merge branch 'master' into knepley/fix-fe-fortran

* master: (56 commits)
PCGAMG: fix spelling in option docs
PetscDS: Fix print
Plex ex19: Output fixes for proper destruction
TSMonitorDrawSolutionPhase: Add missing draw pause
fix warning with gcc-7.1
Plex ex7: Added 600 cell test
Plex: Added 600-cell for S^3
Math: Add the golden ratio
Plex: Handle unused variable
Plex ex7: Meshing a sphere
Plex: Changed to DMPlexCreateSphereMesh() - Now can do quads or triangles
Share: Added large 3D mesh file - Its the only one I have gotten to work
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
PetscDS: Fix update fund for new API
Plex: Fixing insertion of boundary values - If we have a dmAux that is a boundary, we need minHeight = 1, since faces can be mapped from submesh to mesh, but not cells (since they have incomplete cones) - However, in this case the effectiveHeight should be h, whereas when using a boundary as the mesh as well, the effectiveHeight should be h - minHeight - Obviously, we need more thinking here to make the iteration make sense in general
PetscFE: Derivatives for Point class are all zero
Small fixes
TS ex46: Fix sin and cos
TS ex18: Fix sin and cos
...

show more ...

Merge branch 'master' into knepley/feature-ex11-simplex

* master:
Plex ex19: Output fixes for proper destruction
Plex: Handle unused variable
Plex ex19: Now working in parallel - Added partiti

Merge branch 'master' into knepley/feature-ex11-simplex

* master:
Plex ex19: Output fixes for proper destruction
Plex: Handle unused variable
Plex ex19: Now working in parallel - Added partition visualization field - Added some tests
Configure: Bump Pragmatic version
Plex: Update interface and fix boundary labels for Pragmatic
Plex ex19: Create mesh overlap before adaptation - Put coords in metric for checking
Plex: Take overlap generation out of Adapt() for now
Plex ex19: Add another Pragmatic test
Plex: Use new parallel Pragmatic interface
Plex ex19: Added a parallel test
Plex ex19: Make metric tensor a local vector
Plex: Overlap distribution should work in serial
Plex: Fixing Pragmatic interface
Plex: Fix adaptation for serial case - Pragmatic need overlapped mesh
Plex: Starting parallel adaptation interface
Plex: Update to parallel Pragmatic interface
Plex: Use comm variable
Plex ex19: Added parallel Pragmatic test
Plex ex19: Added mesh distribution, and fixed parallel metric calculation

Conflicts:
src/dm/impls/plex/examples/tests/output/ex19_0.out
src/dm/impls/plex/examples/tests/output/ex19_1.out
src/dm/impls/plex/examples/tests/output/ex19_2.out
src/dm/impls/plex/examples/tests/output/ex19_3.out

show more ...

Merge branch 'knepley/feature-plex-pragmatic-parallel'

* knepley/feature-plex-pragmatic-parallel:
Plex: Handle unused variable
Plex ex19: Now working in parallel - Added partition visualization

Merge branch 'knepley/feature-plex-pragmatic-parallel'

* knepley/feature-plex-pragmatic-parallel:
Plex: Handle unused variable
Plex ex19: Now working in parallel - Added partition visualization field - Added some tests
Configure: Bump Pragmatic version
Plex: Update interface and fix boundary labels for Pragmatic
Plex ex19: Create mesh overlap before adaptation - Put coords in metric for checking
Plex: Take overlap generation out of Adapt() for now
Plex ex19: Add another Pragmatic test
Plex: Use new parallel Pragmatic interface
Plex ex19: Added a parallel test
Plex ex19: Make metric tensor a local vector
Plex: Overlap distribution should work in serial
Plex: Fixing Pragmatic interface
Plex: Fix adaptation for serial case - Pragmatic need overlapped mesh
Plex: Starting parallel adaptation interface
Plex: Update to parallel Pragmatic interface
Plex: Use comm variable
Plex ex19: Added parallel Pragmatic test
Plex ex19: Added mesh distribution, and fixed parallel metric calculation

show more ...

Merge branch 'master' into knepley/feature-plex-sphere-mesh

* master: (2633 commits)
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables
ini

Merge branch 'master' into knepley/feature-plex-sphere-mesh

* master: (2633 commits)
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables
initialize ap as complained by the nightly tests
Plex ex3: Update to new point function API
initialize some variables as complained by nightly tests
SNES ex12: Mix misnaming Nf --> Nc
update tests
add tests
add MatView_SeqXAIJ_ASCII_structonly(); enable A->structure_only in MatSetValuesBlocked_SeqBAIJ
typo
fill-in some references
Bib: Fixed ref
Plex: Added PLY to file support
Plex: Adding PLY support
PetscDS: Turn off Fortran bindings for Constants
Plex tests: Updated to new point function API
SNES ex56: Moved tests to new format
DM: Allow field projection functions to take an array of scalars - This will allow us to pass in user data to kernels without void *, more efficiency than aux fields
...

Conflicts:
config/builder.py
include/petscdmplex.h
src/dm/impls/plex/examples/tests/ex1.c
src/dm/impls/plex/plexcreate.c
src/dm/impls/plex/plexgmsh.c
src/snes/examples/tutorials/ex63.c

show more ...

DMPlex: allow uniform refinement for periodic meshes without maxCell/L

Needed an API change for DMSet/GetPeriodicity

Plex: Overlap distribution should work in serial

Merge branch 'master' into knepley/feature-plex-pragmatic-parallel

* master: (1379 commits)
Minor fixes for tchem examples due to changes in PETSc infrastructure
update webpage about API changes

Merge branch 'master' into knepley/feature-plex-pragmatic-parallel

* master: (1379 commits)
Minor fixes for tchem examples due to changes in PETSc infrastructure
update webpage about API changes
add help menu for MatCreateSubMatricesMPI()
rename struct Mat_SubMat -> Mat_SubSppt
add a little information to the manual about using tchem
rename MatCreateSubMatrixComposite -> MatCreateSubMatrixVirtual; Mat_SubMatrix -> Mat_SubVirtual
need to protect the \ passed on the makefile on my Mac
Bib: Updated ref
update src/mat/examples/tests/output/ex23_#.out
rename MatGetSubMat -> MatCreateSubMat
Mark MatCreateMPIAIJWithSplitArrays() and MatMPIAIJSetPreallocationCSR() with MAT_NO_OFF_PROC_ENTRIES
superlu: libray is installed in PREFIX/lib64 - fix this to use PREFIX/lib
backport fixes cc105cc3 2b35c90a - required for 3bf428d - that was backported from master - 80c96bb
'echo -n' fails on opensolaris/osx - so avoid using it.
fix diffs in nightlybuilds due to options in ~/.petscrc
add Kelley03 to src/docs/tex/petsc.bib
build: determine moudles, mpi.h used during buildtime and print in 'info'
fix analyzer warnings
fix build failure on opensolaris - where -lm from flibs is removed due to -lm in cxxlibs. But then - only flib gets used for some configure tests [blaslapack, parms]
config: prefer __restrict over C99 restrict keyword
...

show more ...