[petsc-users] SuperLU convergence problem

[Danyang Su]

Hi Hong,

The binary format of matrix, rhs and solution can be downloaded via the 
link below.

https://www.dropbox.com/s/cl3gfi0s0kjlktf/matrix_and_rhs_bin.tar.gz?dl=0

Thanks,

Danyang

On 15-12-03 10:50 AM, Hong wrote:
> Danyang:
>
>
>
>     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.
>
>     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.
>     _https://www.dropbox.com/s/rr8ueysgflmxs7h/results-check.tar.gz?dl=0_
>
>
> Can you send us matrix in petsc binary format?
>
> e.g., call MatView(M, PETSC_VIEWER_BINARY_(PETSC_COMM_WORLD))
> or '-ksp_view_mat binary'
>
> Hong
>
>
>
>     Below is a summary of the norm from the three solvers at timestep
>     29, newton iteration 1 to 5.
>
>     Timestep 29
>     Norm of residual seq 1.661321e-09, superlu 1.657103e+04, mumps
>     3.731225e-11
>     Norm of residual seq 1.753079e-09, superlu 6.675467e+02, mumps
>     1.509919e-13
>     Norm of residual seq 4.914971e-10, superlu 1.236362e-01, mumps
>     2.139303e-17
>     Norm of residual seq 3.532769e-10, superlu 1.304670e-04, mumps
>     5.387000e-20
>     Norm of residual seq 3.885629e-10, superlu 2.754876e-07, mumps
>     4.108675e-21
>
>     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.
>
>     Please let me know if I did not write down the problem clearly.
>
>     Thanks,
>
>     Danyang
>
>
>
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151203/301bbcef/attachment.html>