30 char etoile_bin_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/etoile_bin.C,v 1.13 2014/10/13 08:52:58 j_novak Exp $" ;
199 Cmp raccord_c1(
const Cmp& uu,
int l1) ;
208 bool irrot,
const Base_vect& ref_triad_i)
209 :
Etoile(mpi, nzet_i, relat, eos_i),
211 ref_triad(ref_triad_i),
213 d_psi(mpi, 1, COV, ref_triad),
214 wit_w(mpi, 1, CON, ref_triad),
217 d_logn_auto(mpi, 1, COV, ref_triad),
218 d_logn_auto_regu(mpi, 1, COV, ref_triad),
219 d_logn_comp(mpi, 1, COV, ref_triad),
221 d_beta_auto(mpi, 1, COV, ref_triad),
222 d_beta_comp(mpi, 1, COV, ref_triad),
223 shift_auto(mpi, 1, CON, ref_triad),
224 shift_comp(mpi, 1, CON, ref_triad),
225 w_shift(mpi, 1, CON, mp.get_bvect_cart()),
227 tkij_auto(mpi, 2, CON, ref_triad),
228 tkij_comp(mpi, 2, CON, ref_triad),
231 bsn(mpi, 1, CON, ref_triad),
236 ssjm1_wshift(mpi, 1, CON, mp.get_bvect_cart()),
264 for (
int i=0; i<3; i++) {
283 for (
int i=0; i<3; i++) {
293 irrotational(et.irrotational),
294 ref_triad(et.ref_triad),
299 logn_comp(et.logn_comp),
300 d_logn_auto(et.d_logn_auto),
301 d_logn_auto_regu(et.d_logn_auto_regu),
302 d_logn_comp(et.d_logn_comp),
303 beta_comp(et.beta_comp),
304 d_beta_auto(et.d_beta_auto),
305 d_beta_comp(et.d_beta_comp),
306 shift_auto(et.shift_auto),
307 shift_comp(et.shift_comp),
309 khi_shift(et.khi_shift),
310 tkij_auto(et.tkij_auto),
311 tkij_comp(et.tkij_comp),
312 akcar_auto(et.akcar_auto),
313 akcar_comp(et.akcar_comp),
315 pot_centri(et.pot_centri),
316 ssjm1_logn(et.ssjm1_logn),
317 ssjm1_beta(et.ssjm1_beta),
318 ssjm1_khi(et.ssjm1_khi),
319 ssjm1_wshift(et.ssjm1_wshift),
320 ssjm1_psi(et.ssjm1_khi),
321 decouple(et.decouple)
331 :
Etoile(mpi, eos_i, fich),
332 ref_triad(ref_triad_i),
334 d_psi(mpi, 1, COV, ref_triad),
335 wit_w(mpi, 1, CON, ref_triad),
338 d_logn_auto(mpi, 1, COV, ref_triad),
339 d_logn_auto_regu(mpi, 1, COV, ref_triad),
340 d_logn_comp(mpi, 1, COV, ref_triad),
342 d_beta_auto(mpi, 1, COV, ref_triad),
343 d_beta_comp(mpi, 1, COV, ref_triad),
344 shift_auto(mpi, 1, CON, ref_triad),
345 shift_comp(mpi, 1, CON, ref_triad),
346 w_shift(mpi, 1, CON, mp.get_bvect_cart()),
348 tkij_auto(mpi, 2, CON, ref_triad),
349 tkij_comp(mpi, 2, CON, ref_triad),
352 bsn(mpi, 1, CON, ref_triad),
357 ssjm1_wshift(mpi, 1, CON, mp.get_bvect_cart()),
589 ost <<
"Star in a binary system" << endl ;
590 ost <<
"-----------------------" << endl ;
593 ost <<
"irrotational configuration" << endl ;
596 ost <<
"corotating configuration" << endl ;
599 ost <<
"Absolute abscidia of the stellar center: " <<
602 ost <<
"Absolute abscidia of the barycenter of the baryon density : " <<
607 double d_tilde = 2 * d_ns / r_0 ;
609 ost <<
"d_tilde : " << d_tilde << endl ;
611 ost <<
"Orientation with respect to the absolute frame : " <<
614 ost <<
"Central value of gam_euler : "
617 ost <<
"Central u_euler (U^X, U^Y, U^Z) [c] : "
618 <<
u_euler(0)(0, 0, 0, 0) <<
" "
619 <<
u_euler(1)(0, 0, 0, 0) <<
" "
620 <<
u_euler(2)(0, 0, 0, 0) << endl ;
623 ost <<
"Central d_psi (X, Y, Z) [c] : "
624 <<
d_psi(0)(0, 0, 0, 0) <<
" "
625 <<
d_psi(1)(0, 0, 0, 0) <<
" "
626 <<
d_psi(2)(0, 0, 0, 0) << endl ;
628 ost <<
"Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : "
629 <<
wit_w(0)(0, 0, 0, 0) <<
" "
630 <<
wit_w(1)(0, 0, 0, 0) <<
" "
631 <<
wit_w(2)(0, 0, 0, 0) << endl ;
633 ost <<
"Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : "
638 ost <<
"Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : "
643 double r_surf =
mp.
val_r(0,1.,M_PI/4,M_PI/4) ;
645 ost <<
"Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : "
646 <<
wit_w(0).val_point(r_surf,M_PI/4,M_PI/4) <<
" "
647 <<
wit_w(1).val_point(r_surf,M_PI/4,M_PI/4) <<
" "
648 <<
wit_w(2).val_point(r_surf,M_PI/4,M_PI/4) << endl ;
650 ost <<
"Central value of loggam : "
651 <<
loggam()(0, 0, 0, 0) << endl ;
655 ost <<
"Central value of log(N) auto, comp : "
659 ost <<
"Central value of beta=log(AN) auto, comp : "
663 ost <<
"Central value of shift (N^X, N^Y, N^Z) [c] : "
664 <<
shift(0)(0, 0, 0, 0) <<
" "
665 <<
shift(1)(0, 0, 0, 0) <<
" "
666 <<
shift(2)(0, 0, 0, 0) << endl ;
668 ost <<
" ... shift_auto part of it [c] : "
673 ost <<
" ... shift_comp part of it [c] : "
678 ost <<
" ... w_shift (NB: components in the star Cartesian frame) [c] : "
680 <<
w_shift(0)(0, 0, 0, 0) <<
" "
681 <<
w_shift(1)(0, 0, 0, 0) <<
" "
682 <<
w_shift(2)(0, 0, 0, 0) << endl ;
684 ost <<
"Central value of khi_shift [km c] : "
685 <<
khi_shift()(0, 0, 0, 0) / km << endl ;
687 ost << endl <<
"Central value of (B^X, B^Y, B^Z)/N [c] : "
688 <<
bsn(0)(0, 0, 0, 0) <<
" "
689 <<
bsn(1)(0, 0, 0, 0) <<
" "
690 <<
bsn(2)(0, 0, 0, 0) << endl ;
693 "Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : "
698 ost <<
"Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : "
704 "Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : "
709 ost <<
"Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : "
715 ost << endl <<
"Central A^2 K^{ij} [c/km] : " << endl ;
716 ost <<
" A^2 K^{xx} auto, comp : "
717 <<
tkij_auto(0, 0)(0, 0, 0, 0) * km <<
" "
718 <<
tkij_comp(0, 0)(0, 0, 0, 0) * km << endl ;
719 ost <<
" A^2 K^{xy} auto, comp : "
720 <<
tkij_auto(0, 1)(0, 0, 0, 0) * km <<
" "
721 <<
tkij_comp(0, 1)(0, 0, 0, 0) * km << endl ;
722 ost <<
" A^2 K^{xz} auto, comp : "
723 <<
tkij_auto(0, 2)(0, 0, 0, 0) * km <<
" "
724 <<
tkij_comp(0, 2)(0, 0, 0, 0) * km << endl ;
725 ost <<
" A^2 K^{yy} auto, comp : "
726 <<
tkij_auto(1, 1)(0, 0, 0, 0) * km <<
" "
727 <<
tkij_comp(1, 1)(0, 0, 0, 0) * km << endl ;
728 ost <<
" A^2 K^{yz} auto, comp : "
729 <<
tkij_auto(1, 2)(0, 0, 0, 0) * km <<
" "
730 <<
tkij_comp(1, 2)(0, 0, 0, 0) * km << endl ;
731 ost <<
" A^2 K^{zz} auto, comp : "
732 <<
tkij_auto(2, 2)(0, 0, 0, 0) * km <<
" "
733 <<
tkij_comp(2, 2)(0, 0, 0, 0) * km << endl ;
735 ost << endl <<
"Central A^2 K_{ij} K^{ij} [c^2/km^2] : " << endl ;
736 ost <<
" A^2 K_{ij} K^{ij} auto, comp : "
788 for (
int i=0; i<3; i++) {
789 v_orb.
set(i) = www(i)(0, 0, 0, 0) ;
794 for (
int l=
nzet; l<=nzm1; l++)
795 for (
int i=0; i<=2; i++)
806 d_psi = d_psi0 + v_orb ;
813 if (d_psi0.
get_etat() == ETATQCQ ) {
815 for (
int i=0; i<3; i++) {
840 double lambda = double(1) / double(3) ;
848 for (
int i=0; i<3; i++) {
850 - (lambda/2./(lambda+1)) * (d_khi(i) + x_d_w(i)) ;
864 double relax_met,
int mer,
int fmer_met) {
866 double relax_ent_jm1 = 1. - relax_ent ;
867 double relax_met_jm1 = 1. - relax_met ;
869 ent = relax_ent *
ent + relax_ent_jm1 * star_jm1.
ent ;
871 if ( (mer != 0) && (mer % fmer_met == 0)) {