Merge remote-tracking branch 'origin/hongzh/copy_l2g_stencil'

Merge branch 'maint'

fortran: fix macro redefinition error

Copy matrix block sizes in MatConvert

The local to global mapping block size must match the matrix block size.

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

Merge remote-tracking branch 'origin/tristan/mat-invert-block-diagonal-mat'

Merge remote-tracking branch 'origin/barry/use-mpi-datatype-F90Array'

Fix MatConvert error when used with hypre

SNES ex56 gives an error when ran with -pc_type hypre -mat_block_size 3. It is because MatSetLocalToGlobalMapping checks the block size and errors out when

Fix MatConvert error when used with hypre

SNES ex56 gives an error when ran with -pc_type hypre -mat_block_size 3. It is because MatSetLocalToGlobalMapping checks the block size and errors out when *M and mat have different values of bs. However it is unnecessary to impose same bs for MatConvert.

Funded-by:
Project: PETSc for Theta
Time: 0.5 hour
Reported-by:
Thanks-to:

show more ...

Merge remote-tracking branch 'origin/barry/fix-matsettype-forsubclasses'

Merge branch 'maint'

Merge remote-tracking branch 'origin/stefano_zampini/cherrypick-commits-for-maint-from-continuousadjoint' into maint

Rename ELL to SELL

SELL is more accurate. It stands for sliced ELLPACK.

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

Merge branch 'maint'

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

Merge remote-tracking branch 'origin/maint'

MatMFFD: remove unneeded MatSetUp_MFFD

remove unneeded guards in MatSetUp and always call PetscLayoutSetUp

Thanks-to: Jed Brown jed@jedbrown.org

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

MatCopy: return if A == B

Mat: implement a fallback when MatFindNonzeroRows_XXX is missing

Fix typos in man pages

F90Array routines now take a MPI_Datatype argument instead of PetscDataType

This is part of the refactorization to limit and maybe eventually eliminate PetscDataType since
it duplicates MPI_Datatype

F90Array routines now take a MPI_Datatype argument instead of PetscDataType

This is part of the refactorization to limit and maybe eventually eliminate PetscDataType since
it duplicates MPI_Datatype functionality (poorly).

Commit-type: style-fix
Thanks-to: Jed Brown <jed@jedbrown.org>

show more ...

MatInvertBlockDiagonalMat: error checking and typo fix

Copy stencil info and localtoglobalmapping in MatConvert

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

Merge remote-tracking branch 'origin/hongzh/add-ell-format'

fixed bug in MatGetRowSum

MatSetType() for matrix subclasses should call MatConvert() instead of recreating matrix

otherwise all the preallocation and other information provided to the matrix is lost.
For example for aijperm

MatSetType() for matrix subclasses should call MatConvert() instead of recreating matrix

otherwise all the preallocation and other information provided to the matrix is lost.
For example for aijperm or aijmkl

Commit-type: feature
Funded-by: ECP
Time: 2 hours
Reported-by: Richard Tran Mills <rtmills@anl.gov>

show more ...