Bugfix in utility function barycentric_coords. Added two utility functions.
Made variables related to time step in VelocityFunc protected instead of private.