double
Gaussian::Cdf(size_t dim_, double val_) const{
double nDevs = (val_ - GetMean(dim_))/GetStDev(dim_);
if (nDevs > fCdfCutOff)
return 1;
if(nDevs < -1 * fCdfCutOff)
return 0;
return gsl_cdf_gaussian_P(val_ - GetMean(dim_), GetStDev(dim_));
}
int main(int argc, const char * argv[]) {
float Array[ARRAYSIZE] = {0};
for (int i = 0; i<ARRAYSIZE; i++) {
Array[i] = i;
}
double result ;
gettimeofday(&start, NULL);
result = GetStDev(Array, ARRAYSIZE);
gettimeofday(&stop, NULL);
result_time = timedifference_msec(stop, start);
printf("C:STDEV RESULT : %f, code executed in %f microsecond \n", result,result_time);
gettimeofday(&start, NULL);
result = GetStDevIntrinsic(Array, ARRAYSIZE);
gettimeofday(&stop, NULL);
result_time = timedifference_msec(stop, start);
printf("C_intric:STDEV RESULT : %f, code executed in %f microsecond \n", result,result_time);
return 0;
}
