Plex: added DMPlexGetRedundantDM

Merge branch 'master' into tisaac/dmforest

Get 3.6.2 patches

* master: (1360 commits)
remove nasty memory leak with CUDA code when using block Jacobi with multiple blocks per process
print abs

Merge branch 'master' into tisaac/dmforest

Get 3.6.2 patches

* master: (1360 commits)
remove nasty memory leak with CUDA code when using block Jacobi with multiple blocks per process
print abs in Fortran test example in order to have same result with real and complex builds
changed outdate mpirun and -np to mpiexec and -n
force sqrt to be taken in the precision of PetscReal in fortran examples
continuation line too soon for some compilers
stop configure if user selected C++ clanguage, single precision, complex combination since we do not support it
fixes for Fortran examples to compile cleanly in single precision that do not involve complex numbers
fixes to Fortran examples to compile cleanly in single precision complex
fix previous mistake
must cast 1. to PetscScalar before dividing by a PetscScalar (C++ complex class is terrible)
must cast 0.0 to PetscScalar before comparison with a PetscScalar
changed order of suggested locations of downloadable PETSc to Barry's preference
Increase patchlevel to 3.6.2
remove debug print statement from two commits before
remove debug print statements I accidently left in the last commit
PetscDraw: Fix dummy definition of PetscDrawCollective{Begin|End}()
remove cleanupMumps flag from MatDestroy_Mumps() it should always free the memory it allocated Bug detected while adding features to libMesh
improved viewing of fractional counts in xml nested log viewing fixed bug in output of doubles to xml nested log files
chombo: build only libs - not examples
fix typo. 'gcc-4.8 --with-cxx-dialect=C++11' gives a warning with this line - even though its in a #ifdef thats never enabled.
...

show more ...

Merge branch 'knepley/feature-label-empty-stratum'

* knepley/feature-label-empty-stratum:
Plex ex11: Added code to test empty strata in DMLabel
Plex: Check the depth label state - If the depth l

Merge branch 'knepley/feature-label-empty-stratum'

* knepley/feature-label-empty-stratum:
Plex ex11: Added code to test empty strata in DMLabel
Plex: Check the depth label state - If the depth label is user managed, then send it when distributing labels
DMLabel: Now has a state - Why are these not PetscObjects?
Plex: If we remove the depth label, clear the optimized access
DMLabel: Check for empty stratum
DMLabel: Make DMLabelAddStratum() public

show more ...

To refactor the comments to follow petsc style.

Merge remote-tracking branch 'petsc/master' into fixed-comment-space

Conflicts:
include/petsc/finclude/petscmat.h
src/mat/examples/tests/makefile

To refactor the comments to follow petsc style.

Merge remote-tracking branch 'petsc/master' into fixed-comment-space

Conflicts:
include/petsc/finclude/petscmat.h
src/mat/examples/tests/makefile
src/mat/impls/aij/mpi/mpiov.c

show more ...

To add a space to all comments I have made before. Barry pointed out
that a comment without a space right after '/*' may cause bfort to crash

Merge remote-tracking branch 'petsc/master' into fixed-c

To add a space to all comments I have made before. Barry pointed out
that a comment without a space right after '/*' may cause bfort to crash

Merge remote-tracking branch 'petsc/master' into fixed-comment-space

Conflicts:
include/petsc/finclude/petscmat.h
include/petscmat.h
src/mat/examples/tests/makefile
src/mat/impls/aij/mpi/mpiaij.h

show more ...

Plex: Check the depth label state
- If the depth label is user managed, then send it when distributing labels

Merge branch 'pr329/master/Fande-Kong/matpartitioning-hierarch' into pr361/master/Fande-Kong/pcgasm-increaseoverlap: need hierarch to build GASM subdomains.

Conflicts:
src/mat/examples/tests/makefi

Merge branch 'pr329/master/Fande-Kong/matpartitioning-hierarch' into pr361/master/Fande-Kong/pcgasm-increaseoverlap: need hierarch to build GASM subdomains.

Conflicts:
src/mat/examples/tests/makefile
src/mat/impls/aij/mpi/mpiaij.h

show more ...

Merge branch 'master' into pr352/dmay/pctelescope/master

Merge branch 'master' into pr358/buptphyphy/addmatcreatenormalhermitian/master

