Question

我尝试开发Android应用。我已经设法使用NDK实现了一个相当小的c ++函数。但现在我想使用我在Google Docs上找到的相当大的c ++库。当我现在尝试使用ndk-build工具构建它时，我在cygwin控制台中收到此错误：

Compile++ thumb  : ndkfoo <= subspace.cpp
Compile++ thumb  : ndkfoo <= classifier.cpp
Compile++ thumb  : ndkfoo <= eigen.cpp
Compile++ thumb  : ndkfoo <= image.cpp
Compile++ thumb  : ndkfoo <= imageio.cpp
Compile++ thumb  : ndkfoo <= local.cpp
Compile++ thumb  : ndkfoo <= matrix.cpp
Compile++ thumb  : ndkfoo <= sample.cpp
SharedLibrary  : libndkfoo.so
./obj/local/armeabi/objs/ndkfoo/eigen.o: In function `LibSubspace::geneigen(doub
le*, double*, int, double*, int, bool)':
U:\workspace\test/jni/eigen.cpp:91: undefined reference to `ilaenv_'
U:\workspace\test/jni/eigen.cpp:96: undefined reference to `dsygv_'
U:\workspace\test/jni/eigen.cpp:128: undefined reference to `dggev_'
./obj/local/armeabi/objs/ndkfoo/eigen.o: In function `LibSubspace::eigen(double*
, double*, double*, int, bool)':
U:\workspace\test/jni/eigen.cpp:55: undefined reference to `ilaenv_'
U:\workspace\test/jni/eigen.cpp:58: undefined reference to `dsyev_'
collect2: ld returned 1 exit status
make: *** [obj/local/armeabi/libndkfoo.so] Error 1

此功能的代码如下：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <float.h>
#include <assert.h>
#include <string.h>

#include "eigen.h"

extern "C" int ilaenv_(int *ispec, const char *name__, const char *opts,
    int *n1, int *n2, int *n3, int *n4, int name_len, int opts_len);

extern "C" int dsyev_(const char *jobz, const char *uplo, int *n, double *a,
    int *lda, double *w, double *work, int *lwork, 
    int *info);

extern "C" int dsygv_(int *itype, const char *jobz, const char *uplo, int *
    n, double *a, int *lda, double *b, int *ldb, 
    double *w, double *work, int *lwork, int *info);

extern "C" int dggev_(const char *jobvl, const char *jobvr, int *n, double *
    a, int *lda, double *b, int *ldb, double *alphar, 
    double *alphai, double *beta, double *vl, int *ldvl, 
    double *vr, int *ldvr, double *work, int *lwork, 
    int *info);

namespace LibSubspace {

int eigen(double* A, double* V, double* E, int n, bool verbose) {
    int info;
    int ispec = 1;
    int lwork;

    lwork = (ilaenv_(&ispec,"DSYEV","U",&n,&n,&n,&n,5,1)+2)*n;
    double *work = (double *)malloc(lwork*sizeof(double));
    memcpy(V,A,n*n*sizeof(double));
    dsyev_("V","U",&n,V,&n,E,work,&lwork,&info);
    free(work);

    //check the return value
    if(info!=0) {
        if(verbose) {
            printf("Error computing eigenvectors: ");
            if(info>0) {
                printf("Algorithm failed to converge\n");
            } else if(info<0) {
                printf("Illegal argument\n");
            } else {
                printf("Unknown error\n");
            }
        }
        return 0;
    }

    return 1;
}


int geneigen(double *A, double *B, int n, double *W, int algorithm,
             bool verbose) {
    //getting the optimal lwork
    int ispec = 1;
    int lwork;
    int info;

    switch(algorithm) {

        case EIGEN_CHOL: 
        {
            //computing eigenvectors
            lwork = (ilaenv_(&ispec,"DSYGV","U",&n,&n,&n,&n,5,1)+2)*n;
            double *work = (double *)malloc(lwork*sizeof(double));
            int problemType = 1;
            char job = 'V';
            char uplo = 'U';
            dsygv_(&problemType,&job,&uplo,&n,A,&n,B,&n,W,work,&lwork,&info);

            free(work);

            //check the return value
            if(info!=0) {
                if(verbose) {
                    printf("Error computing eigenvectors: ");
                    if(info>n) {
                        printf("Matrix B is not positive definite\n");
                    } else if(info<=n) {
                        printf("The problem failed to converge\n");
                    } else if(info<0) {
                        printf("Illegal argument\n");
                    } else {
                        printf("Unknown error\n");
                    }
                }
                return 0;
            }
        }
        break;

        case EIGEN_QZ:
        {
            //more general algorithm
            double *alphar = (double *)malloc(n*sizeof(double));
            double *alphai = (double *)malloc(n*sizeof(double));
            double *beta = (double *)malloc(n*sizeof(double));
            double *VR = (double *)malloc(n*n*sizeof(double));
            lwork = 8*n;
            double *work = (double *)malloc(lwork*sizeof(double));
            dggev_("N","V",&n,A,&n,B,&n,alphar,alphai,beta,
                                   NULL,&n,VR,&n,work,&lwork,&info);
            //eigenvalues
            for(long i=0;i<n;i++) {
                if(beta[i]!=0) {
                    W[i] = alphar[i]/beta[i];
                } else W[i] = 0;
            }
            //eigenvectors
            for(long i=0;i<n;i++) {
                for(long j=0;j<n;j++) {
                    A[i*n+j] = VR[i*n+j];
                }
            }
            free(alphar);
            free(alphai);
            free(beta);
            free(VR);   
            free(work);

            if(info!=0) {
                printf("Error computing eigenvectors: ");
                if(info<0) {
                    printf("Illegal argument\n");
                } else if(info<=n) {
                    printf("QZ iteration failed\n");
                } else {
                    printf("Unknown error\n");
                }
                return 0;
            }
        }
        break;

    }

    return 1;
}

} //namespace

不幸的是，我绝对不明白为什么会出现这些错误。我对C ++并不熟悉，这就是为什么我想使用NDK，希望不必改变现有库的任何内容。

Answer 1

谷歌告诉我这些函数来自Fortran-to-C编译器的运行时。因此，您需要链接到名为libf2c的库，或者找到源并使它们成为NDK项目的一部分。后者听起来像是一种更有前途的方法 - 我认为不存在现成的NDK构建。

NDK未定义对{extern“C”函数}的引用？

1 个答案: