update #define conditions used to turn on Intel intrinsics

Commit-type: portability-fix, style-fix

patches to get it to work with gcc and clang and any non-intel compiler

Commit-type: bug-fix, portability-fix

Remove unused mdiag[] array from BAIJ that was filled but never used

Commit-type: style-fix
Reported-by: Xiangdong <epscodes@gmail.com>

Added case specific intrinsic kernels for VecMdot_Seq, VecMAXPY_Seq, and BAIJ bocksize 9 by 9

Fix several minor code issues in AIJMKL and BAIJMKL code.

Remove improper returning of error codes from PetscFunctionReturn();
remove unused variable; fix some code formatting.

Commit-type: style-f

Fix several minor code issues in AIJMKL and BAIJMKL code.

Remove improper returning of error codes from PetscFunctionReturn();
remove unused variable; fix some code formatting.

Commit-type: style-fix
Reported-by: Barry Smith <bsmith@mcs.anl.gov>

show more ...

Merge remote-tracking branch 'origin/balay/fix-mkl-sparse/maint' into maint

mkl: add mkl_sparse.py as a separate package [as the interface code is organized this way]
and use PETSC_HAVE_MKL_SPARSE to enable/disable this interface code [so that this flag
represents functional

mkl: add mkl_sparse.py as a separate package [as the interface code is organized this way]
and use PETSC_HAVE_MKL_SPARSE to enable/disable this interface code [so that this flag
represents functionality that is more than PETSC_HAVE_MKL flag].

Also disable mkl_sparse, mkl_sparse_optimize for 64-bit-incides

Reported-by: Victor Eijkhout <eijkhout@tacc.utexas.edu>

show more ...

MatConvert: Fix for {Seq|MPI}BAIJ -> {Seq|MPI}AIJ with MAT_REUSE_MATRIX

Bug fixes
Fix for test 107_2. To prevent MKL routine returning error in case when MatSetOption(mat,MAT_STRUCTURE_ONLY,PETSC_TRUE) is called. In this case matrix will be optimized before call to mkl_s

Bug fixes
Fix for test 107_2. To prevent MKL routine returning error in case when MatSetOption(mat,MAT_STRUCTURE_ONLY,PETSC_TRUE) is called. In this case matrix will be optimized before call to mkl_sparse_x_mv.
The same problem will appear in AIJMKl case.
Also added MatCreateBAIJMKL and MatCreateSeqBAIJMKL to external include to call them from examples

show more ...

Minor changes for better code style

Bug fix for MatDuplicate_SeqBAIJMKL
Null pointer was passed to PetscMemcpy

Bug fixes
Assign function poiners on AssemblyEnd so they are not overridden in MatSeqBAIJSetPreallocation_SeqBAIJ()

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 ...

Implementation of baijmkl class that supports mult
and multadd operations using MKL Sparse BLAS IE.
If used version of MKL does not support new API (version 11.3 and
earlier) default PETSc kernels ar

Implementation of baijmkl class that supports mult
and multadd operations using MKL Sparse BLAS IE.
If used version of MKL does not support new API (version 11.3 and
earlier) default PETSc kernels are called

show more ...

for baij and dummy matrices

rename