真的在努力解决这个问题。
我正在编写一段代码来使用前向欧拉方法来解决线性ODE的一般系统。
我相信我的代码在类中运行良好但我在初始化时遇到问题。
这是我的ForwardEulerSolver类的构造函数:
ForwardEulerSolver(ODEInterface& anODESystem, const Vector& initialState,
const double initialTime, const double finalTime, const double stepSize,
Vector (*pFunction)(double time, Vector u), const std::string outputFileName
= "output.dat", const int saveGap = 1, const int printGap = 1);
我认为问题在于Vector(* pFunction)(双倍时间,Vector u)。有问题的函数应该是一个需要时间“t”和向量u的函数,并返回向量u ^ {dot}(即:du / dx)。
然后在我的驱动程序中我定义:
Vector FunctionToSolve(double time, Vector U)
{
Vector Output(1);
Output[0] = 1/(1 + pow(time,2));
return Output;
}
作为一个单独的函数(我知道向量U在这里未被使用,问题是1D,但我希望选项以超过1维和更复杂的f执行此操作)。然后我有:
Vector (*pfunction)(double, Vector);
pfunction = &FunctionToSolve;
在int main()。
然后我尝试按以下方式初始化ForwardEulerClass:
ForwardEulerSolver* EulerSystem(*pODESystem, *initialVector, 0, 10, 0.01, *pfunction);
并在该行上收到以下错误:
error: cannot convert 'Vector (*)(double, Vector)' to 'ForwardEulerSolver*' in initialization
我也收到了大量警告(全部在同一条线上):
warning: right operand of comma operator has no effect [-Wunused-value]
warning: value computed is not used [-Wunused-value]
和
error: expression list treated as compound expression in initializer [-fpermissive]
任何帮助将不胜感激!谢谢!
答案 0 :(得分:1)
首先,您尝试初始化ForwardEulerSolver
的实例并将其分配给指针类型变量。这在语法上是不正确的。你可能想做这样的事情:
ForwardEulerSolver eulerSystem(*pODESystem, *initialVector, 0, 10, 0.01, pfunction);
或许这个:
ForwardEulerSolver* eulerSystem = new ForwardEulerSolver(*pODESystem, *initialVector, 0, 10, 0.01, pfunction);
我改变了EulerSystem
的字母大小写,只是想确保你的意思是变量名,而不是某些派生类名。
另请注意,您不希望derefenrece函数指针。你实际上想要传递指针,而不是当场评估函数。所以我在*
之前删除了pfunction
。