28char bin_bhns_rotaxis_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Bin_bhns/bin_bhns_rotaxis.C,v 1.3 2014/10/13 08:52:41 j_novak Exp $" ;
67 cout <<
"Bin_bhns::rotation_axis:" << endl ;
68 cout <<
" !!! WARNING : error_y is larger than +1 !!!" << endl ;
75 double gg =
pow( (2.-error_y)/(2.-2.*error_y), rot_exp_x) ;
79 cout <<
"Bin_bhns::rotation_axis:" << endl ;
80 cout <<
" error_y : " << error_y <<
" gg : " << gg << endl ;
85 double xbh_new = xbh_old * gg ;
86 double xns_new = xbh_new +
separ ;
88 cout <<
" new X_BH : " << xbh_new / km <<
" [km]"
89 <<
" new X_NS : " << xns_new / km <<
" [km]"
95 (
star.
set_mp()).set_ori(xns_new, yns_old, 0.) ;
110 if (error_x <= -1.) {
111 cout <<
"Bin_bhns::rotation_axis:" << endl ;
112 cout <<
" !!! WARNING : error_x is smaller than -1 !!!" << endl ;
119 double ff =
pow( (2.+error_x)/(2.+2.*error_x), rot_exp_y) ;
123 if ( fabs(error_x) < thres_rot ) {
124 cout <<
"Bin_bhns::rotation_axis:" << endl ;
125 cout <<
" ff is set to 1 because error_x is smaller than" << endl ;
126 cout <<
" the threshold value (" << thres_rot <<
")" << endl ;
131 cout <<
"Local center of mass of NS:" << endl ;
135 cout <<
"Bin_bhns::rotation_axis:" << endl ;
136 cout <<
" error_x : " << error_x <<
" ff : " << ff << endl ;
141 double aa = fact *
separ ;
142 double yns_new = yns_old + aa * (1. - ff) ;
144 cout <<
" new Y_BH : " << 0. / km <<
" [km]"
145 <<
" new Y_NS : " << yns_new / km <<
" [km]"
152 (
star.
set_mp()).set_ori(xns_old, yns_new, 0.) ;
Hole_bhns hole
Black hole.
double ya_barycenter() const
Absolute coordinate Y of the barycenter of the baryon density.
double omega
Angular velocity with respect to an asymptotically inertial observer.
double & set_y_rot()
Sets the absolute coordinate Y of the rotation axis [{\tt r_unit}].
const Tbl & line_mom_bhns() const
Total linear momentum.
double xa_barycenter() const
Absolute coordinate X of the barycenter of the baryon density.
double separ
Absolute orbital separation between two centers of BH and NS.
void rotation_axis_y(double thres_rot, double rot_exp_y, double fact)
Computes the position of the rotation axis Y.
double & set_x_rot()
Sets the absolute coordinate X of the rotation axis [{\tt r_unit}].
Star_bhns star
Neutron star.
void rotation_axis_x(double rot_exp_x)
Computes the position of the rotation axis X.
Map & set_mp()
Read/write of the mapping.
const Map & get_mp() const
Returns the mapping.
double get_mass_bh() const
Returns the gravitational mass of BH [{\tt m_unit}].
double get_ori_y() const
Returns the y coordinate of the origin.
double get_ori_x() const
Returns the x coordinate of the origin.
const Map & get_mp() const
Returns the mapping.
Map & set_mp()
Read/write of the mapping.
Cmp pow(const Cmp &, int)
Power .
Standard units of space, time and mass.