矢量和阵列通过误差

vector and array passing error

本文关键字:误差 阵列      更新时间:2023-10-16

我正在编写一个基于多链Metropolis算法的程序,以找到给定函数的最小能量参数。这个程序基本上运行两次:

  1. 第一次多链老化运行,扫描参数空间,并在满足给定收敛标准时停止
  2. 探索性运行,对第一步中计算的收敛点周围的各种参数可能性进行采样,但仅使用一条链

对于每条链,我保存它通过的所有参数及其相应的能量。每条链的每一步都运行(除其他外(以下代码:

int Step(double * current_energy, vector<double> & energy){
    energy.push_back(current_energy);
    return 0;
}

问题是,在程序的第一部分,我将能量存储在中

vector< vector<double> > burn_energy;

让我可以在一个地方保存多个系列的能量。在第二部分中,我将能量存储在中

vector<double> explore_energy;

因为我只需要节省一系列的能量。

当我不得不调用这些步骤时,问题就出现了。对于第一部分,我会这样称呼他们:

Step(& current_energy[j], burn_energy[j]);

其中j表示链的数目。在第二部分中,由于只有一条链,调用将类似于:

Step(& current_energy, explore_energy);

我看不出我的实现有什么问题,但XCode在第一部分中给了我以下错误:

error: no matching function for call to 'std::vector<double, std::allocator<double> >::push_back(double*&)'

这是矢量能量或阵列current_energy的元素通过的误差吗?

int Step(double * current_energy, vector<double> & energy)

应该是

int Step(double current_energy, vector<double> & energy)

编辑:

在OP发表评论后,我建议这样做;

int Step(double& current_energy, vector<double> & energy)
             ^^^^   

在这种情况下,您可以在不使用任何指针语法的情况下更改当前能量。它被称为参考。

energydoubles的向量,但您正试图将double*推到它上面。不要这样做。