JumpMap with input

Hi,

I think JumpMap with input, i,e,. the one documented in SystemDynamicsBase

$ x^+ = G \delta x + H \delta u \f$

is not yet implemented since the interface of systemDynamics supports only

  virtual VectorFunctionLinearApproximation jumpMapLinearApproximation(scalar_t t, const vector_t& x, const PreComputation& preComp);