This algorithm is implemented in the files vfluid/vfluid.(cc|hh).
Using the time step parameter
an adaptive control of the integration
time step is achieved.
should be chosen to permit a smooth
but steady deformation, and
is used to re-adjust
during the registration.