<div dir="auto">If you do this, please please please also update NGMRES since it looks like you would do a similar linear least squares thing with the directions. It should be factored out. SLEPc has a nice TSQR that we could pull into PETSc for this.<div dir="auto"><br></div><div dir="auto">   Maty</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Feb 28, 2018 14:36, "Jed Brown" <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Nicole has expressed interested in helping with a PETSc implementation,<br>
she just doesn't have much experience with PETSc yet and I haven't had<br>
time to prioritize doing it myself.  If you want to develop a PETSc<br>
implementation, I would suggest reaching out to her and Cc'ing me.<br>
<br>
"Smith, Barry F." <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> writes:<br>
<br>
>       Tried to read the papers, couldn't follow, them but my guess is you need to copy the PETSc CG KSP routines and rework them as a new KSP type for these algorithms.<br>
><br>
>    Barry<br>
><br>
><br>
>> On Feb 28, 2018, at 3:39 AM, Karin&NiKo <<a href="mailto:niko.karin@gmail.com">niko.karin@gmail.com</a>> wrote:<br>
>><br>
>> Dear PETSc team,<br>
>><br>
>> I would like to experiment multi-preconditioned Krylov methods, as presented in the paper or Bridson and Greif (<a href="https://www.cs.ubc.ca/~rbridson/mpcg/" rel="noreferrer" target="_blank">https://www.cs.ubc.ca/~<wbr>rbridson/mpcg/</a>) and more specifically in a context of DD like in the paper of Spillane (<a href="https://hal.archives-ouvertes.fr/hal-01170059/document" rel="noreferrer" target="_blank">https://hal.archives-<wbr>ouvertes.fr/hal-01170059/<wbr>document</a>) or Gosselet (<a href="https://hal.archives-ouvertes.fr/hal-01056928/document" rel="noreferrer" target="_blank">https://hal.archives-<wbr>ouvertes.fr/hal-01056928/<wbr>document</a>).<br>
>><br>
>> It seems to me this is not a basically supported feature of PETSc.<br>
>><br>
>> My question :<br>
>> Is there a way of doing it by playing with the runtime options of PETSc or do I have to implement it in the core of the KSP object? If so, what would be the correct design?<br>
>><br>
>> Thanks,<br>
>> Nicolas<br>
>><br>
</blockquote></div></div>