<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p>I've attached the B matrix (binary) in a zip file. To give more context to the problem, the previous installation of petsc/slecp doesn't give this error. So the matrix being singular shouldn't be the issue since my matrix assembly routines are independent
of the petsc/slepc installation I use. <span style="font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols; font-size:16px">Also, diagonalization happens for smaller
matrices.</span></p>
<p><br>
</p>
<p>Unfortunately, since EPSsolve never got completed, -eps_viewer didn't return anything. Is there any other way I can show the settings used?</p>
<p><br>
</p>
<p>Sidd</p>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Jose E. Roman <jroman@dsic.upv.es><br>
<b>Sent:</b> Wednesday, July 12, 2017 2:02:55 PM<br>
<b>To:</b> Hong Zhang<br>
<b>Cc:</b> petsc-users@mcs.anl.gov; Pandey, Siddhant<br>
<b>Subject:</b> Re: [petsc-users] SLEPC solvers not converging in latest version</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Let me correct my last message: we set this option in the development version (master), so yes, in 3.7 it may be good to try -ksp_error_if_not_converged<br>
Sorry for the noise.<br>
Jose<br>
<br>
<br>
> El 12 jul 2017, a las 19:58, Jose E. Roman <jroman@dsic.upv.es> escribió:<br>
> <br>
> Hong: in the latests versions of SLEPc we turn on this option by default (except in some solvers such as JD), so I don't think the problem is in the KSP solve.<br>
> Jose<br>
> <br>
> <br>
>> El 12 jul 2017, a las 19:53, Hong <hzhang@mcs.anl.gov> escribió:<br>
>> <br>
>> Add option '-ksp_error_if_not_converged'. I suspect the matrix has zero pivot.<br>
>> Hong<br>
>> <br>
>> On Wed, Jul 12, 2017 at 12:45 PM, Jose E. Roman <jroman@dsic.upv.es> wrote:<br>
>> <br>
>>> El 12 jul 2017, a las 18:56, Pandey, Siddhant <spandey2@wpi.edu> escribió:<br>
>>> <br>
>>> We recently updated our installation of PETSC and SLEPC from an older version, and eigenvalue problems converging in the previous installation is not converging anymore (diagonalization never comes to an end). On using -info, the following message keeps
repeating during diagonalization:<br>
>>> <br>
>>> [0] PCSetUp(): Leaving PC with identical preconditioner since operator is unchanged<br>
>>> [0] STMatSolve(): Linear solve iterations=1<br>
>>> [0] PCSetUp(): Leaving PC with identical preconditioner since operator is unchanged<br>
>>> [0] STMatSolve(): Linear solve iterations=1<br>
>>> [0] BV_SafeSqrt(): Zero norm, either the vector is zero or a semi-inner product is being used<br>
>>> [0] PCSetUp(): Leaving PC with identical preconditioner since operator is unchanged<br>
>>> [0] STMatSolve(): Linear solve iterations=1<br>
>>> [0] PCSetUp(): Leaving PC with identical preconditioner since operator is unchanged<br>
>>> [0] STMatSolve(): Linear solve iterations=1<br>
>>> <br>
>>> As can be seen, iteration number remains 1. The issue gets resolved with a smaller sized matrix. Currently I am using a 10648 x 10648 matrix, and PETSC is configured for complex values.<br>
>>> <br>
>>> What changes do I have to make to make my codes work in the new installation? This seems to be a version specific issue.<br>
>>> <br>
>>> <br>
>>> Thanks,<br>
>>> Sidd<br>
>>> WPI<br>
>> <br>
>> I don't understand why are you getting the BV message. Maybe your B-matrix is singular, but it is not sufficient to explain this. If possible, send me your matrices and I will try to reproduce the problem (I would also need to know all the settings, which
can be displayed with -eps_view or -eps_view_pre).<br>
>> <br>
>> Jose<br>
>> <br>
>> <br>
> <br>
<br>
</div>
</span></font>
</body>
</html>