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

[Alexander Grayver]

On 09.12.2011 18:17, Matthew Knepley wrote:
> On Fri, Dec 9, 2011 at 10:31 AM, Alexander Grayver 
> <agrayver at gfz-potsdam.de <mailto:agrayver at gfz-potsdam.de>> wrote:
>
>     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:
>
>
> I believe that Cholesky only works with SBAIJ not AIJ.

Well, according to the documentation 
(http://www.mcs.anl.gov/petsc/petsc-3.2/docs/manual.pdf -- p.89) it 
works with both.
Moreover, this code works well with petsc-3.2-p5.

Regards,
Alexander

>
>    Matt
>
>     [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
>
>
>
>
> -- 
> What most experimenters take for granted before they begin their 
> experiments is infinitely more interesting than any results to which 
> their experiments lead.
> -- Norbert Wiener

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111209/305a9f32/attachment.html>