rm -matptap_viatranspose and -matptap_viamatmatmatmult from MatPtAP()

MatResetPreallocation: use "PetscUseMethod"

instead of "PetscTryMethod" so that users can
be informed when using an unsupported mat type.

check mat->hermitian for MAT_FACTOR_CHOLESKY and ftype == MAT_FACTOR_ICC in MatGetFactor()
cleanup src/mat/examples/tests/ex127.c

MatResetPreallocation: added "()" in comments

MatResetPreallocation: Reset mat to use the original nonzero pattern provided by users.

MPIELL seems to work

Funded-by:
Project: PETSc for Theta
Time: 14 hours
Reported-by:
Thanks-to:

Add seqell format and its matmul kernel

-Make memkind error out when HBW is not usable

Funded-by:
Project: PETSc for Theta
Time: 52 hours
Reported-by:
Thanks-to:

Mat: MatInvertBlockDiagonalMat. KSP: Block diagonal Schur inverse.

Added new matrix classes to include files
Added define MKL_SUPPORTS_BAIJ_ZERO_BASED to suuport different versions
of MKL (older verions support only one based indexing for column major
blocks in baij

Added new matrix classes to include files
Added define MKL_SUPPORTS_BAIJ_ZERO_BASED to suuport different versions
of MKL (older verions support only one based indexing for column major
blocks in baij format)

show more ...

Merge branch 'rmills/add-aijmkl'

* rmills/add-aijmkl:
Adds AIJMKL matrix class to facilitate use of MKL sparse BLAS routines
BitBucket Pull Request 737
This is the initial version of the AIJMKL matr

Merge branch 'rmills/add-aijmkl'

* rmills/add-aijmkl:
Adds AIJMKL matrix class to facilitate use of MKL sparse BLAS routines
BitBucket Pull Request 737
This is the initial version of the AIJMKL matrix class,
a subclass of AIJ that has the goal of using MKL sparse BLAS routines whenever practicable.
If the MKL version is new enough, the sparse inspector-executor
("SpMV2" -- terrible name, but that's what the MKL team called it) routines are used by default.
Currently only MatMult, MatMultAdd, MatMultTranspose, and MatMultTransposeAdd operations are supported,
but more will be added once this initial attempt is at a point
that it provides an acceptable framework for adding additional further functionality.

show more ...

Merge branch 'maint'

Fix error message for MatCholeskyFactor()
Reported-by: Issue #177

Clarified documentation for MatDuplicate() and MatDuplicateOption.

Added clarifications cross links for MatFactorSetSchurIS() and its related routines

Commit-type: documentation
Funded-by: ECP
Time: .45 hours

Development Tools: Vim, Emacs, Eclipse

Merge remote-tracking branch 'origin/jed/matgetrowsum-ones'

fix mem leak as reported by nightly tests

for baij and dummy matrices

add MATOP_DESTROYSUBMATRICES to cleanup MatDestroySubMatrices()
Reported-by: Stefano Zampini <stefano.zampini@gmail.com>

MatGetRowSum: use A*ones to avoid depending on MatGetRow()

This allows use of Jacobi with rowsum (usually desirable for solving
with high-order mass matrices) without needing to implement anything.

MatGetRowSum: use A*ones to avoid depending on MatGetRow()

This allows use of Jacobi with rowsum (usually desirable for solving
with high-order mass matrices) without needing to implement anything.
This can and probably should be extended to allow a MatOp so that
communication can be avoided.

show more ...

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

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

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

Registered PETSc built-in factorization solvers for MATSEQAIJMKL.

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

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