Remove unused code.
Add convergence check, local whitespace cleanup.
Cite CPR article by Gries, Stuben et al.
Put creation of elliptic system in separate function.
Remove unused data members and parameters.