fixed many examples that worked incorrectly for 64 bit indices

PCBDDC: move creation of temporary MATIS matrix before coarse ksp is created.

This way, we can temporary set the coarse MATIS into coarse ksp.
This allows to call pc_bddc_coarse_ksp_fischer_guess <>

PCBDDC: move creation of temporary MATIS matrix before coarse ksp is created.

This way, we can temporary set the coarse MATIS into coarse ksp.
This allows to call pc_bddc_coarse_ksp_fischer_guess <> from option

show more ...

PCBDDC: remove trailing backspace

PCBDDC: stefano_zampini/apply-transpose branch was still using pcbddc->vec4_D. Now the code compiles.

Merge branch 'stefano_zampini/pcbddc-applytranspose' into stefano_zampini/pcbddc-multilevel-improvement

PCBDDC: demand wheter or not to rebuild local scatters to the right routine

PCBDDC: remove prints. Added coarse_subassembling to debugging output

PCBDDC: added reuse of coarse subassembling information

In order to do do, I had to change the MATIS class (I dont know why there was a guard for is->mapping in matis code)

Solved a couple of bugs

PCBDDC: added reuse of coarse subassembling information

In order to do do, I had to change the MATIS class (I dont know why there was a guard for is->mapping in matis code)

Solved a couple of bugs when doing multilevel and a dbg_viewer on a single level only

show more ...

PCBDDC: coarse ranks are no longer contiguous (better memory management)

PCBDDC: remove prints

PCBDDC: KSPRICHARDON with PCBDDC is the default for BDDC coarse solver when doing multilevel

Extreme eigenvalues are computed if the coarse KSP is not preonly to accelerate convergence

PCBDDC: added FETIDPMatMultTranspose and PCApplyTranspose (for Dirichlet FETIDP preconditioner)

PCBDDC: support for PCApplyTranspose completed by solving the transpose of the coarse operator in PCBDDCApplyInterfacePreconditioner.

Function pointer is now set during object creation

PCBDDC: remove unneded guards

these guards were needed in the previous version of the code when PCBDDC coarse problem was handling through the enum PCBDDCCoarseProblemType

PCBDDC: added support for solving the transposed of the substructure correction

rhs, solution and work vectors for the correction are now parameters of the function

PCBDDC: added second parameter to PCBDDCApplyInterfacePreconditioner to determine if the interface preconditioner should be applied tranposed or not

The main reason for this is to avoid code duplica

PCBDDC: added second parameter to PCBDDCApplyInterfacePreconditioner to determine if the interface preconditioner should be applied tranposed or not

The main reason for this is to avoid code duplication

show more ...

PCBDDC: simplify PCBDDCApplyInterfacePreconditioner

coarse_psi_* now always exist. For symmetric problems, they point to the same object of coarse_phi_*

PCBDDC: rename private PCBDDCSolveSaddlePoint to a more explicative PCBDDCSolveSubstructureCorrection

PCBDDC: added initial implementation of PCApplyTranspose_BDDC

PCBDDC: propagate field split and neumann boundaries to coarse problem

PCBDDC: coarse matrix is structurally symmetric

PCBDDC: set compute topography to true with new field split

PCBDDC: set compute topography to true also when Local version of the boundaries API is invoked

Merge branch 'stefano_zampini/pcbddc-globalordering-api' into stefano_zampini/pcbddc-multilevel-improvement

Conflicts:
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcprivate.c

PCBDDC: provide APIs for setting field split in global ordering

previous API moved to PCBDDCSetDofsSplittingLocal

For consistency with updated man pages, which_dof values should default to -1;
this

PCBDDC: provide APIs for setting field split in global ordering

previous API moved to PCBDDCSetDofsSplittingLocal

For consistency with updated man pages, which_dof values should default to -1;
this way, unlisted nodes in field split will belong to a different subset

show more ...