В Tutorial мы узнали, как можно использовать функции генерации make_controlled и make_dense_output для создания контролируемой и плотной выходной степпер из простой степпер или степпер ошибки. Синтаксис этих двух функций очень прост:

auto stepper1 = make_controlled( 1.0e-6 , 1.0e-6 , stepper_type() );
// or with max step size limit:
// auto stepper1 = make_controlled( 1.0e-6 , 1.0e-6 , 0.01, stepper_type() );
auto stepper2 = make_dense_output( 1.0e-6 , 1.0e-6 , stepper_type() );

Первые два параметра являются абсолютными и относительными допусками ошибок, а третий параметр - степпер. Кроме того, существует вторая версия, в которой дополнительно предоставляется максимальный размер шага, который гарантирует, что размер шага не будет увеличен выше этого значения. В C++03 можно вывести тип из механизма result_of:

boost::numeric::odeint::result_of::make_controlled< stepper_type >::type stepper3 = make_controlled( 1.0e-6 , 1.0e-6 , stepper_type() );
(void)stepper3;
boost::numeric::odeint::result_of::make_dense_output< stepper_type >::type stepper4 = make_dense_output( 1.0e-6 , 1.0e-6 , stepper_type() );
(void)stepper4;

Для использования собственных степперов с make_controlled или make_dense_output необходимо специализироваться на двух шаблонах классов. Предположим, что ваши степперы называются custom_stepper, custom_controller и custom_dense_output. Затем, первый класс, который вам нужно специализировать, это boost::numeric::get_controller, метафункция, возвращающая тип контроллера:

namespace boost { namespace numeric { namespace odeint {
template<>
struct get_controller< custom_stepper >
{
    typedef custom_controller type;
};
} } }

Второй — фабричный класс boost::numeric::odeint::controller_factory, который конструирует контроллер из допусков и степпера. В нашей фиктивной реализации этот класс является

namespace boost { namespace numeric { namespace odeint {
template<>
struct controller_factory< custom_stepper , custom_controller >
{
    custom_controller operator()( double abs_tol , double rel_tol , const custom_stepper & ) const
    {
        return custom_controller();
    }
    custom_controller operator()( double abs_tol , double rel_tol , double max_dt ,
                                  const custom_stepper & ) const
    {
        // version with maximal allowed step size max_dt
        return custom_controller();
    }
};
} } }

Это все для использования механизма make_controlled. Вы можете использовать контроллер через

auto stepper5 = make_controlled( 1.0e-6 , 1.0e-6 , custom_stepper() );

Для плотного_output_stepper все работает одинаково. Здесь вы должны специализироваться boost::odeint::get_dense_output и boost::odeint::dense_output_factory. Эти два класса имеют тот же синтаксис, что и их родственники get_controller и controller_factory.

С этими механизмами можно использовать все контроллеры и шаговые устройства с плотным выходом в одейте. В таблице ниже вы найдете, какие степперы построены из make_control или make_dense_output при нанесении на степпер из odeint: