我正在进行模拟,我必须计算很多平均值,我认为使用boost :: accumulators是一个好主意。问题是我想要平均的量之一是对称矩阵,其对角线事先已知。因此,如果i
起初我觉得我可以使用
using namespace boost::accumulators;
using namespace boost::numeric::ublas;
typedef accumulator_set<double, stats<tag::mean> > accumulator;
symmetric_matrix<accumulator, lower> foo; // a symmetric matrix of accumulators
拿着我的蓄电池。但后来我发现这个symmetric_matrix结构可能只适用于保存数值(它们定义了算术运算)或者以某种方式针对这种数据进行了优化。这是对的吗?
如果来自boost的symmetric_matrix不合适,我需要一个数据结构,它可以保持对称矩阵的下三角形而没有对角线,并且它必须适合保持累加器并具有类似矩阵的语法。这是从一些图书馆随处可用的吗?如果没有,这种结构是否容易实现?