 		d_meas.resize(ndata * mfeat);
 		h_meas_unc.resize(ndata * mfeat);
 		d_meas_unc.resize(ndata * mfeat);
-		h_logdens_meas.resize(ndata);
-		d_logdens_meas.resize(ndata);
-		h_logdens_pop.resize(ndata);
-		d_logdens_pop.resize(ndata);
+		h_logdens.resize(ndata);
+		d_logdens.resize(ndata);
 		h_chi.resize(ndata * pchi);
 		d_chi.resize(ndata * pchi);
 		int dim_cholfact = pchi * pchi - ((pchi - 1) * pchi) / 2;

		// get current conditional log-posterior of population
		dvector d_logdensity = Daug.GetDevLogDensPop();
		double logdens_current = thrust::reduce(d_logdensity.begin(), d_logdensity.end());

		// propose new value of population parameter
		hvector h_proposed_theta = Propose();

        CUDA_CHECK_RETURN(cudaDeviceSynchronize());
        // transfer values back to host
        h_chi = d_chi;
        h_logdens = d_logdens;
        h_naccept = d_naccept;

        current_iter++;
