[petsc-users] SuperLU convergence problem

[Hong]

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
> <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/19dbd815/attachment.html>