 
 }
 
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
 BevParams::BevParams(float bev_res, tuple_int bev_xLimits, tuple_int bev_zLimits, tuple_int imSize)
 {
 	

		index++;
	}

	//print1dvector(z_vec_length);
	//print1dvector(x_vec_length);
	
	if(debug)
	{

		}
	}
	
	int numCol = this->numBColumns;

	row_t xi_1(numCol), yi_1(numCol);

	float* result_row_0 = result;
	float* result_row_1 = result + 1*numCol;
	
	vector<int> x_index_vec_copy = this->x_index_vec;
	vector<int> z_index_vec_copy = this->z_index_vec;
	
	vector<int> z_vec_sel(numCol), x_vec_sel(numCol);

	/*
	double e3 = getTickCount();
	
	device_vector<float> result_row_0(result, result + numCol);
	device_vector<float> xi_1(numCol);

	auto it = thrust::copy_if(result_row_0.begin(),result_row_0.end(),xi_1.begin(), checkSanity());
			
	host_vector<float> h_xi_1 (xi_1.begin(),it);
	/

	//double e4 = getTickCount();
	//double time = (e4-e3)/getTickFrequency();
//	cout<<"Time for comapction"<<time<<endl;



	device_vector<float> result_row_0(result, result + numCol);
	device_vector<float> result_row_1(result + numCol,result + numCol);

