<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>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>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>Thanks,<br><br><div><br></div></div>