 
 	double sigma_msmt[5] = {2.2e-4, 3.3e-4, 5.2e-4, 3.7e-4, 2.2e-4};  // Standard deviation for the measurement errors
 
+	for (int j = 0; j < mfeat; ++j) {
+		sigma_msmt[j] *= 1.0;
+	}
+
 	for (int i=0; i<ndata; i++) {
         // Loop over the data indices
 		for (int j=0; j<p; j++) {

        for (int j=0; j<dim_theta; j++) {
            thetafile << h_theta[j] << " ";
        }
        thetafile << std::endl;
        // Save the chi values
        for (int j=0; j<p; j++) {
        	chifile << h_chi[ndata * j] << " ";
        }
        chifile << std::endl;
    }
    thetafile.close();


    thrust::device_vector<double> d_logdens_pop(ndata); // log densities of theta values for each chi

   	// Load the initial guess for theta into constant memory
    double* p_theta = thrust::raw_pointer_cast(&h_theta[0]);
	CUDA_CALL(cudaMemcpyToSymbol(c_theta, p_theta, h_theta.size() * sizeof(*p_theta)));

    // Cuda grid launch
