Date: Fri Feb  4 21:29:00 2011

===================================================================

Relativistic computation                           Irrotational
Star 1  EOS
EOS of class Eos_AkmalPR (Akmal, Pandharipande & Ravenhall 1998) : 
  composition :  n,p,e,mu 
  model : A18+dv+UIX* 

Omega = 1181.18 rad/s                Orbital frequency f = 187.99 Hz
Omega_kepler = 1268.72 rad/s
Coordinate separation : 60 km
1/2 ADM mass :        1.34008 Mo
Total angular momentum : 8.04887 G M_sol^2 / c

Number of steps : 106

===================================================================
       Star no. 1
===================================================================
Grid : 
------ 
Number of domains: 4
  Domain #0: nr = 33, RARE; nt = 25, SYM; np = 24, NONSYM
  Domain #1: nr = 33, FIN; nt = 25, SYM; np = 24, NONSYM
  Domain #2: nr = 33, FIN; nt = 25, SYM; np = 24, NONSYM
  Domain #3: nr = 33, UNSURR; nt = 25, SYM; np = 24, NONSYM



Physical characteristics : 
-------------------------

Relativistic star
-----------------
Number of domains occupied by the star : 1
Equation of state : 
Star 1  EOS
EOS of class Eos_AkmalPR (Akmal, Pandharipande & Ravenhall 1998) : 
  composition :  n,p,e,mu 
  model : A18+dv+UIX* 


Central enthalpy : 0.234475 c^2
Central proper baryon density : 5.3984 x 0.1 fm^-3
Central proper energy density : 5.90629 rho_nuc c^2
Central pressure : 0.918605 rho_nuc c^2

Regularization index of the gravitational potential : k_div = 0
Central lapse N :      0.611883
Central value of A^2 : 2.34256

Coordinate equatorial radius (phi=0) a1 =    8.95868 km
Coordinate equatorial radius (phi=pi/2) a2 = 8.87232 km
Coordinate equatorial radius (phi=pi):       8.99179 km
Coordinate polar radius a3 =                 8.90535 km
Axis ratio a2/a1 = 0.99036  a3/a1 = 0.994047

Baryon mass :        1.4911 M_sol
Gravitational mass : 1.30987 M_sol

Star in a binary system
-----------------------
irrotational configuration
Absolute abscidia of the stellar center: -30 km
Absolute abscidia of the barycenter of the baryon density : -29.9973 km
d_tilde : 6.68875
Orientation with respect to the absolute frame : 0 rad
Central value of gam_euler : 1.00259
Central u_euler (U^X, U^Y, U^Z) [c] : 2.89454e-18  -0.0469609  0
Central d_psi (X, Y, Z) [c] :         8.59461e-18  -0.139439  0
Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : -9.65142e-18  -3.28349e-05  0
Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.316919  0.461017  0.00267077
Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.316919  -0.153913  -0.00267077
Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : 0.0215471  -0.0212666  4.53126e-06
Central value of loggam : 1.23373e-09
Central value of log(N) auto, comp :         -0.458238  -0.0329769
Central value of beta=log(AN) auto, comp :   -0.0650537  -0.000538967
Central value of shift (N^X, N^Y, N^Z) [c] : -1.10971e-17  -0.0894848  0
  ... shift_auto part of it [c] :            -2.07723e-19  -0.104556  0
  ... shift_comp part of it [c] :            0  0.0150716  0
  ... w_shift (NB: components in the star Cartesian frame) [c] :  
7.94197e-18  -0.113324  0
Central value of khi_shift [km c] : 9.01645e-16

Central value of (B^X, B^Y, B^Z)/N [c] : -1.8136e-17  0.0469283  0

Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : -9.59181e-05  2.58684e-18  0
Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : -0.00056837  4.35075e-15  0

Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : -5.87797e-05  1.34007e-18  0
Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : -6.24489e-05  3.63863e-17  0

