[petsc-users] about solving linear equation

[Likun Tan]

Hello,

I am solving a linear equation Ax=0 with Petsc.
Previously i was following Numerical Recipe and using LU decomposition,
for a non-singular matrix A, x should be zero. By applying LU in Recipe, i
got x=0.
When i solve it with ksp, x is some small value to the power of -9.

The condition number of A is about 52, is there any bug in my code or i
should find more stable linear solver?

Thanks,
Likun