<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">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><div>use TS</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><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><div>No.  Certainly not by default.</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></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><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>