[T/RM] Add liakopoulos test using assemble().

Using central differences for Jacobian computation.
Still the HydraulicFlow and Nodal forces are not correct
and don't include the 'M \dot{x} + K x' part.