<div dir="ltr">Hello,<div><br></div><div>I am a new user to Petsc + Slepc. I am trying to extract natural frequency of a finite element model using Slepc. The way I do is </div><div><br></div><div>1. Use other software (Abaqus) to assembly the stiffness and mass matrix.</div>
<div>2. Use Slepc to solve a generalized eigenvalue problem. K x = lamda M x<div>with K, M being stiffness and mass matrix.</div><div><br></div><div>I wrote my petsc/slepc code based on the examples on slepc web. They all compiled and working correctly.</div>
<div><br></div><div>The question I am raising here is what solvers (solver combinations) should I use to be most efficient?   </div><div><br></div><div>Right now I am using "bcgsl" (BiCGSTAB) solver for KSP linear solvers, "JD" jacobian-davison for eigen solver, and "bjacobi" (block jacobian) for my preconditioner. It works, but I need it to be more efficient to solver big problem (millions of degrees of freedom). I am not an expert on knowing how these solvers are different at all!</div>
<div><br></div><div>Is there anybody who has extracted eigenvalues of a Finite element model using Slepc? How can I possibly improve the performance? </div><div><br></div><div>Thank you!</div><div><br></div><div>PS: my running command reads like:</div>
<div>







<p class="">./eigen_solver <span class="">-f1</span> petsc_stiff1.dat <span class="">-f2</span> petsc_mass1.dat -eps_nev <span class="">40</span> -eps_target <span class="">0</span>.<span class="">0</span> -eps_type jd -st_type precond -st_ksp_type bcgsl -st_pc_type bjacobi -st_ksp_rtol <span class="">0</span>.<span class="">001</span> -eps_tol 1e<span class="">-5</span> -eps_harmonic</p>
</div>-- <br><div dir="ltr">Jifeng Zhao<div>PhD candidate at Northwestern University, US</div><div>Theoretical and Applied Mechanics Program</div></div>
</div></div>