
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=koi8-r">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">

<p>Dear PETSc community,</p>

<p><br>

</p>

<p></p>

<span>I wonder if you could give me a hint on how to compute the condition number of a preconditioned matrix in a proper way.</span>

<p>I have a <i style="">MatSchurComplement</i> matrix S and a preconditioner P of the type

<i>PCSHELL</i> (P is a diffusion matrix, which itself is inverted by <i>KSPCG</i>).</p>

<p>I tried to compute the condition number of P^{-1}S "for free" during the outer PCG procedure using <span><i>KSPComputeExtremeSingularValues()</i> routine.</span></p>

<p>Unfortunately, \sigma_min does not converge even if the solution is computed with very high precision.</p>

<p><span>I also looked at SLEPc interface, but did not realised how PC should be included.</span><br>

</p>

<p>Thanks!</p>

<p><br>

</p>

<p>Sincerely,</p>

<p>Vladislav Pimanov</p>

<p><br>

</p>

</div>

</body>

</html>