如何使用动态绑定覆盖纯虚拟 function

Question

我正在尝试解决 Pitt-Francis 文本https://link.springer.com/book/10.1007/978-3-319-73132-2中关于抽象 ODE class 的练习 7.3。

抽象 class 给出为

class AbstractODESolver
{
private:
    ...
    
public:
    ....
    virtual double RightHandSide(double y, double t) = 0;
    virtual double SolveEquation() = 0;
    
};

派生的 class 给出为

class ForwardEuler: public AbstractODESolver
{
public:
     double RightHandSide(double y, double t);
     double SolveEquation();
};

如何在 ForwardEuler 中转换 RightHandSide 以便在运行时定义它？

特别是主文件可以采用以下形式

double RHS(double y, double t);

int main(int argc, char* argv[])
{
    
    AbstractODESolver* p_ODEs = new ForwardEuler;
        
    return 0;
}

double RHS(double y, double t)
{
    return 1+t;
}

我想写一些类似的东西

p_ODEs->RightHandSide = &RHS;` 

指向 function RHS 。

Answer 1

给ForwardEuler class 一个数据成员来存储 function，并在构造函数中传递它。 RightHandSide function 只会调用传入的 function。

class ForwardEuler: public AbstractODESolver
{
public:
    ForwardEuler(std::function<double(double, double)> rhs): rhs(rhs) {}
    double RightHandSide(double y, double t) { return rhs(y, t); }
    double SolveEquation();

private:
    std::function<double(double, double)> rhs;
};

...

double RHS(double y, double t);
AbstractODESolver* p_ODEs = new ForwardEuler(RHS);

Answer 2

这是我不需要std::function的解决方案，我认为这是本章教授的技术所需的解决方案。

（根据@Paul Sanders）

class ForwardEuler: public AbstractODESolver
{
public:
     ForwardEuler(double (*pRhs)(double, double));
     double RightHandSide(double y, double t);
     double SolveEquation();

private:
     double (*mRhs)(double y, double t);     
};

ForwardEuler::ForwardEuler(double (*pRhs)(double, double))
{
    mRhs = pRhs;
}

double ForwardEuler::RightHandSide(double y,double t)
{
    return (*mRhs)(y,t);
}

主function

double RHS(double y, double t);

int main(int argc, char* argv[])
{
    
    AbstractODESolver* p_ODEs = new ForwardEuler(RHS);
    std::cout<<(*p_ODEs).RightHandSide(1.0,1.0)<<std::endl;

    return 0;
}

double RHS(double y, double t)
{
    return 1+t;
}