<div dir="ltr">Actually krylovschur it is giving excellent results when my B matrix is not singular. <div><br></div><div>This time my B matrix is large with very little non-zero entries. There are entire rows and columns with zeros. </div><div>So other than MUMPS, which I will try now, is there another way ?</div><div><br></div><div>Venkatesh</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 14, 2015 at 11:01 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">This has nothing to do with ARPACK.<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
El 14/05/2015, a las 19:29, venkatesh g escribió:<br>
<br>
> I am trying to use ARPACK. will it help?<br>
><br>
> It is not linking Arpack however, when I give ./configure --with-arpack-dir=..... command from manual.<br>
><br>
> Venkatesh<br>
><br>
> On Thu, May 14, 2015 at 10:54 PM, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>> wrote:<br>
><br>
> El 14/05/2015, a las 19:13, venkatesh g escribió:<br>
><br>
> > Hi,<br>
> > 1. I am not calling EPSSetFromOptions in source code. I am using the ex7.c given in the SLEPc webpage.<br>
> > 2. A-matrix is not singular. Only B is singular.<br>
> > 3. I will try with MUMPS.<br>
> > 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>
> > So is there a way ?<br>
> ><br>
><br>
> The ex7.c example calls EPSSetFromOptions(). If it is not called, then command-line options are ignored.<br>
><br>
> It is strange that you get zero pivots when factoring A-sigma*B for different values of sigma. Try with MUMPS.<br>
><br>
> Jose<br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div>