arpackのためにcからfortran 77を呼び出す

Question

ARPACKからFortran関数dsaupdを呼び出そうとしています。私は、コードがコンパイルされ、最初の「ここ」プリントアウトにします

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

でdsaupd呼び出す前に、私は、以前のプログラムのintとしてNUMCOLSを定義しnetlib-java

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

からC宣言を使用しましたしかし、それ以降はストールします。どのようなアイデアを私は間違ってやっているか、またはこの呼び出しをデバッグする方法？

Answer 1

私のシステムでarpackがコンパイルされていても、間違ってリンクしていたと思います。 arpack-ngに切り替えて問題を修正しました。