Added visualization for max density cells.
Tried debugging project() to no avail.
Added a Gauss-Seidel solver object (not yet used).