

def _do_problem(self, problem, integrator, method='adams'):
    f = (lambda t, z: problem.f(z, t))
    jac = None
    if hasattr(problem, 'jac'):
        jac = (lambda t, z: problem.jac(z, t))
    integrator_params = {
        
    }
    if ((problem.lband is not None) or (problem.uband is not None)):
        integrator_params['uband'] = problem.uband
        integrator_params['lband'] = problem.lband
    ig = self.ode_class(f, jac)
    ig.set_integrator(integrator, atol=(problem.atol / 10), rtol=(problem.rtol / 10), method=method, **integrator_params)
    ig.set_initial_value(problem.z0, t=0.0)
    z = ig.integrate(problem.stop_t)
    assert_array_equal(z, ig.y)
    assert_(ig.successful(), (problem, method))
    assert_(problem.verify(array([z]), problem.stop_t), (problem, method))

dataset/ETHPy150Open scipy/scipy/scipy/integrate/tests/test_integrate.py TestODEClass._do_problem/original