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 > > > >