test_clustering_kmeans.cpp

 
#include <iostream>
#include <Pds/Ra>
#include <Pds/Ml>
    
int main(void)
{
    unsigned int L=5000000;//5000000;
    unsigned int K=32;
    Pds::IterationConf Conf;
    
    Pds::Matrix X;
    Pds::Matrix C;
    double Cost;
    Pds::Array<unsigned int> ID;
 
    Pds::Ra::Tic();
    Pds::DataSet::Clustering::Kcumulus(L/K,X,K);
    Pds::Ra::Toc();
 
    std::cout<<std::endl;
    
    Pds::Ra::Tic();
    Conf.SetMinError(0.001);
    C=Pds::Clustering::Kmeans(Conf,X,K,ID);
    Pds::Ra::Toc();
    Conf.Print("Kmeans:\n");
    
    C.Round(3).Print("C:\n");
    
    std::cout<<std::endl;
 
    Pds::Ra::Tic();
    Cost=Pds::Clustering::CostKmeans(X,C,ID);
    Pds::Ra::Toc();
    std::cout<<"Cost:\t"<<Cost<<std::endl;
    
    Pds::Ra::Tic();
    Cost=Pds::Clustering::CostKmeans(X,C);
    Pds::Ra::Toc();
    std::cout<<"Cost:\t"<<Cost<<std::endl;
    
    return 0;
}