Hey Tim,<br> There is no parallel ILU support in petsc.<br>In parallel, the closest thing is block jacobi with ilu on each block.<br>This is the default parallel preconditioner in petsc. <br><br>If you have installed hypre, it provide a parallel ILU, via pilut.<br>
<br>Cheers,<br> Dave<br><br><br><div class="gmail_quote">On Thu, Nov 13, 2008 at 10:22 PM, Tim Kroeger <span dir="ltr"><<a href="mailto:tim.kroeger@cevis.uni-bremen.de">tim.kroeger@cevis.uni-bremen.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear PETSc team,<br>
<br>
In my application, I got this (on 4 CPUs):<br>
<br>
[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[0]PETSC ERROR: No support for this operation for this object type!<br>
[0]PETSC ERROR: Matrix type mpiaij symbolic ILU!<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 11, Thu Mar 27 09:16:38 CDT 2008 HG revision: 0a79e1469a317b0468b3b3310524bf13867c66f5<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: ./checkConvergence on a linux-gnu named inf056 by tkroger Thu Nov 13 04:06:32 2008<br>
[0]PETSC ERROR: Libraries linked from /home/collab/tkroger/packages/petsc-2.3.3-p11/lib/linux-gnu-c-debug<br>
[0]PETSC ERROR: Configure run at Fri Oct 17 06:57:08 2008<br>
[0]PETSC ERROR: Configure options --with-mpi-dir=/usr/local/lam-mpi/ --with-mpi-compilers=1 --download-f-blas-lapack=1 --with-shared=0<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: MatILUFactorSymbolic() line 4788 in src/mat/interface/matrix.c<br>
[0]PETSC ERROR: PCSetUp_ILU() line 541 in src/ksp/pc/impls/factor/ilu/ilu.c<br>
[0]PETSC ERROR: PCSetUp() line 787 in src/ksp/pc/interface/precon.c<br>
[0]PETSC ERROR: KSPSetUp() line 234 in src/ksp/ksp/interface/itfunc.c<br>
[0]PETSC ERROR: KSPSolve() line 347 in src/ksp/ksp/interface/itfunc.c<br>
[0]PETSC ERROR: User provided function() line 687 in unknowndirectory/src/numerics/petsc_linear_solver.C<br>
<br>
<br>
I wonder what that means because the AIJ matrix type should be able to do ILU preconditioning, shouldn't it? In other words: If not, which matrix type will then be able to do that?<br>
<br>
Best Regards,<br>
<br>
Tim<br><font color="#888888">
<br>
-- <br>
Dr. Tim Kroeger Phone +49-421-218-7710<br>
<a href="mailto:tim.kroeger@mevis.de" target="_blank">tim.kroeger@mevis.de</a>, <a href="mailto:tim.kroeger@cevis.uni-bremen.de" target="_blank">tim.kroeger@cevis.uni-bremen.de</a> Fax +49-421-218-4236<br>
<br>
MeVis Research GmbH, Universitaetsallee 29, 28359 Bremen, Germany<br>
<br>
Amtsgericht Bremen HRB 16222<br>
Geschaeftsfuehrer: Prof. Dr. H.-O. Peitgen<br>
<br>
</font></blockquote></div><br>