<div dir="ltr"><div><div>Thanks for the answers,<br><br></div>Jed, PCPG stands for Preconditioned Conjugate Projected Gradient. Since any FETI literature suggests PCPG, I am planning to go stick with it. <br>My plan is to extend your KSPCG algorithm with optional application of projection space and re-orthogonalizations.<br><br></div>Best regards,<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Dec 29, 2014 at 2:33 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  Take a look at KSPFischerGuessCreate() and the material it points to. From the command line you can run for example<br>
<br>
-ksp_fischer_guess 1,20<br>
<br>
This method "works" by saving information about Krylov directions and then projecting those directions out of the NEXT linear solve at the beginning of the new linear solve (constructing a "better" initial guess), hence it does not remove these directions out each KSP iteration, just for each new linear solve. It can be used with the preconditioned conjugate gradient method.  There is a tiny community of people who claim this helps significantly on their problems, we'd love to hear your experience.<br>
<span class="HOEnZb"><font color="#888888"><br>
  Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
> On Dec 28, 2014, at 11:24 AM, Alp Kalpalp <<a href="mailto:alpkalpalp@gmail.com">alpkalpalp@gmail.com</a>> wrote:<br>
><br>
> Thanks for the answers,<br>
><br>
> Please forgive me, I forgot to say that my stiffness matrix is not changing during time steps. I could not remember directly but just after a google search..I just hit this<br>
><br>
> <a href="http://web.stanford.edu/group/frg/publications/recent/FETI-stoch.pdf" target="_blank">http://web.stanford.edu/group/frg/publications/recent/FETI-stoch.pdf</a><br>
><br>
> please look around eq37<br>
><br>
> My problem is not related to this random paper I found. But, I think I can find several others that shows the enhancing power of orthogonalization with successive directions when the system's behaviour is not changing rapidly. In my current sample case a gradually increasing force is applied to a linear system.<br>
><br>
> Since I use FETIDP, preconditioned conjugate projected gradient (PCPG) is crucial in order to select any generalized inverse for the system.<br>
><br>
> So, any suggestions on how to complete these tasks?<br>
><br>
> For example anyway of obtaining search direction from KSPCG?<br>
><br>
> or<br>
><br>
> how to implement a projection space?<br>
><br>
> Is it posible or too difficuly to code a variant of a KSPCG that meets my requirements?<br>
><br>
> On Sun, Dec 28, 2014 at 7:08 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
> On Sun, Dec 28, 2014 at 11:02 AM, Umut Tabak <<a href="mailto:u.tabak@tudelft.nl">u.tabak@tudelft.nl</a>> wrote:<br>
> On 12/28/2014 05:54 PM, Alp Kalpalp wrote:<br>
>> Hi,<br>
>><br>
>> Thank you Mark.<br>
>><br>
>> Let me clarify my questions;<br>
>><br>
>> 1-)How to implement or activate a Reorthogonalization procedure for KSPCG..<br>
>> As you know, search directions can be found more rapidly (with less numer of iterations) by using previous successive directions<br>
> Without answering the PETSc related questions, interesting discussion,<br>
><br>
> indeed, but at the cost of purging the previous directions(which means explicit orthogonalizations with respect to these vectors also), so I am not sure if you can gain something with this, cost wise...<br>
><br>
> This has been proposed many times, but it has never been shown to work. I have tried every variant I could<br>
> find and it did not work. You can try LGMRES, which is the closest one to working in my opinion. There is<br>
> definitely no theoretical relation between Krylov directions from subsequent solves unless the operator is<br>
> identical.<br>
><br>
>    Matt<br>
><br>
>> 2-) How to implement or activate a projection space over CG. A sample projection can be;<br>
>> P = I - G*((G'*G)\G').<br>
>> I need to insert project,scale,precondition,re-scae,re-project steps during each KSPCG iteration. How can I utilize this?<br>
>><br>
> Just a side note, I had previous experience on this that these kinds of practice increase the cost more...<br>
> BR,<br>
> Umut<br>
><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<br>
><br>
<br>
</div></div></blockquote></div><br></div>