On Tue, Aug 21, 2012 at 11:25 AM, Alexander Grayver <span dir="ltr"><<a href="mailto:agrayver@gfz-potsdam.de" target="_blank">agrayver@gfz-potsdam.de</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    <div>On <a href="tel:21.08.2012%2018" value="+12108201218" target="_blank">21.08.2012 18</a>:14, Hong Zhang wrote:<br>
    </div>
    <blockquote type="cite">Alexander :
      <div><br>
        <div>Testing your matrix A.dat
          using petsc-3.3/src/ksp/ksp/examples/tutorials/ex10.c,</div>
        <div>I find that A is numerically singular:</div>
      </div>
    </blockquote>
    <br>
    Hong,<br>
    <br>
    Sorry, but I forgot to mention that this is double complex symmetric
    matrix. Did you load and treat it like that?<br>
    <br>
    <blockquote type="cite">
      <div>
        <div>I'm convinced that matrix A is corrupted or obtained from
          an incorrect numerical model.</div>
        <div>Also note, your matrix A has size of approx. 60k, which is
          too large for sequential direct solver.</div>
      </div>
    </blockquote>
    <br>
    MUMPS takes only several minutes and 6 GB of memory to factorize it.
    <br>
    This factorization gives residual on the order of 10e-12 and
    solution is indeed correct.<br>
    <br>
    Nevertheless, you're right, there is numerical null-space in this
    matrix since it comes <br>
    from the discretization of equation that contains curl curl
    operator, but practically this <br>
    case is not really the worst one. <br></div></blockquote><div><br></div><div>This makes no sense whatsoever. How can you LU factor a matrix that has a null space?</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">
<div bgcolor="#FFFFFF" text="#000000">
    <blockquote type="cite">
      <div>
        <div>
          <br>
        </div>
        <div>The error on MatMatSolve() when calling </div>
        <div><span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">' -pc_type
            cholesky -pc_factor_mat_solver_package mumps'</span></div>
        <div><span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">is
            a bug in our petsc-mumps interface. I've pushed a fix to
            petsc-3.3</span></div>
        <div><a href="http://petsc.cs.iit.edu/petsc/releases/petsc-3.3/rev/8badc49a596e" target="_blank">http://petsc.cs.iit.edu/petsc/releases/petsc-3.3/rev/8badc49a596e</a></div>
        <div><br>
        </div>
        <div>Thanks for your report.</div>
      </div>
    </blockquote>
    <br>
    Thank you for the patch.<span class="HOEnZb"><font color="#888888"><br>
    <br>
    <pre cols="72">-- 
Regards,
Alexander</pre>
  </font></span></div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>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<br>