ColumnVector find_top_n(ColumnVector W, int n, double cindex){
	ColumnVector out(n);
	int count=0;
	double max=0.0;
	out(0)=cindex;
	for(int i=0;i<W.rows();i++){
		if(W(i)!=0.0) count +=1;
		if(W(i)<0.0) W(i)=(-1)*W(i);
	}
	out(1)=count;
	int k=2;
	while(k<n){
		out(k)=0;
		max=W(0);
		for(int i=1;i<W.rows();i++){
			if(W(i)>max){
				max=W(i);
				out(k)=i;
			}
		}
		W(out(k))=0.0;
		k +=1;
	}
	return out;
}
 Uniform::Uniform (const ColumnVector& center, const ColumnVector& width)
   : Pdf<ColumnVector> ( center.rows() )
   , _samples(DimensionGet())
 {
   // check if inputs are consistent
   assert (center.rows() == width.rows());
   _Lower = center - width/2.0;
   _Higher = center + width/2.0;
   _Height = 1;
   for (int i=1 ; i < width.rows()+1 ; i++ )
   {
       _Height = _Height / width(i);
   }
 }
 void
 Uniform::UniformSet (const ColumnVector& center,const ColumnVector& width)
 {
   assert(center.rows() == width.rows());
   _Lower = center - width/2.0;
   _Higher = center + width/2.0;
   _Height = 1;
   for (int i=1 ; i < width.rows()+1 ; i++ )
   {
       _Height = _Height / width(i);
   }
   if (this->DimensionGet() == 0) this->DimensionSet(center.rows());
   assert(this->DimensionGet() == center.rows());
 }
 Probability Uniform::ProbabilityGet(const ColumnVector& input) const
 {
   // test if input is located in area of Uniform distribution
   for (int i=1; i<input.rows()+1; i++)
   {
       if ( ( input(i)>_Higher(i) ) || ( input(i) < _Lower(i) ) ) return 0;
   }
   return _Height;
 }
Example #5
0
double pall_dot(RowVector feature, ColumnVector beta){
	double sum=0.0; int nf=beta.rows();
	omp_set_num_threads(36);	
	#pragma omp parallel
	{
		#pragma omp for reduction(+:sum)
		for(int i=0;i<nf;i++)
			sum +=feature[i]*beta[i];
	}
		return sum;
}