<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Hi Hong,<br>
<br>
The binary format of matrix, rhs and solution can be downloaded via
the link below.<br>
<br>
<a class="moz-txt-link-freetext" href="https://www.dropbox.com/s/cl3gfi0s0kjlktf/matrix_and_rhs_bin.tar.gz?dl=0">https://www.dropbox.com/s/cl3gfi0s0kjlktf/matrix_and_rhs_bin.tar.gz?dl=0</a><br>
<br>
Thanks,<br>
<br>
Danyang<br>
<br>
On 15-12-03 10:50 AM, Hong wrote:<br>
<blockquote
cite="mid:CAGCphBu5FOoAgwqhDFDzxEr=Up25WAu5+Gy8dYnDte+_GQC=Kg@mail.gmail.com"
type="cite">
<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 moz-do-not-send="true"
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>
</blockquote>
<br>
</body>
</html>