Hello Hong, 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.
Best regards, Jifeng Zhao On Wed, Jul 30, 2014 at 1:16 PM, [email protected] <[email protected]> wrote: > Jifeng, > Using petsc-3.5 and slepc-3.5 (maint branch), I tested your matrices > > mpiexec -n 32 ./ex7 -f1 $data/zhao/petsc_stiff20.dat -f2 > $data/zhao/petsc_mass20.dat -mat_view ascii::ascii_info -mat_type aij > -eps_gen_hermitian -eps_type krylovschur -eps_target -0.1 -st_ksp_type > minres -st_pc_type bjacobi -st_sub_pc_type icc -st_ksp_rtol 1.e-2 > -eps_tol 1.e-2 -eps_nev 40 -st_type sinvert -log_summary |& tee > log.mat20.ex7.np32.tol-2 > > The output is attached. > > When using '-mat_type sbaij', I got same error as reported, which > might be a bug in our libraries. I'll investigate it. > Switching to '-mat_type aij' works (I'm not sure if the computed > eigenvalues are meaningful though), but the linear solver > minres/bjacobi/icc converges very slow, takes around 500 - 1000 linear > iterations for each eps iteration. You may experiment with stronger > pc_type, e.g., '-pc_type asm' etc. > > I'll let you know why sbaij crashes. > > Hong > > On Mon, Jul 28, 2014 at 4:52 PM, Jose E. Roman <[email protected]> wrote: > > > > El 28/07/2014, a las 23:48, jifeng zhao escribió: > > > >> Hi Barry, > >> > >> Thanks for the reply. I see. I am still testing it and looking for the > bug. > >> > >> It is pretty weird that the error appears, since the code I am running > is ex7 and I didn't change anything. > >> > >> 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. > >> > >> 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(); > >> > >> Any ideas? > >> > >> Best regards, > >> Jifeng Zhao > >> > > > > 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. > > > > Jose > > > -- Jifeng Zhao PhD candidate at Northwestern University, US Theoretical and Applied Mechanics Program
