 	std::string chifile("dusthm_chi_summary.dat");
 	write_chis(chifile, chi_samples);
 
-}
-
+	time(&timer2);
+	double seconds = difftime(timer2, timer1);
+
+	std::cout << "Program took " << seconds << " seconds." << std::endl;
+
+}
+

	}
	return chisqr;
}

//__device__ __host__
//double tanh(double x) {
//	return (exp(2.0 * x) - 1.0) / (exp(2.0 * x) + 1.0);
//}

/*
 * This function returns the logarithm of the conditional density of the characteristics given the


int main(int argc, char** argv)
{
	time_t timer1, timer2;  // keep track of how long the program takes to run
	time(&timer1);
	/*
	 * Read in the data for the measurements, meas, and their standard deviations, meas_unc.
	 */
	std::string datafile = "../data/cbt_sed_1000.dat";
	int ndata = get_file_lines(datafile) - 1;  // subtract off one line for the header
	std::cout << "Loaded " << ndata << " data points." << std::endl;
