<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Feb 7, 2017 at 7:53 AM, Michel Kern <span dir="ltr"><<a href="mailto:michel.kern@inria.fr" target="_blank">michel.kern@inria.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:#000000"><div>Thanks Mat,</div><div><br></div><div>according tot he hypre documentation, pilut is "no longer supported by the hypre team", as less efficient than Euclid, so I wouldn't consider an option. I wsa looking for a "natural" upgrade path.</div></div></div></blockquote><div><br></div><div>I don't know of another ILU.</div><div><br></div><div>I have never understood the use of ILU. It is a fragile, not very scalable algorithm which is soundly beaten in every</div><div>arena where we understand the problem. I also do not understand why it is used as part of CPR. I can't believe that</div><div>is the best way to solve the problem. But maybe I do not understand.</div><div><br></div><div>  Thanks,</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"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:#000000"><div>The problems we've looked at for the moment can be up to a few millions cells (that's the size of the pressure block, the saturation block is larger, but should be easier to solve). We're aiming for 10^6 cells, and 10^3 cores.</div><div>Yes, direct methods deserve a second look !</div><div><br></div><div>Michel</div><div><br></div><hr id="m_-8590045050658432464zwchr"><blockquote style="border-left:2px solid #1010ff;margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Feb 7, 2017 at 2:49 AM, Michel Kern <span dir="ltr"><<a href="mailto:michel.kern@inria.fr" target="_blank">michel.kern@inria.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:#000000"><div>I have a code that uses euclid for ILU preconditioning (as part of a CPR-AMG preconditioner for a mulitphase flow solver). As noted in the release notes, euclid is no longer supported starting from 3.6. What is the recommended replacement ? Is block Jacobi plus a local solver an alternative ? </div></div></div></blockquote><div><br></div><div>There is pilut from Hypre.</div><div><br></div><div>How big is the problem? Did you try sparse direct on it (SuperLU, MUMPS)?</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:#000000"><div>Thanks</div><div>Michel</div><div><br></div><div>-- <br></div><div><span></span><div><span style="font-family:arial,helvetica,sans-serif">Michel Kern</span></div><div><span style="font-size:12pt"><span style="font-family:arial,helvetica,sans-serif">Centre de Recherche Inria de Paris,    Maison de la Simulation</span></span></div><div><span style="font-family:arial,helvetica,sans-serif">2 rue Simone Iff, bureau A413,            Digiteo Labs, bât 565, <span style="font-size:12pt">CEA Saclay</span><span style="font-size:12pt"> </span></span></div><div><span style="font-family:arial,helvetica,sans-serif">75012 Paris                                          91191 Gif-sur-Yvette cedex </span></div><div><span style="font-family:arial,helvetica,sans-serif">tel. : <a href="tel:+33%201%2080%2049%2042%2036" target="_blank">+33 (0)1 80 49 42 36</a>, <span class="m_-8590045050658432464m_-8652995892514175014Apple-tab-span" style="white-space:pre-wrap"> </span>                tel. : <a href="tel:+33%201%2069%2008%2058%2032" target="_blank">+33 (0)1 69 08 58 32</a>            </span></div><div style="font-family:arial,helvetica,sans-serif;color:#000000;font-size:12pt"><br></div><span></span><br></div></div></div></blockquote></div><br><br clear="all"><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_-8590045050658432464gmail_signature">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</div></font></span></div></div></blockquote><span class="HOEnZb"><font color="#888888"><div><br><br></div><div><br></div><div>-- <br></div><div><span name="x"></span>Michel Kern<br>Centre de recherche Inria de Paris,  Maison de la Simulation, <br>2 rue Simone Iff, bureau A413,       Digiteo Labs, bât. 565, CEA Saclay   <br>75012 Paris                          91191 Gif-sur-Yvette Cedex <br>tel. : <a href="tel:+33%201%2080%2049%2042%2036" value="+33180494236" target="_blank">+33 (0)1 80 49 42 36</a>,         tel. : <a href="tel:+33%201%2069%2008%2058%2032" value="+33169085832" target="_blank">+33 (0)1 69 08 58 32</a>            <span name="x"></span><br></div></font></span></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">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</div>
</div></div>