int genotypeConfidence(const double *prob);

 int genotypeConfidence2(double probability);

 int intInSet(const int *x, const int *set, const int *n);
 int genotypeCall(const double *pAA, const double *pAB, const double *pBB);
 double sdCorrection(const int *n);
 int sort_double(const double *a1, const double *a2);
 void trimmed_mean(double *datavec, int *classvec, int class, double trim, int cols, int rows, double *m1, double *m2, double *m3, int i_ext);
 void trimmed_stats(double *data, double *m1, double *m2, double *m3, int *class, int rows, int cols, double *trim);

 long CMatrixElementPosition(long row, long col, long totrow);