Add a few comments to indicate what some #endif directives match.

Use PETSC_ERR_PLIB, not PETSC_ERR_LIB, in MatPtAPNumeric_SeqAIJMKL_SeqAIJMKL_SymmetricReal() when input matrix not marked as symmetric.

Thanks-to: Stefano Zampini <stefano.zampini@gmail.com>

Define MatSeqAIJMKL_setup_structure_from_mkl_handle() and MatSeqAIJMKL_update_from_mkl_handle() only if PETSC_HAVE_MKL_SPARSE_SP2M_FEATURE is defined.

This change prevents warnings about functions t

Define MatSeqAIJMKL_setup_structure_from_mkl_handle() and MatSeqAIJMKL_update_from_mkl_handle() only if PETSC_HAVE_MKL_SPARSE_SP2M_FEATURE is defined.

This change prevents warnings about functions that are defined but unused
when building against an older version of MKL that has mkl_sparse_optimize()
but does not have the two stage matrix-multiplication API.

show more ...

Make error messages for MKL failures more consistent in MATSEQAIJMKL.

Rename MatPtAPNumeric_SeqAIJMKL_SeqAIJMKL() to MatPtAPNumeric_SeqAIJMKL_SeqAIJMKL_SymmetricReal() for clarity and consistency.

Style and formatting changes suggested by Barry, plus a few other minor edits.

Thanks-to: Barry Smith <bsmith@mcs.anl.gov>

Add MatSeqAIJMKL_view_mkl_handle().

Add support for new MatProductXXX() API in MATSEQAIJMKL.

This commit also:
Removes MatSeqAIJMKL_create_from_mkl_handle(), as it is no longer needed.
Adds code to handle empty sparse matrices in AIJM

Add support for new MatProductXXX() API in MATSEQAIJMKL.

This commit also:
Removes MatSeqAIJMKL_create_from_mkl_handle(), as it is no longer needed.
Adds code to handle empty sparse matrices in AIJMKL; mostly needed for MatProductXXX() routines, but also for MatMatMultNumeric.
Makes various cleanup, formatting, and typo fixes.

For helpful comments and suggestions,
Thanks-to: Stefano Zampini <stefano.zampini@gmail.com>
Thanks-to: Barry Smith <bsmith@mcs.anl.gov>
Thanks-to: Hong Zhang <hzhang@mcs.anl.gov>

show more ...

Add MatMatMultSymbolic_SeqAIJMKL_SeqAIJMKL_SpMV2().

Merge remote-tracking branch 'origin/maint'

Merge branch 'rmills/fix-matptapnumeric-mpiaij-objectstate/maint' into 'maint'

Bug fix: MatPtAPNumeric_MPIAIJ_MPIAIJ() must increase ObjectState of AP_loc.

See merge request petsc/petsc!2943

Fix logic error in MatSeqAIJMKL_create_mkl_handle().

Bug fix: MatPtAPNumeric_MPIAIJ_MPIAIJ() must increase ObjectState of AP_loc.

MUMPS: set default ICNTL(7) per mumps manual

Most external factorization packages do not use the ordering provided by PETSc, therefor only provide it when needed.

Adds MatOrderingType external to indicate not to generate an ordering and use wh

Most external factorization packages do not use the ordering provided by PETSc, therefor only provide it when needed.

Adds MatOrderingType external to indicate not to generate an ordering and use what the package needs
Now -pc_view will not print the wrong PETSc ordering when the ordering is done externally

Saves some memory and compute time.

Commit-type: optimization
Time: 1.2 hours
Reported-by: Junchao Zhang <junchao.zhang@gmail.com>
Thanks-to: Stefano Zampini <stefano.zampini@gmail.com>

show more ...

Merge branch 'stefanozampini/feature-cusparse-notransgen-herm' into 'master'

feature cusparse notransgen herm

See merge request petsc/petsc!2915

PCHPDDM: add ops->matapply

MATSEQAIJCUSPARSE: test multiple transpose generation flags

MATSEQAIJCUSPARSE: log transpose generation

MATSEQAIJCUSPARSE: remove wrong CPUtoGPU memory logging

MatMultHermitianTranspose: use MatTranspose directly for real builds

MATSEQAIJCUSPARSE: implement hermitian transpose operation

Simplify the class, using a single routine for all the MatMult variants
Make transpose generation optional (not default)

MUMPS: support distributed rhs

MATSCALAPACK: new Mat type that wraps a ScaLAPACK matrix

Merge remote-tracking branch 'origin/maint'