<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Dec 4, 2014 at 1:00 PM, Raphaël Couturier <span dir="ltr"><<a href="mailto:raphael.couturier@univ-fcomte.fr" target="_blank">raphael.couturier@univ-fcomte.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Matt,<br>
<br>
Thank you for all your work on PETSc.<br>
I would like to know if you or your colleagues have some examples of code with linear systems in PETSc that are difficult to solve with more of less all preconditioners and GMRES?<br></blockquote><div><br></div><div>I am Cc'ing petsc-dev since better feedback might come from there. ex15 is the p-Laplacian, ex16</div><div>is large deformation elasticity (ex54 is infinitesimal). These are somewhat difficult, but most of the</div><div>difficulty is getting MG to work right. These are likely candidates to be solved with FMG, neglecting</div><div>Krylov.</div><div><br></div><div>In the past, hard things to solve have generally been very non-normal matrices. We do not maintain</div><div>a catalogue since we (I at least) have almost given up on black box solvers in favor of focusing</div><div>on preconditioners for certain problem types. However, other people on here might have specific</div><div>problems that they need solved.</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">
In fact, I am working on a new method which solves, for example ex15 or ex54 of ksp/ksp/examples/tutorials with large instances (for example 200,000,000 components or more) with the mg preconditioner and for example 8192 cores, faster than GMRES.<br>
So I would like to know if it is possible to test other difficult examples?<br>
<br>
I have submitted my work to a conference, and as soon as this work is accepted, I will be able to give more details on the method.<br>
<br>
Best wishes,<br>
<br>
Raphaël Couturier<br>
<a href="http://members.femto-st.fr/raphael-couturier/" target="_blank">http://members.femto-st.fr/<u></u>raphael-couturier/</a><br>
<br>
<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="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>