LORENE: eos_glendnh3.C Source File

00001 /*
00002  *  Methods of class Eos_GlendNH3
00003  *
00004  *  (see file eos_tabul.h for documentation).
00005  *
00006  */
00007 
00008 // Headers Lorene
00009 #include "headcpp.h"
00010 #include "eos.h"
00011 
00012             //----------------------------//
00013             //      Constructors          //
00014             //----------------------------//
00015 
00016 // Standard constructor
00017 // --------------------         
00018 Eos_GlendNH3::Eos_GlendNH3(const char* path)
00019         : Eos_tabul("EOS GlendNH3",
00020                     "eos_glendnh3.d", path)
00021 {}
00022 
00023 
00024 // Constructor from binary file
00025 // ----------------------------
00026 Eos_GlendNH3::Eos_GlendNH3(FILE* fich) : Eos_tabul(fich) {}
00027 
00028 
00029 
00030 // Constructor from a formatted file
00031 // ---------------------------------
00032 Eos_GlendNH3::Eos_GlendNH3(ifstream& fich) :
00033             Eos_tabul(fich, "eos_glendnh3.d") {}
00034 
00035 
00036 
00037             //--------------//
00038             //  Destructor  //
00039             //--------------//
00040 
00041 Eos_GlendNH3::~Eos_GlendNH3(){
00042 
00043     // does nothing
00044 
00045 }
00046 
00047 
00048             //------------------------//
00049             //  Comparison operators  //
00050             //------------------------//
00051 
00052 
00053 bool Eos_GlendNH3::operator==(const Eos& eos_i) const {
00054 
00055     bool resu = true ;
00056 
00057     if ( eos_i.identify() != identify() ) {
00058     cout << "The second EOS is not of type Eos_GlendNH3 !" << endl ;
00059     resu = false ;
00060     }
00061 
00062     return resu ;
00063 
00064 }
00065 
00066 bool Eos_GlendNH3::operator!=(const Eos& eos_i) const {
00067 
00068     return !(operator==(eos_i)) ;
00069 
00070 }
00071 
00072             //------------//
00073             //  Outputs   //
00074             //------------//
00075 
00076 
00077 ostream& Eos_GlendNH3::operator>>(ostream & ost) const {
00078 
00079     ost <<
00080     "EOS of class Eos_GlendNH3 : Glendenning, N.K. 1985 ApJ293, 470, case 3 "
00081         << endl ;
00082         
00083     ost << "  composition :  n,p,H" << endl ;
00084     ost << "  model : Lagrangian field theory; MFT"
00085         << endl ;
00086     ost << "  BPS EOS below neutron drip point" << endl ;
00087     ost << "  Sly4 EOS up to the liquid core" << endl ;
00088     ost << "  Crust bottom at n = 0.1 fm^{-3}" << endl ;
00089 
00090     return ost ;
00091 
00092 }