Add -ts_monitor_extreme

This is useful for monitoring ODEs where the solutions should not grow (for example advection)

Commit-type: feature

refactor bsi and add the stormer verlet method

wikipedia is wrong about the order of updating p and p

Add the first symplectic integrator SIEuler

The interface is designed particularly for separable Hamiltonian systems

Merge branch 'master' into knepley/fix-plex-project-tests

* master: (477 commits)
Added correct report number for PETSc Developers Guide
Update PETSc Developers Guide title and final pages to ne

Merge branch 'master' into knepley/fix-plex-project-tests

* master: (477 commits)
Added correct report number for PETSc Developers Guide
Update PETSc Developers Guide title and final pages to new ANL style
Hash: Quick documentation fix
PetscBinaryViewerRead -> PetscViewerBinaryRead
cast in GAMG_HASH so that it works for PetscInt32 near PETSC_INT_MAX without overflow to negative numbers
Add link for DrMemory which may serve as an alternative to valgrind on Windows
Finish Gail's copyedits for developers manual.
openmpi: fix namespace conflict [with static build on ppc64le]
two pages of Gail's edits to developers guide
Another chapter of Gail's copyedits for the developers guide
Gail's copyedits on two more chapters of the developers guide
Fix more %d to %D in GAMG code. Fix code style in GAMG code
DMPlexVTK: vertices is int! need to use %d instead of %D
TS ex18: update outputs (reported by nightly tests)
Use %D in error messages for GAMG hash; fix style issues with GAMG hash
config: Update petsc4py git commit
Fix bad plain extern declaration.
Add Gail's copyedits for the initial portion of the developers guide
Add missing common block to subroutine
fix mpiuni tests with xl compilers
...

show more ...

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

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 remote-tracking branch 'origin/jed/const-XType'

Remove "const MatType" and similar (const is already inside the typedef)

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 'psanan/manpages-notes-newlines' [PR #932]

* psanan/manpages-notes-newlines:
A large number of man pages are improperly formatted because of the lack of a newline after "Notes:", "Fortr

Merge branch 'psanan/manpages-notes-newlines' [PR #932]

* psanan/manpages-notes-newlines:
A large number of man pages are improperly formatted because of the lack of a newline after "Notes:", "Fortran Notes:", or "Developer Notes:".
This batch adds those newlines, with 4 spaces of indentation afterwards.
The uniform indentation arguably makes the man pages a little less readable in the source itself,
but overall this hopefully has a large positive impact on the readability of the HTML pages.

show more ...

Man pages: add newlines after "Notes:"

This allows for proper formatting from sowing.

On OS X (using gsed, not the default BSD sed), from the PETSc root directory:

find src include -type f \(

Man pages: add newlines after "Notes:"

This allows for proper formatting from sowing.

On OS X (using gsed, not the default BSD sed), from the PETSc root directory:

find src include -type f \( -name "*.c" -or -name "*.h" -or -name "*.cxx" \) | xargs gsed -i 's/Notes\s*:\s*\(\w.*\)/Notes:\n \1/'

This adds a newline and 4 spaces whenever "Notes:" is followed by any "word" character, in any .c, .h, or .cxx file in src/ or include/

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 remote-tracking branch 'origin/knepley/fix-snes-ex69' into tisaac/feature-dmfield

* origin/knepley/fix-snes-ex69: (311 commits)
PetscFECreateDefault: fix custom Fortran interface to use MPI_

Merge remote-tracking branch 'origin/knepley/fix-snes-ex69' into tisaac/feature-dmfield

* origin/knepley/fix-snes-ex69: (311 commits)
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
cleanup very ugly PCGAMG manual page
minor clarifications for VecScatterCreate() manual page
gcov: PETSC_ARCH is required onlyfor -run_gcov - and not for -merge_gcov [where it doesn't exist]
make clearer in help statements what --with-external-packages-dir is for.
Website: Added news entry for PETSc User Meeting 2018 on main page.
CUDA: Further test output fixes for mat_tests-ex2
Turn off elemental test for parmetis running since this does not work on batch systems.
3.8 Changes: KSPCHEBYSHEV
CUDA: Fixed most test errors due to MatMultTranspose*
website: update urls
website: fix Jason's name to match pic. Also add Todd's pic
Dev manual: copyedit
DMPlexDistributeLabels: Propagate the output flag
cuda/cusp: update minimum supported versions of cuda to 7.5 and cusp to 0.5.0
...

show more ...

Merge branch 'hongzh/update-ts-adjoint-api' [PR #870]

* hongzh/update-ts-adjoint-api:
Reorganizes TSAdjoint code.
Deprecates TSAdjointSetRHSJacobian() TSAdjointComputeDRDP() TSAdjointComputeDRDY()

Merge branch 'hongzh/change-vecarray-to-mat-tlm'

* hongzh/change-vecarray-to-mat-tlm:
Uses a dense array instead of multiple vectors to store TLM sensitivity
and removes TSForwardSetRHSJacobianP (us

Merge branch 'hongzh/change-vecarray-to-mat-tlm'

* hongzh/change-vecarray-to-mat-tlm:
Uses a dense array instead of multiple vectors to store TLM sensitivity
and removes TSForwardSetRHSJacobianP (use TSAdjointSetRHSJacobian instead).

show more ...

Reorganize TSAdjoint and deprecate some functions

The derivative routines are independent of integration methods, thus can be shared by adjoint and TLM.

Funded-by:
Project: FASTMATH
Time: 1.0 hours

Reorganize TSAdjoint and deprecate some functions

The derivative routines are independent of integration methods, thus can be shared by adjoint and TLM.

Funded-by:
Project: FASTMATH
Time: 1.0 hours
Reported-by:
Thanks-to:

show more ...

Merge branch 'barry/feature-ts-burgers-spectral'

Add code for automatically checking if shell matrix vector and transpose matrix vector products work with TS

Added MatShellTestMult(), MatShellTestMultTranspose(), TSRHSJacobianTest(), and TSRHSJaco

Add code for automatically checking if shell matrix vector and transpose matrix vector products work with TS

Added MatShellTestMult(), MatShellTestMultTranspose(), TSRHSJacobianTest(), and TSRHSJacobianTestTranspose()

Commit-type: feature

show more ...

-ts_monitor_error can now draw error as well as print the 2-norm

Commit-type: feature

Refactor the TLM model

- The forward sensitivities are now stored in dense matrices. ex20fwd is modified to test the change.
- TSForwardSetRHSJacobianP has been removed. TSForward uses the JacobianP

Refactor the TLM model

- The forward sensitivities are now stored in dense matrices. ex20fwd is modified to test the change.
- TSForwardSetRHSJacobianP has been removed. TSForward uses the JacobianP matrix that one provides via the TSAdjoint interface.

Funded-by:
Project: FASTMath
Time: 8.0 hours
Reported-by:
Thanks-to:

show more ...

Merge branch 'barry/correct-radau5'

Start implementation of Radau5 interface

Reproduce Radau5 result for simple test case

Commit-type: feature

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

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