26 char star_bin_xcts_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/star_bin_xcts.C,v 1.9 2014/10/13 08:53:39 j_novak Exp $" ;
73 Cmp raccord_c1(
const Cmp& uu,
int l1) ;
85 :
Star(mpi, nzet_i, eos_i),
88 d_psi(mpi, COV, mpi.get_bvect_cart()),
89 wit_w(mpi, CON, mpi.get_bvect_cart()),
91 bsn(mpi, CON, mpi.get_bvect_cart()),
100 w_beta(mpi, CON, mpi.get_bvect_cart()),
102 dcov_Psi(mpi, COV, mpi.get_bvect_cart()),
103 dcov_chi(mpi, COV, mpi.get_bvect_cart()),
104 flat(mpi, mpi.get_bvect_cart()),
105 beta_auto(mpi, CON, mpi.get_bvect_cart()),
106 beta_comp(mpi, CON, mpi.get_bvect_cart()),
107 haij_auto(mpi, CON, mpi.get_bvect_cart()),
108 haij_comp(mpi, CON, mpi.get_bvect_cart()),
114 ssjm1_wbeta(mpi, CON, mpi.get_bvect_cart()) {
167 irrotational(star.irrotational),
173 pot_centri(star.pot_centri),
174 chi_auto(star.chi_auto),
175 chi_comp(star.chi_comp),
176 Psi_auto(star.Psi_auto),
177 Psi_comp(star.Psi_comp),
183 dcov_Psi(star.dcov_Psi),
184 dcov_chi(star.dcov_chi),
186 beta_auto(star.beta_auto),
187 beta_comp(star.beta_comp),
188 haij_auto(star.haij_auto),
189 haij_comp(star.haij_comp),
190 hacar_auto(star.hacar_auto),
191 hacar_comp(star.hacar_comp),
192 ssjm1_chi(star.ssjm1_chi),
193 ssjm1_psi(star.ssjm1_psi),
194 ssjm1_khi(star.ssjm1_khi),
195 ssjm1_wbeta(star.ssjm1_wbeta) {
206 :
Star(mpi, eos_i, fich),
208 d_psi(mpi, COV, mpi.get_bvect_cart()),
209 wit_w(mpi, CON, mpi.get_bvect_cart()),
211 bsn(mpi, CON, mpi.get_bvect_cart()),
213 chi_auto(mpi, *(mpi.get_mg()), fich),
215 Psi_auto(mpi, *(mpi.get_mg()), fich),
220 w_beta(mpi, mpi.get_bvect_cart(), fich),
221 khi(mpi, *(mpi.get_mg()), fich),
222 dcov_Psi(mpi, COV, mpi.get_bvect_cart()),
223 dcov_chi(mpi, COV, mpi.get_bvect_cart()),
224 flat(mpi, mpi.get_bvect_cart()),
225 beta_auto(mpi, mpi.get_bvect_cart(), fich),
226 beta_comp(mpi, CON, mpi.get_bvect_cart()),
227 haij_auto(mpi, CON, mpi.get_bvect_cart()),
228 haij_comp(mpi, CON, mpi.get_bvect_cart()),
231 ssjm1_chi(mpi, *(mpi.get_mg()), fich),
232 ssjm1_psi(mpi, *(mpi.get_mg()), fich),
233 ssjm1_khi(mpi, *(mpi.get_mg()), fich),
234 ssjm1_wbeta(mpi, mpi.get_bvect_cart(), fich) {
240 bool status = fread(&
irrotational,
sizeof(
bool), 1, fich) ;
242 cout <<
"Star_bin_xcts::Constructor from a file: Problem with reading ! " << endl ;
471 ost <<
"Number of domains occupied by the star : " <<
nzet << endl ;
473 ost <<
"Equation of state : " << endl ;
476 ost << endl <<
"Central enthalpy : " <<
ent.
val_grid_point(0,0,0,0) <<
" c^2" << endl ;
478 <<
" x 0.1 fm^-3" << endl ;
480 <<
" rho_nuc c^2" << endl ;
482 <<
" rho_nuc c^2" << endl ;
490 ost <<
"Central value of Psi^4 : " << psi4_local.
val_grid_point(0,0,0,0) << endl ;
493 <<
"Coordinate equatorial radius (phi=0) a1 = "
494 <<
ray_eq()/km <<
" km" << endl ;
495 ost <<
"Coordinate equatorial radius (phi=pi/2) a2 = "
497 ost <<
"Coordinate equatorial radius (phi=pi): "
499 ost <<
"Coordinate polar radius a3 = "
506 double mass_shedd_chi = fabs( dent_eq / dent_pole ) ;
508 ost <<
"Mass-shedding estimator = " << mass_shedd_chi << endl ;
510 ost << endl <<
"Baryon mass : " <<
mass_b() / msol <<
" M_sol" << endl ;
511 ost <<
"Gravitational mass : " <<
mass_g() / msol <<
" M_sol" << endl ;
516 ost <<
"Star in a binary system" << endl ;
517 ost <<
"-----------------------" << endl ;
520 ost <<
"irrotational configuration" << endl ;
523 ost <<
"corotating configuration" << endl ;
526 ost <<
"Absolute abscidia of the stellar center: " <<
529 ost <<
"Absolute abscidia of the barycenter of the baryon density : " <<
534 double d_tilde = 2 * d_ns / r_0 ;
536 ost <<
"d_tilde : " << d_tilde << endl ;
538 ost <<
"Central value of gam_euler : "
541 ost <<
"Central u_euler (U^r, U^t, U^p) [c] : "
542 <<
u_euler(1).val_grid_point(0, 0, 0, 0) <<
" "
543 <<
u_euler(2).val_grid_point(0, 0, 0, 0) <<
" "
544 <<
u_euler(3).val_grid_point(0, 0, 0, 0) << endl ;
547 ost <<
"Central d_psi (r, t, p) [c] : "
548 <<
d_psi(1).val_grid_point(0, 0, 0, 0) <<
" "
549 <<
d_psi(2).val_grid_point(0, 0, 0, 0) <<
" "
550 <<
d_psi(3).val_grid_point(0, 0, 0, 0) << endl ;
552 ost <<
"Central vel. / co-orb. (W^r, W^t, W^p) [c] : "
553 <<
wit_w(1).val_grid_point(0, 0, 0, 0) <<
" "
554 <<
wit_w(2).val_grid_point(0, 0, 0, 0) <<
" "
555 <<
wit_w(3).val_grid_point(0, 0, 0, 0) << endl ;
557 ost <<
"Max vel. / co-orb. (W^r, W^t, W^p) [c] : "
562 ost <<
"Min vel. / co-orb. (W^r, W^t, W^p) [c] : "
567 double r_surf =
mp.
val_r(0,1.,M_PI/4,M_PI/4) ;
569 ost <<
"Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : "
570 <<
wit_w(1).val_point(r_surf,M_PI/4,M_PI/4) <<
" "
571 <<
wit_w(2).val_point(r_surf,M_PI/4,M_PI/4) <<
" "
572 <<
wit_w(3).val_point(r_surf,M_PI/4,M_PI/4) << endl ;
574 ost <<
"Central value of loggam : "
578 ost <<
"Central value of Psi auto, comp : "
582 ost <<
"Central value of beta (N^r, N^t, N^p) [c] : "
583 <<
beta(1).val_grid_point(0, 0, 0, 0) <<
" "
584 <<
beta(2).val_grid_point(0, 0, 0, 0) <<
" "
585 <<
beta(3).val_grid_point(0, 0, 0, 0) << endl ;
587 ost <<
" ... beta_auto part of it [c] : "
588 <<
beta_auto(1).val_grid_point(0, 0, 0, 0) <<
" "
589 <<
beta_auto(2).val_grid_point(0, 0, 0, 0) <<
" "
590 <<
beta_auto(3).val_grid_point(0, 0, 0, 0) << endl ;
592 ost << endl <<
"Central value of (B^r, B^t, B^p)/N [c] : "
593 <<
bsn(1).val_grid_point(0, 0, 0, 0) <<
" "
594 <<
bsn(2).val_grid_point(0, 0, 0, 0) <<
" "
595 <<
bsn(3).val_grid_point(0, 0, 0, 0) << endl ;
598 ost << endl <<
"Central \\hat{A}^{ij} [c/km] : " << endl ;
599 ost <<
" \\hat{A}^{xx} auto, comp : "
600 <<
haij_auto(1, 1).val_grid_point(0, 0, 0, 0) * km <<
" "
601 <<
haij_comp(1, 1).val_grid_point(0, 0, 0, 0) * km << endl ;
602 ost <<
" A^{xy} auto, comp : "
603 <<
haij_auto(1, 2).val_grid_point(0, 0, 0, 0) * km <<
" "
604 <<
haij_comp(1, 2).val_grid_point(0, 0, 0, 0) * km << endl ;
605 ost <<
" A^{xz} auto, comp : "
606 <<
haij_auto(1, 3).val_grid_point(0, 0, 0, 0) * km <<
" "
607 <<
haij_comp(1, 3).val_grid_point(0, 0, 0, 0) * km << endl ;
608 ost <<
" A^{yy} auto, comp : "
609 <<
haij_auto(2, 2).val_grid_point(0, 0, 0, 0) * km <<
" "
610 <<
haij_comp(2, 2).val_grid_point(0, 0, 0, 0) * km << endl ;
611 ost <<
" A^{yz} auto, comp : "
612 <<
haij_auto(2, 3).val_grid_point(0, 0, 0, 0) * km <<
" "
613 <<
haij_comp(2, 3).val_grid_point(0, 0, 0, 0) * km << endl ;
614 ost <<
" A^{zz} auto, comp : "
615 <<
haij_auto(3, 3).val_grid_point(0, 0, 0, 0) * km <<
" "
616 <<
haij_comp(3, 3).val_grid_point(0, 0, 0, 0) * km << endl ;
618 ost << endl <<
"Central \\hat{A}_{ij}\\hat{A}^{ij} [c^2/km^2] : "
620 ost <<
" \\hat{A}_{ij}\\hat{A}^{ij} auto, comp : "
653 for (
int i=1; i<=3; i++)
654 v_orb.
set(i) = www(i).val_grid_point(0, 0, 0, 0) ;
664 d_psi = d_psi0 + v_orb ;
666 for (
int i=1; i<=3; i++) {
667 if (
d_psi(i).get_etat() == ETATZERO)
678 for (
int i=1; i<=3; i++) {
680 d_psi_i.
va.
set_base( d_psi0(i).get_spectral_va().base ) ;
681 d_psi_i = raccord_c1(d_psi_i,
nzet) ;
695 double relax_ent_jm1 = 1. - relax_ent ;
696 double relax_met_jm1 = 1. - relax_met ;
698 ent = relax_ent *
ent + relax_ent_jm1 * star_jm1.
ent ;
700 if ( (mer != 0) && (mer % fmer_met == 0)) {
703 + relax_met_jm1 * star_jm1.
Psi_auto ;
706 + relax_met_jm1 * star_jm1.
chi_auto ;