> 1) The solution norm you report is NaN, which indicates memory corruption > in your code I guess the reason may be a division by zero. The diagonals of the matrix have some zeros. >2) The number of allocated nonzeros is really suspicious. There are only 2 nonzeros per row? At least 5. Hui