index_tnb_neighbors=index_t(num_points);// since there are only seven points right now, this number will probably be more like 50 later

index_tnb_neighbors=50;//index_t(num_points); // since there are only seven points right now, this number will probably be more like 50 later

std::vector<index_t>neighbors(nb_neighbors);

std::vector<index_t>total_neighbors(nb_neighbors*num_points);// i don't know if this is the correct number if, for example, nb points and nb neighbors are not the same (which they usually will not be)