 	int cell = blockIdx.x;
 	for(int j=0; j<cell; j++)
 	{
-		sample_covariance[cell*N]
+		sample_covariance[cell*N] = 0.0;
 	}
 }
 

	curandState *d_state;
	cudaMalloc(&d_state, nblocks);





}

	int N = 20; /* Number of Neurons */
	int samps_per_block = 1024;
	int nblocks = 1024
	int nsamps_tot = samps_per_block*nblocks
	int nsweeps = 10;

	/* 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 samples_size = N*nsamps_tot*sizeof(float);

	/* allocate result memory */
	float *alpha_res = (float*)malloc(mean_size);
	float *beta_res = (float*)malloc(cov_size);
	float *sample_mean_res = (float*)malloc(mean_size);

	/* Allocate device memory */
	float *sample_mean;
