

def _get_solver(self, f, jac):
    solver = ode(f, jac)
    if self.solver_uses_jac:
        solver.set_integrator(self.solver_name, atol=1e-09, rtol=1e-07, with_jacobian=self.solver_uses_jac)
    else:
        solver.set_integrator(self.solver_name, atol=1e-09, rtol=1e-07)
    return solver

dataset/ETHPy150Open scipy/scipy/scipy/integrate/tests/test_integrate.py ODECheckParameterUse._get_solver/original