是否有gsl_vector的count_if函数？ C / C ++

Question

我正在使用gnu科学图书馆（GSL）。 假设我有一个像这样的gsl_vector ：7 0 -6 5 8 0 10 -2

这是一个包含正数，负数和零作为元素的向量。

我想计算这个gsl_vector非零元素或零元素的gsl_vector 。 我知道有一个名为count_if的函数用于C ++向量。 但我搜索gsl_vector.h和gsl_blas.h ，没有相同的功能。 我可以通过gsl_vector_get()评估所有元素，然后问问题。

int counter = 0;

for(int i = 0;i<length_of_the_gsl_vector;++i){ 
  if(fabs(gsl_vector_get(y,i))<0.5) ++counter;
}

return counter;

但是我几乎花了一天时间想知道GSL中是否已经存在这样一种效率更高的功能。

或者gsl_array可能有一个count_if函数？

Answer 1

您可以使用gsl_vector_ptr获取数据指针，然后在指针上使用std::count_if ：

struct Predicate{
    inline bool operator()(double x) const {
        return fabs(x) < 0.5 ;
    }
} ;

int res = std::count_if( 
    gsl_vector_ptr(y,0), gsl_vector_ptr(y,0) + size, 
    Predicate()
) ;

Answer 2

他们陈述的是GSL实现的std :: valarray的模型，我读到这样的语句作为授权我使用的东西，如：

gsl_block vec;
/* initialize ... */
std::valarray<double> a(vec.data, vec.size);
/* use std:: galore */