bugfix of MatCreateRedundantMatrix() for rectangular matrices
Reported-by: Jose E. Roman <jroman@dsic.upv.es>

Registered MATSEQAIJMKL as a subtype of MATSEQAIJ.

Removed DEBUG statements.

Fixes to AIJMKL MatDuplicate and MatAssemblyEnd.

Added another #ifdef to allow use of MKL that does not have
mkl_sparse_optimize().

Code cleanup to eliminate some compiler warnings.

Make sure that mkl_sparse_destroy() is only called if
mkl_sparse_optimize() has been called; more work to allow AIJMKL code to
work with older versions of MKL.

Added code to handle case in AIJMKL in which the MKL library version is
older and does not have the sparse inspector-executor routines.

Fixed AIJMKL code breaking things when PETSc is built without Intel MKL.

Use of SpMV2 inspector-executor model can now be turned off via command
line options.

Added support for the MKL SpMV2 inspector-executor routines for
matrix-vector multiplication.

Added MATMPIAIJMKL and fixed errors in MatMultAdd_SeqAIJMKL() and
MatMultTransposeAdd_SeqAIJMKL().

Added support for MatMultTranspose() and MatMultTransposeAdd() in
SEQAIJMKL.

Added some missing function pointer resests in
MatConvert_SeqAIJMKL_SeqAIJ.

Initial skeleton for MATSEQAIJMKL matrix type that uses MKL sparse BLAS
routines. Currently only implements mult and multadd operations, but
ready to be expanded.

Reduce Vec type checks on CUDA/CUSP/VIENNACL

Add several Vec type checks on CUDA/CUSP/VIENNACL

Add macros PetscCheckTypeName(s)

Merge remote-tracking branch 'origin/jose/cusparse-sethandle'

Merge branch 'jose/aijcusparse-fix'

* jose/aijcusparse-fix:
Fixes problems in MatMult and MatTransposeMult for matrices with empty rows when using -mat_type aijcusparse.
The fix is to initialize the

Merge branch 'jose/aijcusparse-fix'

* jose/aijcusparse-fix:
Fixes problems in MatMult and MatTransposeMult for matrices with empty rows when using -mat_type aijcusparse.
The fix is to initialize the result vector with zero, which has some performance overhead.
Arguably, this is a bug in cuSparse, which we have to circumvent here.

show more ...

Merge branch 'jose/fix-compile-cxx-cuda'

* jose/fix-compile-cxx-cuda:
Fixes compilation problems encountered when CUDA is used and --with-clanguage=cxx is set.

MatCUSPARSESetHandle: do not destroy handle to set

Merge remote-tracking branch 'origin/barry/feature-aij-subtypes'

Fix MatMult/MatMultTranspose SeqAIJCUSPARSE

Cannot call MatSetFromOptions() on inner matrices of MPIAIJ because it breaks with -mat_type mpiaij

This undoes the addition in the previous commit.
Also change some -mat_type xxx values in tests to

Cannot call MatSetFromOptions() on inner matrices of MPIAIJ because it breaks with -mat_type mpiaij

This undoes the addition in the previous commit.
Also change some -mat_type xxx values in tests to -dm_mat_type since the matrices arise from DMs

Commit-type: bug-fix
Reported-by: nightly builds
Thanks-to: Satish Balay <balay@mcs.anl.gov>
Development Tools: Vim, Emacs, Eclipse

show more ...