16 double rp0 = 13.955599;
17 double e0 = 0.3749609;
22 getEBBH(m1, m2, rp0, e0, hp, hx);
26 cout <<
"distance source Kpc : " << distance_source_Kpc <<
" (Kpc) " << endl;
33 for(
int i=0;
i<hp.
size();
i++) hrss_hp+=hp[
i]*hp[
i];
34 hrss_hp*=
dt; hrss_hp=sqrt(hrss_hp);
38 for(
int i=0;
i<hx.
size();
i++) hrss_hx+=hx[
i]*hx[
i];
39 hrss_hx*=
dt; hrss_hx=sqrt(hrss_hx);
42 hrss_hp*=distance_source_Kpc/10.;
43 hrss_hx*=distance_source_Kpc/10.;
45 cout <<
"hrss_hp @ 10Kpc : " << hrss_hp << endl;
46 cout <<
"hrss_hx @ 10Kpc : " << hrss_hx << endl;
49 double strain_10kpc=sqrt(hrss_hp*hrss_hp+hrss_hx*hrss_hx);
50 cout <<
"strain @ 10Kpc : " << strain_10kpc << endl;
53 double strain = 10*strain_10kpc/
distance;
54 cout <<
"strain @ 4509.8798 Kpc : " << strain << endl;
double distance_source_Kpc
virtual void rate(double r)
virtual size_t size() const
double GravitationalConstant()
getEBBH(m1, m2, rp0, e0, hp, hx)
double SpeedOfLightInVacuo()