67 double compact = ggrav * mtot / rr ;
78 omega2 = ggrav * mtot /
pow(rr, 3)
79 * (1. - 2.75 * compact + 8.625 * compact*compact ) ;
91 double ired = double(5)/double(3) * ( 1. - double(6) / M_PI / M_PI ) ;
92 omega2 = ggrav * mtot /
pow(rr, 3)
93 * (1. - compact * ( 2.75 + 2.*a0sr*a0sr * ired
94 - 0.48*
pow(a0sr, 4) * ired*ired )
95 + compact*compact * ( 8.625 + 2.75*a0sr*a0sr * ired
96 + 2.*
pow(a0sr, 4) * ired*ired ) ) ;
Cmp sqrt(const Cmp &)
Square root.
Star_bin star1
First star of the system.
Standard units of space, time and mass.
double omega
Angular velocity with respect to an asymptotically inertial observer.
double separation() const
Returns the coordinate separation of the two stellar centers [r_unit].
void del_deriv() const
Deletes all the derived quantities.
Cmp pow(const Cmp &, int)
Power .
bool is_irrotational() const
Returns true for an irrotational motion, false for a corotating one.
double ray_eq() const
Coordinate radius at , [r_unit].
virtual double mass_g() const
Gravitational mass.
Star_bin star2
Second star of the system.
void analytical_omega()
Sets the orbital angular velocity to some 2-PN analytical value (Keplerian value in the Newtonian cas...