 // 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 const int mfeat = 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__ double c_nu[m] = {4.286e12, 1.765e12, 1.200e12, 8.571e11, 6.000e11};
+static const double nu[mfeat] = {4.286e12, 1.765e12, 1.200e12, 8.571e11, 6.000e11};
+static __constant__ double c_nu[mfeat] = {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;
 static __constant__ double c_2h_over_csqr = 1.4745002e-47;  // 2 * hplanck / clight^2, used in the Planck function
 static __constant__ double c_h_over_k = 4.7992157e-11;  // hplanck / kboltz
+
+// Global random number generator and distributions for generating random numbers on the host. The random number generator used
+// is the Mersenne Twister mt19937 from the BOOST library.
+boost::random::mt19937 rng;
+boost::random::normal_distribution<> snorm(0.0, 1.0); // Standard normal distribution
+boost::random::uniform_real_distribution<> uniform(0.0, 1.0); // Uniform distribution from 0.0 to 1.0
 
 // Function to return the model SED at a input frequency, given the characteristics chi = (norm, beta, temp)
 __device__

        for (int j=0; j<c_p; j++) {
            snorm_deviate[j] = curand_normal_double(&localState);
        }
        
        // propose a new chi value
        double new_chi[p];
        double scaled_proposal[p];

/*
 * dust_seds.cu
 *
 *  Created on: Jun 2, 2013
 *      Author: brandonkelly
 */

// Cuda Includes
#include "cuda_runtime.h"
#include "device_launch_parameters.h"
