【发布时间】:2011-10-07 06:23:34
【问题描述】:
我正在编写一个基于多链 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的错误吗?
【问题讨论】:
标签: c++ vector multidimensional-array parameter-passing