[petsc-dev] trying to make MKL work with GAMG

[Mark Adams]

I have having to fix AIJ methods that don't get the type from the args to
set created matrix type, so as to keep the MKL typing.

I am not sure how to fix this because I don't know to say 'make MPI from
SEQ' in this method

PetscErrorCode MatCreateMPIAIJWithSeqAIJ(MPI_Comm comm,Mat A,Mat B,const
PetscInt garray[],Mat *mat)
{
 ....
  ierr = MatCreate(comm,mat);CHKERRQ(ierr);
  ierr = MatGetSize(A,&m,&n);CHKERRQ(ierr);
....
  ierr = MatSetSizes(*mat,m,n,PETSC_DECIDE,N);CHKERRQ(ierr);
  ierr = MatSetType(*mat,MATMPIAIJ);CHKERRQ(ierr);

How should I do this?  Other places in the code I get the type from input
(MPI) args and set created matrix types accordingly. I could just put a
switch with an error for any new types that might come along in the future
...

Mark