Merge branch 'tisaac/plex-feature-tree-fv'

Make the PetscFV gradient reconstruction method work for the
non-conforming meshes introduced by the Tree extension to DMPlex.

* tisaac/plex-feature-tree-

Merge branch 'tisaac/plex-feature-tree-fv'

Make the PetscFV gradient reconstruction method work for the
non-conforming meshes introduced by the Tree extension to DMPlex.

* tisaac/plex-feature-tree-fv: (52 commits)
Plex ex3: make test that works for single precision
PlexSNES: address compiler warning
DMPlex ex3: make test output PASS/FAIL rather than precision specific
Plex: fix type punning problem stuffing (PetscFECellGeom *) in a Vec
Plex: allow ex3 testFVgrad examples in complex
PetscFV: allow SVD pseudoinverse least squares gradient calculation for complex.
Plex ex3: put fv tensor examples in C_NoComplex tests
Plex ex3: update regression output for nonconforming examples
Plex ex3: make separate fe and fv tests
Plex ex3: use MPIU_MAX, not MPI_MAX
DMPlex: added DMPlexReconstructGradientsFVM()
Plex: fix typo in change to cohesive cell creation
Plex: add ex4 and ex5 to makefile so they can be used for regression tests
DMPlex: remove overzealous error condition from FV gradient reconstruction
DMPlex: correct pointSF size in new point shift approach
DMPlexMigrate(): use correct global indices when migrating tree
DMPlexMigrate(): Don't invalidate old DM by overwriting cones
Plex ex3: use DMPlexDistributeOverlap()
DMPlexPartitionLabelClosure: traverse the tree up and down
PetscPartitionerSimple: work in parallel
...

show more ...

Merge branch 'tisaac/dt-fix-geom-punning' into tisaac/plex-feature-tree-fv

Pull in a fix to FE geometry that affects the complex,64bit case

* tisaac/dt-fix-geom-punning: (240 commits)
Plex: fix t

Merge branch 'tisaac/dt-fix-geom-punning' into tisaac/plex-feature-tree-fv

Pull in a fix to FE geometry that affects the complex,64bit case

