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 branch 'master' into knepley/feature-pc-patch

* master: (440 commits)
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html
release: set v3.9 strings
fix grep us

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

* master: (440 commits)
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html
release: set v3.9 strings
fix grep usage on solaris [use egrep]
opensolaris: fix errors
Update timer reporting and minor reporting fix
Remove redundant if statement. Found by @jedbrown
Fixed problem with order of developer photos on website.
superlu_dist: use correct flag in example [PETSC_HAVE_SUPERLU -> PETSC_HAVE_SUPERLU_DIST]
Bib: PERMON ref
fix for osx
doc: more author fixes
doc: more author fixes
Italicized "orderings" in MPI process placement discussion.
Added Richard to the user manual author list.
Added Richard to the list of active developers on the website.
add string check for MS-MPI and IMPI timeout errors
opensolaris does not have multiple grep utils - but default one does not support -F. So use fgrep instead
Standardizing timeout checks
Report wall clock time of total tests
Use time instead of date %s
...

Conflicts:
include/petscpc.h
src/ksp/ksp/interface/itfunc.c

show more ...

Merge branch 'karlrupp/feature-remove-cusp'

* karlrupp/feature-remove-cusp:
Removes VECCUSP vectors and AIJCUSP matrices, because they provide no added benefit over VECCUDA and AIJCUSPARSE.
Also rem

Merge branch 'karlrupp/feature-remove-cusp'

* karlrupp/feature-remove-cusp:
Removes VECCUSP vectors and AIJCUSP matrices, because they provide no added benefit over VECCUDA and AIJCUSPARSE.
Also removes the CUSP-based preconditioners for now.
SACUSP will be re-added later, after the new release is out.

show more ...

CUSP: Removed VECCUSP and AIJCUSP as well as preconditioners.

SA-AMG from CUSP will be re-added in a follow-up commit,
made to work with AIJCUSPARSE and VECCUDA.

Merge branch 'rmills/fix-rename-csrperm-files'

Renamed source files from "csrperm" to "aijperm" to be consistent with AIJPERM matrix class name.

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

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

Merge branch 'rmills/add-aijmkl'

* rmills/add-aijmkl:
Adds AIJMKL matrix class to facilitate use of MKL sparse BLAS routines
BitBucket Pull Request 737
This is the initial version of the AIJMKL matr

Merge branch 'rmills/add-aijmkl'

