<div dir="ltr">Hello Hong,<div><br></div><div>Thanks a lot. The result makes sense, but I think I will need higher accuracy (smaller tolerance). I will take a try myself at the meanwhile. </div><div><br></div><div>Best regards,</div>
<div>Jifeng Zhao</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jul 30, 2014 at 1:16 PM, <a href="mailto:hong@aspiritech.org">hong@aspiritech.org</a> <span dir="ltr"><<a href="mailto:hong@aspiritech.org" target="_blank">hong@aspiritech.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Jifeng,<br>
Using petsc-3.5 and slepc-3.5 (maint branch), I tested your matrices<br>
<br>
mpiexec -n 32 ./ex7 -f1 $data/zhao/petsc_stiff20.dat -f2<br>
$data/zhao/petsc_mass20.dat -mat_view ascii::ascii_info -mat_type aij<br>
-eps_gen_hermitian -eps_type krylovschur -eps_target -0.1 -st_ksp_type<br>
minres -st_pc_type bjacobi -st_sub_pc_type icc -st_ksp_rtol 1.e-2<br>
-eps_tol 1.e-2 -eps_nev 40 -st_type sinvert -log_summary |& tee<br>
log.mat20.ex7.np32.tol-2<br>
<br>
The output is attached.<br>
<br>
When using '-mat_type sbaij', I got same error as reported, which<br>
might be a bug in our libraries. I'll investigate it.<br>
Switching to '-mat_type aij' works (I'm not sure if the computed<br>
eigenvalues are meaningful though), but the linear solver<br>
minres/bjacobi/icc converges very slow, takes around 500 - 1000 linear<br>
iterations for each eps iteration. You may experiment with stronger<br>
pc_type, e.g., '-pc_type asm' etc.<br>
<br>
I'll let you know why sbaij crashes.<br>
<span class="HOEnZb"><font color="#888888"><br>
Hong<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
On Mon, Jul 28, 2014 at 4:52 PM, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>> wrote:<br>
><br>
> El 28/07/2014, a las 23:48, jifeng zhao escribió:<br>
><br>
>> Hi Barry,<br>
>><br>
>> Thanks for the reply. I see. I am still testing it and looking for the bug.<br>
>><br>
>> It is pretty weird that the error appears, since the code I am running is ex7 and I didn't change anything.<br>
>><br>
>> My matrices are assembled in binary files separately. As I explained earlier, I have a series of matrices to solve, each of them have the same patterns but with different values. The way I assembly the matrix has proven to be correct for all the other cases. Only for this one, when the matrix has an zero eigenvalue, ex7 failed to solve it.<br>

>><br>
>> In fact, I don't understand why this "Out of range" error could possibly appear? 1. my matrices have the correct parrellel lay up. 2. In ex7, I am not explicitly accessing to any elements at all! All I did is EPSSolve();<br>

>><br>
>> Any ideas?<br>
>><br>
>> Best regards,<br>
>> Jifeng Zhao<br>
>><br>
><br>
> Tomorrow we will release slepc-3.5. I would suggest trying with the new version. If the problem persist we can help to track down the problem.<br>
><br>
> Jose<br>
><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Jifeng Zhao<div>PhD candidate at Northwestern University, US</div><div>Theoretical and Applied Mechanics Program</div></div>
</div>