 public:
 	// constructor
 	PopulationPar(double rate, DataAugmentation<ChiType>& D, dim3& nB, dim3& nT) :
-		target_rate(rate), daug(D), nBlocks(nB), nThreads(nT)
+		target_rate(rate), Daug(D), nBlocks(nB), nThreads(nT)
 	{
 		h_theta.resize(dim_theta);
 		d_theta = h_theta;

		logdensity_pop<ChiType><<<nBlocks,nThreads>>>(p_theta, p_chi, p_logdens);
        CUDA_CHECK_RETURN(cudaDeviceSynchronize());

		// reset the number of MCMC iterations
		current_iter = 1;
		naccept = 0;


	// compute the conditional log-posterior density of the characteristics given the population parameter
	template<class ChiType> __global__
	void logdensity_pop(double* theta, double* chi, double* logdens)
	{
		int idata = blockDim.x * blockIdx.x + threadIdx.x;
		if (idata < c_ndata)
		{
			ChiType Chi(1, idata);
			logdens[idata] = Chi.logdensity_pop(chi, theta);
		}
	}

	// propose a new value of the population parameters
