Added interface to pyOptSparse with automatic total derivatives; linear Gauss-Seidel tested