[petsc-users] SuperLU convergence problem

[Hong]

Danyang:
Using petsc/src/ksp/ksp/examples/tutorials/ex10.c,
I tested a_flow_check_168.bin

mpiexec -n 4 ./ex10 -f0
/Users/Hong/Downloads/matrix_and_rhs_bin/a_flow_check_168.bin -rhs
/Users/Hong/Downloads/matrix_and_rhs_bin/x_flow_check_168.bin -pc_type lu
-pc_factor_mat_solver_package superlu_dist
Number of iterations =   1
Residual norm 1.38302e-05

petsc (sequential) and mumps all give residual norm = 1.e-5.

For 169,
mpiexec -n 4 ./ex10 -f0
/Users/Hong/Downloads/matrix_and_rhs_bin/a_flow_check_169.bin -rhs
/Users/Hong/Downloads/matrix_and_rhs_bin/x_flow_check_169.bin -pc_type lu
-pc_factor_mat_solver_package superlu_dist
Number of iterations =   1
Residual norm 6.83854e-10

For 170, 171, 172, I get
mpiexec -n 4 ./ex10 -f0
/Users/Hong/Downloads/matrix_and_rhs_bin/a_flow_check_172.bin -rhs
/Users/Hong/Downloads/matrix_and_rhs_bin/x_flow_check_172.bin -pc_type lu
-pc_factor_mat_solver_package superlu_dist
Number of iterations =   1
  Residual norm < 1.e-12

Can you use this example test these matrices?
I use superlu_dist-v4.1.

Hong

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
>> <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/63eba366/attachment-0001.html>