Rewritten from scratch: 3 of 5 initialization steps contained in base class; eliminated unnecessary duplication of data in each system; nonlinear GS and Jacobi, linear GS and Jacobi done; not yet tested