29 char star_bin_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/star_bin.C,v 1.19 2014/10/13 08:53:37 j_novak Exp $" ;
107 Cmp raccord_c1(
const Cmp& uu,
int l1) ;
116 bool irrot,
bool conf_flat0)
117 :
Star(mpi, nzet_i, eos_i),
120 d_psi(mpi, COV, mpi.get_bvect_cart()),
121 wit_w(mpi, CON, mpi.get_bvect_cart()),
123 bsn(mpi, CON, mpi.get_bvect_cart()),
127 dcov_logn(mpi, COV, mpi.get_bvect_cart()),
128 dcon_logn(mpi, CON, mpi.get_bvect_cart()),
132 dcov_phi(mpi, COV, mpi.get_bvect_cart()),
133 dcon_phi(mpi, CON, mpi.get_bvect_cart()),
134 flat(mpi, mpi.get_bvect_cart()),
136 beta_auto(mpi, CON, mpi.get_bvect_cart()),
137 beta_comp(mpi, CON, mpi.get_bvect_cart()),
138 hij(mpi, CON, mpi.get_bvect_cart()),
139 hij_auto(mpi, CON, mpi.get_bvect_cart()),
140 hij_comp(mpi, CON, mpi.get_bvect_cart()),
141 tkij_auto(mpi, CON, mpi.get_bvect_cart()),
142 tkij_comp(mpi, CON, mpi.get_bvect_cart()),
148 ssjm1_wbeta(mpi, CON, mpi.get_bvect_cart()),
156 conf_flat(conf_flat0){
212 irrotational(star.irrotational),
218 pot_centri(star.pot_centri),
219 logn_auto(star.logn_auto),
220 logn_comp(star.logn_comp),
221 dcov_logn(star.dcov_logn),
222 dcon_logn(star.dcon_logn),
223 lnq_auto(star.lnq_auto),
224 lnq_comp(star.lnq_comp),
226 dcov_phi(star.dcov_phi),
227 dcon_phi(star.dcon_phi),
230 beta_auto(star.beta_auto),
231 beta_comp(star.beta_comp),
233 hij_auto(star.hij_auto),
234 hij_comp(star.hij_comp),
235 tkij_auto(star.tkij_auto),
236 tkij_comp(star.tkij_comp),
237 kcar_auto(star.kcar_auto),
238 kcar_comp(star.kcar_comp),
239 ssjm1_logn(star.ssjm1_logn),
240 ssjm1_lnq(star.ssjm1_lnq),
241 ssjm1_khi(star.ssjm1_khi),
242 ssjm1_wbeta(star.ssjm1_wbeta),
243 ssjm1_h11(star.ssjm1_h11),
244 ssjm1_h21(star.ssjm1_h21),
245 ssjm1_h31(star.ssjm1_h31),
246 ssjm1_h22(star.ssjm1_h22),
247 ssjm1_h32(star.ssjm1_h32),
248 ssjm1_h33(star.ssjm1_h33),
249 decouple(star.decouple),
250 conf_flat(star.conf_flat)
259 :
Star(mpi, eos_i, fich),
261 d_psi(mpi, COV, mpi.get_bvect_cart()),
262 wit_w(mpi, CON, mpi.get_bvect_cart()),
264 bsn(mpi, CON, mpi.get_bvect_cart()),
266 logn_auto(mpi, *(mpi.get_mg()), fich),
268 dcov_logn(mpi, COV, mpi.get_bvect_cart()),
269 dcon_logn(mpi, CON, mpi.get_bvect_cart()),
270 lnq_auto(mpi, *(mpi.get_mg()), fich),
273 dcov_phi(mpi, COV, mpi.get_bvect_cart()),
274 dcon_phi(mpi, CON, mpi.get_bvect_cart()),
275 flat(mpi, mpi.get_bvect_cart()),
277 beta_auto(mpi, mpi.get_bvect_cart(), fich),
278 beta_comp(mpi, CON, mpi.get_bvect_cart()),
279 hij(mpi, CON, mpi.get_bvect_cart()),
280 hij_auto(mpi, mpi.get_bvect_cart(), fich),
281 hij_comp(mpi, CON, mpi.get_bvect_cart()),
282 tkij_auto(mpi, CON, mpi.get_bvect_cart()),
283 tkij_comp(mpi, CON, mpi.get_bvect_cart()),
286 ssjm1_logn(mpi, *(mpi.get_mg()), fich),
287 ssjm1_lnq(mpi, *(mpi.get_mg()), fich),
288 ssjm1_khi(mpi, *(mpi.get_mg()), fich),
289 ssjm1_wbeta(mpi, mpi.get_bvect_cart(), fich),
290 ssjm1_h11(mpi, *(mpi.get_mg()), fich),
291 ssjm1_h21(mpi, *(mpi.get_mg()), fich),
292 ssjm1_h31(mpi, *(mpi.get_mg()), fich),
293 ssjm1_h22(mpi, *(mpi.get_mg()), fich),
294 ssjm1_h32(mpi, *(mpi.get_mg()), fich),
295 ssjm1_h33(mpi, *(mpi.get_mg()), fich),
303 fread(&
conf_flat,
sizeof(
bool), 1, fich) ;
440 ssjm1_wbeta = star.ssjm1_wbeta ;
501 ssjm1_wbeta.
sauve(fich) ;
510 fwrite(&
conf_flat,
sizeof(
bool), 1, fich) ;
529 ost <<
"Star in a binary system" << endl ;
530 ost <<
"-----------------------" << endl ;
533 ost <<
"irrotational configuration" << endl ;
536 ost <<
"corotating configuration" << endl ;
539 ost <<
"Absolute abscidia of the stellar center: " <<
542 ost <<
"Absolute abscidia of the barycenter of the baryon density : " <<
547 double d_tilde = 2 * d_ns / r_0 ;
549 ost <<
"d_tilde : " << d_tilde << endl ;
551 ost <<
"Central value of gam_euler : "
554 ost <<
"Central u_euler (U^r, U^t, U^p) [c] : "
555 <<
u_euler(1).val_grid_point(0, 0, 0, 0) <<
" "
556 <<
u_euler(2).val_grid_point(0, 0, 0, 0) <<
" "
557 <<
u_euler(3).val_grid_point(0, 0, 0, 0) << endl ;
560 ost <<
"Central d_psi (r, t, p) [c] : "
561 <<
d_psi(1).val_grid_point(0, 0, 0, 0) <<
" "
562 <<
d_psi(2).val_grid_point(0, 0, 0, 0) <<
" "
563 <<
d_psi(3).val_grid_point(0, 0, 0, 0) << endl ;
565 ost <<
"Central vel. / co-orb. (W^r, W^t, W^p) [c] : "
566 <<
wit_w(1).val_grid_point(0, 0, 0, 0) <<
" "
567 <<
wit_w(2).val_grid_point(0, 0, 0, 0) <<
" "
568 <<
wit_w(3).val_grid_point(0, 0, 0, 0) << endl ;
570 ost <<
"Max vel. / co-orb. (W^r, W^t, W^p) [c] : "
575 ost <<
"Min vel. / co-orb. (W^r, W^t, W^p) [c] : "
580 double r_surf =
mp.
val_r(0,1.,M_PI/4,M_PI/4) ;
582 ost <<
"Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : "
583 <<
wit_w(1).val_point(r_surf,M_PI/4,M_PI/4) <<
" "
584 <<
wit_w(2).val_point(r_surf,M_PI/4,M_PI/4) <<
" "
585 <<
wit_w(3).val_point(r_surf,M_PI/4,M_PI/4) << endl ;
587 ost <<
"Central value of loggam : "
592 ost <<
"Central value of log(N) auto, comp : "
596 ost <<
"Central value of beta (N^r, N^t, N^p) [c] : "
597 <<
beta(1).val_grid_point(0, 0, 0, 0) <<
" "
598 <<
beta(2).val_grid_point(0, 0, 0, 0) <<
" "
599 <<
beta(3).val_grid_point(0, 0, 0, 0) << endl ;
601 ost <<
" ... beta_auto part of it [c] : "
602 <<
beta_auto(1).val_grid_point(0, 0, 0, 0) <<
" "
603 <<
beta_auto(2).val_grid_point(0, 0, 0, 0) <<
" "
604 <<
beta_auto(3).val_grid_point(0, 0, 0, 0) << endl ;
606 ost << endl <<
"Central value of (B^r, B^t, B^p)/N [c] : "
607 <<
bsn(1).val_grid_point(0, 0, 0, 0) <<
" "
608 <<
bsn(2).val_grid_point(0, 0, 0, 0) <<
" "
609 <<
bsn(3).val_grid_point(0, 0, 0, 0) << endl ;
612 ost << endl <<
"Central A^{ij} [c/km] : " << endl ;
613 ost <<
" A^{xx} auto, comp : "
614 <<
tkij_auto(1, 1).val_grid_point(0, 0, 0, 0) * km <<
" "
615 <<
tkij_comp(1, 1).val_grid_point(0, 0, 0, 0) * km << endl ;
616 ost <<
" A^{xy} auto, comp : "
617 <<
tkij_auto(1, 2).val_grid_point(0, 0, 0, 0) * km <<
" "
618 <<
tkij_comp(1, 2).val_grid_point(0, 0, 0, 0) * km << endl ;
619 ost <<
" A^{xz} auto, comp : "
620 <<
tkij_auto(1, 3).val_grid_point(0, 0, 0, 0) * km <<
" "
621 <<
tkij_comp(1, 3).val_grid_point(0, 0, 0, 0) * km << endl ;
622 ost <<
" A^{yy} auto, comp : "
623 <<
tkij_auto(2, 2).val_grid_point(0, 0, 0, 0) * km <<
" "
624 <<
tkij_comp(2, 2).val_grid_point(0, 0, 0, 0) * km << endl ;
625 ost <<
" A^{yz} auto, comp : "
626 <<
tkij_auto(2, 3).val_grid_point(0, 0, 0, 0) * km <<
" "
627 <<
tkij_comp(2, 3).val_grid_point(0, 0, 0, 0) * km << endl ;
628 ost <<
" A^{zz} auto, comp : "
629 <<
tkij_auto(3, 3).val_grid_point(0, 0, 0, 0) * km <<
" "
630 <<
tkij_comp(3, 3).val_grid_point(0, 0, 0, 0) * km << endl ;
632 ost << endl <<
"Central A_{ij} A^{ij} [c^2/km^2] : " << endl ;
633 ost <<
" A_{ij} A^{ij} auto, comp : "
667 for (
int i=1; i<=3; i++) {
668 v_orb.
set(i) = www(i).val_grid_point(0, 0, 0, 0) ;
679 d_psi = d_psi0 + v_orb ;
680 for (
int i=1; i<=3; i++) {
681 if (
d_psi(i).get_etat() == ETATZERO)
691 for (
int i=1; i<=3; i++) {
693 d_psi_i.
va.
set_base( d_psi0(i).get_spectral_va().base ) ;
694 d_psi_i = raccord_c1(d_psi_i,
nzet) ;
702 double relax_met,
int mer,
int fmer_met) {
704 double relax_ent_jm1 = 1. - relax_ent ;
705 double relax_met_jm1 = 1. - relax_met ;
707 ent = relax_ent *
ent + relax_ent_jm1 * star_jm1.
ent ;
709 if ( (mer != 0) && (mer % fmer_met == 0)) {
733 cout <<
"La jauge de Dirac est elle bien satisfaite ??" << endl ;
734 cout <<
"Vector Hi" << endl ;
735 for (
int i=1; i<=3; i++)
740 cout <<
"Pour comparaison valeur de D_i(g^1i)" << endl ;
741 for (
int i=1; i<=3; i++)
743 (1, i, i)/(nr*nt*np)) << endl ;