<div dir="ltr"><div><div><div><div><div>Hi, <br><br></div>Thank you Mark. <br><br></div>Let me clarify my questions;<br><br></div>1-)How to implement or activate a Reorthogonalization procedure for KSPCG..<br></div><div>As you know, search directions can be found more rapidly (with less numer of iterations) by using previous successive directions<br></div><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></div>I need to insert project,scale,precondition,re-scae,re-project steps during each KSPCG iteration. How can I utilize this?<br><div><br></div><div>Thanks again and merry christmas to all<br></div><div><br><div><div><div><br></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Dec 28, 2014 at 4:28 PM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Sat, Dec 27, 2014 at 5:00 AM, Alp Kalpalp <span dir="ltr"><<a href="mailto:alpkalpalp@gmail.com" target="_blank">alpkalpalp@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div>Hi,<br><br></div>I implemented a newmark time stepping algorithm without using TS structure. I am following ex59 about PCBDDC.<br><br></div></div></div></div></div></div></div></blockquote><div><br></div></span><div>use TS</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div></div><div>ComputeMatrix<br></div>ComputeKSPBDDC <br><br>for // a time loop<br>{<br></div><div>Compute RHS<br></div><div>KSPSolve<br></div><div>GatherResults<br>MoveToNextTS<br></div><div>}<br></div><br><br></div>However, when I watch the iteration counts of KSPSolve they do not decrease signicantly..Decrease was around 5%. </div><div><br>So I have some problems;<br><br></div>1-) I guess currently, factorization is not taking place for each time step in my code. Ok this is expected. But I wonder whether Kspsolve stores the Krylov subspace vectors and reuse them for the next time step.<br><br></div></div></div></blockquote><div><br></div></span><div>No.  Certainly not by default.</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div></div>2-) PCBDDC uses KSPCG and AFAIK petsc doesnot have preconditioned conjugate projected gradient (PCPG). Is it possible to simulate PCPG iteration in some way?<br><br></div></div></blockquote><div><br></div></span><div>KSP has a PC object and all methods use it AFAIK.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div>Thanks,<br><br><div><br></div></div>
</blockquote></div><br></div></div>
</blockquote></div><br></div>