     // Initialize the random number generator states on the GPU
     setup_kernel<<<nBlocks,nThreads>>>(devStates);
 
-    // Wait until everything is done running on the GPU, make sure everything went OK
+    // Wait until RNG stuff is done running on the GPU, make sure everything went OK
     CUDA_CALL(cudaDeviceSynchronize());
 
-    // Now run the MCMC sampler
+    /*******************************
+    	Now run the MCMC sampler
+	*******************************/
 
     //std::ofstream chifile("chis.dat");
     std::ofstream thetafile("thetas.dat");

	}
}

int main(void)
{
    /* TODO: Should use the boost random number generator libraries for doing this on the host */

        return 2;
    }

    /*
    Set up parallel random number generations on the GPU
    */
    curandState* devStates;  // Create state object for random number generators on the GPU
    // Allocate memory on GPU for RNG states
    CUDA_CALL(cudaMalloc((void **)&devStates, nThreads.x * nBlocks.x *
