我需要做些什么才能使odeint集成函数在另一个类中编译

What do I need to do to make the odeint integrate function compile within another class?

本文关键字：函数 另一个 编译 集成 odeint 做些什么      更新时间：2023-10-16

几天来，我一直在尝试在类中编译ODEINT的integrate_adapive（）函数。如果我把函数放在类之外，并有通用的函数名，那么程序就会根据需要编译和工作。（out只是一个文件的输出，并且是以前定义的）

void waterbodies::writeInfo(){
    bulirsch_stoer_dense_out< state_type > stepper( 1E-9 , 1E-9 , 1.0 , 0.0 );
    state_type x1 = { 0.0024 , 0.00001 , 0.0008,0.6,100,5,5,4,4,4,4,20000,0,10,10,100 }; // initial conditions
    out.open( Name.c_str() );
    out.precision(16);
    double t = 0.0;
    double max = 100;
    double dt = 0.01;
    boost::numeric::odeint::integrate_adaptive(
            rhs,
            x1 , t , max , dt, write_rhs);
    out.close();
}

我有一个叫做rhs的函数，还有一个叫做write_rhs 的函数

void waterbodies::write_rhs( const state_type &x , const double t )
{
    out << t << 't' << x[0] << 't' << x[1] << 't' << x[2] << 't' << x[3] << 't' << x[4] << 't' << x[5] << 't' << x[6] << 't' << x[7] << 't' << x[8] << 't' << x[9] << 't' << x[10] << 't' << x[11] << 't' << x[12] << 't' << x[13] << 't' << x[14] << 't' << x[15] << std::endl;
}

void waterbodies::rhs( const state_type &x , state_type &dxdt , const double t  )
{}

state_type在标头中定义如下

typedef boost::array< double , 16 > state_type;

错误输出为：

    ..srcwaterbodies.cpp: In member function 'void waterbodies::writeInfo()':
..srcwaterbodies.cpp:84:32: error: no matching function for call to 'integrate_adaptive(<unresolved overloaded function type>, waterbodies::state_type&, double&, double&, double&, <unresolved overloaded function type>)'
    x1 , t , max , dt, write_rhs);
                                ^
..srcwaterbodies.cpp:84:32: note: candidates are:
In file included from C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate.hpp:26:0,
                 from C:softwareboost_1_59_0/boost/numeric/odeint.hpp:67,
                 from ..srcwaterbodies.cpp:4:
C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate_adaptive.hpp:36:8: note: template<class Stepper, class System, class State, class Time, class Observer> size_t boost::numeric::odeint::integrate_adaptive(Stepper, System, State&, Time, Time, Time, Observer)
 size_t integrate_adaptive(
        ^
C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate_adaptive.hpp:36:8: note:   template argument deduction/substitution failed:
..srcwaterbodies.cpp:84:32: note:   could not resolve address from overloaded function '((waterbodies*)this)->waterbodies::write_rhs'
    x1 , t , max , dt, write_rhs);
                                ^
In file included from C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate.hpp:26:0,
                 from C:softwareboost_1_59_0/boost/numeric/odeint.hpp:67,
                 from ..srcwaterbodies.cpp:4:
C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate_adaptive.hpp:60:8: note: template<class Stepper, class System, class State, class Time, class Observer> size_t boost::numeric::odeint::integrate_adaptive(Stepper, System, const State&, Time, Time, Time, Observer)
 size_t integrate_adaptive(
        ^
C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate_adaptive.hpp:60:8: note:   template argument deduction/substitution failed:
..srcwaterbodies.cpp:84:32: note:   could not resolve address from overloaded function '((waterbodies*)this)->waterbodies::write_rhs'
    x1 , t , max , dt, write_rhs);
                                ^
In file included from C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate.hpp:26:0,
                 from C:softwareboost_1_59_0/boost/numeric/odeint.hpp:67,
                 from ..srcwaterbodies.cpp:4:
C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate_adaptive.hpp:79:8: note: template<class Stepper, class System, class State, class Time> size_t boost::numeric::odeint::integrate_adaptive(Stepper, System, State&, Time, Time, Time)
 size_t integrate_adaptive(
        ^
C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate_adaptive.hpp:79:8: note:   template argument deduction/substitution failed:
..srcwaterbodies.cpp:84:32: note:   could not resolve address from overloaded function '((waterbodies*)this)->waterbodies::write_rhs'
    x1 , t , max , dt, write_rhs);
                                ^
In file included from C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate.hpp:26:0,
                 from C:softwareboost_1_59_0/boost/numeric/odeint.hpp:67,
                 from ..srcwaterbodies.cpp:4:
C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate_adaptive.hpp:91:8: note: template<class Stepper, class System, class State, class Time> size_t boost::numeric::odeint::integrate_adaptive(Stepper, System, const State&, Time, Time, Time)
 size_t integrate_adaptive(
        ^
C:softwareboost_1_59_0/boost/numeric/odeint/integrate/integrate_adaptive.hpp:91:8: note:   template argument deduction/substitution failed:
..srcwaterbodies.cpp:84:32: note:   could not resolve address from overloaded function '((waterbodies*)this)->waterbodies::write_rhs'
    x1 , t , max , dt, write_rhs);

如果你需要我解释什么，请告诉我，我会尽量让我说的更清楚。

感谢

Simon

这个问题已经回答了好几次：

在类中使用boost:：numeric:：odeint

您需要将rhs和write_rhs绑定到当前类：

// not tested
// the code also requires a C++14 compiler
auto rhs_wrapper = [this]( auto const& x , auto& dxdt , auto t ) {
    this->rhs( x , dxdt , t ); };
auto obs = [this]( auto const& x , auto t ) {
    this->write_rhs( x , t );
integrate_adaptive( stepper , rhs_wrapper , x , max , dt , obs );
相关文章：