/* --------------- */
/* --- mutil.h --- */
/* --------------- */

#ifndef __MUTIL_H__
#define __MUTIL_H__

#ifdef __cplusplus
#pragma message ("C++")
extern "C" {
#endif

#define RUN 5
#define FREQ 2.4e9

//#define CHRONO(X,t) tmin = 1e10; for(r=0; r<run; r++) { t0 = chrono(); for(i=0;i<iter;i++) X;  t1 = chrono(); dt = t1-t0; if(dt<tmin) tmin = dt; } t = tmin*FREQ/(iter*n*n)
    #define CHRONO(X,t) tmin = 1e10; for(r=0; r<run; r++) { t0 = chrono(); for(i=0;i<iter;i++) X;  t1 = chrono(); dt = t1-t0; if(dt<tmin) tmin = dt; } t = tmin/iter
//#define CHRONO(X,t) tmin = 1e10; for(r=0; r<run; r++) { t0 = chrono(); for(i=0;i<iter;i++) X;  t1 = chrono(); dt = t1-t0; if(dt<tmin) tmin = dt; } t=tmin/iter

float32 forceCalc(float32 **X, int i0, int i1, int j0, int j1);

double cpp(double t, int n);

double cp0DotProduct  (double t, int n);
double cp0Euclide     (double t, int n);
double cp0Manhattan   (double t, int n);
double cp0Mahalanobis (double t, int n);

double mflopsDotProduct (double t, int n);
double mflopsEuclide    (double t, int n);
double mflopsManhattan  (double t, int n);
double mflopsMahalanobis(double t, int n);
double mflopsMul        (double t, int n);

uint8 ui8rand(void);
uint32 ui32rand(void);
uint32 f32rand(void);

void rand_ui8vector(uint8 *X, int i0, int i1);
void rand_ui32vector(uint32 *X, int i0, int i1);
void rand_f32vector(float32 *X, int i0, int i1);

int getIter(int size);
int getIterAV(int size);

float32 gauss(float32 sigma, float32 x);


#ifdef __cplusplus
}
#endif

#endif /* __MUTIL_H__ */