Question

我是 c ++ 的新手，我正在尝试修改代码并在名为 Cosmology 的类中使用{c} c库。为了分配成员函数来形成gsl的指针，我使用了通过在网上查找找到的回调程序

更新：Cosmology.h

gsl_integration

Cosmology.cpp

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


struct CCallbackHolder
{
  Cosmology* cls;
  void* data;
};

class Cosmology {
private:
    static const 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);
    static double CCallback(double z,void* param)
    {
      CCallbackHolder* h = static_cast<CCallbackHolder*>(param);
      return h->cls->Hz(h->data);
    }
public:
    double distH, timeH;
    Cosmology();
    Cosmology(double);
    Cosmology(double , double );
    Cosmology(double , double , double );
    Cosmology(double , double , double , double );
    Cosmology(double , double , double , double , double );
    double distC(double);
  } cosmo;

当我尝试编译代码时，我收到以下错误：

#include <cmath>
#include <gsl/gsl_integration.h>
#include "Cosmology.h"
#include<iostream>
using namespace std;
double Cosmology::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 Cosmology::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 = &Cosmology::CCallback;
    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;
    }
using namespace std;    
int main () {
  Cosmology cosmo;
  cout << "Comoving Distance: " << cosmo.distC (0.3);
  return 0;
}

通过以下行编译：{{1}}。我该如何修复代码？

Answer 1

您的程序的主要问题是您在声明它们时初始化类成员的非静态变量，这不是合法的C ++操作。具体来说，问题是这些问题：

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;

c，Mpc2Km，Yrs2Sec，H0，OmegaM，OmegaL，OmegaG，OmegaK 是非静态的，非常数和双重的。如果你想初始化它们（即为它们赋值），你必须在Cosmology类构造函数中执行它，或者将它们移到类声明之外，从而使它们成为全局变量，但请注意，除非你知道什么，否则应该避免全局变量你做到了，别无他法。

如果你想像你想要的那样在头文件中初始化它们，它们必须是整数类型 - 不适用于双打。

很简单：第一行不起作用，因为类型是double（只能在类声明中初始化整数静态const类型）;第二个不起作用，因为类型是非const，非静态和double。

一旦这些线被修复，其余的错误可能会消失。

应用回调来调用gsl_integration

1 个答案: