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

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

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

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.

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

PCBDDC: adaptive selection of constraints now works without explicitly changing the basis on the velocities

PCBDDC: change the basis during scaling/average operations

PCBDDC: start supporting no-net-flux condition as a constraint

PCBDDC: construct a fake change of basis that will be needed if we want to use deluxe scaling withouth change the basis on the pressures

The added code also goes in the direction of supporting a mix

PCBDDC: construct a fake change of basis that will be needed if we want to use deluxe scaling withouth change the basis on the pressures

The added code also goes in the direction of supporting a mix of user-defined constraints and adaptively generated primal spaces.

show more ...

PCBDDC: add runtime error. I still have to implement this.

PCBDDC: assign enough space for nonzeros in B_0

PCBDDC: if zerodiag_subs is not available from the reusage of the local solvers, compute it

PCBDDC: include primal vertices (if any) in the computation of the local Schur complements even if we are not doing adaption

When converting from PetscViewerSetFormat() to PetscViewerPushFormat() I missed adding the PetscViewerPopFormat() in several places

Reported-by: Jed Brown <jed@jedbrown.org>

Deprecate PetscViewerSetFormat() since it better to use PetscViewerPush/PopFormat() and using
PetscViewerSetFormat() can lead to confusing situations for users when later viewers do not recognize the

Deprecate PetscViewerSetFormat() since it better to use PetscViewerPush/PopFormat() and using
PetscViewerSetFormat() can lead to confusing situations for users when later viewers do not recognize the set
format and thus do not produce any information.

Reported-by: Jed Brown <jed@jedbrown.org>

show more ...

PCBDDC: prevent from unneeded computations in PCPreSolve_BDDC with benign trick