* rmills/add-aijmkl:
Adds AIJMKL matrix class to facilitate use of MKL sparse BLAS routines
BitBucket Pull Request 737
This is the initial version of the AIJMKL matrix class,
a subclass of AIJ that has the goal of using MKL sparse BLAS routines whenever practicable.
If the MKL version is new enough, the sparse inspector-executor
("SpMV2" -- terrible name, but that's what the MKL team called it) routines are used by default.
Currently only MatMult, MatMultAdd, MatMultTranspose, and MatMultTransposeAdd operations are supported,
but more will be added once this initial attempt is at a point
that it provides an acceptable framework for adding additional further functionality.

show more ...

Initial skeleton for MATSEQAIJMKL matrix type that uses MKL sparse BLAS
routines. Currently only implements mult and multadd operations, but
ready to be expanded.

Merge branch 'stefano_zampini/feature-convert-matnest-matis' into stefano_zampini/feature-pcbddc-saddlepoint

I have to merge it here, since there's an API addition for MatIS in stefano_zampini/featu

Merge branch 'stefano_zampini/feature-convert-matnest-matis' into stefano_zampini/feature-pcbddc-saddlepoint

I have to merge it here, since there's an API addition for MatIS in stefano_zampini/feature-pcbddc-saddlepoint
that prevents me to smoothly merge stefano_zampini/feature-convert-matnest-matis into next.

This will also facilitate merging stefano_zampini/feature-pcbddc-saddlepoint into master

Conflicts:
src/docs/website/documentation/changes/dev.html

show more ...

Merge branch 'master' into barry/fix-successful-numerical-factorization-after-failure/master

The maint patch for the given problem cannot merge cleanly into next due to changes
in the same code in m

Merge branch 'master' into barry/fix-successful-numerical-factorization-after-failure/master

The maint patch for the given problem cannot merge cleanly into next due to changes
in the same code in master (the introduction of MatFactorGetError()) therefor need
to merge in master and use the new construct before merging to next

show more ...

Merge branch 'master' into psanan/pc-telescope-nearnullspace

Merge branch 'master' into mark/snes-ex56

* master: (93 commits)
Web: Added a tutorial
typo
rename X -> U
sowing: update to version that creates manpage URLs for 'enum' types/values
Increa

Merge branch 'master' into mark/snes-ex56

* master: (93 commits)
Web: Added a tutorial
typo
rename X -> U
sowing: update to version that creates manpage URLs for 'enum' types/values
Increase patchlevel to 3.7.3
move runex3_3 to NOTSINGLE [ partial backport 9ed0bd76d5c329eb91401ebabfbd851446214e20 ]
DMPlex: local coordinate vector should be a sequential vector
DM: Name coordinate vectors in subdomain and restrict hooks
tests: disable a couple of TS examples from single precision tests
reduce TS ex12_2 output
fix PetscFormatRealArray() to cast PetscReal in print statement to double so works for __float128 precision
Vec: enable Fortran interfaces for VecGetLocalVectorRead and VecRestoreLocalVectorRead
Bib: Added refs
PCShell: man page fixes
update handling of latex in example code to use https
fix typo in isltog comment
test: Avoid conversion of object <-> function pointer type [-Wpedantic]
fortran: enable interface for VecGetLocalVector(), VecRestoreLocalVector()
update citation for cmz2016
PCBDDC: silent compilation warnings from nightlybuilds
...

# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.

show more ...

Merge branch 'karlrupp/feature-viennacl-hybrid'

Conflicts:
src/mat/impls/aij/mpi/makefile

* karlrupp/feature-viennacl-hybrid:
Added hybrid capabilities (CUDA, OpenCL, OpenMP) for ViennaCL bindi

Merge branch 'karlrupp/feature-viennacl-hybrid'

Conflicts:
src/mat/impls/aij/mpi/makefile

* karlrupp/feature-viennacl-hybrid:
Added hybrid capabilities (CUDA, OpenCL, OpenMP) for ViennaCL bindings.
Depending on whether the user configures with CUDA, OpenCL, or OpenMP,
the respective backends will be available in ViennaCL.
Available options:
- viennacl_view: Lists available and selected backends
- viennacl_backed: One out of CUDA, OpenCL, OpenMP (if enabled in configure)
- viennacl_opencl_device_type: cpu, gpu, accelerator (OpenCL device type selection)

show more ...

Remove aij/seq/strumpack.
Rename aij/mpi/strumpack_mpi to aij/mpi/strumpack.
Make aij/mpi/strumpack work for MPISEQAIJ.
Rename CENTRALIZED to REPLICATED.

Add interface to STRUMPACK-sparse.

STRUMPACK-sparse is a sparse direct solver and
preconditioner (using Hierarchically Semi-Separable matrices),
developed at LBNL, released under a BSD license:
http

Add interface to STRUMPACK-sparse.

STRUMPACK-sparse is a sparse direct solver and
preconditioner (using Hierarchically Semi-Separable matrices),
developed at LBNL, released under a BSD license:
http://portal.nersc.gov/project/sparse/strumpack/

This commit adds 2 interfaces to STRUMPACK-sparse: one for
aij/seq and one for aij/mpi matrices. It also adds a script to
the buildsystem to download/compile/install STRUMPACK-sparse.

I can get everything compiled with
./configure --download-strumpack --with-cxx-dialect=C++11 \
--with-openmp --download-scalapack --download-parmetis
--download-metis --download-ptscotch
(There was an issue compiling STRUMPACK-sparse as a shared
library, so it is now always static)

I tested the new code with the src/snes/examples/tutorials/ex5.c
example:
OMP_NUM_THREADS=2 mpirun -n 2 ./ex5 -da_grid_x 100 \
-da_grid_y 100 \
-pc_type lu -pc_factor_mat_solver_package strumpack_mpi \
-ksp_monitor -mat_strumpack_mpi_matinput DISTRIBUTED \
--sp_mc64job 0 --sp_reordering_method scotch --sp_verbose

The aij/mpi solver is called 'strumpack_mpi', the aij/seq solver
is 'strumpack'. All command line options are passed to the
STRUMPACK solver, including -h, which prints out the STRUMPACK
options.

STRUMPACK-sparse has 3 interfaces:
* A sequential/nulti-threaded interface. This one is used with
aij/seq matrices by selecting the 'strumpack' solver.
* A distributed memory solver with replicated interface. This
interface has the same input on every MPI proc (entire matrix
and rhs and solution vectors). Select this for aij/mpi, with
"-pc_type lu -pc_factor_mat_solver_package strumpack_mpi -mat_strumpack_mpi_matinput GLOBAL"
* A fully distributed solver, works with the aij/mpi block-row
distributed input matrix directly. Use with:
"-pc_type lu -pc_factor_mat_solver_package strumpack_mpi -mat_strumpack_mpi_matinput DISTRIBUTED"

Pieter

show more ...

ViennaCL: Added hybrid capabilities (CUDA, OpenCL, OpenMP)

Depending on whether the user configures with CUDA, OpenCL, or OpenMP,
the respective backends will be available in ViennaCL.
CUDA has prec

ViennaCL: Added hybrid capabilities (CUDA, OpenCL, OpenMP)

Depending on whether the user configures with CUDA, OpenCL, or OpenMP,
the respective backends will be available in ViennaCL.
CUDA has precedence over OpenCL, which has precedence over OpenMP.
If none of the backends is selected, ViennaCL is built in sequential mode.

Remaining tasks to better leverage these hybrid capabilities:
- Select backend at runtime (API + options)
- Eliminate copy-overhead for OpenMP backend
- Documentation, examples, and tests to demonstrate usage

show more ...

Merge branch 'master' into pr263/PierreBdR/ts-domain-error-in-stage/master

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

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

Merge branch 'master' into tisaac/dmforest

directory layout changed again: merging to switching branches back and forth
doesn't cause reconfigure to have a hard time

* master: (357 commits)
remov

Merge branch 'master' into tisaac/dmforest

directory layout changed again: merging to switching branches back and forth
doesn't cause reconfigure to have a hard time

* master: (357 commits)
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
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
add a test for MatPtAP()
small cleanup; add tests
cleanup MatPtAPSymbolic_MPIAIJ_MPIAIJ()
Fix nonzerostate tracking in all MATMPI types.
update location of test matrices in faq.html file
unused variable still declared
Fix for MatLoad_MPIDense() from sparsely stored matrix also fix for MatLoad_MPIDense_DenseInFile() when matrix has not had its local sizes preset by user
Preserve the B nonzerostate during disassembly.
TSEvent: Fix arch-mswin nightly errors due to declaring variables in the middle of a function.
Fix error in parallel loading of dense format matrix
removed KSPSetNullSpace() since mat already carries that information and having two locations for it is confusing, unnecessary and prone to error added MatSetTransposeNullSpace() for providing the null space of the transpose of the matrix and added support in KSPSolve() to remove this null space from the right handside of the linear system thus making the system consistent.
...

show more ...

Merge branch 'master' into knepley/feature-composable-callbacks

* master: (151 commits)
PetscMalloc: Fix handling of zero-sized allocations
add note to DMGetLocal/GlobalVector() about intended s

Merge branch 'master' into knepley/feature-composable-callbacks

* master: (151 commits)
PetscMalloc: Fix handling of zero-sized allocations
add note to DMGetLocal/GlobalVector() about intended short-time usage of the vectors.
sf: use patched openmpi tarball that has a potential fix for sf example hang
SNES problems that set VI bounds must use a VI solver since the other solvers ignore the bounds
pounders viewer did not cast PetscReal to double in printf nor use %D instead of %d
fix ts ex20adj
fix ts theta
fixed another variable declaration in the middle of a code block
example main() had PetscInt return type and argument that should be int
DMDAVTKWriteAll_VTS: Fix handling of higher-dimension coordinates
fixed example that declared variables in the middle of a function
C programs do not return PetscInt nor do they take an argument that is PetscInt, fixed example that did
You cannot declare variables in the middle of functions in PETSc.
add ierr checking to PetscOptionsBegin/End() in example
configure: add error check if user specifies PETSC_ARCH=-dbg [i.e starts with a '-']
scale norm of vector difference by norm of vector to insure it is small enough to consider the two vectors equal
missing PetscViewerDestroy() added
MatNestGetISs: use custom fortranstub to handle PETSC_NULL_OBJECT
update location of open-mx website on FAQ page
If the user sets bounds and calls a solver that does not support them, throw an error.
...

show more ...

Merge branch 'pr192/Fande-Kong/scalable-matincreaseoverlap/master' of https://bitbucket.org/petsc/petsc

Merged petsc/petsc into master