Central A^2 K^{ij} [c/km] : 
  A^2 K^{xx} auto, comp : -5.63863e-15  -2.23328e-16
  A^2 K^{xy} auto, comp : 4.82403e-06  -5.92126e-05
  A^2 K^{xz} auto, comp : 0  0
  A^2 K^{yy} auto, comp : 8.54601e-15  4.46023e-16
  A^2 K^{yz} auto, comp : 0  0
  A^2 K^{zz} auto, comp : -2.90738e-15  -2.22695e-16

Central A^2 K_{ij} K^{ij} [c^2/km^2] : 
   A^2 K_{ij} K^{ij}  auto, comp : 2.88788e-10  -1.33833e-09


===================================================================
       Star no. 2
===================================================================
Grid : 
------ 
Number of domains: 4
  Domain #0: nr = 33, RARE; nt = 25, SYM; np = 24, NONSYM
  Domain #1: nr = 33, FIN; nt = 25, SYM; np = 24, NONSYM
  Domain #2: nr = 33, FIN; nt = 25, SYM; np = 24, NONSYM
  Domain #3: nr = 33, UNSURR; nt = 25, SYM; np = 24, NONSYM



Physical characteristics : 
-------------------------

Relativistic star
-----------------
Number of domains occupied by the star : 1
Equation of state : 
Star 2  EOS
EOS of class Eos_AkmalPR (Akmal, Pandharipande & Ravenhall 1998) : 
  composition :  n,p,e,mu 
  model : A18+dv+UIX* 


Central enthalpy : 0.234475 c^2
Central proper baryon density : 5.3984 x 0.1 fm^-3
Central proper energy density : 5.90629 rho_nuc c^2
Central pressure : 0.918605 rho_nuc c^2

Regularization index of the gravitational potential : k_div = 0
Central lapse N :      0.611883
Central value of A^2 : 2.34256

Coordinate equatorial radius (phi=0) a1 =    8.95868 km
Coordinate equatorial radius (phi=pi/2) a2 = 8.87232 km
Coordinate equatorial radius (phi=pi):       8.99179 km
Coordinate polar radius a3 =                 8.90535 km
Axis ratio a2/a1 = 0.99036  a3/a1 = 0.994047

Baryon mass :        1.4911 M_sol
Gravitational mass : 1.30987 M_sol

Star in a binary system
-----------------------
irrotational configuration
Absolute abscidia of the stellar center: 30 km
Absolute abscidia of the barycenter of the baryon density : 29.9973 km
d_tilde : 6.68875
Orientation with respect to the absolute frame : 3.14159 rad
Central value of gam_euler : 1.00259
Central u_euler (U^X, U^Y, U^Z) [c] : -2.89454e-18  0.0469609  0
Central d_psi (X, Y, Z) [c] :         -8.59461e-18  0.139439  0
Central vel. / co-orb. (W^X, W^Y, W^Z) [c] : 9.65142e-18  3.28349e-05  0
Max vel. / co-orb. (W^X, W^Y, W^Z) [c] : 0.316919  0.153913  0.00267077
Min vel. / co-orb. (W^X, W^Y, W^Z) [c] : -0.316919  -0.461017  -0.00267077
Velocity at (r_surf,pi/4,pi/4) / co-orb. [c] : -0.0215471  0.0212666  4.53126e-06
Central value of loggam : 1.23373e-09
Central value of log(N) auto, comp :         -0.458238  -0.0329769
Central value of beta=log(AN) auto, comp :   -0.0650537  -0.000538967
Central value of shift (N^X, N^Y, N^Z) [c] : 1.10971e-17  0.0894848  0
  ... shift_auto part of it [c] :            2.07723e-19  0.104556  0
  ... shift_comp part of it [c] :            0  -0.0150716  0
  ... w_shift (NB: components in the star Cartesian frame) [c] :  
7.94197e-18  -0.113324  0
Central value of khi_shift [km c] : 9.01645e-16

