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