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

Not sure if this has done the right thing because there are no tests
yet. Also, the fortran stubs break (it doesn't know what to do with
SNESCo

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

Not sure if this has done the right thing because there are no tests
yet. Also, the fortran stubs break (it doesn't know what to do with
SNESCompositeType), and I don't know how to fix it.

show more ...

Merged petsc/petsc into master

Merged petsc/petsc into master

Merge branch 'stefano_zampini/cherrypick-from-continuous-adjoint' [PR #1250]

* stefano_zampini/cherrypick-from-continuous-adjoint:
Few improvements from the continuous adjoint branch
This branch ext

Merge branch 'stefano_zampini/cherrypick-from-continuous-adjoint' [PR #1250]

* stefano_zampini/cherrypick-from-continuous-adjoint:
Few improvements from the continuous adjoint branch
This branch extracts few things that are in the continuous adjoint branch that should go into PETSc

show more ...

TS: added TSHISTORY object that allows to inquire a TSTrajectory by time (as key) and not by time step number only

Merge from master

Merged petsc/petsc into master

Merge branch 'master' into arcowie-rem/feature-error-logging

* master: (345 commits)
Update dev notes for the coming release
Increase patchlevel to 3.9.4
doc fixes
Added some AIJSELL tests t

Merge branch 'master' into arcowie-rem/feature-error-logging

* master: (345 commits)
Update dev notes for the coming release
Increase patchlevel to 3.9.4
doc fixes
Added some AIJSELL tests that mirror the SELL tests.
Bib: More refs
Removed use of hardwired -1 to represent PETSC_NULL_XXX object, since the null value is now 0.
fix memleaks flagged by valgrind
ex62: disable long running test and workaround differences in iteration counts on various builds.
Add FAQ we can point to when users get error messages such as: "PETSc was configured with one OpenMPI mpi.h version but now appears to be compiling using a different OpenMPI mpi.h version"
swith test from mpich to openmpi as mpich ch3:sock is unuserable per PetscShmCommGet()
DM: Forgot to move include statement
fix nightlybuild issues - disable hypre [for openmp] and fix uninitialized variables
tests: enable openmp and hwloc with one of the mumps tests
Default BoomerAMG smoothers
Configure: Revert changes from debb90c, so that we are back to the original Libmesh
Configure: Reverted change for Libmesh from 83cd5d7
Split vecimpl.h into vecimpl.h and vecscatterimpl.h
SNES ex62: Turned off some tests for other precisions and shortened output of sensitive tests
Use mmap to allocate the shared pthread barrier if available
Add configure test for functional mmap
...

Conflicts:
include/petscconvest.h
src/snes/examples/tutorials/ex13.c
src/snes/interface/snes.c
src/snes/utils/convest.c

show more ...

Merge branch 'master' into knepley/fix-superdm-is

* master: (123 commits)
Fix a few more Fortran objects to have their NULL value match the C NULL value
remove dead code
Add variable point blo

Merge branch 'master' into knepley/fix-superdm-is

* master: (123 commits)
Fix a few more Fortran objects to have their NULL value match the C NULL value
remove dead code
Add variable point block Jacobi for MPIAIJ matrix.
should not pass user contexts through CHKFORTRANNULLOBJECT in Fortran stubs
fix mpiuni test
Remove CHKFORTRANNULLOBJECTDEREFERENCE since it is no longer needed with Fortran null equal to C null
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Basic as external
Change Fortran null pointer to match C null pointer
one more fix
jenkins: run 'cleantest' after examples are run - to recover diskspace
logical order bug in diagbrdn update fixed, tests updated again
replace new -petscspace_order instances with -petscspace_degree
PETSCSPACETENSOR: clear up naming confusion of subspaces
regression.py: remove configure code related to old test suite - now that SLEPc no longer needs it
PetscSpaceSetup_Poly(): make idempotent
Plex ex3: update tests
Lift changes to defunct dtfe.c from [PR #1074]
PetscSpace_Tensor: clarification of vector-valued
PetscSpace_Tensor: allow height subspaces, fixes tabulation initialization
...

show more ...

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

* master: (131 commits)
fix mpiuni test
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Ba

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

* master: (131 commits)
fix mpiuni test
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Basic as external
one more fix
jenkins: run 'cleantest' after examples are run - to recover diskspace
logical order bug in diagbrdn update fixed, tests updated again
replace new -petscspace_order instances with -petscspace_degree
PETSCSPACETENSOR: clear up naming confusion of subspaces
regression.py: remove configure code related to old test suite - now that SLEPc no longer needs it
PetscSpaceSetup_Poly(): make idempotent
Plex ex3: update tests
Lift changes to defunct dtfe.c from [PR #1074]
PetscSpace_Tensor: clarification of vector-valued
PetscSpace_Tensor: allow height subspaces, fixes tabulation initialization
PetscSpace: add a dimension member
PetscSpaceSetDegree(): now take maxDegree argument
-petscspacepoly_tensor -> -petscspace_type tensor
PetscSpaceEvaluate_Polynomial(): fix temp array dependencies
PETSCSPACETENSOR: A generic tensor product of spaces.
PetscSpace: -petscspace_order -> -petscspace_degree
...

Conflicts:
src/dm/dt/interface/dtfe.c
src/snes/examples/tutorials/ex13.c
src/snes/examples/tutorials/ex17.c
src/snes/examples/tutorials/ex62.c
src/snes/examples/tutorials/output/ex62_bddc_quad.out
src/snes/examples/tutorials/output/ex62_fetidp_2d_quad.out
src/snes/examples/tutorials/output/ex62_fetidp_2d_tri.out
src/snes/examples/tutorials/output/ex62_fetidp_3d_hex.out
src/snes/examples/tutorials/output/ex62_fetidp_3d_tet.out

show more ...

Merged petsc/petsc into master

Merge branch 'hongzh/symplectic-integrator' [PR #1049]

* hongzh/symplectic-integrator:
This branch implements symplectic time integrators.

We should come up with a unifying framework that works for

Merge branch 'hongzh/symplectic-integrator' [PR #1049]

* hongzh/symplectic-integrator:
This branch implements symplectic time integrators.

We should come up with a unifying framework that works for all multi- methods such as multirate methods (see PR #1043). The idea of this branch is to ask users to provide a two-way splitting of the separable Hamiltonian system.

TSGetRHSFunctionSplit2w(TS ts,Vec *r,TSRHSFunctionSplit2w *fun1,TSRHSFunctionSplit2w *fun2,void **ctx)
TSSetRHSFunctionSplit2w(TS ts,Vec r,TSRHSFunctionSplit2w fun1,TSRHSFunctionSplit2w fun2,void *ctx)

The split/partitioned functions typically have an output vector that has lower DOF. For example, RHSFunction() represents a mapping Rn → Rn , but RHSFunctionSplit2w() provides a mapping Rn → Rn/2 since the subvector is always half of the combined vector. This is also why I use nested vectors for simplicity.

Using IS is a more general way to support multi-way splitting. I am not sure if we should use IS instead of nested vectors for symplectic integrators. Perhaps we can keep both approaches.

show more ...

Allow RHSSplit to return sub TS

Funded-by:
Project: PETSc time integration
Time: 6.0 hours
Reported-by:
Thanks-to:

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

* master: (984 commits)
netcdf: set zlib as optional dependency corretly [wrt --enable-netcdf-4] - otherwise netcdf configure errors out

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

* master: (984 commits)
netcdf: set zlib as optional dependency corretly [wrt --enable-netcdf-4] - otherwise netcdf configure errors out
fix gcov obj path
-options_table was changed to options_view in 2014 but documentation was not updated to reflect it
configure: fix F2003 test
Dev manual: remove duplicate entry
Plex test ex5: Removed "future" test, as it resulted in issues in the nightly tests.
Some minor typos in TS
Add ADIOS2.py which I forgot to git add in the past
Plex ex5: Disable test
Doc: Added note for interface change
KSPLSQRMonitorDefault fixed for quad precision
Tao: Added alternative output for chwirut2 test.
TS ex45,ex46: Removed overly refined tests for faster nightly tests.
KSP ex72: Fixed needlessly long-running PCTFS tests.
SNES ex56: Reduced test size to speed up nightly test suite.
KSPSolve_LSQR: don't compute unused norm of rhs
SNES ex13: Reduced problem sizes for nightly tests.
PetscOptionsDeprecated(): fix for NULL prefix
PetscOptionsDeprecated(): add option to print additional info
SNES ex12: Reduced sizes for BDDC tests.
...

show more ...

Merge branch 'master' into knepley/feature-plex-elasticity-test

* master: (1553 commits)
-options_table was changed to options_view in 2014 but documentation was not updated to reflect it
config

Merge branch 'master' into knepley/feature-plex-elasticity-test

* master: (1553 commits)
-options_table was changed to options_view in 2014 but documentation was not updated to reflect it
configure: fix F2003 test
Plex test ex5: Removed "future" test, as it resulted in issues in the nightly tests.
Some minor typos in TS
Add ADIOS2.py which I forgot to git add in the past
Plex ex5: Disable test
Doc: Added note for interface change
KSPLSQRMonitorDefault fixed for quad precision
Tao: Added alternative output for chwirut2 test.
TS ex45,ex46: Removed overly refined tests for faster nightly tests.
KSP ex72: Fixed needlessly long-running PCTFS tests.
SNES ex56: Reduced test size to speed up nightly test suite.
KSPSolve_LSQR: don't compute unused norm of rhs
SNES ex13: Reduced problem sizes for nightly tests.
PetscOptionsDeprecated(): fix for NULL prefix
PetscOptionsDeprecated(): add option to print additional info
SNES ex12: Reduced sizes for BDDC tests.
PetscFileRetrieve(): Suppress output to STDERR when downloading a file.
MatMatMult: Fixed unused variable pN.
DMDestroy: Guard call to dm->ops->destroy * Similar guards are in place in many other methods (IS,Vec,Mat) and permit the DM to be destroyed even if DMSetType() was never called.
...

show more ...

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

* master: (2931 commits)
remove N in docs for SNESQN
replaced test tolerance with PETSC_SQRT_MACHINE_EPSILON
fortran: PetscComp

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

* master: (2931 commits)
remove N in docs for SNESQN
replaced test tolerance with PETSC_SQRT_MACHINE_EPSILON
fortran: PetscComplex type is missing from fortran
initialize variables so that static analyzer doesn't mistakenly think they are used uninitialized.
Move access of prefix until it is actually used
Remove DMDASplitComm2d(), used in one example that was not tested
Functions that do not return PetscErrorCode should not use PetscFunctionBegin/End paradgm
remove unneeded reassignment of same value
Missing CHKERRQ(); Mostly in Tao code
Remove clang static analyser warning about passing NULL variable to strncat()
missing newline in -ts_monitor_extreme
minor fixes for integer counters
Configure: Deal with '__func__' in C89 builds with newer Clang and GCC.
Jenkinsfile: Increased timeout for tests to 5 minutes.
3.9 changes: further clarify removal of -pc_factor_mat_solver_package
Options deprecations: add "PETSc" to inline warning
removing unused variable
MatShift support added to MATSUBMAT. Mimics MATSHELL functionality. Needed for matrix-free Hessians in TAO Newton-Line-Search methods.
Add -ts_monitor_extreme
Dev manual: add note on deprecating options database keys
...

Conflicts:
src/dm/impls/plex/examples/tests/ex4.c

show more ...

Merge branch 'master' into arcowie-rem/feature-error-logging

* master: (353 commits)
Configure: Pragmatic now requires Cxx11
fix compile warning
mpiuni: build somempifort.o after mpiuni.mod

Merge branch 'master' into arcowie-rem/feature-error-logging

* master: (353 commits)
Configure: Pragmatic now requires Cxx11
fix compile warning
mpiuni: build somempifort.o after mpiuni.mod
Remove use of common block from MPIUNI code
use was in wrong location
Fortran line extends past 72 char
Moved the rest of the PETSc common block variables over to petscsys module
tests: fix breakage due to cd93b0e13c6da70ee1e038f55989cf888bc1d1c5
Need to use petscsys to get PETSC_COMM_SELF defined in example
Mat ex23 and ex40: fix errors reported from nightly tests
Put declaration on too lines because otherwise it generates an error for truncated line
Move PETSC_COMM_WORLD and PETSC_COMM_SELF from common blocks to petscsys module
tests: fix requires - and create out files for seqaijmkl tests
KSP ex72: nested dissection is likely to be different between machines
cmake: downgrade to version 3.9 - as newer versions require cxx11 compiler [which is not universal]. However stick with version 3.11 on opensolaris
Mat ex23: fix memory leaks
Bib: Some new refs
KSP examples: add missing casts to double in printf
updated GAMG KSP tests
PetscRound() should be PetscRoundReal() in __fp16
...

Conflicts:
include/petsclog.h
src/sys/logging/plog.c
src/sys/logging/utils/eventlog.c

show more ...

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 remote-tracking branch 'origin/hongzh/separate-tssen-from-ts'

Move tssen.c into a subsection of TS

Sensitivity related functions should appear in a separate page than TS

Funded-by:
Project: FASTMATH adjoint
Time: 0.5 hours
Reported-by:
Thanks-to:

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()

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