Спасибо, но используя векторное представление, я хотел избежать временного вектора, потому что, если начальный вектор очень большой, я потеряю время, копируя его во временный ... Или это "tempForStats (stl_v.begin (), stl_v.end ());» очень эффективный?
рограмма манипулируетSTL векторы целых чисел но время от времени мне нужно посчитать немного статистики по ним. Поэтому я используюФункции GSL, Чтобы избежать копирования вектора STL в вектор GSL, я создаюGSL векторный види передайте его функциям GSL, как в этом фрагменте кода:
#include <iostream>
#include <vector>
#include <gsl/gsl_vector.h>
#include <gsl/gsl_statistics.h>
using namespace std;
int main( int argc, char* argv[] )
{
vector<int> stl_v;
for( int i=0; i<5; ++i )
stl_v.push_back( i );
gsl_vector_int_const_view gsl_v = gsl_vector_int_const_view_array( &stl_v[0], stl_v.size() );
for( int i=0; i<stl_v.size(); ++i )
cout << "gsl_v_" << i << "=" << gsl_vector_int_get( &gsl_v.vector, i ) << endl;
cout << "mean=" << gsl_stats_mean( (double*) gsl_v.vector.data, 1, stl_v.size() ) << endl;
}
После компиляции (gcc -lstdc ++ -lgsl -lgslcblas test.cpp) этот код выдает следующее:
gsl_v_0=0
gsl_v_1=1
gsl_v_2=2
gsl_v_3=3
gsl_v_4=4
mean=5.73266e-310
Векторное представление создано правильно, но я не понимаю, почему среднее значение неверно (оно должно быть равно 10/5 = 2). Любая идея? Заранее спасибо.