0

我有一个在 Visual Studio 2012 中使用 win 32 平台编译和运行的程序。将平台更改为 X64 后,它给出以下错误消息:

" test_integration.cpp(27): 错误 C2664: 'hcubature_v' : 无法将参数 2 从 'int (__cdecl *)(unsigned int,unsigned int,const double *,void *,unsigned int,double *)' 转换为 'integrand_v '"

有人可以提供一些建议吗?谢谢!

#include <iostream>
#include "cubature.h"
#include <ctime>
#include <limits>
using namespace std;

int f_v(unsigned ndim, unsigned npts, const double *x, void *fdata, unsigned fdim, double *fval) {
    double sigma = *((double *) fdata);
    unsigned i, j;
    for (j = 0; j < npts; ++j) { // evaluate the integrand for npts points
        double sum = 0;
        for (i = 0; i < ndim; ++i) sum += x[j*ndim+i] * x[j*ndim+i];
        fval[j] = exp(-sigma * sum);
    };
    return 0; // success
}

int main(){

    double xmin[3] = {-2,-2,-2}, xmax[3] = {2,2,2}, sigma = 0.5;
    double  val_v,err_v, time_v;
    const clock_t begin_time_v = clock();

    for (int j = 0; j<1000;j++)
    {
    // hcubature_v calculates a multidimensional integration.
    hcubature_v(1, f_v, &sigma, 1, xmin, xmax, 0, 0, 1e-4, ERROR_INDIVIDUAL, &val_v, &err_v);
    }

    cout << float( clock () - begin_time_v ) /  CLOCKS_PER_SEC<<endl;
    printf("Computed integral   = %0.10g +/- %g\n", val_v, err_v);
    cin.get();
}

以下是一些额外的定义:

/* a vector integrand of a vector of npt points: x[i*ndim + j] is the
   j-th coordinate of the i-th point, and the k-th function evaluation
   for the i-th point is returned in fval[i*fdim + k].  Return 0 on success
   or nonzero to terminate the integration. */
typedef int (*integrand_v) (unsigned ndim, size_t npt,
                const double *x, void *,
                unsigned fdim, double *fval);

/* as hcubature, but vectorized integrand */
int hcubature_v(unsigned fdim, integrand_v f, void *fdata,
        unsigned dim, const double *xmin, const double *xmax, 
        size_t maxEval, double reqAbsError, double reqRelError, 
        error_norm norm,
        double *val, double *err);
4

1 回答 1

0

您应该使用用于定义的确切类型integrand_v

代替:

int f_v(unsigned ndim, unsigned npts, const double *x, void *fdata, unsigned fdim, double *fval) {
                       ^^^^^^

采用

int f_v(unsigned ndim, size_t npt,  const double *x, void *fdata, unsigned fdim, double *fval) {
                       ^^^^^^

它可能在 32 位下工作,因为size_t并且unsigned恰好是同一类型。

它不适用于 64 位,因为它们是不同的类型。

于 2015-09-22T20:13:58.960 回答