@@ -7,17 +7,11 @@ template<typename scalar_t>
77int nanoflann_neighbors (vector<scalar_t >& queries, vector<scalar_t >& supports,
88 vector<long >& neighbors_indices, float radius, int dim, int max_num){
99
10- // Initiate variables
11- // ******************
12-
1310 const scalar_t search_radius = static_cast <scalar_t >(radius*radius);
1411
1512 // Counting vector
1613 size_t max_count = 1 ;
1714
18- // Nanoflann related variables
19- // ***************************
20-
2115 // CLoud variable
2216 PointCloud<scalar_t > pcd;
2317 pcd.set (supports, dim);
@@ -54,32 +48,15 @@ int nanoflann_neighbors(vector<scalar_t>& queries, vector<scalar_t>& supports,
5448 scalar_t * query_pt = new scalar_t [dim];
5549 std::copy (p0.begin (), p0.end (), query_pt);
5650
57- // for(int i=0; i < p0.size(); i++)
58- // std::cout << query_pt[i] << '\n';
59-
6051 list_matches[i0].reserve (max_count);
6152 std::vector<std::pair<size_t , scalar_t > > ret_matches;
6253
6354 const size_t nMatches = index->radiusSearch (query_pt, search_radius+eps, ret_matches, search_params);
6455
65- // cout << "radiusSearch(): radius=" << search_radius << " -> " << nMatches << " matches\n";
66- // for (size_t i = 0; i < nMatches; i++)
67- // cout << "idx["<< i << "]=" << ret_matches[i].first << " dist["<< i << "]=" << ret_matches[i].second << endl;
68- // cout << "\n";
69-
7056 list_matches[i0] = ret_matches;
7157 if (max_count < nMatches) max_count = nMatches;
7258 i0++;
7359
74-
75- // Get worst (furthest) point, without sorting:
76-
77- // cout << "\n neighbors: " << nMatches << "\n";
78-
79- // Get worst (furthest) point, without sorting:
80- // for(int i=0; i < ret_matches.size(); i++)
81- // std::cout << ret_matches.at(i) << '\n';
82-
8360 }
8461 // Reserve the memory
8562 if (max_num > 0 ) {
0 commit comments