Central value of (B^X, B^Y, B^Z)/N [c] : 1.8136e-17  -0.0469283  0

Central (d/dX,d/dY,d/dZ)(logn_auto) [km^{-1}] : 9.59181e-05  -2.58684e-18  0
Central (d/dX,d/dY,d/dZ)(logn_comp) [km^{-1}] : 0.00056837  -4.35075e-15  0

Central (d/dX,d/dY,d/dZ)(beta_auto) [km^{-1}] : 5.87797e-05  -1.34007e-18  0
Central (d/dX,d/dY,d/dZ)(beta_comp) [km^{-1}] : 6.24489e-05  -3.63863e-17  0

Central A^2 K^{ij} [c/km] : 
  A^2 K^{xx} auto, comp : -5.63863e-15  -2.23328e-16
  A^2 K^{xy} auto, comp : 4.82403e-06  -5.92126e-05
  A^2 K^{xz} auto, comp : 0  0
  A^2 K^{yy} auto, comp : 8.54601e-15  4.46023e-16
  A^2 K^{yz} auto, comp : 0  0
  A^2 K^{zz} auto, comp : -2.90738e-15  -2.22695e-16

Central A^2 K_{ij} K^{ij} [c^2/km^2] : 
   A^2 K_{ij} K^{ij}  auto, comp : 2.88788e-10  -1.33833e-09



===================================================================
Diff_ent :  star 1 : 9.59568e-08   star 2 : 9.59568e-08
Relative difference between the baryon masses of the two stars : 0
dH/dx at r = 0 :  star 1 : -1.2319e-08   star 2 : -1.2319e-08
Relative error on the virial theorem : 
   VE(M)= -5.2562e-06   VE(GB)= -2.34162e-06   VE(FUS)= 9.37339e-06

================================================================
	    PARAMETERS USED FOR THE COMPUTATION : 
================================================================
# Parameters for the binary equilibrium computation by coal
###################################################################
ini.d
0.6	fact_separ : Factor by which the initial separation is multiplied
1.491102233	   mbar_voulue[0] : Baryon mass required for star 1 [M_sol]
1.491102233	   mbar_voulue[1] : Baryon mass required for star 2 [M_sol]
###################################################################
400	mermax : Maximum number of steps in the main iteration 
0.5 	relax :	 Relaxation factor in the main iteration 
1	mermax_eqb : Maximum number of steps in Etoile_bin::equilibrium 
1	prompt : 1 if no pause during the computation
0	graph : 1 if graphical outputs during the computation 
1.E-7 	seuil : Threshold on the enthalpy relative change for ending the computation 
2	fmer_stop : Step interval between pauses in the main iteration 
5	fmer_save : Step interval between safeguards of the whole configuration 
4	mermax_poisson : Maximum number of steps in Map_et::poisson
1.5	relax_poisson :  Relaxation factor in Map_et::poisson
4 	mermax_potvit : Maximum number of steps in Map_radial::poisson_compact
0.5	relax_potvit :  Relaxation factor in Map_radial::poisson_compact
20	mer_masse : Step from which the baryon mass is forced to converge
0.4	aexp_masse : Exponent for the increase factor of the central enthalpy
8	fmer_udp_met : Step interval between metric updates
1	ind_rel_met : 1 if relaxation of the metric, 0 if not
0.65	relax_met : Relaxation factor of the metric (used only if ind_rel_met=1)
0.75	relax_omeg : Relaxation factor on Omega (orbital angular velocity)
0.7	fact_omeg_min : fact_omeg_min * omega = low bound in the omega search
1.3	fact_omeg_max : fact_omeg_max * omega = high bound in the omega search
0.	thres_adapt1 : threshold on dH/dr for the adaptation of the mapping in star 1
0.	thres_adapt2 : threshold on dH/dr for the adaptation of the mapping in star 2






================================================================
	    IDENTIFICATION OF THE CODE : 
================================================================

================================================================
	    CPU TIME and MEMORY infos : 
================================================================

