0
Я пытаюсь вызвать функцию fortran dsaupd из ARPACK. Я использовал декларацию C от netlib-javaвызов fortran 77 из c для arpack
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);
тогда я определила numcols как ИНТ ранее в программе перед вызовом 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");
}
код компилирует и делает его первым «Здесь» распечатке, но киосков после этого. Любая идея, что я делаю неправильно, или как отладить этот вызов?
1) Не отбрасывайте результат 'malloc' & friends в C. 2) Инициализаторы для составных типов должны быть постоянными. – Olaf
Я дал обе эти попытки, но я все еще получаю стойло. инициализатор имеет смысл (хотя здесь maxiter является константой, поэтому не должно быть проблемой), но почему я не должен приводить malloc к правильным типам? – AatG
'maxiter' окончательно не является константой, а переменной. Даже если он был квалифицирован 'const', он не был бы постоянным. C не является C++, он не имеет символических констант, кроме очень ограниченных констант _enum. – Olaf