PCBDDC: nullspace info for the benign trick

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: make code for checking operator a function

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: unused variable

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: fix memory leak

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

Merging with Barry's refactor of PetscViewerGetSubComm

Merge branch 'barry/fix-petscviewer-attempt-2' into dmay/pcsemiredundant

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of git. I think it is not justified to ever have anything in next for more than a few weeks at most.
Either take it out of next if it is broken or put it in master if it is not broken. Hanging around in next but
not master for ever is not a good policy.

show more ...

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