Add support for using -snes_mf_operator with SNESSetPicard()

Applies Newton's method using matrix-free application of the true Jacobian, preconditioned by user provided Picard matrix

This involved

Add support for using -snes_mf_operator with SNESSetPicard()

Applies Newton's method using matrix-free application of the true Jacobian, preconditioned by user provided Picard matrix

This involved adding (*computemffunction)(SNES,Vec,Vec,void*); to be used by the MatCreateSNESMF() operator when it is set and routines to work with the new function.

The implementation requires maintaining a copy of the matrix that is filled by the users routine that computes the matrix at each matrix-vector product. The original matrix is used for constructing the preconditioner and therefor cannot be overwritten.

Also dropped the requirement that the user provide a right hand side function to SNESSetPicard()

Commit-type: feature, example
/spend 8h
Thanks-to: Saransh Saxena <saransh.saxena5571@gmail.com>

show more ...

Merge remote-tracking branch 'origin/dalcinl/fix-null'

TS/SNES/Tao: Fix 0 ->NULL for pointers

See merge request petsc/petsc!2783

SNES: Replace 0 -> NULL for pointers

Merge branch 'dalcinl/binary-io' into 'master'

Binary IO Enhancements

See merge request petsc/petsc!2594

PetscBinaryWrite: Remove `istemp` boolean flag

Merge branch 'jczhang/fix-dm-originaldm/maint' into maint [PR #1785]

* jczhang/fix-dm-originaldm/maint:
Set DMKSP originaldm when doing copy or getting write access.

Merge branch 'jczhang/fix-dm-originaldm/maint' [PR #1785]

* jczhang/fix-dm-originaldm/maint:
Set DMKSP originaldm when doing copy or getting write access.

Initialize DMSNES originaldm when creating DMSNES

Merge branch 'master' into jpusztay/feature-swarm-symplectic-example

Merged petsc/petsc into master

Merged petsc/petsc into master

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

Merge branch 'master' into knepley/feature-multiple-ds

* master: (156 commits)
Added missing CHKERRQ()
KSP ex56: Added alt output file for ex56_2.
Revert "KSP ex56: Updated test output and dom

Merge branch 'master' into knepley/feature-multiple-ds

* master: (156 commits)
Added missing CHKERRQ()
KSP ex56: Added alt output file for ex56_2.
Revert "KSP ex56: Updated test output and domain size for test _2."
DMStag tutorial ex4: fix stencil errors
Hierarchical partitioner: added an alt output file
Add support for DMGetMatrix() with periodic boundary conditions in one dimension for DMDA.
Fixed incorrect 'ntext' argument in PetscOptionsEList call in SNES ex12.c
ISCreateSubIS: Initialize variables owner and lidx to avoid compiler warnings
metis: provide -download-metis-use-doubleprecision option
KSP ex56: Updated test output and domain size for test _2.
Hierarchical partitioner: add a bit documentation
ISCreateSubIS: Added checkers for parameters
Added a test for hierarchical partitioning vertex weights
Plex+Doc: Document options for DMPlexCreateBoxMesh()
Hierarchical partitioner supports vertex weights
Create an API ISCreateSubIS that extracts
Add an alt output file
Plex: Can now give bounds to DMPlexCreateBoxMesh() from options - Use -dm_plex_box_lower and dm_plex_box_upper
PCMG: These methods should be optional
PCMG: Make PCMGGet/SetLevels() proper member functions - Also need to compose them in GAMG
...

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

show more ...

Merge remote-tracking branch 'origin/pr1206/arcowie-rem/feature-error-logging/master'

DMSNES: Do not copy the data member

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 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 'denera/fix-tablevels-for-views' [PR #869]

* denera/fix-tablevels-for-views:
Instead of manual push-pop tab operations, views now inherit tablevels automatically.
This makes views consi

Merge branch 'denera/fix-tablevels-for-views' [PR #869]

* denera/fix-tablevels-for-views:
Instead of manual push-pop tab operations, views now inherit tablevels automatically.
This makes views consistent with default convergence monitors.

show more ...

All views changed to inherit tab levels from associated objects

Merge commit 'd47bf9aaf1e5266cc3f0ff499f934c85788965a9' into jed/fix-matcreatempibaij/maint

Obtain changed (v3.5) handling of MatMPIBAIJSetPreallocationCSR to match
documentation.

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

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

Merge branch 'barry/remove-funct'

Remove the use and definition of __FUNCT__ throughout the code

Since all modern C/C++ compilers provide this functionality we no longer need to provide it manually in PETSc
Time: 1.5 hours
Thanks-t

Remove the use and definition of __FUNCT__ throughout the code

Since all modern C/C++ compilers provide this functionality we no longer need to provide it manually in PETSc
Time: 1.5 hours
Thanks-to: Andreas Mang <andreas@ices.utexas.edu>

show more ...