Question

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

C：\ Documents \ C code \ RK4 \ addRK4.h | 7 | error：expected'）'before'（'token |

之后还有一堆其他错误消息，但我不认为它们很重要。我无法弄清楚什么是错的，特别是因为当我在main.cpp中定义RK4的原型时，一切都运行得很好。相关代码如下。任何有关此事的帮助（或者，如果您对c ++很新，我有任何其他建议）将不胜感激！

的main.cpp

#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);
}

}

addRK4.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;
}

fitzhough.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;

}

Answer 1

你有一个符号冲突。

您#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

Answer 2

~~types Vector2d和VectorXd似乎有问题。~~

编辑：好抓@Paddy。留下其余的答案，因为它仍然有效。

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

将函数传递给函数并使用头文件时出错

2 个答案: