Fixes memory error on 3d case with US, other
Re-use some factorizations
Fail on nonconvergence of inner iterations
Set fluid storage coefficient to reasonable value for spinal_cord
Add notes about Gauss-Seidel, number equations