<font size=2 face="sans-serif">Thanks</font>
<br><font size=2 face="sans-serif"><br>
------------------------------------------------------------------------------------------------------------------------------------<br>
Mathematical Sciences, TJ Watson Research Center<br>
mhender@watson.ibm.com<br>
</font><a href=http://www.research.ibm.com/people/h/henderson/><font size=2 face="sans-serif">http://www.research.ibm.com/people/h/henderson/</font></a><font size=2 face="sans-serif"><br>
</font><a href=http://multifario.sourceforge.net/><font size=2 face="sans-serif">http://multifario.sourceforge.net/</font></a><font size=2 face="sans-serif"><br>
</font>
<br>
<br>
<br>
<br><font size=1 color=#5f5f5f face="sans-serif">From:
</font><font size=1 face="sans-serif">Matthew Knepley <knepley@gmail.com></font>
<br><font size=1 color=#5f5f5f face="sans-serif">To:
</font><font size=1 face="sans-serif">PETSc users list <petsc-users@mcs.anl.gov></font>
<br><font size=1 color=#5f5f5f face="sans-serif">Date:
</font><font size=1 face="sans-serif">02/01/2011 10:44 AM</font>
<br><font size=1 color=#5f5f5f face="sans-serif">Subject:
</font><font size=1 face="sans-serif">Re: [petsc-users]
options for gmres w/ ilu precondioner for sparse</font>
<br><font size=1 color=#5f5f5f face="sans-serif">Sent by:
</font><font size=1 face="sans-serif">petsc-users-bounces@mcs.anl.gov</font>
<br>
<hr noshade>
<br>
<br>
<br><font size=3>On Tue, Feb 1, 2011 at 9:03 AM, Michael E Henderson <</font><a href=mailto:mhender@us.ibm.com><font size=3 color=blue><u>mhender@us.ibm.com</u></font></a><font size=3>>
wrote:</font>
<br><font size=2 face="sans-serif">Good morning,</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
I'm trying to use gmres with an ilu preconditioner and having trouble getting
the options right. I figure it's got to be something simple, so hope it's
an easy question.</font><font size=3> </font>
<br>
<br><font size=3>This is the joy of using other packages. Looking at the
source, I see 4 packages which can factor a parallel (MPIAIJ) matrix:</font>
<br>
<br><font size=3> MUMPS, SuperLU_dist, Spooles (now unsupported),
Pastix</font>
<br>
<br><font size=3>These are all usable from -pc_factor_mat_solver_package
when using MPIAIJ.</font>
<br>
<br><font size=3>1) We do not consider SPAI a matrix factorization package.
You just use it with --pc_type spai</font>
<br>
<br><font size=3>2) We cannot see inside Hypre, and thus it is hard to
get into this framework. You use Euclid with</font>
<br>
<br><font size=3> -pc_type hypre -pc_hypre_type euclid</font>
<br>
<br><font size=3> Matt</font>
<br><font size=3> </font>
<br><font size=2 face="sans-serif">With options:</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
-ksp_type gmres</font><font size=3> </font><font size=2 face="sans-serif"><br>
-pc_type ilu</font><font size=3> </font><font size=2 face="sans-serif"><br>
-pc_factor_levels 10</font><font size=3> </font><font size=2 face="sans-serif"><br>
-pc_factor_fill 10</font><font size=3> </font><font size=2 face="sans-serif"><br>
-pc_factor_mat_solver_package spai</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
I get the message:</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
unknown: [1|MatGetFactor() line 3646 in src/mat/interface/matrix.c:
Matrix format mpiaij does not have a solver spai. Perhaps you must config/configure.py
with --download-spai</font><font size=3><br>
</font><font size=2 face="sans-serif"><br>
I checked the configuration output and spai was indeed configured and built.
I also tried spooles with a similar result.</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
The table </font><a href="http://www.mcs.anl.gov/petsc/petsc-as/documentation/linearsolvertable.html" target=_blank><font size=2 color=blue face="sans-serif"><u>http://www.mcs.anl.gov/petsc/petsc-as/documentation/linearsolvertable.html</u></font></a><font size=2 face="sans-serif">
seems to be saying that only hypre/euclid can be used for ilu(k) w/ aij.
Is that true?</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
-pc_factor_mat_solver_package hypre <br>
-pc_hypre_type euclid</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
also gives</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
unknown: [1MatGetFactor() line 3646 in src/mat/interface/matrix.c: Matrix
format mpiaij does not have a solver hypre. Perhaps you must config/configure.py
with --download-hypre </font><font size=3><br>
</font><font size=2 face="sans-serif"><br>
I'm using hypre as a preconditioer elsewhere, so I'm sure it's installed.
Am I doing something obviously wrong?</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
Thanks,</font><font size=3> <br>
</font><font size=2 face="sans-serif"><br>
Mike Henderson</font><font size=3> </font><font size=2 face="sans-serif"><br>
------------------------------------------------------------------------------------------------------------------------------------<br>
Mathematical Sciences, TJ Watson Research Center</font><font size=2 color=blue face="sans-serif"><u><br>
</u></font><a href=mailto:mhender@watson.ibm.com target=_blank><font size=2 color=blue face="sans-serif"><u>mhender@watson.ibm.com</u></font></a><font size=3 color=blue><u><br>
</u></font><a href=http://www.research.ibm.com/people/h/henderson/ target=_blank><font size=2 color=blue face="sans-serif"><u>http://www.research.ibm.com/people/h/henderson/</u></font></a><font size=3 color=blue><u><br>
</u></font><a href=http://multifario.sourceforge.net/ target=_blank><font size=2 color=blue face="sans-serif"><u>http://multifario.sourceforge.net/</u></font></a>
<br><font size=3><br>
<br>
<br>
-- <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</font>
<br>