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

Added MATMPIAIJMKL and fixed errors in MatMultAdd_SeqAIJMKL() and
MatMultTransposeAdd_SeqAIJMKL().

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

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bddc.c

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcfetidp.c
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcschurs.c
src/mat/impls/is/matis.c
src/mat/impls/is/matis.h

show more ...

Merge remote-tracking branch 'origin/pr464/pieter_ghysels/pghysels/strumpack-interface-rebased/master'

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

Rename the aij/mpi files from strumpack to strumpack_mpi (to avoid conflicts with aij/seq).

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

Merge branch 'master' into stefano_zampini/pcbddc-adaptive-selection

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c
src/mat/impls/aij/mpi/mumps/mumps.c
s

Merge branch 'master' into stefano_zampini/pcbddc-adaptive-selection

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c
src/mat/impls/aij/mpi/mumps/mumps.c
src/mat/impls/is/matis.h

show more ...

Merged master into semihozmen/ksp-fcg-eigen-comp-support

Merge branch 'master' into pr192/Fande-Kong/scalable-matincreaseoverlap/master

Needed to merge master to update to new calling sequence PetscCommBuildTwoSided

Merged petsc/petsc into master

Merge branch 'barry/fix-directory-layout'