70 cout <<
"Bin_bhns::rotation_axis:" << endl ;
71 cout <<
" !!! WARNING : error_y is larger than +1 !!!" << endl ;
78 double gg =
pow( (2.-error_y)/(2.-2.*error_y), rot_exp_x) ;
82 cout <<
"Bin_bhns::rotation_axis:" << endl ;
83 cout <<
" error_y : " << error_y <<
" gg : " << gg << endl ;
88 double xbh_new = xbh_old * gg ;
89 double xns_new = xbh_new +
separ ;
91 cout <<
" new X_BH : " << xbh_new / km <<
" [km]" 92 <<
" new X_NS : " << xns_new / km <<
" [km]" 98 (
star.
set_mp()).set_ori(xns_new, yns_old, 0.) ;
113 if (error_x <= -1.) {
114 cout <<
"Bin_bhns::rotation_axis:" << endl ;
115 cout <<
" !!! WARNING : error_x is smaller than -1 !!!" << endl ;
122 double ff =
pow( (2.+error_x)/(2.+2.*error_x), rot_exp_y) ;
126 if ( fabs(error_x) < thres_rot ) {
127 cout <<
"Bin_bhns::rotation_axis:" << endl ;
128 cout <<
" ff is set to 1 because error_x is smaller than" << endl ;
129 cout <<
" the threshold value (" << thres_rot <<
")" << endl ;
134 cout <<
"Local center of mass of NS:" << endl ;
138 cout <<
"Bin_bhns::rotation_axis:" << endl ;
139 cout <<
" error_x : " << error_x <<
" ff : " << ff << endl ;
144 double aa = fact *
separ ;
145 double yns_new = yns_old + aa * (1. - ff) ;
147 cout <<
" new Y_BH : " << 0. / km <<
" [km]" 148 <<
" new Y_NS : " << yns_new / km <<
" [km]" 155 (
star.
set_mp()).set_ori(xns_old, yns_new, 0.) ;
double get_mass_bh() const
Returns the gravitational mass of BH [{ m_unit}].
Hole_bhns hole
Black hole.
const Tbl & line_mom_bhns() const
Total linear momentum.
double ya_barycenter() const
Absolute coordinate Y of the barycenter of the baryon density.
double & set_x_rot()
Sets the absolute coordinate X of the rotation axis [{ r_unit}].
double get_ori_y() const
Returns the y coordinate of the origin.
Standard units of space, time and mass.
double xa_barycenter() const
Absolute coordinate X of the barycenter of the baryon density.
double get_ori_x() const
Returns the x coordinate of the origin.
const Map & get_mp() const
Returns the mapping.
double separ
Absolute orbital separation between two centers of BH and NS.
double omega
Angular velocity with respect to an asymptotically inertial observer.
Map & set_mp()
Read/write of the mapping.
Star_bhns star
Neutron star.
const Map & get_mp() const
Returns the mapping.
void rotation_axis_x(double rot_exp_x)
Computes the position of the rotation axis X.
Cmp pow(const Cmp &, int)
Power .
double & set_y_rot()
Sets the absolute coordinate Y of the rotation axis [{ r_unit}].
void rotation_axis_y(double thres_rot, double rot_exp_y, double fact)
Computes the position of the rotation axis Y.
Map & set_mp()
Read/write of the mapping.