[petsc-users] eps_target malloc issue

[jifeng zhao]

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, hong at aspiritech.org <hong at aspiritech.org>
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 <jroman at dsic.upv.es> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140730/740d0001/attachment-0001.html>