Можно ли использовать импульсные аккумуляторы с векторами?

Я хотел использовать импульсные накопители для вычисления статистики переменной, которая является вектором. Есть простой способ сделать это. Я думаю, что невозможно использовать тупицу:

  using namespace boost::accumulators;
  //stuff...

  accumulator_set<vector<double>, stats<tag::mean> > acc;
  vector<double> some_vetor;
  //stuff 
  some_vector = doStuff();
  acc(some_vector);

Возможно, это очевидно, но я все равно пробовал.: P

Я хотел иметь аккумулятор, который бы вычислил вектор, который является средним из компонентов многих векторов. Есть ли простой выход?

EDIT:

Я не знаю, был ли я совершенно ясен. Я не хочу этого:

 for_each(vec.begin(), vec.end(),acc);

Это вычислило бы среднее значение элементов данного вектора. Мне нужно другое. У меня есть функция, которая будет плевать векторы:

 vector<double> doSomething(); 
 // this is a monte carlo simulation;

И мне нужно запустить это много раз и вычислить векторное среднее этих векторов:

  for(int i = 0; i < numberOfMCSteps; i++){
  vec = doSomething();
  acc(vec);
  }
  cout << mean(acc);

И я хочу, чтобы mean (acc) был самим вектором, чья запись [i] была бы средством записей [i] накопленных векторов.

Это намек на это в документах Boost, но ничего явного. И я немного тупой.: P

Ответы

Ответ 1

Я немного рассмотрел ваш вопрос, и мне кажется, что Boost.Accumulators уже поддерживает std::vector. Вот что я могу найти в разделе руководства пользователя:

Другой пример, где число Подсистема операторов полезна когда тип не определяет перегрузки оператора, необходимые для его использования для некоторых статистических расчетов. Например, std::vector<> не перегружает никаких арифметических операторов, но может быть полезно использовать std::vector<>как образец или тип переменной. Подкатегория числовых операторов определяет необходимая перегрузка оператора в boost::numeric::operatorsпространство имен, которое вводится в объем с помощью Системы Аккумуляторов с используя директиву.

В самом деле, после проверки файл boost/accumulators/numeric/functional/vector.hpp содержит необходимые операторы для работы "наивного".

Я считаю, вы должны попробовать:

  • Включая либо
    • boost/accumulators/numeric/functional/vector.hpp перед любым другим заголовком аккумуляторов
    • boost/accumulators/numeric/functional.hpp при определении BOOST_NUMERIC_FUNCTIONAL_STD_VECTOR_SUPPORT
  • Приведение операторов в область видимости с помощью using namespace boost::numeric::operators;.

Осталось только одна последняя деталь: выполнение будет прерываться во время выполнения, потому что начальное накопленное значение построено по умолчанию, и при попытке добавить вектор размера n в пустой вектор произойдет утверждение. Для этого кажется, что вы должны инициализировать аккумулятор с (где n - количество элементов в вашем векторе):

accumulator_set<std::vector<double>, stats<tag::mean> > acc(std::vector<double>(n));

Я попробовал следующий код, mean дает мне std::vector размера 2:

int main()
{
    accumulator_set<std::vector<double>, stats<tag::mean> > acc(std::vector<double>(2));

    const std::vector<double> v1 = boost::assign::list_of(1.)(2.);
    const std::vector<double> v2 = boost::assign::list_of(2.)(3.);
    const std::vector<double> v3 = boost::assign::list_of(3.)(4.);
    acc(v1);
    acc(v2);
    acc(v3);

    const std::vector<double> &meanVector = mean(acc);
}

Я считаю, что это то, что вы хотели?

Ответ 2

У меня нет настройки, чтобы попробовать прямо сейчас, но если все boost:: accumulators нуждаются в корректном определении математических операторов, то вы можете уйти с другим векторным типом: http://www.boost.org/doc/libs/1_37_0/libs/numeric/ublas/doc/vector.htm

Ответ 3

А как насчет документации?

// The data for which we wish to calculate statistical properties:
std::vector< double > data( /* stuff */ );

// The accumulator set which will calculate the properties for us:    
accumulator_set< double, features< tag::min, tag::mean > > acc;

// Use std::for_each to accumulate the statistical properties:
acc = std::for_each( data.begin(), data.end(), acc );