<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Danyang:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
  

    
  
  <div text="#000000" bgcolor="#FFFFFF"><br>
    <br>
    To my surprising, solutions from SuperLU at timestep 29 seems not
    correct for the first 4 Newton iterations, but the solutions from
    iteration solver and MUMPS are correct. <br>
    <br>
    Please find all the matrices, rhs and solutions at timestep 29 via
    the link below. The data is a bit large so that I just share it
    through Dropbox. A piece of matlab code to read these data and then
    computer the norm has also been attached. <br>
    <u><a href="https://www.dropbox.com/s/rr8ueysgflmxs7h/results-check.tar.gz?dl=0" target="_blank">https://www.dropbox.com/s/rr8ueysgflmxs7h/results-check.tar.gz?dl=0</a></u></div></blockquote><div><br></div><div>Can you send us matrix in petsc binary format?</div><div><br></div><div>e.g., call MatView(M, PETSC_VIEWER_BINARY_(PETSC_COMM_WORLD))</div><div>or '-ksp_view_mat binary'</div><div><br></div><div>Hong</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div text="#000000" bgcolor="#FFFFFF"><br>
    <br>
    Below is a summary of the norm from the three solvers at timestep
    29, newton iteration 1 to 5.<br>
    <br>
    Timestep 29<br>
    Norm of residual seq 1.661321e-09, superlu 1.657103e+04, mumps
    3.731225e-11 <br>
    Norm of residual seq 1.753079e-09, superlu 6.675467e+02, mumps
    1.509919e-13 <br>
    Norm of residual seq 4.914971e-10, superlu 1.236362e-01, mumps
    2.139303e-17 <br>
    Norm of residual seq 3.532769e-10, superlu 1.304670e-04, mumps
    5.387000e-20 <br>
    Norm of residual seq 3.885629e-10, superlu 2.754876e-07, mumps
    4.108675e-21 <br>
    <br>
    Would anybody please check if SuperLU can solve these matrices?
    Another possibility is that something is wrong in my own code. But
    so far, I cannot find any problem in my code since the same code
    works fine if I using iterative solver or direct solver MUMPS. But
    for other cases I have tested,  all these solvers work fine.<br>
    <br>
    Please let me know if I did not write down the problem clearly.<br>
    <br>
    Thanks,<br>
    <br>
    Danyang<br>
    <br>
    <br>
    <br>
  </div>

</blockquote></div><br></div></div>