所以我有一个带有以下头文件的基类:
#ifndef ODESolver_H
#define ODESolver_H
#endif
#include "DoublePendulum.h"
class ODESolver
{
public:
ODESolver(DoublePendulum &);
virtual void Predict(const double &)=0;
DoublePendulum DP;
};
执行:
#include "ODESolver.h"
#include <iomanip>
ODESolver::ODESolver(DoublePendulum &DPRef) :DP(DPRef)
{}
还有一个继承自它的类:
#ifndef ODEEuler_H
#define ODEEuler_H
#endif
#include "ODESolver.h"
class ODEEuler : public ODESolver
{
public:
ODEEuler(DoublePendulum &);
void Predict(const double &);
};
执行:
#include "ODEEuler.h"
ODEEuler::ODEEuler (DoublePendulum &DPRef) :ODESolver(DPRef)
{}
void ODEEuler::Predict(const double &dt=0.005)
{
DP=DP+DP.Derivate()*dt;
cout << DP.getUp().getTheta() << endl;
}
问题是在上面的实现中,DP=DP+DP.Derivate()*dt;
不会做任何事情。我不知道这怎么可能。我检查了.Derivate()
成员函数是否有效并且它实际上返回了正确的值/对象。最重要的是,我在那里做什么并不重要,例如DP=DoublePendulum RandomObject(1,2,3,4)
仍然会保持 DP 不变。有任何想法吗?我希望我提供了足够的信息...
还有 DoublePendulum 等的实现和声明: