.. _program_listing_file_src_define.h:

Program Listing for File define.h
=================================

|exhale_lsh| :ref:`Return to documentation for file <file_src_define.h>` (``src/define.h``)

.. |exhale_lsh| unicode:: U+021B0 .. UPWARDS ARROW WITH TIP LEFTWARDS

.. code-block:: cpp

   #ifndef DEFINE_H_
   #define DEFINE_H_
   
   /*---------------
   
     Copyright notices?
     Require gsl >= 2.0
     
     ------------------*/
   
   #include <iostream>
   #include <fstream>
   #include <gsl/gsl_odeiv2.h>
   #include <gsl/gsl_errno.h>  // gsl_success
   #include <gsl/gsl_roots.h>  // root solver
   #include <gsl/gsl_multiroots.h> // multidimensional root solver
   #include <gsl/gsl_sf_debye.h>   // Debye function
   #include <stdexcept>        // runtime_error
   #include <sys/time.h>       // gettimeofday
   #include <cmath>
   #include <vector>
   #include <limits>       // quiet_NaN
   #include <sstream>
   #include <iomanip>
   #include <numeric>
   #include <random>
   
   
   using std::cout;
   using std::cin;
   using std::endl;
   using std::string;
   using std::stringstream;
   using std::ifstream;
   using std::ofstream;
   using std::streampos;
   using std::max;
   using std::min;
   using std::vector;
   using std::abs; //In C++, it's always sufficient to use std::abs; it's overloaded for all the numerical types. 
   using std::runtime_error;
   using std::numeric_limits;
   using std::accumulate;
   using std::default_random_engine;
   using std::normal_distribution;
   using std::uniform_real_distribution;
   using std::lognormal_distribution;
   using std::random_device;
   
   const double pi=3.14159265359;
   const double G=6.67384E-8;
   const double c=2.997924589E10;
   const double kb=1.3806488E-16;
   const double R = 8.31446E7; // gas constant in erg K-1 mol-1
   const double e=4.80320451E-10;          // electron charge in e.s.u unit
   const double me=9.10938291E-28;
   const double mp=1.672621777E-24;
   const double ME=5.9722E27;  // Earth mass defined by IAU
   const double RE=6.3781E8;   // Earth equatorial radius
   const double mFe=55.845;    // iron molar mass in grams
   const double mSi=28.0855;
   const double mMg=24.305;
   const double mO=15.999;
   const double mH=1.00794;
   const double mC=12.011;
   const double mCa=40.078;    // calcium molar mass in grams
   const double mAl=26.9815;   // aluminum molar mass in grams
   const double mNa=22.9898;   // sodium molar mass in grams
   const double NA=6.02214E23;
   const double AAtocm3 = 0.3011;  // 1 \AA^3 equivalents to NA / (2*10^24) = 0.3011 cm^3 / mol.
   const double eV_erg = 1.60217710E-12;   // 1 eV=1.6E-12 erg
   
   template<class T>
   inline double sq(const T &a)
   {return double(a*a);}
   
   struct double_params
   {
     vector<double> x;
   };
   
   #endif // DEFINE_H_