<div dir="ltr"><div style="font-size:12.8000001907349px">Hi, <br></div><div style="font-size:12.8000001907349px">1. I am not calling EPSSetFromOptions in source code. I am using the ex7.c given in the SLEPc webpage. <br></div><div style="font-size:12.8000001907349px">2. A-matrix is not singular. Only B is singular.<br></div><div style="font-size:12.8000001907349px">3. I will try with MUMPS.<br></div><div style="font-size:12.8000001907349px">4. Matrix sizes are same as in Matlab. If B is not singular it is giving same result as Matlab, but if B is singular it gives the zero pivot issue which Matlab doesnt give. <br></div><div style="font-size:12.8000001907349px">So is there a way ?</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 14, 2015 at 2:03 PM, Jose E. Roman <span dir="ltr"><<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Are you calling EPSSetFromOptions in your source code?<br>
Is your A-matrix also singular with a nullspace in common with B?<br>
You can try with MUMPS as explained in section 3.4.1 of the users manual.<br>
Also check that your matrices are the same as in Matlab.<br>
<span class="HOEnZb"><font color="#888888">Jose<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
El 14/05/2015, a las 08:52, venkatesh g escribió:<br>
<br>
> Hi Jose, all,<br>
><br>
> I tried that too, it is telling the same error of "Detected Zero pivot in LU factorization".<br>
> My matrix B in AX=lambda BX has large no. of zeros and is singular.<br>
><br>
> Is there another way ?<br>
><br>
> Venkatesh<br>
><br>
> On Wed, May 13, 2015 at 9:02 PM, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>> wrote:<br>
><br>
> El 13/05/2015, a las 17:28, venkatesh g escribió:<br>
><br>
> > Dear All,<br>
> ><br>
> > I am having a block matrix Eigenvalue problem AX = lambda BX<br>
> ><br>
> > My B is singular, it is a block matrix and only one block has non-zero values.<br>
> ><br>
> > Matlab Eigs is giving the result.<br>
> ><br>
> > I am using option: "./ex7 -f1 A1 -f2 B1 -st_type sinvert"<br>
> ><br>
> > Result is: Detected zero pivot in LU factorization.<br>
> ><br>
> > The ex7 is an example from SLEPc.<br>
> ><br>
> > Pls advice what to do ?<br>
> ><br>
> > cheers,<br>
> > Venkatesh<br>
> ><br>
> ><br>
><br>
> If you know where your wanted eigenvalues are, set the target value, for instance:<br>
> ./ex7 -f1 A1 -f2 B1 -st_type sinvert -eps_target 0.5<br>
><br>
> Jose<br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div>