我正在为一项研究任务编写遗传算法。我不是一个非常有经验的C ++程序员(因为我是一名数学家),但我管理我的程序在Windows环境中使用MSVC 2008和g ++编译器正确编译和运行(我使用CygWin在Windows 7上运行g ++)。现在,问题是该程序必须在一个在Linux环境中使用g ++的集群计算机上运行,但它总是在执行时崩溃(尽管至少它正确编译)。我已经解决了明显的细节(使用/而不是文件名),但无法在Linux上正常运行。
然而,当我评论下面的函数'Cruce'时,程序结束,尽管它显然没有执行所需的任务。这是这样一个函数的代码,以及它调用的其他函数(唯一的全局变量是BITGEN,它是一个值为30的整数,用于浮点数的格雷码转换):
void xPC_BLX(double d, const vector<double>& P1, const vector<double>& P2,
vector<double>& Hijo1, vector<double>& Hijo2, int genes){
double I, A1, C1;
int i;
for (i=0; i < genes; i++)
{
I= d * fabs(P1[i]-P2[i]);
A1=P1[i]-I; if (A1<0) A1=0.0;
C1=P1[i]+I; if (C1>1) C1=1.0;
Hijo1[i]= A1 + Rand()*(C1-A1);
A1=P2[i]-I; if (A1<0) A1=0.0;
C1=P2[i]+I; if (C1>1) C1=1.0;
Hijo2[i]= A1 + Rand()*(C1-A1);
}
}
/**********************************************************/
/* Itoc and Ctoi translate ints to strings and vice versa */
/**********************************************************/
unsigned long int Ctoi(char *Cad_ent, int length){
int i;
unsigned long n;
n = (unsigned long) 0;
for (i=0; i<length; i++)
{
n <<= 1;
n += (*Cad_ent++ - (int) '0');
}
return(n);
}
void Itoc(unsigned long int n, char *Cad_sal, int length){
int i;
for (i=length-1; i>=0; i--)
{
Cad_sal[i] = (char)('0' + (n & 1));
n >>= 1;
}
}
/*****************************************************************/
/* Translations between fixed point ints and reflected Gray code */
/*****************************************************************/
void Gray(char *Cad_ent, char *Cad_sal, int length){
int i;
char last;
last = '0';
for (i=0; i<length; i++)
{
Cad_sal[i] = (char)('0' + (Cad_ent[i] != last));
last = Cad_ent[i];
}
}
/*************************************************************************/
/* Translations between string representation and floating point vectors */
/*************************************************************************/
void StringRep(const vector<double> vect, char *Cad_sal, int genes){
int i;
unsigned long int n;
int pos;
double INCREMENTO;
static char *tmpstring;
static int flag = 1;
if (flag) {
tmpstring = (char*) calloc (genes*BITGEN,sizeof(char));
flag = 0;
}
pos = 0;
for (i=0; i < genes; i++)
{
INCREMENTO=(1-0)/(pow(2.0, (double) BITGEN) - 1.0);
n = (int) ((vect[i] - 0) / INCREMENTO + 0.5);
Itoc(n, tmpstring, BITGEN);
Gray(tmpstring, &Cad_sal[pos], BITGEN);
pos += BITGEN;
}
Cad_sal[pos] = '\0';
}
/*****************************************/
int DistHam(char *Cr_1, char *Cr_2, int genes){
int i, dist;
dist=0;
for (i=0; i<genes*BITGEN; i++) if (Cr_1[i]!=Cr_2[i]) dist++;
return dist ;
}
/**********************/
/* CROSS OPERATOR */
/**********************/
void Cruce(int& fin, int genes, vector<vector<double> >& POPULATION,
vector<vector<double> >& CONTROL, double GA_THR)
{
int i, j, temp, mom, dad;
static char *String1, *String2;
static int flag=1;
vector<double> newind(genes+1,0), newcont(2,0);
if (flag) {
String1 = (char*) calloc (genes*BITGEN,sizeof(char));
String2 = (char*) calloc (genes*BITGEN,sizeof(char));
flag = 0;
}
vector<int> sample(TAMPOP,0);
for (i=0; i < TAMPOP; i++) sample[i] = i;
for (i=0; i < TAMPOP; i++)
{
j = Randint(i,TAMPOP-1);
temp = sample[j];
sample[j] = sample[i];
sample[i] = temp;
}
for (i=0; i < TAMPOP; i++) CONTROL[i][0]=0;
fin=TAMPOP;
for (i=0; i < TAMPOP/2; i++)
{
mom=sample[2*i];
dad=sample[2*i+1];
StringRep(POPULATION[mom], String1, genes);
StringRep(POPULATION[dad], String2, genes);
if (DistHam(String1, String2, genes)/2.0 > GA_THR)
{
POPULATION.push_back(newind);
POPULATION.push_back(newind);
CONTROL.push_back(newcont);
CONTROL.push_back(newcont);
xPC_BLX(1,POPULATION[mom],POPULATION[dad],POPULATION[fin],
POPULATION[fin+1], genes);
CONTROL[fin][0]=1;
CONTROL[fin+1][0]=1;
fin=fin+2;
}
}
}
我已经尽力了好几天才能找到可能的错误,但这是徒劳的。我猜错误必须与内存问题有关(因为我总是得到分段错误错误),但我无法确定它在哪里。 Windows和Linux处理内存的方式之间应该存在一个我不了解的关键区别。你们中的任何一个可以帮助我吗?
提前致谢!
答案 0 :(得分:2)
一些可能对您有帮助的一般提示:
使用-Wall -Wextra
选项编译到g++
并改进源代码,直到所有警告都消失为止
还尝试使用-Wall
编译器与-g
(或类似)选项进行编译
使用-Wall
选项进行编译(除了bt
)以获取调试信息
使用clang++调试您的程序(并在调试器下运行程序,然后使用calloc
进行回溯等...)
使用gdb查找内存泄漏
避免使用malloc
和std::string
,更好地使用标准C ++容器(例如std::vector<>
,std::map<>
,{{1}}等。)。
此外,在您自己的笔记本电脑或台式机上安装Linux发行版,并学习通过命令行使用Linux。
答案 1 :(得分:2)
我发现代码中存在潜在错误。
static char *tmpstring;
static int flag = 1;
if (flag) {
tmpstring = (char*) calloc (genes*BITGEN,sizeof(char));
flag = 0;
}
如果使用不断增长的StringRep()
参数多次调用函数genes
,会发生什么:缓冲区不会增长,因此有人会在分配的第一个缓冲区之外写入。请注意,同样适用于函数Cruce()
。
然后我认为关于fin=TAMPOP;
的第二个错误。什么是TAMOP
没有解释。但是,在执行POPULATION[fin+1]
,CONTROL[fin+1][0]=1;
或fin=fin+2;
时,您是不是在数组之外。
问候。