我正在编写一个基于多链Metropolis算法的程序,以找到给定函数的最小能量参数。这个程序基本上运行了两次:
对于每个链,我保存它通过的所有参数及其相应的能量。每个链的每个步骤(以及其他内容)运行以下代码:
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
的错误吗?
答案 0 :(得分:2)
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)
^^^^
在这种情况下,您可以在没有任何指针语法的情况下更改当前能量它被称为参考。
答案 1 :(得分:1)
energy
是double
s的向量,但您尝试将double*
推到其上。不要那样做。