Enhancements for KSP solve with singular matrices

update KSPSolveTranspose() to utilize the nullspaces attached to the matrix with MatNullSpaceSet()
add test examples that utilizes the singular matr

Enhancements for KSP solve with singular matrices

update KSPSolveTranspose() to utilize the nullspaces attached to the matrix with MatNullSpaceSet()
add test examples that utilizes the singular matrices nullspaces for KSPSolve() and KSPSolveTranspose() with
right and left preconditioning.
Lock vectors in MatNullSpaceCreate() so user doesn't accidently change them when they should not
Add details to KSPSolve() manual page for singular systems.

Thanks-to: Fande Kong <fande.kong@inl.gov>

show more ...

PCBDDC: fix from pipeline build

ViennaCL: Added "row-scaling" preconditioner.

This is a simple Jacobi-like diagonal preconditioner.
The preconditioner computes the norms of each row of the system matrix and uses these values for t

ViennaCL: Added "row-scaling" preconditioner.

This is a simple Jacobi-like diagonal preconditioner.
The preconditioner computes the norms of each row of the system matrix and uses these values for the diagonal of the preconditioner matrix.
In contrast to PETSc's native Jacobi preconditioner, this preconditioner bypasses a couple of unnecessary host-device transfers (via e.g. PCI-Express) for GPUs.

show more ...

PCBDDC: fix usage of nvtxs_csr and output

PCBDDC: Rearrangement of code

added support field=PETSC_DECIDE

singular problems?

When the Schur complement form of PCFIELDSPLIT is called PCSetUp() needs to be called before PCFieldSplitGetSubKSP()

Add very useful error message to this case
Time: .1 hour
Reported-by: David Nolte

When the Schur complement form of PCFIELDSPLIT is called PCSetUp() needs to be called before PCFieldSplitGetSubKSP()

Add very useful error message to this case
Time: .1 hour
Reported-by: David Nolte
Thanks-to: Arthur Soprano, Lisandro Dalcin

show more ...

PCBDDC: fix for missing gesvd

ViennaCL: Fixed parallel runs with block-Jacobi.

Now it is possible to run ViennaCL with multiple MPI ranks again.

PCBDDC: allow for null adjncy array (nvtxs > 0 and ii[0...nvtxs] = 0 is a valid csr)

Added example that uses DM_BOUNDARY_MIRROR for Neumann boundary conditions

Put basic tests for DMDA DM_BOUNDARY_MIRROR in 1d and 2d into makefile
Improved manual page for DM_BOUNDARY_MIRROR
Made DMS

Added example that uses DM_BOUNDARY_MIRROR for Neumann boundary conditions

Put basic tests for DMDA DM_BOUNDARY_MIRROR in 1d and 2d into makefile
Improved manual page for DM_BOUNDARY_MIRROR
Made DMSetUp_DA_3d() generate an error on DM_BOUNDARY_MIRROR since it is not currently supported

Time: 2 hours
Thanks-to: Mohammed Al-Mosallam <Al-MosallamMB@cardiff.ac.uk>

show more ...

PCBDDC: add flag to compute range or kernel

ViennaCL: Added parallel ILU preconditioner as proposed by Chow and Patel.

See paper "Fine-Grained Parallel Incomplete LU Factorization", SIAM SISC 37(2)
http://epubs.siam.org/doi/abs/10.1137/140968

ViennaCL: Added parallel ILU preconditioner as proposed by Chow and Patel.

See paper "Fine-Grained Parallel Incomplete LU Factorization", SIAM SISC 37(2)
http://epubs.siam.org/doi/abs/10.1137/140968896

Example for running with OpenCL with KSP ex2:
./ex2 -mat_type aijviennacl -vec_type viennacl -pc_type chowiluviennacl -viennacl_view -viennacl_backend opencl

show more ...

PCBDDC: add local subs to the output of PCView_BDDC

PCBDDC: fix computation of twodim flag

ViennaCL: Added smoothed aggregation preconditioner (saviennacl) similar to sacusp.

Requires commit 6e09ddc from ViennaCL repository. Release 1.7.1 is insufficient.

Merge branch 'maint'

PCBDDC: fixes for uniprocessor runs

bugfix: initialize pc->failedreason = PC_NOERROR for PCSetUp_xxx()
Reported-by: David Knezevic <david.knezevic@akselos.com>

PCBDDC: use t_coarse_mat_is to compute subassembling pattern if multilevel has not been requested

PCBDDC: need to remove dirichlet dofs shared by more than two subregions

if a nodal dof is connected to a dirichlet dof set it to corner

added test example for ex42.c

Improved manual page
Changed MatSchurComplementComputeExplicitOperator() to actually compute
the Schur complement and not ignore the 1,1 block of the matrix

Time: .1 h

added test example for ex42.c

Improved manual page
Changed MatSchurComplementComputeExplicitOperator() to actually compute
the Schur complement and not ignore the 1,1 block of the matrix

Time: .1 hours
Reported-by: Hoang Giang Bui <hgbk2008@gmail.com>

show more ...

PCBDDC: add extra information to ASCIIView

PCBDDC: fix bug with tetrahedral meshes. A nodal dof is a corner also if the set of sharing subdomains is different from that of the edge dofs

Merge branch 'mark/gamg-aijtest'

PCTelescope: add citation