 				);
 
 			//now we get the clipped polygon stored in "polygon", do something.
-			/*action(
+			action(
 				current_polygon,
 				current_seed,
 				retdata
-				);*/
-			MyAtomicAdd(&retdata[current_seed], 1);
-			MyAtomicAdd(&retdata[2000], 1);
+				);
+			//MyAtomicAdd(&retdata[current_seed], 1);
+			//MyAtomicAdd(&retdata[2000], 1);
 			//store_info(tid, current_seed, current_polygon, &retdata[tid * 400 + counter * 40]);
 			//Propagate to adjacent seeds
 			for (index_t v = 0; v < current_polygon.vertex_nb; ++v)

							has_visited_flag = true;
					}
					//the neighbor seed is new!
					if (!has_visited_flag)
					{
						to_visit[idx][to_visit_pos++] = ns;
						has_visited[has_visited_nb++] = ns;

	}

	void CudaRestrictedVoronoiDiagram::update_points(){
		//x_->clear();
		
		for (int i = 0; i < points_nb_; ++i){
			double d = host_ret_[i * 4 + 3];
			if (fabs(host_ret_[i * 4 + 3]) > 1e-30){
				host_ret_[i * 4 + 0] /= host_ret_[i * 4 + 3];
				host_ret_[i * 4 + 1] /= host_ret_[i * 4 + 3];
				host_ret_[i * 4 + 2] /= host_ret_[i * 4 + 3];

				x_->set_vertex(&host_ret_[i * 4], dimension_, i);
			}
			else
			{
				//std::cout << "point  " << i << " hsa no RVD." << std::endl;
				
			}
			//x_->add_vertexd(&host_ret_[i * 4], dimension_);
		}

		std::vector<int> sample_facet(points_nb_);
