我正在编写一个基于多链 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*&)'
这是传递 vector 能量或数组 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)
^^^^
在这种情况下,您无需任何指针语法即可更改当前能量。这称为引用。
关于c++ - vector 和数组传递错误,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/6726315/