<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
On 09.12.2011 18:17, Matthew Knepley wrote:
<blockquote
cite="mid:CAMYG4Gk_9ru9QF+y5okuN9H_pgT7TWcgYn6DQRwNPwE94Tx4fA@mail.gmail.com"
type="cite">On Fri, Dec 9, 2011 at 10:31 AM, Alexander Grayver <span
dir="ltr"><<a moz-do-not-send="true"
href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>></span>
wrote:<br>
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
0.8ex; border-left: 1px solid rgb(204, 204, 204);
padding-left: 1ex;">
Hi dev-team,<br>
<br>
I have this code:<br>
<br>
call KSPCreate(comm3d,ksp,ierr);CHKERRQ(ierr)<br>
call KSPSetOperators(ksp,A,A,DIFFERENT_NONZERO_PATTERN,ierr);CHKERRQ(ierr)<br>
call KSPSetType(ksp,KSPPREONLY,ierr);CHKERRQ(ierr)<br>
call KSPGetPC(ksp,pc,ierr);CHKERRQ(ierr)<br>
call PCSetType(pc,PCCHOLESKY,ierr);CHKERRQ(ierr)<br>
call PCFactorSetMatSolverPackage(pc,MATSOLVERMUMPS,ierr);CHKERRQ(ierr)<br>
call PCFactorSetUpMatSolverPackage(pc,ierr);CHKERRQ(ierr)<br>
call PCFactorGetMatrix(pc,F,ierr);CHKERRQ(ierr)<br>
<br>
which works well under petsc-3.2-p5, but produces error with
petsc-dev:<br>
</blockquote>
<div><br>
</div>
<div>I believe that Cholesky only works with SBAIJ not AIJ.</div>
</div>
</blockquote>
<br>
Well, according to the documentation
(<a class="moz-txt-link-freetext" href="http://www.mcs.anl.gov/petsc/petsc-3.2/docs/manual.pdf">http://www.mcs.anl.gov/petsc/petsc-3.2/docs/manual.pdf</a> -- p.89) it
works with both. <br>
Moreover, this code works well with petsc-3.2-p5.<br>
<br>
Regards,<br>
Alexander<br>
<br>
<blockquote
cite="mid:CAMYG4Gk_9ru9QF+y5okuN9H_pgT7TWcgYn6DQRwNPwE94Tx4fA@mail.gmail.com"
type="cite">
<div class="gmail_quote">
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
0.8ex; border-left: 1px solid rgb(204, 204, 204);
padding-left: 1ex;">
[0]PETSC ERROR: --------------------- Error Message
------------------------------------<br>
[0]PETSC ERROR: No support for this operation for this object
type!<br>
[0]PETSC ERROR: Matrix format mpiaij does not have a built-in
PETSc CHOLESKY!<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Development HG revision:
a4f22f4fcb371a016e046e537076dcfd2ce5087f HG Date: Fri Dec 09
09:08:30 2011 -0600<br>
[0]PETSC ERROR: See docs/changes/index.html for recent
updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about trouble
shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual pages.<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: /home/model on a openmpi-i named node233 by
agrayver Fri Dec 9 17:23:30 2011<br>
[0]PETSC ERROR: Libraries linked from
/home/lib/petsc-dev/openmpi-intel-complex-release-f-ds/lib<br>
[0]PETSC ERROR: Configure run at Fri Dec 9 16:51:10 2011<br>
[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<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: MatGetFactor() line 3943 in
/home/lib/petsc-dev/src/mat/interface/matrix.c<br>
[0]PETSC ERROR: PCFactorSetUpMatSolverPackage_Factor() line 17
in /home/lib/petsc-dev/src/ksp/pc/impls/factor/factimpl.c<br>
[0]PETSC ERROR: PCFactorSetUpMatSolverPackage() line 26 in
/home/lib/petsc-dev/src/ksp/pc/impls/factor/factor.c<br>
<br>
Any idea?<br>
<br>
Regards,<br>
Alexander<br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to
which their experiments lead.<br>
-- Norbert Wiener<br>
</blockquote>
<br>
</body>
</html>