CUSP,ViennaCL: Fixed ILU and Cholesky preconditioner fallback.

The new registration mechanism for matrix factorizations was not
applied for CUSP and ViennaCL, thus the default ILU preconditioners di

CUSP,ViennaCL: Fixed ILU and Cholesky preconditioner fallback.

The new registration mechanism for matrix factorizations was not
applied for CUSP and ViennaCL, thus the default ILU preconditioners did not work.
This commit adds the missing registrations for these packages.

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

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

Merge branch 'master' into tisaac/ts-ex11-amr

* master: (141 commits)
chombo: fix dirname in tarball and update to use snapshot with gcc-6.1 compile fix
petsc4py: fix tarball url [for arbitary c

Merge branch 'master' into tisaac/ts-ex11-amr

* master: (141 commits)
chombo: fix dirname in tarball and update to use snapshot with gcc-6.1 compile fix
petsc4py: fix tarball url [for arbitary commitid] - and dirname in tarball
Plex: Moved ReconstructGradients_Internal() and ApplyLimiter_Internal() to Plex from SNES
CUDA: Added external declaration for cuBLAS handle.
docs: fix reference to KSP_Create in parent commit
Add information on PETSC_EXTERN and PETSC_INTERN to developers manual
Plex: Symbol is shared between DM and SNES, so must be EXTERN
Plex: Remove unused variables
Configure: Fix check for spaces in PETSC_DIR
Plex: Updated examples to new DMPlexCreateBoxMesh()
sick of failure due to incorrect master branch
SNES ex12: Fixed unstructured GMG test
SNES ex12: Added another unstructured GMG test
petsc4py: use the snapshot with the fix for dmplex change
Examples that require Fortran to compile but were listed in tests that did not require Fortran
Plex ex19: Added DMSetFromOptions()
Plex: Store a flag indicating whether the boundary should change during remeshing - Set using -dm_plex_remesh_bd - Changed Remesh_Internal() signature
Always coarsen boundaries.
Plex: Fixed option name
Plex: Protect unused variables
...

Conflicts:
src/dm/impls/plex/makefile
src/snes/utils/dmplexsnes.c

show more ...

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

Conflicts:
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcprivate.c
src/mat/impls/aij/seq/mkl_pardiso/mkl_pardiso.c

Merge branch 'master' into barral/feature-dmplex-adapt

Merge branch 'barry/housekeeping-1'

removed a bunch of unneeded inclusions of xxximpl.h which exposed private data where it did not have to be exposed

also made more systematic the inclusion of the base xxximpl.h in the derived yyyimp

removed a bunch of unneeded inclusions of xxximpl.h which exposed private data where it did not have to be exposed

also made more systematic the inclusion of the base xxximpl.h in the derived yyyimpl.h
so the inclusion is done in the include file rather than in each source code that includes it

show more ...

Merge branch 'jose/maint/tabs-print-class' into maint

Merge branch 'master' into knepley/feature-parallel-mesh-load

* master: (59 commits)
Increase patchlevel to 3.7.1
Fix manual page for MatGetTransposeNullSpace()
remove unneeded use of PetscCon

Merge branch 'master' into knepley/feature-parallel-mesh-load

* master: (59 commits)
Increase patchlevel to 3.7.1
Fix manual page for MatGetTransposeNullSpace()
remove unneeded use of PetscContainer to hold iscoloring in SNESComputeNGSDefaultSecant()
added SNESView for ngs, added options database to ngs manual page, removed dead routine
gcov: enable in mpiuni test
incorrect capitalization of an even name
force same number of iterations for output of ml example
changed output for ml example to not be machine dependent
DMP4est: access data after setup
SNESLineSearch_Apply was labeled with PETSC_EXTERN in one place but when declared this was not visible hence library build failed under Microsoft.
Example had print statements of complex numbers
Fix more errors in agmres for 64 bit integers
builds: fix nampespase conflict with 'analyzer-buid' and checkBuilds.py
gcov: add more packages to test
add a couple of KSP that were not represented in testing
gcov: enable more externalpackages for this test
fixes for agmres for 64 bit indices, fix for findvisibleprivatefunctions to be portable for Mac and Linux
get GAMG to get null space (more)
fix c++ compile errors
Add test example for kacmarz preconditioner since gcov indicated there was NONE
...

show more ...

Merge branch 'master' of bitbucket.org:petsc/petsc into mark/gamg-nullspacefix

Merge branch 'barry/cleanup-external-mat-solver-packages'

Merge branch 'barry/make-functions-static'

Merge branch 'master' into barral/feature-dmplex-adapt

Merge branch 'master' into knepley/feature-parallel-mesh-load

* master: (180 commits)
Bib: Added ref
-random_view caused error because PetscRandomView() was not provided in PetscRandomCreate()

Merge branch 'master' into knepley/feature-parallel-mesh-load

* master: (180 commits)
Bib: Added ref
-random_view caused error because PetscRandomView() was not provided in PetscRandomCreate()
Update the xSDKTrilinos test message to use interactive session and ctest
gcc6: fix comple warnings with -Wmisleading-indentation
fix compile warning
bugfix for MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable()
minor formatting cleanup of gamg.c
add ml example to tutorials page for xsdk
bugfix of MatDuplicate for mpiaij matrix product Reported-by: Miorelli, Federico <Federico.Miorelli@cgg.com>
The deprecated -log_summary_exclude was still being used when it should have been -log_exclude
cpardiso: requires mkl_blacs so check for that. Also common code to detect MKL for both pardiso and cpardiso
cpardiso: fix segv and memory leaks
resume MatGetSubMatrix_SeqSBAIJ() because it cannot share routine with BAIJ due to cast
fix bad relative links
replace MatSeqBAIJSetPreallocation_SeqBAIJ with MatSeqBAIJSetPreallocation to skip this call when input mat is sbaij type
fix a variety of set but not used variables in examples detected by clang static analyzer
removed declared but not used variable
ANL Report formatting requires listing authors by affiliation
DMP4est: uncrustify
DMP4est: remove memset()
...

show more ...

Updated script for detecting visible symbols

Removed unneeded visibility of some symbols
Fixed some namespacing
Removed the use of unneeded implementation includes in the TAO code

Merge branch 'master' into barry/cleanup-external-mat-solver-packages

removed bstream code since it was not in nightly tests and broken

Merge branch 'maint'

The deprecated -log_summary_exclude was still being used when it should have been -log_exclude

Reported-by: Jose E. Roman <jroman@dsic.upv.es>

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

Conflicts:
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcgraph.c
src/ksp/pc/impls/bddc/bddcprivate.c

Merge branch 'tisaac/sf-fix-multi-sf-leaves' into tisaac/dmp4est-feature-injection

* tisaac/sf-fix-multi-sf-leaves: (174 commits)
PetscSF: fix PetscSFGetMultiSF() for sparse leaves
PetscSF: make

Merge branch 'tisaac/sf-fix-multi-sf-leaves' into tisaac/dmp4est-feature-injection

* tisaac/sf-fix-multi-sf-leaves: (174 commits)
PetscSF: fix PetscSFGetMultiSF() for sparse leaves
PetscSF: make regression test that fails
few more clang static analyzer fixes
many issues detected by clang static analyzer
Update hypre to its release version
Updates for xSDKTrilinos builds
Added --download-xsdktrilinos
Communicator passed to VecCreateSeq should be PETSC_COMM_SELF
PCMG man page: fix typo in options name
Bib: Added refs
DMLabel: Removed unused variable
Metis tries to use backtrace which requires -lexecinfo on some systems such as freebsd Remove -lexecinfo from freebsd configure files since it is only needed by metis
add alternative output that occurs on some systems due to use of random number generator in partitioner
Plex: When creating cohesive submeshes, fix memory leak
Plex ex11: Updated output for new distribution algorithm
Plex: When explicitly distributing the depth label, we must recreate empty strata
Plex: When partitioning, the cell numbering should include hybrid points
Plex: Allow cell and vertex numberings to include hybrid points
need to show diffs of ex22 if it fails
Support two different output for test example, since due to roundoff this example can produce slightly different convergence history
...

show more ...