<t>I'm finding that the jacobi, Gauss-Sediel and SOR methods don't perform well when the matrix is not srtrictly diagonal. These were the solvers I first tried when implimenting the Galekin method using cubic polynomials. That's when I moved to LU.Maybe there are "fixes" to improve the performance -...