PCBDDC: optimize refinement of connected components

The previous code used a lot of memory constructing an ad-hoc csr information to refine the connected components
Changed the way the connected com

PCBDDC: optimize refinement of connected components

The previous code used a lot of memory constructing an ad-hoc csr information to refine the connected components
Changed the way the connected component are searched in PCBDDCComputeConnectedComponentsLocal: now the search is symmetric

show more ...

PCBDDC: output the global id of the wrongly computed coarse index

Merge branch 'barry/fix-private-functions'

Merge branch 'hzhang/fix-pcredundant'

PCBDDC: use BLAS rank updates routines

PCBDDC: less stringent check on benign trick applicability

Added DMShellSetRestriction() and ability of PCMG to use DMCreateRestriction()

Time: .2 hours
Reported-by: anton <popov@uni-mainz.de>

fix MatSolverPackageGet() so it will never return a different one than what was requested (because the requested one did not exist). This is done by setting the initial package to null and then order

fix MatSolverPackageGet() so it will never return a different one than what was requested (because the requested one did not exist). This is done by setting the initial package to null and then ordering the petsc factorization routines before any other package

Time: .5 hours
Reported-by: Hong Zhang <hzhang@mcs.anl.gov>

show more ...

PCBDDC: remove faster deluxe switch

Fix -ksp_monitor_singular_value
Fix doc for failed factorization in PC

PCBDDC: use MPIU_REPLACE to for SFReduce ops when setting the p0 values (avoid VecAssemblyBegin/End)

PCBDDC: remove the horrible logic around the remap of the coarse problem

the code now is more readable, maintanable. In case of multilevel, it is also faster, since one shift is avoided.
when the co

PCBDDC: remove the horrible logic around the remap of the coarse problem

the code now is more readable, maintanable. In case of multilevel, it is also faster, since one shift is avoided.
when the coarse matrix has to be assembled, it still uses MatISSubassemble followed by MatISGetMPIXAIJ

fixed some bugs in MatISSubassemble and MatISGetSubassemblingPattern, with preliminary checks added

show more ...

MatSolverPackageGet() search for a working package
if input package is not found or input package = NULL;
PCSetUp_Redundant() set gmres/bjacobi when matsolver package is not found

PCBDDC: multilevel fixes

allow simple redistribution on void procs (if any) at level 0

PCBDDC: allow processes to pass 0 nvtxs and void arrays

PCBDDC: make change_primal_sub a member of sub_schurs

we need it during adaptive selection
sub_schurs->change and sub_schurs->chane_primal_sub are created when needed

Merge branch 'master' into barry/fix-private-functions

PCBDDC: Added support for adaptive selection of constraints with pointwise scalings

PCBDDC: use PtAP for SeqDense from master

PCBDDC: use MAT_INPLACE_MATRIX

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
src/ksp/pc/imp

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
src/ksp/pc/impls/bddc/bddcscalingbasic.c
src/ksp/pc/impls/bddc/bddcschurs.c
src/mat/impls/is/matis.c

show more ...

small cleanup

add PetscSubcommSetOptionsPrefix()

labeled a variety of additional internal functions as static
to prevent them from being publically available in the library

Merge branch 'stefano_zampini/fix-for-release'