#include "daal.h"
#include "service.h"
#include <iostream>
using namespace std;
using namespace daal;
using namespace daal::algorithms;
using namespace daal::algorithms::pca::quality_metric;
using namespace daal::algorithms::pca::quality_metric_set;
const string dataFileName = "../data/batch/pca_normalized.csv";
const size_t nVectors = 1000;
const size_t nComponents = 5;
int main(int argc, char *argv[])
{
checkArguments(argc, argv, 1, &dataFileName);
FileDataSource<CSVFeatureManager> dataSource(dataFileName, DataSource::doAllocateNumericTable,
DataSource::doDictionaryFromContext);
dataSource.loadDataBlock(nVectors);
pca::Batch<float, pca::svdDense> algorithm;
algorithm.input.set(pca::data, dataSource.getNumericTable());
algorithm.compute();
pca::quality_metric_set::Batch qms(nComponents);
services::SharedPtr<algorithms::Input> algInput =
qms.getInputDataCollection()->getInput(explainedVariancesMetrics);
explained_variance::InputPtr varianceMetrics = explained_variance::Input::cast(algInput);
varianceMetrics->set(explained_variance::eigenvalues, algorithm.getResult()->get(pca::eigenvalues));
qms.compute();
explained_variance::ResultPtr qmsResult = explained_variance::Result::cast
(qms.getResultCollection()->getResult(explainedVariancesMetrics));
printNumericTable(qmsResult->get(explained_variance::explainedVariances),
"Explained variances:");
printNumericTable(qmsResult->get(explained_variance::explainedVariancesRatios),
"Explained variance ratios:");
printNumericTable(qmsResult->get(explained_variance::noiseVariance),
"Noise variance:");
return 0;
}