Made paired down constructor for DiscreteLagrangian to do the tracing
Variational Midpoint nearly works. But calling member function Integrator.model_->rhs not allowed in static function.
Test-Passing State Restored