Runs way too long to be a good test. Maybe also a good test case for solvers though. And looks fairly cool
Update: I got a big speedup by using the SIMPLE-preconditioned BiCGSTABSolver of !1989