Close to first release. Testing on 1D and 2D complete. Needs splitting, implicity, parallelization of element boundaries, and boundary conditions.