<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, May 12, 2014 at 12:27 PM,  <span dir="ltr"><<a href="mailto:likunt@caltech.edu" target="_blank">likunt@caltech.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Petsc developers,<br>
<br>
I am solving a linear system Ax=b, while A is constant and b is changing<br>
in each time step. Here is the code I wrote:<br>
<br>
/****************************************************************<br>
...compute matrix A...<br>
KSPCreate(PETSC_COMM_WORLD, &ksp);  CHKERRQ(ierr);<br>
KSPSetOperators(ksp, A, A, SAME_PRECONDITIONER);  CHKERRQ(ierr);<br>
KSPSetTolerances(ksp, 1.e-5, 1.E-50, PETSC_DEFAULT, PETSC_DEFAULT);<br>
KSPSetFromOptions(ksp);<br>
for(int step=0; step<STEP; step++)<br>
{<br>
    ... compute vector b ...<br>
    KSPSolve(ksp, b, x);<br>
}<br>
*****************************************************************/<br>
<br>
I tested a system with size 1725*1725, on 4 processors, it takes 0.06s.<br>
Would you please let me know if there is a way to improve its efficiency?<br></blockquote><div><br></div><div>It would be amazing if we could do that given the description. First, we do not know</div><div>exactly what solver is being used (-ksp_view), but lets assume its GMRES/ILU(0)</div>
<div>which is the default. Second, we have no idea what the convergence was like</div><div>(-ksp_monitor_true_residual -ksp_converged_reason), so we do not know what</div><div>the bottleneck is, and have no performance monitoring (-log_summary). Lastly,</div>
<div>even if we had that we have no idea what the operator is so that we could make</div><div>intelligent suggestions for other preconditioners.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Thanks.<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <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
</div></div>