                 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;


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