<div dir="ltr"><div><div><div><div>Hello,<br></div>          I see an interesting behavior when i call EPSSolve in a loop. Can you help me figure out whats going on ?<br><br></div>I have a setup like this to solve a A x = lambda B x   (Generalized Hermitian problem)<br>
<br></div>type eigenSolver_t<br></div><div>  EPS :: eps<br></div><div>  <br></div>end type eigenSolver_t<br clear="all"><div><div><div><div><div><br>the type has bound procedure which calls  EPSSetOperators and EPSSolve when eigenSolver%solve() is called<br>
<br></div><div>Now i run a for loop <br></div><div><br></div><div>do ic=1,111<br></div><div>  call eigenSolver(ic)%solve()  <br></div><div>end do<br></div><div><br></div><div>I print the time for each solve.  The operators A,B =A1,B1 for ic=1 to 50<br>
</div><div>A2,B2 for ic=51 to 80  and A3,B3 for ic=81 to 111<br><br></div><div>Now i see that time per solve per ic is almost constant when i use -eps_type lapack.<br></div><div>But for, defualt solver, time per solve seem to increase with increasing ic. Have look at the attached timing information.<br>
<br></div><div>Also, Time per solve using lapack is lower than any of the iterative solvers i have tried. Problem size is about 100 x 100, operators are tri-diagonal. <br><br><br></div><div>Thanks<br>Reddy<br></div><div>-- <br>
-----------------------------------------------------<br>Dharmendar Reddy Palle<br>Graduate Student<br>Microelectronics Research center,<br>University of Texas at Austin,<br>10100 Burnet Road, Bldg. 160<br>MER 2.608F, TX 78758-4445<br>
e-mail: <a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a><br>Phone: +1-512-350-9082<br>United States of America.<br>Homepage: <a href="https://webspace.utexas.edu/~dpr342" target="_blank">https://webspace.utexas.edu/~dpr342</a><br>

</div></div></div></div></div></div>