* tisaac/dt-fix-geom-punning: (240 commits)
Plex: fix type punning problem stuffing (PetscFECellGeom *) in a Vec
PGI compiler cannot handle [restrict] usage therefor disable restrict for PGI compiler Fix ex48.c to use PETSc standard way of marking restrict
PETSC_UNUSED must occur before the typename. With PGI compilers they error if it is after the type name
on solaris in TSAdaptChoose_Basic() negative root for zero could produce -infinity thus causing time step to be decreased even though the error estimate was zero
MatGetSchurComplement: revert test from LUMP back to DIAG and update output
there are no tests for opencl or viennacl so prevent warning message: make[2]: *** No rule to make target `testexamples_OPENCL', needed by `tree'. Stop. make[2]: *** No rule to make target `testexamples_VIENNACL', needed by `tree'. Stop.
incorrect initial conditions function that could produce floating point exception value (and does on windows) due to squareroot of non-negative number
forgot to fix rule for building ex194
move ex193.c to ex194 because someone else created a 193
remove outdate manual page comment
add missing VecRestoreArray()
add missing VecRestoreArray()
missing restore array in example
formating for print statements of test example with __float128 precision
fix example for __float128 precision by using PETSc math macros and casting PetscReal in print statements
remove Tao example custom monitor that produced different results for __float128 precision
all PETSc built in monitors require a valid PetscViewer, NULL is not acceptable
PETSc built in monitors always require viewers to be passed in, no longer support a NULL viewer
use tigher haptol for minres for quad precision because otherwise false positive of indefinite matrix appeared in example
check for PetscMPIInt overflow in ISAllGather()
...

show more ...

Merge branch 'master' into jed/process-options-handler

Merge branch 'tisaac/dmplex-stable-tet-refinement' into tisaac/plex-feature-tree-fv

This feature branch also changed the topology of nonconforming examples:
I am merging so that I can make the regre

Merge branch 'tisaac/dmplex-stable-tet-refinement' into tisaac/plex-feature-tree-fv

This feature branch also changed the topology of nonconforming examples:
I am merging so that I can make the regression output that should be
correct in master and next

* tisaac/dmplex-stable-tet-refinement: (594 commits)
plex: added tet shape regression test
plexrefine: remove replaced code
plex: update test output for new tet refinement order
plex: fix tet refinement for stable element shape
DMPlex: added test of shape quality to ex1
Better to use the proper error handling version of PetscViewerASCIIGetStdout() instead of PETSC_VIEWER_STDOUT_()
Generate a very useful error message if PetscLogView(), PetscLogEventGetPerfInfo() or PetscLogEventGetFlops() is called with a preceeding call to PetscLogBegin() and friends (also -log_summary which triggers a PetscLogBegin())
fix compile warnings
Plex: Fix bug with bcComps initialization
Turn off elemental shared libraries if --with-shared-libraries=0 is used
Bib: Added refs
Bib: Fixed a bunch of errors
add warning/error flag for '-mat_superlu_dist_parsymbfact' when MatInputMode=GLOBAL
makefile: fix format
tests: add back ex6f, runex6f targets - that got removed in 2f1d3d565a6a43067e6df75b043e6343b05fcd1d
accessing option that may never have been set
typo in DMNetworkSetNumVariables() reported by Sean Shahkarami
superlu_dist: update to using latest 4.1 version
fix compile warnings
Generate an error message when trying to overwrite default options of CN and BEULER
...

show more ...

DMPlexMigrate(): use correct global indices when migrating tree

DMPlexMigrate(): Don't invalidate old DM by overwriting cones

Merging with Barry's refactor of PetscViewerGetSubComm

Merge branch 'barry/fix-petscviewer-attempt-2' into dmay/pcsemiredundant

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of git. I think it is not justified to ever have anything in next for more than a few weeks at most.
Either take it out of next if it is broken or put it in master if it is not broken. Hanging around in next but
not master for ever is not a good policy.

show more ...

Replaced PetscViewerASCIISynchronizedAllow() with PetscViewerASCIIPushSynchronized() PetscViewerASCIIPopSynchronized()

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

Merge branch 'master' into knepley/feature-snes-deflation

* master: (206 commits)
Bib: Added ref
fix the output file for the ts example ex20adj
Bugfix: setting post-step function in ztsf.c
f

Merge branch 'master' into knepley/feature-snes-deflation

* master: (206 commits)
Bib: Added ref
fix the output file for the ts example ex20adj
Bugfix: setting post-step function in ztsf.c
fix makefile format error with 'all-legacy'
doc: handle missing dev.html in a release tarball
doc: add in petscistypes.h to doc sourcefile list
Increase patchlevel to 3.6.1
update some example output files
Doc+Vec: Fixed VecSetValuesBlocked() Reported by: Shao-Ching Huang <huangsc@gmail.com>
clarify that PETSc's "thread safe" code only works with OpenMP
allow adjoint integration over split time intervals
DT: Goddman it. Why are people stepping on my commits
Bib: Added ref
fix the bug in the endpoint version of discrete adjoint theta and correct the affected example outputs
SNES: Added SNESQN types to Fortran Reported-by: Chris Bradley <c.bradley@auckland.ac.nz>
fix: buggy change in 33da457b864643cfcac2768e47197c286471a5a4
updated petscthreads.html with information on thread per object model
fix for ex56
removed the dead links to threads.html
send petsc repository messages to Karl no longer from his mcs address
...

Conflicts:
config/builder.py
src/dm/impls/plex/examples/tests/ex3.c
src/dm/impls/plex/examples/tests/makefile
src/dm/impls/plex/makefile

show more ...

Merge branch 'master' into knepley/feature-plex-point-location

* master: (202 commits)
Bib: Added ref
fix makefile format error with 'all-legacy'
Increase patchlevel to 3.6.1
update some exa

Merge branch 'master' into knepley/feature-plex-point-location

* master: (202 commits)
Bib: Added ref
fix makefile format error with 'all-legacy'
Increase patchlevel to 3.6.1
update some example output files
Doc+Vec: Fixed VecSetValuesBlocked() Reported by: Shao-Ching Huang <huangsc@gmail.com>
clarify that PETSc's "thread safe" code only works with OpenMP
allow adjoint integration over split time intervals
DT: Goddman it. Why are people stepping on my commits
Bib: Added ref
fix the bug in the endpoint version of discrete adjoint theta and correct the affected example outputs
SNES: Added SNESQN types to Fortran Reported-by: Chris Bradley <c.bradley@auckland.ac.nz>
fix: buggy change in 33da457b864643cfcac2768e47197c286471a5a4
updated petscthreads.html with information on thread per object model
fix for ex56
removed the dead links to threads.html
send petsc repository messages to Karl no longer from his mcs address
Plex: Fix for all the hidebound morons living with c89. Good luck doing computational science in the 80s!
DT: Fix tanh-sinh - Removed wrong var
SNES: Added SNESQN types to Fortran Reported-by: Chris Bradley <c.bradley@auckland.ac.nz>
Bib: More refs
...

Conflicts:
config/builder.py

show more ...

Merge branch 'mesgarnejad/feature-plex-globaltonatural'

* mesgarnejad/feature-plex-globaltonatural:
Plex ex15: This requires HDF5
Plex ex15: Fixed variable type
Plex: Fixed variable init
Ple

Merge branch 'mesgarnejad/feature-plex-globaltonatural'

* mesgarnejad/feature-plex-globaltonatural:
Plex ex15: This requires HDF5
Plex ex15: Fixed variable type
Plex: Fixed variable init
Plex ex15: Added test for reading and writing HDF5 in native format
Vec: Straighten out defines for VECOP
DM: Reminder to handle g2n in DMCreateSubDM()
Plex ex15: One test of g2n working
Plex: Fixed global to natural mapping creation
Plex: When creating g2n map, must use parallel DM
PetscSF: Corrected bug in PetscSFCreateInverseSF() for sparse leaves - Have to index over entire leaf space
PetscSection: PetscSectionCreateGlobalSection() now allows the user to keep local offsets - This is useful is all you want is to filter out the constrained points
Plex: Create the global to natural mapping during distribution if the flag is set
Plex: Added viewing in native format using the global to natural mapping
Plex: Added global to natural ordering support - The natural ordering is the one you have before distribution - Added DMPlexCreateGlobalToNaturalSF() - Added DMPlexGlobalToNaturalBegin/End() - Added DMPlexNaturaToGloballBegin/End()
SF: Added PetscSFCreateEmbeddedLeafSF() - This name sucks, but it subsets on leaves instead of roots
DM: Added an SF for the global to local mapping and a flag to enable its construction - Added DMGet/SetUseNatural() - Should hook this up to DMDA
SF: Added Fortran binding for SFView()
DM: Automatic Fortran binding should work for DMCreateSubDM()
SF: Added PetscSFCompose()
SF: Added PetscSFCreateFromZero()

show more ...

Plex: When creating g2n map, must use parallel DM

Plex: Create the global to natural mapping during distribution if the flag is set

Merge branch 'master' into knepley/solkx

* master: (6933 commits)
Bib: Added reference
Includegraph: Updated for compatibility with new folder layout.
Add SNESLineSearchReason to fortran inclu

Merge branch 'master' into knepley/solkx

* master: (6933 commits)
Bib: Added reference
Includegraph: Updated for compatibility with new folder layout.
Add SNESLineSearchReason to fortran includes Also truncate names to fit 32 char fortran limit
fix compile warnings /usr/home/balay/petsc.clone-2/src/vec/vec/interface/vector.c:1944: warning: division by zero in '1.0e+0 / 0.' /usr/home/balay/petsc.clone-2/src/ksp/ksp/interface/itfunc.c:508: warning: 'vec_rhs' may be used uninitialized in this function
updated output for new trust region initial size
initial size of the trust region is set as a percent of the norm of the initial guess, it should not be related to the initial norm of the function (which is kind of nuts).
The default bounds for SNESVISetVariableBounds() in the manual page were reversed from correct values
Bib: Added refs
Bib: Update MPICH webpage
Plex: Forgot to close file
removed nonexistent ex46f from makefile test rule
Revert "fixed bad merge into master"
fixed bad merge into master
fix memory leak in mkl_pardiso fix some formatting in mkl_pardiso code
mv runex111 from TESTEXAMPLES_C to TESTEXAMPLES_DATAFILESPATH
Mat+Doc: More Fortran docs
Plex ex3: Run tests with the correct number of field components
Bib: Added TetGen cite
parmetis: make sure 'ldd libparmetis.so' points to the correct libmetis.so that its linked with.
small fix as reported by the nightly test
...

Conflicts:
config/builder.py

show more ...