PCBDDC: Solved compilation errors and warnings for complex builds

Example code changed by using LAPACK steqr method.

rm white spaces

provide informative error message when MatGetMultiPorcBlock_xxx() is not supported

PCBDDC: add wrappers for PetscStackCallBLAS

PCBDDC: Final fix for locally periodic problems

Added mirrors to graph structure

PCNNCreateCoarseMatrix: fix allocation error

Estimate the number of off-diagonal nonzeros using the number of local
neighbors, but be permissive because this is not sufficient in general.
Exact prea

PCNNCreateCoarseMatrix: fix allocation error

Estimate the number of off-diagonal nonzeros using the number of local
neighbors, but be permissive because this is not sufficient in general.
Exact preallocation would require the connectivity graph for subdomains,
but is not critical for PCNN performance.

show more ...

PCIS: consider myself as a neighbor when classifying interface nodes

Local periodic nodes should be considered has interface nodes.

PCBDDC: Remove misuses of neighbours_set in FETI-DP code

PCBDDC: Coarse solution lives in the span of phi (not psi)

PCBDDC: More robust method to compute local boundary nodes

Using counter vector to infer local boundary nodes is wrong in case a
subdomain has local periodicity.

PCBDDC: Add partial support for periodic problems replacing Sort with RemoveDups

* Local-to-global maps which are not one-to-one if restricted to any
subdomain will still be a problem.
* Graph vie

PCBDDC: Add partial support for periodic problems replacing Sort with RemoveDups

* Local-to-global maps which are not one-to-one if restricted to any
subdomain will still be a problem.
* Graph view now prints also the global number of the queued indices.
* Optimal allocation for queue_in_global_numbering.

show more ...

PCBDDC: Add support for non-symmetric linear systems

PCBDDC: Stiffness scaling should be rebuilt in any case (just to be safe)

PCBDDC: Remove trailing blanks

PCBDDC: Consistency of CSR graph passed in postponed to PCSetUp

CSR graph is freed if the PC type is not PCBDDC and PETSC_OWN_POINTER
has been requested.

(Thanks to Lisandro Dalcin to report on thi

PCBDDC: Consistency of CSR graph passed in postponed to PCSetUp

CSR graph is freed if the PC type is not PCBDDC and PETSC_OWN_POINTER
has been requested.

(Thanks to Lisandro Dalcin to report on this)

show more ...

PCBDDC: Preliminary work to split PCBDDCCoarseSetUp in simpler pieces

PCBDDC: allow compilation without metis header

Coarse problem type default to PARALLEL_BDDC
Some irrelevant warning will show up without metis headers (will be solved soon)

PCBDDC: Optimal allocation in PCBDDCGetPrinalVerticesLocalIdx

PCBDDC: add verbosity level to graph view

PCBDDC: change boolean debug flag to integer -pc_bddc_check_level

Distribution of local primal indices requires level > 1.

PCBDDC: Workaround for BAIJ local matrices: convert to AIJ for coarse matrix

PCIS: Use PETSC_COMM_SELF instead of MPI_COMM_SELF

PCBDDC: many changes

- Source code split
- Added Reset procedures
- Support for different MatStructure flags

Complex fixes for classical

Mat: remove third argument to MatNullSpaceRemove

Out of place removals were dangerous because the MatNullSpace held the
only reference to the resulting Vec, which is dangerous. This use case
is rar

Mat: remove third argument to MatNullSpaceRemove

Out of place removals were dangerous because the MatNullSpace held the
only reference to the resulting Vec, which is dangerous. This use case
is rare (never used in PETSc itself) and can easily be accomplished:

VecDuplicate(x,&y);
VecCopy(x,y);
MatNullSpaceRemove(nullsp,y);

Suggested-by: Tobin Isaac <tisaac@ices.utexas.edu>

show more ...