AGMG运行时出现分段故障

Question

我已经将一个C++程序与AGMG库中的fortran dagmg.f90例程链接起来。最初，稀疏矩阵是由特征库以CRS格式构建的，然后我将其转换为fortran例程所需的格式。AGMG运行并计算非零的数量，并显示未知的数量。但随后它突然显示分段错误。我不知道背后的原因是什么。

    int jatest[NNZ];
int iatest[nodes+1];
double bftest[nodes];
double VXtest[nodes];
// code to convert Af into a, ja, ia vectors for AGMG
double *aptr; int* japtr; int* iaptr;
aptr = Af.valuePtr();
japtr = Af.innerIndexPtr();
iaptr = Af.outerIndexPtr();
for (i = 0; i < NNZ ; i++ ){
  atest[i] = aptr[i];
  jatest[i] = japtr[i] + 1;
 // cout << atest[i] << "t" << japtr[i] << endl;
}
for ( i = 0; i <= nodes; i++){
  iatest[i] = iaptr[i] + 1;
   }
 dagmg_(nodes,atest,jatest,iatest,bftest,VXtest,ijob,iprint,nrest,iter,tol);

输出如下：

*进入AGMG*

未知数：17**非零：51（每行：3.00）

分段故障

为什么我会犯这个错？？？

Answer 1

问题已解决：

事实上，问题在于AGMG调用LAPACK函数以直接解决最粗略的级别。在改变了链接到AGMG的lapack库之后，解决了分割故障的问题。

仅供参考：

g++ *.o ~/lapack-3.4.1/liblapack.a ~/CLAPACK-3.2.1/lapack_LINUX.a ~/mylibs/blas_LINUX.a ~/CLAPACK-3.2.1/blas_LINUX.a ~/CLAPACK-3.2.1/F2CLIBS/libf2c.a ~/lapack-3.4.1/liblapack.a -o myexe -lgfortran