add ode_mat_exp for ODE solver based on matrix exponential