                 scaled_proposal[j] += h_theta_cholfact[cholfact_index] * snorm_deviate[k];
                 cholfact_index++;
             }
-            scaled_proposal[j] = 0.0;
             proposed_theta[j] = h_theta[j] + scaled_proposal[j];
         }
 

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

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


     collection of the mean and variances of the normal distributions.
    */

    int n = 1000; // # of data points

    double mu_norm = 8.5 * log(10.0);  // Average value of the natural logarithm of the SED normalization
    double sig_norm = 0.5 * log(10.0);  // Standard deviation in the SED normalization
