00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030 char eos_fps_C[] = "$Header: /cvsroot/Lorene/C++/Source/Eos/eos_fps.C,v 1.2 2002/10/16 14:36:35 j_novak Exp $" ;
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051 #include "headcpp.h"
00052 #include "eos.h"
00053
00054
00055
00056
00057
00058
00059
00060 Eos_FPS::Eos_FPS(const char* path)
00061 : Eos_tabul("EOS FPS",
00062 "eos_fps.d", path)
00063 {}
00064
00065
00066
00067
00068 Eos_FPS::Eos_FPS(FILE* fich) : Eos_tabul(fich) {}
00069
00070
00071
00072
00073
00074 Eos_FPS::Eos_FPS(ifstream& fich) : Eos_tabul(fich, "eos_fps.d") {}
00075
00076
00077
00078
00079
00080
00081
00082 Eos_FPS::~Eos_FPS(){
00083
00084
00085
00086 }
00087
00088
00089
00090
00091
00092
00093
00094 bool Eos_FPS::operator==(const Eos& eos_i) const {
00095
00096 bool resu = true ;
00097
00098 if ( eos_i.identify() != identify() ) {
00099 cout << "The second EOS is not of type Eos_FPS !" << endl ;
00100 resu = false ;
00101 }
00102
00103 return resu ;
00104
00105 }
00106
00107 bool Eos_FPS::operator!=(const Eos& eos_i) const {
00108
00109 return !(operator==(eos_i)) ;
00110
00111 }
00112
00113
00114
00115
00116
00117
00118 ostream& Eos_FPS::operator>>(ostream & ost) const {
00119
00120 ost <<
00121 "EOS of class Eos_FPS : Friedman-Pandharipande + Skyrme "
00122 << endl ;
00123 ost << " (Lorenz, Ravenhall, Petick)" << endl ;
00124
00125 ost << " composition : n,p,e,mu" << endl ;
00126 ost << " model : effective nucleon energy functional, FPS" << endl ;
00127 ost << " BPS EOS below neutron drip point" << endl ;
00128 ost << " Crust bottom at n = 0.0957 fm^{-3}" << endl ;
00129
00130 return ost ;
00131
00132 }
00133
00134