simo-removero as matrix + bcs updated to t^{n+1/2}
gearing up for poiseuille