     thrust::device_vector<double> centroids(k * d);
     
     
-    kmeans::detail::find_centroids(n, d, k, data, labels, centroids);
+    kmeans::kmeans(iterations, n, d, k, data, labels, centroids);
 
 
 }

#include <thrust/device_vector.h>
#include "centroids.h"
#include <iostream>

template<typename T>

}

int main() {
    int iterations = 50;
    int n = 1e6;
    int d = 50;
    int k = 100;
