unites.h
/*
* Lorene's units
*
*/
/*
* Copyright (c) 1999-2001 Eric Gourgoulhon
*
* This file is part of LORENE.
*
* LORENE is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* LORENE is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with LORENE; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*/
/*
* $Id: unites.h,v 1.8 2015/03/17 14:20:00 j_novak Exp $
* $Log: unites.h,v $
* Revision 1.8 2015/03/17 14:20:00 j_novak
* New class Hot_eos to deal with temperature-dependent EOSs.
*
* Revision 1.7 2014/10/13 08:52:37 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.6 2014/06/30 14:34:57 j_novak
* Update of the values of some constants (G, M_sol and MeV).
*
* Revision 1.5 2014/05/13 10:06:12 j_novak
* Change of magnetic units, to make the Lorene unit system coherent. Magnetic field is now expressed in Lorene units. Improvement on the comments on units.
*
* Revision 1.4 2004/12/01 12:28:32 p_grandclement
* Include math.h in unite.h
*
* Revision 1.3 2004/03/25 10:28:56 j_novak
* All LORENE's units are now defined in the namespace Unites (in file unites.h).
*
* Revision 1.2 2004/03/22 13:12:43 j_novak
* Modification of comments to use doxygen instead of doc++
*
* Revision 1.1.1.1 2001/11/20 15:19:27 e_gourgoulhon
* LORENE
*
* Revision 1.2 2000/03/17 15:28:29 eric
* Ajout de ggrav (G).
*
* Revision 1.1 1999/12/06 13:35:10 eric
* Initial revision
*
*
* $Header: /cvsroot/Lorene/C++/Include/unites.h,v 1.8 2015/03/17 14:20:00 j_novak Exp $
*
*/
#include <cmath>
//--------------------------//
// Standard LORENE units //
//--------------------------//
namespace Lorene {
namespace Unites {
const double g_si = 6.6738E-11 ;
const double c_si = 2.99792458E+8 ;
const double kB_si = 1.3806488E-23 ;
const double rhonuc_si = 1.66E+17 ;
const double km_si = 1.E+3 ;
const double msol_si = 1.9885E+30 ;
const double mev_si = 1.602176565E-13 ;
const double r_unit = 1.e4 ;
const double v_unit = c_si ;
const double rho_unit = rhonuc_si ;
const double t_unit = r_unit/v_unit ;
const double m_unit = rho_unit * pow(r_unit, 3.) ;
const double g_unit = 1./(rho_unit*t_unit*t_unit) ;
const double f_unit = 1./t_unit ;
const double ggrav = g_si / g_unit ;
const double qpig = 4 * M_PI * ggrav ;
const double msol = msol_si/m_unit ;
const double km = km_si/r_unit ;
const double mevpfm3 = mev_si/( rho_unit * v_unit *v_unit) *1.e45 ;
const double m_u_mev = rho_unit / 1.e44 *c_si*c_si / mev_si ;
}
//----------------------------------//
// Electro-magnetic LORENE units //
//----------------------------------//
namespace Unites_mag {
using namespace Unites ;
const double mu_si = 1.2566370614359173e-6 ;
const double j_unit = 1e11 ;
const double mag_unit = rho_unit*v_unit*v_unit/ (r_unit * j_unit) ;
const double elec_unit = mag_unit * v_unit ;
const double mu0_unit = rho_unit*v_unit*v_unit / (pow(j_unit ,2)*pow(r_unit,2));
const double mu0 = mu_si / mu0_unit ;
}
}
Lorene::pow
Cmp pow(const Cmp &, int)
Power .
Definition cmp_math.C:348
Lorene
Lorene prototypes.
Definition app_hor.h:64
Unites_mag
Standard electro-magnetic units.
Unites
Standard units of space, time and mass.