83 int ii = (
star.
get_mp()).get_mg()->get_nr(0) - 1 ;
84 int jj = (
star.
get_mp()).get_mg()->get_nt(0) - 1 ;
86 int kb = (
star.
get_mp()).get_mg()->get_np(0) / 2 ;
102 double shiftx_a = shift(1).val_grid_point(0,ka,jj,ii) ;
103 double shiftx_b = shift(1).val_grid_point(0,kb,jj,ii) ;
104 double shifty_a = shift(2).val_grid_point(0,ka,jj,ii) ;
105 double shifty_b = shift(2).val_grid_point(0,kb,jj,ii) ;
106 double shiftz_a = shift(3).val_grid_point(0,ka,jj,ii) ;
107 double shiftz_b = shift(3).val_grid_point(0,kb,jj,ii) ;
117 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
189 double sa = shiftx_a*shiftx_a+shifty_a*shifty_a+shiftz_a*shiftz_a ;
190 double sb = shiftx_b*shiftx_b+shifty_b*shifty_b+shiftz_b*shiftz_b ;
192 double ta = -shiftx_a*yns_rot + shifty_a*(ra+xns_rot) ;
193 double tb = -shiftx_b*yns_rot + shifty_b*(-rb+xns_rot) ;
195 double ua = yns_rot*yns_rot + (ra+xns_rot)*(ra+xns_rot) ;
196 double ub = yns_rot*yns_rot + (-rb+xns_rot)*(-rb+xns_rot) ;
201 double aaa = psi4_a * gam2_a * ua - psi4_b * gam2_b * ub ;
202 double bbb = psi4_a * gam2_a * ta - psi4_b * gam2_b * tb ;
203 double ccc = psi4_a * gam2_a * sa - psi4_b * gam2_b * sb
204 - lap2_a * gam2_a / con2_a + lap2_b * gam2_b / con2_b ;
209 double ddd = bbb*bbb - aaa*ccc ;
213 "!!! WARNING : Omega (from two points) does not exist !!!"
220 double omega_1 = (-bbb +
sqrt(ddd)) / aaa ;
221 double omega_2 = (-bbb -
sqrt(ddd)) / aaa ;
223 cout <<
"Bin_bhns::omega_two_points:" << endl ;
224 cout <<
" omega_1 : " << omega_1 * f_unit <<
" [rad/s]"
226 cout <<
" omega_2 : " << omega_2 * f_unit <<
" [rad/s]"
229 omega_two = omega_1 ;