Added convergence criterion when calling Gauss_Seidel subroutine