convert private functions to static to decrease the unneccessary function in
the public interface in the libraries

PCFactorSetUpMatSolverPackage() should not crash but instead produce a useful
error message if called before it is in the proper state.

Reported-by: David Knezevic <david.knezevic@akselos.com>

recover PCRedundantGetKSP; call it inside PCSetUp_Redundant() instead

PCBDDC: graph rebuild needs to use the same custom vertex size as the original graph

PCBDDC: fix for static switch

The previous code was not applying M_3 preconditioner as outlined in Dohrmann's reference article in case the preconditioner was freezed

PCBDDC: fix for upcoming release.

When the user provides a change of basis matrix, pcis data structure may not be initialized, and we need to use pc->pmat->rmap->mapping

replace if (!red->psubcomm) { ...} block with PCSetUp() for simplification
as suggested by Jakub Kruzik

PCBDDC: get size of R nodes after creating msolv_ctx->is_R

win32: compile fix

libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_0
libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_1
libpets

win32: compile fix

libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_0
libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_1
libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_2
libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_3
libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_4
libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_S
libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_L
libpetsc.lib(dlregisksp.o) : error LNK2001: unresolved external symbol KSP_Solve_FS_U

show more ...

PCBDDC: Finally, I was able to get rid of changing the iteration matrix during the KSP, which I've never liked from the very beginning

The code is way simpler now. The change of basis (if any) is pe

PCBDDC: Finally, I was able to get rid of changing the iteration matrix during the KSP, which I've never liked from the very beginning

The code is way simpler now. The change of basis (if any) is performed during PCApply and PCApplyTranspose.
Added interior flag to PCBDDCSetChangeOfBasisMat in order to save some work, if possible

show more ...

PCBDDC: remove unneeded code. benign_vec should be set to zero

PCBDDC: added test for constraints with benign trick

A_RR matrix is not extracted unless it is needed (dbg_flag > 0 or reuse_solver not present)

PCBDDC: Finally I found the error causing dbg output to report errors in the Neumann solver and in the coarse submatrix

Schur rhs and interior sol need to be modified.
Now it works great.
Tested wit

PCBDDC: Finally I found the error causing dbg output to report errors in the Neumann solver and in the coarse submatrix

Schur rhs and interior sol need to be modified.
Now it works great.
Tested with RT0 and BDM2

show more ...

Merge branch 'master' into barry/add-pcfieldsplit-logevent

PCBDDC: the change of basis context in sub_schurs is no longer needed if we have changed the basis

still does not give the same condition number as without automatic change of basis (this is why it

PCBDDC: the change of basis context in sub_schurs is no longer needed if we have changed the basis

still does not give the same condition number as without automatic change of basis (this is why it is still disabled)

show more ...

PCBDDC: propagate info to coarser levels

PCBDDC: compute the quadrature weights for the no-net-flux condition if requested

Added logging of individual KSPSolves within PCFIELDSPLIT

Thanks-to: Xujun Zhao <xzhao99@gmail.com>

PCBDDC: add vertex size customization by command line

This overcomes the issues when some subdomain faces are made up by a single element face by treating them as constraints instead of primal verti

PCBDDC: add vertex size customization by command line

This overcomes the issues when some subdomain faces are made up by a single element face by treating them as constraints instead of primal vertices
When the benign subspace technique is used without any change of variables, the S_CV/S_VC blocks of the coarse matrix are not properly computed (still need to understand why)

show more ...

PCBDDC: Separate deluxe subproblems

Recompute_topography now drives the setup/reuse of deluxe
There's still something which does not work properly when some subdomain faces are made up by just one e

PCBDDC: Separate deluxe subproblems

Recompute_topography now drives the setup/reuse of deluxe
There's still something which does not work properly when some subdomain faces are made up by just one element face
There's no need to extract B0_R and B0_V separately
Solved a couple of bugs with pcbddc->switch_static = PETSC_TRUE

show more ...

PCBDDC: add PCView_BDDC for ASCII

Adds a braindead PCView_BDDC which simply prints out all the parameters
which are exposed via the options database. Only implemented for the
ASCII viewer.

Fix error message from Must call PCMGSetInterpolation() or PCMGSetInterpolation()

Reported-by: Boris Boutkov <borisbou@buffalo.edu>

PCBDDC: if we are not setting up the sub_schurs data structures, we have to change explicitly the basis on the pressures

PCBDDC: fix for recompute_topography flag

Objects reference now are taken first

PCBDDC: add two runtime options