<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Oct 28, 2017 at 7:11 AM, Franck Houssen <span dir="ltr"><<a href="mailto:franck.houssen@inria.fr" target="_blank">franck.houssen@inria.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I just added that before EPSSetOperators:<br>
PetscViewer viewerA; PetscViewerBinaryOpen(PETSC_<wbr>COMM_WORLD,"Atau.out",FILE_<wbr>MODE_WRITE,&viewerA); MatView(A,viewerA);<br>
PetscViewer viewerB; PetscViewerBinaryOpen(PETSC_<wbr>COMM_WORLD,"Btau.out",FILE_<wbr>MODE_WRITE,&viewerB); MatView(B,viewerB);<br>
<br>
At first, I avoided binary as I didn't know if the format handles big/little endianness... So I prefered ASCII.<br></blockquote><div><br></div><div>Everything is big endian. We convert it from that to what we need on the machine.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Binary data are attached.<br>
<br>
Franck<br>
<br>
PS : running debian with little endian.<br>
>> python -c "import sys;print(0 if sys.byteorder=='big' else 1)"<br>
1<br>
<br>
<br>
----- Mail original -----<br>
> De: "Jose E. Roman" <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>><br>
> À: "Franck Houssen" <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>><br>
> Cc: "For users of the development version of PETSc" <<a href="mailto:petsc-dev@mcs.anl.gov">petsc-dev@mcs.anl.gov</a>><br>
> Envoyé: Vendredi 27 Octobre 2017 18:52:56<br>
> Objet: Re: [petsc-dev] SLEPc failure<br>
><br>
> I cannot load the files you sent. Please send the matrices in binary format.<br>
> The easiest way is to run your program with -eps_view_mat0 binary:Atau.bin<br>
> -eps_view_mat1 binary:Btau.bin<br>
><br>
> However, the files are written at the end of EPSSolve() so if the solve fails<br>
> then it will not create the files. You can try running with -eps_max_it 1<br>
> or add code in your main program to write the matrices.<br>
><br>
> Jose<br>
><br>
><br>
> > El 27 oct 2017, a las 12:28, Franck Houssen <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>><br>
> > escribió:<br>
> ><br>
> > Maybe could be convenient for the users to have an option (or an EPSSetXXX)<br>
> > to relax that check ?<br>
> > Data are attached.<br>
> ><br>
> > Franck<br>
> ><br>
> > ----- Mail original -----<br>
> >> De: "Jose E. Roman" <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>><br>
> >> À: "Franck Houssen" <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>><br>
> >> Cc: "For users of the development version of PETSc"<br>
> >> <<a href="mailto:petsc-dev@mcs.anl.gov">petsc-dev@mcs.anl.gov</a>><br>
> >> Envoyé: Vendredi 27 Octobre 2017 10:15:44<br>
> >> Objet: Re: [petsc-dev] SLEPc failure<br>
> >><br>
> >> There is no new option. What I mean is that from 3.7 to 3.8 we changed the<br>
> >> line that produces this error. But it seems that it is still failing in<br>
> >> your<br>
> >> problem. Maybe your B matrix is indefinite or not exactly symmetric. Can<br>
> >> you<br>
> >> send me the matrices?<br>
> >> Jose<br>
> >><br>
> >>> El 27 oct 2017, a las 9:57, Franck Houssen <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>><br>
> >>> escribió:<br>
> >>><br>
> >>> I use the development version (bitbucket clone). How to relax the check ?<br>
> >>> At command line option ?<br>
> >>><br>
> >>> Franck<br>
> >>><br>
> >>> ----- Mail original -----<br>
> >>>> De: "Jose E. Roman" <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>><br>
> >>>> À: "Franck Houssen" <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>><br>
> >>>> Cc: "For users of the development version of PETSc"<br>
> >>>> <<a href="mailto:petsc-dev@mcs.anl.gov">petsc-dev@mcs.anl.gov</a>><br>
> >>>> Envoyé: Jeudi 26 Octobre 2017 18:49:22<br>
> >>>> Objet: Re: [petsc-dev] SLEPc failure<br>
> >>>><br>
> >>>><br>
> >>>>> El 26 oct 2017, a las 18:36, Franck Houssen <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>><br>
> >>>>> escribió:<br>
> >>>>><br>
> >>>>> Here is a stack I end up with when trying to solve an eigen problem<br>
> >>>>> (real,<br>
> >>>>> sym, generalized) with SLEPc. My understanding is that, during the Gram<br>
> >>>>> Schmidt orthogonalisation, the projection of one basis vector turns out<br>
> >>>>> to<br>
> >>>>> be null.<br>
> >>>>> First, is this correct ? Second, in such cases, are there some<br>
> >>>>> recommended<br>
> >>>>> "recipe" to test/try (options) to get a clue on the problem ? (I would<br>
> >>>>> unfortunately perfectly understand the answer could be no !... As this<br>
> >>>>> totally depends on A/B).<br>
> >>>>><br>
> >>>>> With arpack, the eigen problem is solved (so the matrix A and B I use<br>
> >>>>> seems<br>
> >>>>> to be relevant). But, when I change from arpack to<br>
> >>>>> krylovschur/ciss/arnoldi, I get the stack below.<br>
> >>>>><br>
> >>>>> Franck<br>
> >>>>><br>
> >>>>> [0]PETSC ERROR: #1 BV_SafeSqrt()<br>
> >>>>> [0]PETSC ERROR: #2 BVNorm_Private()<br>
> >>>>> [0]PETSC ERROR: #3 BVNormColumn()<br>
> >>>>> [0]PETSC ERROR: #4 BV_NormVecOrColumn()<br>
> >>>>> [0]PETSC ERROR: #5 BVOrthogonalizeCGS1()<br>
> >>>>> [0]PETSC ERROR: #6 BVOrthogonalizeGS()<br>
> >>>>> [0]PETSC ERROR: #7 BVOrthonormalizeColumn()<br>
> >>>>> [0]PETSC ERROR: #8 EPSFullLanczos()<br>
> >>>>> [0]PETSC ERROR: #9 EPSSolve_KrylovSchur_Symm()<br>
> >>>>> [0]PETSC ERROR: #10 EPSSolve()<br>
> >>>><br>
> >>>> Is this with SLEPc 3.8? In SLEPc 3.8 we relaxed this check so I would<br>
> >>>> suggest<br>
> >>>> trying with it.<br>
> >>>> Jose<br>
> >>>><br>
> >>>><br>
> >><br>
> >><br>
> > <ABtau.tar.gz><br>
><br>
><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>