solving a particular linear equation system
cmay at phys.ethz.ch
Fri May 15 14:29:09 CDT 2009
I have trouble solving a linear equation system with PETSc 3.0.0.
I tried the following combinations of solvers and preconditioners and list
cg+sor: Linear solve did not converge due to DIVERGED_NAN
cg+bjacobi: Linear solve did not converge due to DIVERGED_INDEFINITE_PC
cg+ilu: Linear solve did not converge due to DIVERGED_INDEFINITE_MAT
gmres+sor: Linear solve did not converge due to DIVERGED_NAN
gmres+bjacobi:Linear solve did not converge due to DIVERGED_ITS
gmres+ilu: Linear solve did not converge due to DIVERGED_ITS
It might be that the matrix is not positive definite, which explains
some of the problems. In the last two cases it kind of works, but convergence
is horribly slow and it aborts after 10000 iterations.
If somebody could please have a look at the matrix and the right hand side
vector, I've put them here:
The archive contains a file named A.txt, consisting of three columns
i j Aij and a file b.txt with two columns i and bi.
Any hint on how to solve this linear equation system (in reasonable time)
Thanks a lot in advance!
More information about the petsc-users