PCBDDC: no need to change the iteration matrix when the user specified the change of basis with the benign subspace trick

PCBDDC: avoid PtAP during presolve using the benign subspace trick

With the trick, there's no reason to use pmat != amat, so avoid projecting amat local since pmat local has been already projected

PCBDDC: bugfix

PCBDDC: use benign_original_mat instead of M

PCBDDC: allow adaptive selection of constraints using the exact Schur complement instead of the (possibly) incorrect one

In order to implement this efficiently and to reuse the solvers, some modific

PCBDDC: allow adaptive selection of constraints using the exact Schur complement instead of the (possibly) incorrect one

In order to implement this efficiently and to reuse the solvers, some modifications has been made to the code

show more ...

PCBDDC: bugfix in when computing cc when local subdomain info is present

PCBDDC: added support to local primal computations with factored matrix whose size is not equal to n_R

PCBDDC: added API for the specification of primal vertices in global ordering

propagate info on primal vertices to coarser levels when using the benign subspace trick (still need to be improved)

PCBDDC: fixes for multiple PCSetUp calls using the benign trick

MatMUMPS: added API to set info on Schur complement using MUMPS convention

When A is just symmetric and we eliminate all the dofs accounting for the negative eigenvalues, S is positive
Like P2/P0 St

MatMUMPS: added API to set info on Schur complement using MUMPS convention

When A is just symmetric and we eliminate all the dofs accounting for the negative eigenvalues, S is positive
Like P2/P0 Stokes eliminating the pressures

show more ...

PCBDDC: added MatSeqAIJCompress to remove zero entries from SeqAIJ matrices

this is very very useful when changing the local basis with a PtAP operation, which does not care about possible zero entr

PCBDDC: added MatSeqAIJCompress to remove zero entries from SeqAIJ matrices

this is very very useful when changing the local basis with a PtAP operation, which does not care about possible zero entries in the result
MatDetectSubdomains now implements a threshold check for square dense matrices

show more ...

PCBDDC: call restore indices after having used idxs

PCBDDC: disable local subdomain info before refining connected components

PCBDDC: major improvement in benign code

Allows for multiple local benign dofs (in case any of the local matrices represents disconnected subdomains)
Extended accordingly BDDC class and BDDCGraph cl

PCBDDC: major improvement in benign code

Allows for multiple local benign dofs (in case any of the local matrices represents disconnected subdomains)
Extended accordingly BDDC class and BDDCGraph class
Added MatDetectDisconnectedComponent to identify them

show more ...

PCBDDC: avoid checks if the option has been set in the matrices and spd check if the benign trick has been requested

PCBDDC: usage of benign trick even without change of basis on velocities

this could happen if the user has already changed the basis on the velocities for local matrices before calling KSPSetUp

PCBDDC: complex build fixes

PCBDDC: Workaround for rhs in benign subspace

When the rhs is already in the benign space for Darcy, CG suffers a lot (still need to write down some formulas to understand why); instead GMRES conver

PCBDDC: Workaround for rhs in benign subspace

When the rhs is already in the benign space for Darcy, CG suffers a lot (still need to write down some formulas to understand why); instead GMRES converges in 1 (!) iteration.
With this, CG converges without iterating

show more ...

PCBDDC: nullspace info for the benign trick

PCBDDC: remove unneeded lines of code

PCBDDC: handle possible non-zero diagonal terms in pressure blocks (wells)

By default, the code uses the zero diagonal dofs as pressure dofs.
If a local dof split is present, assumes that the last f

PCBDDC: handle possible non-zero diagonal terms in pressure blocks (wells)

By default, the code uses the zero diagonal dofs as pressure dofs.
If a local dof split is present, assumes that the last field represents pressures
Then checks the two fields and avoid change of basis on pressures if the two fields are not the same

show more ...

PCBDDC: avoid possible deadlock

PCBDDCGraphGetDirichletDofs is collective whereas the pressure dofs (if any) are local

PCBDDC: raise error if CSR graph is not consistent

PCBDDC: extend benign trick to PCApplyTranspose_BDDC

PCBDDC: make code for checking operator a function