 		cudaMalloc((void**)&dev_facets_, sizeof(index_t) * facet_nb_ * dimension_);
 		cudaMalloc((void**)&dev_points_nn_, sizeof(index_t) * points_nb_ * k_);
 		cudaMalloc((void**)&dev_facets_nn_, sizeof(index_t) * facet_nb_ * 1);
-
+		cudaMalloc((void**)&dev_ret_, sizeof(double) * points_nb_ * 4);
 		//Output result.
 		//cudaMalloc((void**)&dev_ret_, sizeof(double) *  facet_nb_ * 10 * 40);
 		//Copy


		
		for (index_t t = 0; t < iter_nb_; ++t){
			cudaMalloc((void**)&dev_ret_, sizeof(double) * points_nb_ * 4);
			
			knn_search(); 
				CudaStopWatcher iter_watcher("iteration");
				iter_watcher.start();				
				cudaMemcpy(dev_points_, points_, DOUBLE_SIZE * points_nb_ * dimension_, cudaMemcpyHostToDevice);
				cudaMemcpy(dev_points_nn_, points_nn_, sizeof(index_t) * points_nb_ * k_, cudaMemcpyHostToDevice);
				cudaMemcpy(dev_facets_nn_, facets_nn_, sizeof(index_t) * facet_nb_ * 1, cudaMemcpyHostToDevice);

			*/
		//}
		//__syncthreads();
	}

	void CudaRestrictedVoronoiDiagram::knn_search(){
