28char HiggsMonopole_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Compobj/higgs_monopole.C,v 1.3 2014/10/13 08:52:49 j_novak Exp $" ;
69 ifstream file(file_name) ;
71 cerr <<
"Problem in opening the file " << file_name << endl ;
82 file.ignore(1000,
'\n') ;
87 cout <<
"nz : " << nz << endl ;
92 double explamb_last, nn_last, hh_last, press_last ;
93 for (
int l=0; l<nz; l++) {
94 cout <<
"l = " << l << endl ;
99 double* explamb_tab =
new double[nr] ;
100 double* nn_tab =
new double[nr] ;
101 double* hh_tab =
new double[nr] ;
102 double* press_tab =
new double[nr] ;
106 explamb_tab[0] = explamb_last ;
107 nn_tab[0] = nn_last ;
108 hh_tab[0] = hh_last ;
109 press_tab[0] = press_last ;
111 for (
int i=i_min; i<nr; i++) {
113 file >> explamb_tab[i] ;
116 file >> press_tab[i] ;
118 explamb_last = explamb_tab[nr-1] ;
119 nn_last = nn_tab[nr-1] ;
120 hh_last = hh_tab[nr-1] ;
121 press_last = press_tab[nr-1] ;
123 for (
int i=0; i<nr; i++) {
124 cout <<
" explamb, nn, hh : " << explamb_tab[i] <<
" "
125 << nn_tab[i] <<
" " << hh_tab[i] << endl ;
128 for (
int k=0; k<np; k++) {
129 for (
int j=0; j<nt; j++) {
130 for (
int i=0; i<nr; i++) {
139 delete[] explamb_tab ;
185 ost << endl <<
"Higgs monopole" << endl ;
Base class for stationary compact objects (under development).
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
Scalar nn
Lapse function N .
Map & mp
Mapping describing the coordinate system (r,theta,phi)
Higgs monopole (under development).
HiggsMonopole(Map &mp_i, const char *file_name)
Standard constructor.
Scalar hh
Higgs scalar field.
virtual ~HiggsMonopole()
Destructor.
Scalar press
Fluid pressure.
char description2[256]
String describing the model.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
char description1[256]
String describing the model.
Scalar grr
Metric coefficient g_rr.
Base class for coordinate mappings.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
virtual void allocate_all()
Sets the logical state to ETATQCQ (ordinary state) and performs the memory allocation of all the elem...
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
Standard units of space, time and mass.