def eqom_jacobian_ekf(x,t,u): return cp_dynamics.eqom_det_jac(x,t)
def jac_gmm(x, t, u=None): return cp_dynamics.eqom_det_jac(x, t)