Question

我有一个旧的C代码，我正在更新到C ++。我已将所有数组（指向数组的指针）更改为向量。此代码涉及多个文件。相关部分如下。

当我运行它时，我在fn.cpp中计算m [i]的第一次迭代时得到分段错误错误。如果我在语句之外定义m，在减速时，我得到w [i]上的错误。我不确定如何解决这个问题

在main.cpp文件中：

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <iostream>
#include <complex>
#include "fn.h"
using namespace std;

#include <gsl/gsl_rng.h>


// ===========================
// Variables
// ===========================
double eta = 0.13;
double w_max = 3;
int N_bath = 60;
vector<double> m(N_bath);
vector<double> c(N_bath);
vector<double> w(N_bath);

void main(){

bath_para(eta, w_max);

}

在fn.h文件中：

#ifndef FUNCTIONS
#define FUNCTIONS

extern int N_bath;
extern vector<double> c;
extern vector<double> m;
extern vector<double> w;

#endif

在fn.cpp文件中：

#include   <stdlib.h>
#include   <stdio.h>
#include   <math.h>
#include   <iostream>
#include   <complex>
#include   <vector>
#include   "fn.h"
using namespace std;

void bath_para(double eta, double w_max){

  double w_0;
  w_0 = (1 - exp(-w_max))/N_bath;

    for (int i = 0; i < N_bath; ++i){
        m[i] = 1.0
        w[i] = -log(1-(i+1)*w_0);
        c[i] = sqrt(eta*w_0*m[i])*w[i];
    }
}

Answer 1

你可能在其他地方有错误。使用以下fn.h，代码不会产生错误（添加了both_para和namespace的声明，我假设您错误地将其删除为无关紧要）：

#ifndef FUNCTIONS
#define FUNCTIONS
using namespace std;
extern int N_bath;
extern vector<double> c;
extern vector<double> m;
extern vector<double> w;
void bath_para(double eta, double w_max);
#endif

for循环的分段错误

1 个答案: