矢量和阵列通过误差

我正在编写一个基于多链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)

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