Question

我试图从ARPACK调用fortran函数dsaupd。我使用了来自netlib-java

的C声明

  extern void dsaupd_(int *ido, char *bmat, int *n, char *which,
                       int *nev, double *tol, double *resid,
                       int *ncv, double *V, int *ldv,
                       int *iparam, int *ipntr, double *workd,
                       double *workl, int *lworkl, int *info);

然后我在使用

调用dsaupd之前将numcols定义为程序中的int

int ido = 0;
int ncv = 2*numeigs;
int maxiter = 30;
double tol = 1e-13;
double * v = (double *) malloc(numcols * ncv *sizeof(double));
int iparam[11] = {1, 0, maxiter, 1, 0, 0, 1, 0, 0, 0, 0};
int ipntr[11];
double * workd = (double *) malloc(3*numcols*sizeof(double));
int lworkl = ncv*(ncv + 8);
double * workl = (double *) malloc(lworkl*sizeof(double));
int arpack_info = 0;

char bmat = 'I';
char which[2] = {'L', 'M'};
MPI_Barrier(comm);
if (mpi_rank == 0) {
    printf("Here!\n");
    dsaupd_(&ido, &bmat, &numcols, which,
            &numeigs, &tol, vector, 
            &ncv, v, &numcols,
            iparam, ipntr, workd, 
            workl, &lworkl, &arpack_info);
    printf("Here!\n");
}

代码编译并转到第一个＆＃34; Here＆＃34;打印输出，但在那之后停止。知道我做错了什么，或者如何调试这个电话？

Answer 1

我认为即使arpack正在我的系统上编译，它也是错误的链接。切换到arpack-ng修复了问题。

从c调用fortran 77 for arpack

1 个答案: