added storage calc in Gauss-Seidel