     */
     
     int n = 2000; // # of data points
-	int m = 5; // # of features per data point (i.e., # of points in the SED for the i^th source)
     
     double mu_norm = 8.5 * log(10);  // Average value of the natural logarithm of the SED normalization
     double sig_norm = 0.5 * log(10);  // Standard deviation in the SED normalization


// Perform the MHA update on the n characteristics, done in parallel on the GPU.
__global__
void update_chi(double* chi, double* meas, int m, double* meas_unc, int n, double* logdens, double* prop_cholfact,
                curandState* state, int current_iter)
{
	int i = blockDim.x * blockIdx.x + threadIdx.x;

static __constant__ double c_theta[dim_theta];

// Constants used for the model SED
static __constant__ double c_nu0 = 2.3e11;  // reference frequency, nu0 = 230 GHz
static const double nu0 = 2.3e11;
static const int m = 5;  // the number of measured features per data point
static __constant__ int c_m = 5;
    // observational frequencies, corresponding to the Herschel PACS and SPIRES bands
static const double nu[m] = {4.286e12, 1.765e12, 1.200e12, 8.571e11, 6.000e11};
static __constant__ c_nu[m] = {4.286e12, 1.765e12, 1.200e12, 8.571e11, 6.000e11};
static const double hplanck = 6.6260755e-27;  // Planck's constant, in cgs
static const double clight = 2.997925e10;
static const double kboltz = 1.380658e-16;
