PCBDDC: provide API for global to local translation (more general than ISGlobalToLocal and memory scalable)

PCBDDC: provide API for setting Dirichlet/Neumann boundaries in global ordering

PCBDDC: replace SetXXXBoundaries with SetXXXBoundariesLocal (better logic)

PCBDDC: fix for new reuse logic in the multilevel case

PCBDDC: remove unneeded code to check nearnullspace

PCBDDC: remove trailing whitespaces from source

PCBDDC: solved compilation error with complex build

PCBDDC: introduce a new variable in PCBDDC data stucture to store the actual number of current vertices (corners with constraints with basis changed)

removed calls to PCBDDCGetPrimalVerticesLocalIdx

PCBDDC: introduce a new variable in PCBDDC data stucture to store the actual number of current vertices (corners with constraints with basis changed)

removed calls to PCBDDCGetPrimalVerticesLocalIdx which are no longer needed

show more ...

PCBDDC: make BtoNmap a member of bddc data struct (created once after PCISSetUp)

Use GlobalToLocalMapping to map local to boundary indices
Inserted some memcpy to speed up memory movement

PCBDDC: check of dirichlet and neumann solvers is not needed if pcbddc->NullSpace is not present

Added some PetscSynchronizedAllow since on some machines the verbose output is broken otherwise

PCBDDC: improved reuse of coarse matrix

MatCovert_IS_AIJ renamed in MatISGetMPIXAIJ (MatConvert does not support out-of-place reuse of destination matrix)
Added some lines in MATIS code (commented).

PCBDDC: better management of primal space. Now a coarse numbering is computed if and only if the primal space has been really changed

PCBDDC: allow different customization of dirichlet and neumann boundary indices among successive linear solves

PCBDDC: now PCBDDC reuses part of its private data if NearNullSpace customization is changed between different linear solves

I know it is a large commit, but this is the minimal set of changes to ge

PCBDDC: now PCBDDC reuses part of its private data if NearNullSpace customization is changed between different linear solves

I know it is a large commit, but this is the minimal set of changes to get it working

show more ...

PCBDDC: local work vectors should be created only if new constraints are created and the relevant sizes are different

PCBDDC/PCIS: make vec4_D a member of pcis, not pcbddc

PCBDDC: Prepare for different primal spaces customization among different linear solves

now bddc primal space is built using a NullSpace attached to pmat. It only remains to determine whether the Ne

PCBDDC: Prepare for different primal spaces customization among different linear solves

now bddc primal space is built using a NullSpace attached to pmat. It only remains to determine whether the NearNullSpace object has changed

show more ...

PCBDDC: demand object destruction to setup routines for BDDCGraph and constraints

PCBDDC: updated example using local to global mapping generated by DMDA

PCBDDC: solved a bug when multiple constraints are present and some are eliminated by SVD

change_of_basis flag should be set AFTER SVD has recognized independent components
Some minor modications to

PCBDDC: solved a bug when multiple constraints are present and some are eliminated by SVD

change_of_basis flag should be set AFTER SVD has recognized independent components
Some minor modications to the code (due to broken MatMult_BAIJ with compressed rows)

show more ...

PCBDDC: index set for R nodes now stores block size information

When IS block size is equal to Mat block size, there's no need to convert local matrices anymore

PCBDDC: updated todo list

added a further guard before calling SVDs

PCBDDC: updated man pages

PCBDDC: bddc now works with sbaij matrices.

This is a partial solution, since in practice the Neumann problem matrix will always be seqaij

Implemented user switch for field split

PCBDDC: added field in bddc data struct to determine if fields splitting has been provided by the user

Initialization of number of field splits was missing. Now fixed