 	}
 	__syncthreads();
 	
-	for (int i = 0; i < seeds_nb * 4; ++i)
+	/*for (int i = 0; i < seeds_nb * 4; ++i)
 	{
 		ret_seeds[i] = SeedsInformation[i];
-	}
+	}*/
 	return;
 }
 

		
		intersection_clip_facet_with_knn(polygon, vertex_nb, current_seed, seeds_pointer, seeds_nb, seeds_neighbor_index, k, tid, i);

		/*if (tid == 0 & i == 2)
		{
			for (int m = 0; m < seeds_nb * 4; ++m)
				ret_seeds[m] = m;
			ret_seeds[0] = v1.x;
			ret_seeds[1] = v1.y;
			ret_seeds[2] = v1.z;

			for (int m = 0; m < vertex_nb * 4; ++m)
			{
				ret_seeds[4 + m] = polygon[m];
			}
			ret_seeds[21] = vertex_nb;
		}
		if (i == 5)
			return;*/
		//now we get the clipped polygon stored in "polygon"
		//take care of the sychonize
		//change the polygon data into "weight" and "position"

void compute_RVD_with_knn(double* seeds_pointer, int seeds_nb,
double* mesh_vertex, int mesh_vertex_nb,
int* mesh_facet, int mesh_facet_nb,
int* facet_center_neighbor_index, int* seeds_neighbor_index, int k, double* ret_seeds, double* test_seeds)
{
	int tid = blockIdx.x * blockDim.x + threadIdx.x;
	if (tid >= mesh_facet_nb) return;
