<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, May 20, 2016 at 9:57 AM, Alejandro Cosimo <span dir="ltr"><<a href="mailto:alecosimo@gmail.com" target="_blank">alecosimo@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>Hi,<br><br></div>I would like to know if I'm using EPSSolve() correctly in the following situation. I must solve a series of eigenvalue problems where the size of the matrix and the EPS options remain the same, but the entries of the matrix change for each problem. In pseudo-code, something like this:<br><br>Mat A;<br></div><div>/* create and set matrix */<br></div><div>EPSCreate(comm,&eps);<br></div>/* here set EPS options */ <br></div>for (i=0;i<n;i++)  {<br></div><div>/* set entries of matrix A which vary for each index i of the for-loop */<br>EPSSetOperators(eps,A,NULL);<br>EPSSolve(eps);<br></div><div>}<br><br></div><div>My questions: <br>1* Is it correct to call at each iteration of the loop to EPSSetOperators() before proceeding to solve the eigenvalue problem? <br></div></div></blockquote><div><br></div><div>Yes</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>2* Would be wrong to call EPSSetOperators() before entering the loop and to call to EPSReset() at every iteration instead of calling to EPSSetOperators()?<br></div></div></blockquote><div><br></div><div>Yes, reset frees all memory and destroys the implementation object</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>3* A little different question, but related. Suppose a similar situation where you have a KSP solver instead of an EPS solver. In this context, do the answers to the previous questions apply similarly? That is, is it correct to call at every iteration to KSPSetOperators() before  calling to KSPSolve() or KSPSetUp()?<br></div></div></blockquote><div><br></div><div>Yes.</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"><div dir="ltr"><div></div><div>Thanks,<br></div><div>Alejandro<br></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>