 __constant__ int c_mfeat; // # of measured features per data point
 __constant__ int c_pchi; // # of characteristics per data point
 __constant__ int c_dim_theta; // # dimension of population parameters
-__constant__ double c_target_rate; // MCMC sampler target acceptance rate
-__constant__ double c_decay_rate; // decay rate of robust adaptive metropolis algorithm
 
 // Function to compute the rank-1 Cholesky update/downdate. Note that this is done in place.
 __device__ __host__

		exit(1);															\
	} }

typedef std::vector<std::vector<double> > vecvec;
typedef hvector hvector;
typedef dvector dvector;


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