PCBDDC: make code more clear

PCBDDC: move check into separate function

PCBDDC: move code to a separate function

PCBDDC: use consistent name for setting or retrieving P0 values

PCBDDC: remove stupid requirement of nonzero near null space on vertices

PCBDDC: Enable benign trick even with constraints

The user is responsible to either

1) use an iteration matrix whose basis is already changed to make \int_F v*n a primal dof for all F
or
2) attach

PCBDDC: Enable benign trick even with constraints

The user is responsible to either

1) use an iteration matrix whose basis is already changed to make \int_F v*n a primal dof for all F
or
2) attach a proper MatNullSpace object which guaranties that B0_R is null
or
3) provide a proper change of basis

In cases 1) and 3), adaptive selection of constraints could be used.

show more ...

PCBDDC: added TODO

PCBDDC: avoid playing with pointers. Use API instead

PCBDDC: unused variable

PCBDDC: Dirichlet dofs on internal nodes should be imported

this is needed by the benign trick because
b(v_I,p_0) should be zero for all v_I = 0 \Partial\Omega (including Essential boundary conditio

PCBDDC: Dirichlet dofs on internal nodes should be imported

this is needed by the benign trick because
b(v_I,p_0) should be zero for all v_I = 0 \Partial\Omega (including Essential boundary conditions)

show more ...

PCBDDC: remove primal vertices from face/edge deluxe scaling

this does not change the behaviour of standard BDDC.
It is required by the benign trick in order to
keep the iterates in the benign space

PCBDDC: remove primal vertices from face/edge deluxe scaling

this does not change the behaviour of standard BDDC.
It is required by the benign trick in order to
keep the iterates in the benign space

remove extra code to check benign-ness inside PCBDDCMatMult_Private

show more ...

PCBDDC: raise error if benign trick and user defined change of basis are both present

PCBDDC: Wrong update during postsolve and forgot to restore the original local matrix

Now the benign trick works
Still need to do something on the deluxe side

PCBDDC: fix memory leak

PCBDDC: default values should be set during creation

PCBDDC: added benign correction during presolve and check on the BDDC operator

still does not work, need to understand what is missing

there is some extra code that should be removed from PCBDDCMat

PCBDDC: added benign correction during presolve and check on the BDDC operator

still does not work, need to understand what is missing

there is some extra code that should be removed from PCBDDCMatMult_Private

show more ...

PCBDDC: MUMPS Cholesky solver could be used for symmetric matrices

This has also the advantage that POTRF/S/I will be called inside BDDC
SYTR/S/I will be called inside MUMPS (can I fix this?)

PCBDDC: remove zerodiag from bddc data structure

PCBDDC: forget to protect line by n_benign guard

PCBDDC: fix error in coarse submatrix computation for the benign trick

debug code now is properly working

PCBDDC: forgot to destroy objects

PCBDDC: added mechanism to pass p0 from residual vector to coarse space

still does not work completely
need to add more checks

PCBDDC: move change of basis on pressure during PreSolve

Added check for b(v_I,p_0) = 0 for all v_I (o on Gamma and \partial\Omega)

PCBDDC: extend computation of coarse basis functions and coarse matrix for the benign trick

debug output temporary suppressed,
the code is still not working.
- GMRES with BDDC setup and dummy applic

PCBDDC: extend computation of coarse basis functions and coarse matrix for the benign trick

debug output temporary suppressed,
the code is still not working.
- GMRES with BDDC setup and dummy application (PCApply as VecCopy) does not converge
- still need to extract p0 from residual

show more ...

PCBDDC: forgot to set columns values