#include "KMEANS.h"
//using   namespace   std;
class Gaussian
{
public:
	double *devia;
	double *mean;
	double *a;

	double **samplelist;
	int NumPatterns;
	int SizeVector;
	int NumMixture;

	int *clu_idx;
	int *num_samples_clu;

	CKmeans m_kmeans;
	double result;
	void erhand(char err_msg[]);
	double ** matrix(int n,int m);
	double * vect(int n);
	void Classifier(double *vector);
	void Initialize(char *fname);
	void GaussEM(double **samples, int num_samples, int vec_size, int num_cluster, int *cluster_index,int* num_samples_in_clusters,double *means,double* delta,double* weight,int maxIter);
	void Cal_mean_delta(double *data,aCluster *Clus,double *Uni_mean,double *Uni_devia);

	double ComputeUniModeGauss(double* vect, double* mu, double* inv_var, double log_var_val, int vect_size);
	double ComputeGaussMixture(double* vect, double* mu, double* inv_var, double* log_var_val, int vect_size, double* weight, int num_mix );
};