[petsc-dev] CHOLESKY with MUMPS doesn't work properly

[Alexander Grayver]

Hi dev-team,

I have this code:

   call KSPCreate(comm3d,ksp,ierr);CHKERRQ(ierr)
   call 
KSPSetOperators(ksp,A,A,DIFFERENT_NONZERO_PATTERN,ierr);CHKERRQ(ierr)
   call KSPSetType(ksp,KSPPREONLY,ierr);CHKERRQ(ierr)
   call KSPGetPC(ksp,pc,ierr);CHKERRQ(ierr)
   call PCSetType(pc,PCCHOLESKY,ierr);CHKERRQ(ierr)
   call PCFactorSetMatSolverPackage(pc,MATSOLVERMUMPS,ierr);CHKERRQ(ierr)
   call PCFactorSetUpMatSolverPackage(pc,ierr);CHKERRQ(ierr)
   call PCFactorGetMatrix(pc,F,ierr);CHKERRQ(ierr)

which works well under petsc-3.2-p5, but produces error with petsc-dev:

[0]PETSC ERROR: --------------------- Error Message 
------------------------------------
[0]PETSC ERROR: No support for this operation for this object type!
[0]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc 
CHOLESKY!
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Development HG revision: 
a4f22f4fcb371a016e046e537076dcfd2ce5087f  HG Date: Fri Dec 09 09:08:30 
2011 -0600
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: /home/model on a openmpi-i named node233 by agrayver Fri 
Dec  9 17:23:30 2011
[0]PETSC ERROR: Libraries linked from 
/home/lib/petsc-dev/openmpi-intel-complex-release-f-ds/lib
[0]PETSC ERROR: Configure run at Fri Dec  9 16:51:10 2011
[0]PETSC ERROR: Configure options --download-metis --download-mumps 
--download-parmetis --download-superlu_dist 
--with-blacs-include=/opt/intel/Compiler/11.1/072/mkl/include 
--with-blacs-lib=/opt/intel/Compiler/11.1/072/mkl/lib/em64t/libmkl_blacs_openmpi_lp64.a 
--with-blas-lapack-dir=/opt/intel/Compiler/11.1/072/mkl/lib/em64t 
--with-debugging=0 --with-fortran-interfaces=1 --with-fortran-kernels=1 
--with-mpi-dir=/opt/mpi/intel/openmpi-1.4.2 
--with-petsc-arch=openmpi-intel-complex-release-f-ds 
--with-precision=double 
--with-scalapack-include=/opt/intel/Compiler/11.1/072/mkl/include 
--with-scalapack-lib=/opt/intel/Compiler/11.1/072/mkl/lib/em64t/libmkl_scalapack_lp64.a 
--with-scalar-type=complex --with-x=0 
PETSC_ARCH=openmpi-intel-complex-release-f-ds
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: MatGetFactor() line 3943 in 
/home/lib/petsc-dev/src/mat/interface/matrix.c
[0]PETSC ERROR: PCFactorSetUpMatSolverPackage_Factor() line 17 in 
/home/lib/petsc-dev/src/ksp/pc/impls/factor/factimpl.c
[0]PETSC ERROR: PCFactorSetUpMatSolverPackage() line 26 in 
/home/lib/petsc-dev/src/ksp/pc/impls/factor/factor.c

Any idea?

Regards,
Alexander