     }
 
     /*
-    Set up parallel random number generations on the GPU
+    Set up parallel random number generators on the GPU
     */
     curandState* devStates;  // Create state object for random number generators on the GPU
     // Allocate memory on GPU for RNG states

    //std::ofstream chifile("chis.dat");
    std::ofstream thetafile("thetas.dat");

    int mcmc_iter = 10000;
    int naccept_theta = 0;
    std::cout << "Running MCMC Sampler...." << std::endl;


                scaled_proposal_j += cholfactor[cholfact_index] * snorm_deviate[k];
                cholfact_index++;
            }
            scaled_proposal_j = 0.0;
            new_chi[j] = chi[n * j + i] + scaled_proposal_j;
            scaled_proposal[j] = scaled_proposal_j;
        }
