0

当我尝试使用头文件添加函数 RK4 时尝试运行以下代码时收到以下错误消息。

C:\Documents\C 代码\RK4\addRK4.h|7|错误：在 '(' 标记 | 之前应为 ')'

之后还有很多其他错误消息，但我认为它们并不重要。我不知道出了什么问题，尤其是当我在 main.cpp 中定义 RK4 的原型时，一切运行良好。相关代码如下。对此问题的任何帮助（或者如果您有任何其他建议，因为我对 c++ 很陌生）将不胜感激！

主文件

#include <iostream>
#include <fstream>
#include <Eigen/Dense>
#include "gnuplot.h"
#include "addfitzhough.h"
#include "addRK4.h"

using namespace std;
using namespace Eigen;

int main()
{

//int mydims = 2;

double u = 0;
double *Iion;
double h = .5;

double y1ans[800];
double y2ans[800];
double tans[800];


Vector2d ycurr;

Vector2d Ynot, yplus;

Ynot << .2,
        .1;

y1ans[0] = Ynot(0);
y2ans[0] = Ynot(1);
tans[0] = 0.0;

for(int i = 1;i<800;i++){
tans[i] = tans[i-1] + h;
ycurr << y1ans[i-1],
         y2ans[i-1];
yplus = RK4(fitzhough,tans[i],ycurr,h,u,Iion,2);
y1ans[i] = yplus(0);
y2ans[i] = yplus(1);
}

}

添加RK4.h

#ifndef RK4
#define RK4

using namespace Eigen;

VectorXd RK4(VectorXd (*f) (double t, Vector2d Y, double u, double * Iion), double t, VectorXd z, double h, double u, double *Iion, int d);


#endif // RK4

RK4.cpp

#include <Eigen/Dense>

using namespace std;
using namespace Eigen;

Vector2d RK4(Vector2d (*f)(double, Vector2d, double, double*), double t, VectorXd z, double h, double u, double *Iion, int d){

VectorXd Y1(d), Y2(d), Y3(d), Y4(d), Y1buf(d), Y2buf(d), Y3buf(d);

Y1 = z;
Y1buf = (*f)(t,Y1,u, Iion);
Y2 = z + 0.5*h*Y1buf;
Y2buf = (*f)(t+.5*h,Y2,u, Iion);
Y3 = z + 0.5*h*Y2buf;
Y3buf = (*f)(t+.5*h,Y3,u, Iion);
Y4 = z + h*Y3buf;


Vector2d yn = z + (h/6.0)*(Y1buf + 2.0*Y2buf + 2.0*Y3buf + (*f)(t+h,Y4,u, Iion));

return yn;
}

菲州.cpp

#include <Eigen/Dense>

using namespace std;
using namespace Eigen;
Vector2d fitzhough(double t, Vector2d Y, double u, double * Iion){

Vector2d dy;

double v = Y(0);
double w = Y(1);

double a = .13;
double b = .013;
double c1 = .26;
double c2 = .1;
double d = 1.0;

dy(0) = c1*v*(v-a)*(1-v)-c2*w*v + u;
dy(1) = b*(v-d*w);

*Iion = dy(0)-u;

return dy;

}
4

2 回答 2

2

你有一个符号冲突。

您#define是符号RK4，然后您尝试使用该名称创建一个函数。因为您已将其定义为空宏，所以它将被替换为空。然后编译器将此视为您的函数声明：

VectorXd (VectorXd (*f) (double t, Vector2d Y, double u, double * Iion), double t, VectorXd z, double h, double u, double *Iion, int d);

为您的头笼添加额外的字符是个好主意。就像是：

#ifndef RK4__H
#define RK4__H
于 2013-01-11T03:07:57.610 回答
1

类型 Vector2d和似乎有问题VectorXd。

编辑：很好的赶上@Paddy。将其余答案留在这里，因为它仍然有效。

此外，您在中的声明与中addRK4.h 的定义不匹配RK4.cpp。这将是您要修复的下一个错误。

于 2013-01-11T03:07:14.240 回答