I'm trying to run the following code:
#include <iostream>
#include "mkl.h"
int main () {
/* Define vector of obseravations (5 2D observation points) */
float pObservations [] = {1., 2., 3., 4.2, 5, 9., 10., 8., 7., 6., 5., 9.};
/* Creates and initializes a new summary statistics task descriptor */
VSLSSTaskPtr task;
const int p = 2;
const int n = 5;
const int xstorage = VSL_SS_MATRIX_STORAGE_ROWS;
int status = 0;
status = vslsSSNewTask (&task, &p, &n, &xstorage, pObservations, NULL, NULL);
if (status != VSL_STATUS_OK) {
std::cout << "Failed to create a new summary statistics task descriptor"<< std::endl;
throw false;
}
/* Modifies array pointers related to multivariate mean calculation */
float* pMean = new float [p];
status = vslsSSEditTask (task, VSL_SS_ED_MEAN, pMean);
if (status != VSL_STATUS_OK) {
std::cout << "Failed to modifies array pointers related to multivariate mean calculation"<< std::endl;
throw false;
}
/* Computes Summary Statistics estimates - mean calculation */
status = vslsSSCompute(task, VSL_SS_MEAN, VSL_SS_METHOD_FAST);
if (status != VSL_STATUS_OK) {
std::cout << "Failed to compute summary statistics estimates with error code "<< status << std::endl;
throw false;
}
// Print mean values
for (int ip = 0; ip < p; ip++)
std::cout << pMean [ip] << std::endl;
/* Modifies array pointers related to multivariate outliers detection */
const int nParams = 0;
float* pWeights = new float [n];
status = vslsSSEditOutliersDetection (task, &nParams, NULL, pWeights);
if (status != VSL_STATUS_OK) {
std::cout << "Failed to modifies array pointers related to multivariate outliers detection"<< std::endl;
throw false;
}
/* Computes Summary Statistics estimates - outlier detection */
status = vslsSSCompute(task, VSL_SS_OUTLIERS, VSL_SS_METHOD_BACON);
if (status != VSL_STATUS_OK) {
std::cout << "Failed to compute summary statistics estimates with error code "<< status << std::endl;
throw false;
}
return (0);
}
for a 2D data set consisting of 5 pairs of observation. The output of the program reads:
3.04
8
Failed to compute summary statistics estimates with error code -4002
terminate called after throwing an instance of 'bool'
Abort (core dumped)
The first two numbers are the mean values of the observations in each dimension (2), and the result is accurate.
I'm using the same dataset for the outlier detection BACON algorithm, but can an error -4002, which means that the number of input observation (5 in my case) is either 0 or negative.
Is this a bug in MKL, or something wrong on my side.
Thanks,
Yaniv
