Temporal Gauss-Seidel solver