matlab代码优化直方图C++

Question

LIBIQTOOL_API void Hist(std::vector<double>input, std::vector<double> bins, std::vector<double>& histogram) 
{ 
    double minY = *std::min_element(std::begin(input), std::end(input)); 
    double maxY = *std::max_element(std::begin(input), std::end(input)); 
    std::vector<double> edges;  
    edges.push_back(-1 * std::numeric_limits<double>::infinity()); 
    for (int i = 0; i < bins.size() - 1; i++) 
    { 
     edges.push_back(bins[i] + 0.0100/2); 
    } 
    edges.push_back(std::numeric_limits<double>::infinity()); 


    //histC 

    histogram.resize(edges.size() - 1); 
#pragma omp parallel for 
    for (int i = 0; i < input.size(); i++) 
    { 
     for (int j = 0; j < edges.size() - 1; j++) 
     { 
      if ((edges[j] < input[i]) && (input[i] <= edges[j + 1])) 
      { 
       histogram[j] = histogram[j] + 1; 
       break; 
      } 
     } 
    } 

    histogram[histogram.size() - 1] = histogram[histogram.size() - 1] + histogram[histogram.size() - 2]; 
    histogram.pop_back(); 
} 

输入向量大小为3000000 ++，bin数为7000。
我已经采用了Matlab的Hist（）函数并在C++中创建了我需要的代码。但是它需要很长的时间才能运行，你能看到更多的运行时优化，可以在这里完成吗？
我做过：
a。 b。当你发现箱子放置当前编号
b。使用OpenMP的

Answer 1

可能的优化：

• 不按值传递你的输入数据，但常引用
• 不检查下界，做了线性搜索时只上限为每个仓正确的bin。
• 或者：由于您的垃圾箱单调排列且没有间隙，请对正确的垃圾箱进行二进制搜索，而不是线性搜索。

最后一个应该给你最大的收益，其他的更实施。

顺便说一句，你填充边缘向量看起来很奇怪。