#include <iostream>
#include <vector>
#include <algorithm>

typedef std::vector<double> vec;

void rand_vector(vec& v){
  generate(v.begin(),v.end(), [](){return rand()/((double)RAND_MAX);});  
}

int vector_chop(vec& v, const double lim){
  int count = 0;
  for_each(v.begin(),v.end(), [&count,lim](double& x){if(x<lim) {x=0.0; count++;}});  
  return count;
}

int main()
{
   vec v(7);
   rand_vector(v);
   for(auto x:v) std::cout<<x<<"\n";
   const double limit=0.4;
   std::cout<<"chopped "<<vector_chop(v,limit)<<" values below "<<limit<<"\n";
   for(auto x:v) std::cout<<x<<"\n";
}