Merge branch 'maint'

Fixed mkl_sparse_sypr() calls that were using an unsupported value in the matrix descriptor.

Commit-type: bug-fix
Reported-by: Maria Zhukova <maria.zhukova@intel.com>

Merge remote-tracking branch 'origin/rmills/fix-aijmkl-replace-ifdefs'

Merge remote-tracking branch 'origin/jed/const-XType'

Removed use of #ifdef/#ifndef in favor of #if defined()/#if !defined().

Commit-type: style-fix

Remove "const MatType" and similar (const is already inside the typedef)

MatIS: added implementation for PtAP operation

supported P are {Seq|MPI}(B)AIJ and HYPRE
Interface for {MPI|Seq}BAIJ and HYPRE code works, however:
- local PtAP ops are missing for BAIJ
- MatCreateS

MatIS: added implementation for PtAP operation

supported P are {Seq|MPI}(B)AIJ and HYPRE
Interface for {MPI|Seq}BAIJ and HYPRE code works, however:
- local PtAP ops are missing for BAIJ
- MatCreateSubMatr* ops are missing for HYPRE (will be added in a separate branch)

If requested, carry over the l2l representation of P by attaching it to the resulting mat
This is needed by the MFEM constructor for BDDC

added few more tests to ex23.c
L2L storage is disabled because of a bug in MatDestroySubMatrices that will be fixed in a separate branch

show more ...

Old sparse BLAS interfaces are now not referenced by MKL versions that deprecate them.

The old interfaces were deprecated in MKL version 18, update 2.
The inspector-executor routines can still be "t

Old sparse BLAS interfaces are now not referenced by MKL versions that deprecate them.

The old interfaces were deprecated in MKL version 18, update 2.
The inspector-executor routines can still be "turned off" by the user, but
in version 18, update 2 and onwards, the new inspector-executor interface is
used, with mkl_sparse_optimize() not called if the user does not want the
inspector-executor functionality.

show more ...

Old sparse BLAS interfaces are now not referenced by MKL versions that deprecate them.

The old interfaces were deprecated in MKL version 18, update 2.
The inspector-executor routines can still be "t

Old sparse BLAS interfaces are now not referenced by MKL versions that deprecate them.

The old interfaces were deprecated in MKL version 18, update 2.
The inspector-executor routines can still be "turned off" by the user, but
in version 18, update 2 and onwards, the new inspector-executor interface is
used, with mkl_sparse_optimize() not called if the user does not want the
inspector-executor functionality.

show more ...

Man pages: add newlines after "Notes:"

This allows for proper formatting from sowing.

On OS X (using gsed, not the default BSD sed), from the PETSc root directory:

find src include -type f \(

Man pages: add newlines after "Notes:"

This allows for proper formatting from sowing.

On OS X (using gsed, not the default BSD sed), from the PETSc root directory:

find src include -type f \( -name "*.c" -or -name "*.h" -or -name "*.cxx" \) | xargs gsed -i 's/Notes\s*:\s*\(\w.*\)/Notes:\n \1/'

This adds a newline and 4 spaces whenever "Notes:" is followed by any "word" character, in any .c, .h, or .cxx file in src/ or include/

show more ...

Add more protection to intrinsics

Funded-by:
Project: PETSc for Theta
Time: 0.5 hours
Reported-by: "Nourgaliev, Robert Nr" <nourgaliev1@llnl.gov>
Thanks-to:

Merge remote-tracking branch 'origin/hzhang/superlu_dist-matinertia-hermitian'

Updated manual page for MatCreateSeqAIJMKL().

Merge remote-tracking branch 'origin/rmills/feature-aijmkl-add-matptap-symmetric'

add MatGetInertia_SuperLU_DIST for complex Hermitian matrix

Improved MKL error checking in MatSeqAIJMKL_create_from_mkl_handle().

Removed functions such as MatScale_SeqAIJMKL(); no longer needed because we track ObjectState.

Merge branch 'rmills/vectorize-aij-avx512-only'

Merge remote-tracking branch 'origin/hzhang/vec_node/master-rebased-v2'

Fixed bugs preventing proper tracking of object state.

Added MatPtAP() implementations that use MKL routines optimized for symmetric matrices.

Also fixed bugs affecting MatMatMultNumeric() and MatPtAPNumeric().

Merge branch 'karlrupp/feature-remove-cusp'

* karlrupp/feature-remove-cusp:
Removes VECCUSP vectors and AIJCUSP matrices, because they provide no added benefit over VECCUDA and AIJCUSPARSE.
Also rem

Merge branch 'karlrupp/feature-remove-cusp'

* karlrupp/feature-remove-cusp:
Removes VECCUSP vectors and AIJCUSP matrices, because they provide no added benefit over VECCUDA and AIJCUSPARSE.
Also removes the CUSP-based preconditioners for now.
SACUSP will be re-added later, after the new release is out.

show more ...

Merge branch 'hzhang/superlu_dist-matinertia-fix'

Adjusted whitespace in some variable declarations.

CUDA: Fixed MatMultTranspose_MPIAIJCUSPARSE()

This was a broken copy&paste from MatMult_MPIAIJCUSPARSE()
without adjusting the VecScatters.