28 char 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 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.
double get_ori_y() const
Returns the y coordinate of the origin.
const Tbl & line_mom_bhns() const
Total linear momentum.
Cmp pow(const Cmp &, int)
Power .
void rotation_axis_x(double rot_exp_x)
Computes the position of the rotation axis X.
double get_mass_bh() const
Returns the gravitational mass of BH [{\tt m_unit}].
Map & set_mp()
Read/write of the mapping.
Standard units of space, time and mass.
void rotation_axis_y(double thres_rot, double rot_exp_y, double fact)
Computes the position of the rotation axis Y.
double separ
Absolute orbital separation between two centers of BH and NS.
double omega
Angular velocity with respect to an asymptotically inertial observer.
double xa_barycenter() const
Absolute coordinate X of the barycenter of the baryon density.
double & set_x_rot()
Sets the absolute coordinate X of the rotation axis [{\tt r_unit}].
const Map & get_mp() const
Returns the mapping.
Star_bhns star
Neutron star.
double ya_barycenter() const
Absolute coordinate Y of the barycenter of the baryon density.
double & set_y_rot()
Sets the absolute coordinate Y of the rotation axis [{\tt r_unit}].