On Tue, Feb 1, 2011 at 9:03 AM, Michael E Henderson <span dir="ltr"><<a href="mailto:mhender@us.ibm.com">mhender@us.ibm.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<font size="2" face="sans-serif">Good morning,</font>
<br>
<br><font size="2" face="sans-serif">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>
<br></blockquote><div><br></div><div>This is the joy of using other packages. Looking at the source, I see 4 packages which can factor a parallel (MPIAIJ) matrix:</div><div><br></div><div> MUMPS, SuperLU_dist, Spooles (now unsupported), Pastix</div>
<div><br></div><div>These are all usable from -pc_factor_mat_solver_package when using MPIAIJ.</div><div><br></div><div>1) We do not consider SPAI a matrix factorization package. You just use it with --pc_type spai</div><div>
<br></div><div>2) We cannot see inside Hypre, and thus it is hard to get into this framework. You use Euclid with</div><div><br></div><div> -pc_type hypre -pc_hypre_type euclid</div><div><br></div><div> Matt</div><div>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><font size="2" face="sans-serif">With options:</font>
<br>
<br><font size="2" face="sans-serif">-ksp_type gmres</font>
<br><font size="2" face="sans-serif">-pc_type ilu</font>
<br><font size="2" face="sans-serif">-pc_factor_levels 10</font>
<br><font size="2" face="sans-serif">-pc_factor_fill 10</font>
<br><font size="2" face="sans-serif">-pc_factor_mat_solver_package spai</font>
<br>
<br><font size="2" face="sans-serif">I get the message:</font>
<br>
<br><font size="2" face="sans-serif"> 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<br>
</font>
<br><font size="2" face="sans-serif">I checked the configuration output and
spai was indeed configured and built. I also tried spooles with a similar
result.</font>
<br>
<br><font size="2" face="sans-serif">The table </font><a href="http://www.mcs.anl.gov/petsc/petsc-as/documentation/linearsolvertable.html" target="_blank"><font size="2" face="sans-serif">http://www.mcs.anl.gov/petsc/petsc-as/documentation/linearsolvertable.html</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>
<br>
<br><font size="2" face="sans-serif">-pc_factor_mat_solver_package hypre
</font>
<br><font size="2" face="sans-serif">-pc_hypre_type euclid</font>
<br>
<br><font size="2" face="sans-serif">also gives</font>
<br>
<br><font size="2" face="sans-serif">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>
<br>
<br><font size="2" face="sans-serif">I'm using hypre as a preconditioer elsewhere,
so I'm sure it's installed. Am I doing something obviously wrong?</font>
<br>
<br><font size="2" face="sans-serif">Thanks,</font>
<br>
<br><font size="2" face="sans-serif">Mike Henderson</font>
<br><font size="2" face="sans-serif">------------------------------------------------------------------------------------------------------------------------------------<br>
Mathematical Sciences, TJ Watson Research Center<br>
<a href="mailto:mhender@watson.ibm.com" target="_blank">mhender@watson.ibm.com</a><br>
</font><a href="http://www.research.ibm.com/people/h/henderson/" target="_blank"><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/" target="_blank"><font size="2" face="sans-serif">http://multifario.sourceforge.net/</font></a><font size="2" face="sans-serif"><br>
</font></blockquote></div><br><br clear="all"><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<br>