类成员函数的 GSL 集成答案

【问题标题】：GSL integration of a class member function类成员函数的 GSL 集成
【发布时间】：2013-06-12 21:26:54
【问题描述】：

我有一个名为 Universe 的课程。该类包括一个计算距离的成员函数，并且需要对一个丑陋的函数进行数值积分。我试图使用 GSL 执行集成，但在尝试编译库时出现以下错误 -

$ g++ -c -O3 -std=c++11 Universe.cpp -o Universe.o
$ error: cannot convert ‘Universe::Hz’ from type ‘double (Universe::)(double, void*)’ to type ‘double (*)(double, void*)’

这是没有构造函数的 Universe 类（为简洁起见）：

宇宙.h

#ifndef UNIVERSE_H
#define UNIVERSE_H

#include <cmath>
#include <gsl/gsl_integration.h>

using namespace std;

class Universe {
private:
    static constexpr double c = 299792458.0, Mpc2Km = 3.08567758e+19, Yrs2Sec = 3.15569e7;
    double H0 = 67.77, OmegaM = (0.022161+0.11889)/(H0*H0), OmegaL = 0.6914, OmegaG = 8.24e-5, OmegaK = 0.0009;
    double Ez(double z);
    double Hz(double z, void* params);
public:
    double distH, timeH;
    Universe() = default;
    Universe(double h0);
    Universe(double omegaM, double omegaL);
    Universe(double h0, double omegaM, double omegaL);
    Universe(double omegaM, double omegaL, double omegaG, double omegaK);
    Universe(double h0, double omegaM, double omegaL, double omegaG, double omegaK);
    //double radius();
    //double age();
    double distC(double z);
    };

#endif

Universe.cpp

#include <cmath>
#include <gsl/gsl_integration.h>
#include "Universe.h"

using namespace std;

double Universe::Hz(double z, void* params)  {
    double result = 1.0/pow(OmegaL + pow(1.0+z,3.0)*OmegaM + pow(1.0+z,4.0)*OmegaG + pow(1.0+z,2.0)*OmegaK, 0.5);
    return result;
    }

double Universe::distC(double z) { 
    double lower_limit = 0.0, abs_error = 1.0e-8, rel_error = 1.0e-8, alpha = 0.0, result, error;
    gsl_integration_workspace *work_ptr = gsl_integration_workspace_alloc(1000);
    gsl_function Hz_function;
    void* params_ptr = &alpha;
    Hz_function.function = Universe::Hz;
    Hz_function.params = params_ptr;
    gsl_integration_qags(&Hz_function, lower_limit, z, abs_error, rel_error, 1000, work_ptr, &result, &error);
    return distH*result;
    }

我不太清楚如何解决这个问题，我是第一次使用 GSL，基于以下文档： http://www.gnu.org/software/gsl/manual/html_node/Numerical-integration-examples.html 和以下指南： http://www.physics.ohio-state.edu/~ntg/780/gsl_examples/qags_test.cpp 感谢您的关注和任何答案！

【问题讨论】：

在 [1] [1] 有一个众所周知的解决方案：stackoverflow.com/questions/13074756/…

标签： c++ class integration gsl

【解决方案1】：

尝试以下操作：将您的 Hz 函数设为静态，如下所示：

static double Hz(double z, void* params)

我已验证这适用于您的代码。

我不是专家，但我相信（以下是手动解释）基本问题是，由于gsl_function 结构需要一个指向函数的指针（即double *），编译器不满意使用尚未实例化的“对象”中的类方法。然而，在类中将函数设为静态意味着编译器甚至在类的实例存在之前就知道可以在哪里找到该函数，并且可以使用指向该函数的指针来处理。（也许有人可以给出比这更好的解释，但希望我在这里不会偏离轨道太远。）

希望这会有所帮助。

【讨论】：

【解决方案2】：

当你更改 HZ 的方法名称时会发生什么？

【讨论】：