 	/* Compute sizes of various data structures */
 	size_t mean_size = N*sizeof(float);
 	size_t cov_size = (N*(N+1)/2)*sizeof(float);
+	size_t full_cov_size = N*N*sizof(float);
 	size_t samples_size = N*nsamps_tot*sizeof(float);
 	
 	/* allocate result memory */

		printf("%f\n", sample_mean_res[i]);
	}
	write_paths_to_csv("maxent_samples.csv", samples_res, nsamps_tot, N, "w");
	write_paths_to_csv("maxent_cov.csv", sample_covariance_res, N, N, "w");


	/* free memory */

__global__ void compute_sample_covariance(float* samples, float* sample_covariance, int nsamps, int N)
{
	int i = blockIdx.x;
	
	sample_covariance[i] = 0.0;
	for(int samp = 0; samp<nsamps; samp++)
	{
		sample_covariance[i] += samples[samp*N+i]*samples[samp*N+j];
	}
	sample_covariance[i] = sample_covariance[i]/nsamps;

}

