88 #include "graphique.h" 89 #include "utilitaires.h" 108 if ( (comp.
logn_auto).get_etat() == ETATZERO ) {
133 if ( (comp.
lnq_auto).get_etat() == ETATZERO ) {
150 for(
int i=1; i<=3; i++)
151 for(
int j=i; j<=3; j++) {
186 for(
int i=1; i<=3; i++)
187 for(
int j=i; j<=3; j++) {
190 gtilde_con.set(i,j) =
hij(i,j) +
flat.
con()(i,j) ;
216 double* max_det =
new double[nz] ;
217 double* min_det =
new double[nz] ;
218 double* moy_det =
new double[nz] ;
220 for (
int i=0; i<nz; i++){
226 for (
int l=0; l<nz; l++)
227 for (
int k=0; k<np; k++)
228 for (
int j=0; j<nt; j++)
229 for (
int i=0; i<nr; i++){
232 if (det_gtilde.val_grid_point(l,k,j,i) > max_det[l]){
233 max_det[l] = det_gtilde.val_grid_point(l,k,j,i) ;
235 if (det_gtilde.val_grid_point(l,k,j,i) < min_det[l]){
236 min_det[l] = det_gtilde.val_grid_point(l,k,j,i) ;
240 cout <<
"average determinant of gtilde in each zone : " << endl ;
241 for (
int l=0; l<nz; l++){
242 cout << moy_det[l]/(nr*nt*np) <<
" " ;
247 cout <<
"maximum of the determinant of gtilde in each zone : " << endl ;
248 for (
int l=0; l<nz; l++){
249 cout << max_det[l] <<
" " ;
253 cout <<
"minimum of the determinant of gtilde in each zone : " << endl ;
254 for (
int l=0; l<nz; l++){
255 cout << min_det[l] <<
" " ;
281 double relax,
double om) {
294 if ( (comp.
logn_auto).get_etat() == ETATZERO ) {
320 if ( (comp.
lnq_auto).get_etat() == ETATZERO ) {
338 for(
int i=1; i<=3; i++)
339 for(
int j=i; j<=3; j++) {
349 double relaxjm1 = 1. - relax ;
359 for(
int i=1; i<=3; i++)
360 for(
int j=i; j<=3; j++) {
363 + relaxjm1 * (star_jm1.
hij_comp)(i,j) ;
395 for(
int i=1; i<=3; i++)
396 for(
int j=i; j<=3; j++) {
399 gtilde_con.set(i,j) =
hij(i,j) +
flat.
con()(i,j) ;
423 double* max_det =
new double[nz] ;
424 double* min_det =
new double[nz] ;
425 double* moy_det =
new double[nz] ;
427 for (
int i=0; i<nz; i++){
433 for (
int l=0; l<nz; l++)
434 for (
int k=0; k<np; k++)
435 for (
int j=0; j<nt; j++)
436 for (
int i=0; i<nr; i++){
439 if (det_gtilde.val_grid_point(l,k,j,i) > max_det[l]){
440 max_det[l] = det_gtilde.val_grid_point(l,k,j,i) ;
442 if (det_gtilde.val_grid_point(l,k,j,i) < min_det[l]){
443 min_det[l] = det_gtilde.val_grid_point(l,k,j,i) ;
447 cout <<
"average determinant of gtilde in each zone : " << endl ;
448 for (
int l=0; l<nz; l++){
449 cout << moy_det[l]/(nr*nt*np) <<
" " ;
454 cout <<
"maximum of the determinant of gtilde in each zone : " << endl ;
455 for (
int l=0; l<nz; l++){
456 cout << max_det[l] <<
" " ;
460 cout <<
"minimum of the determinant of gtilde in each zone : " << endl ;
461 for (
int l=0; l<nz; l++){
462 cout << min_det[l] <<
" " ;
Sym_tensor hij_comp
Deviation of the inverse conformal metric from the inverse flat metric generated principally by the ...
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Cmp exp(const Cmp &)
Exponential.
Map & mp
Mapping associated with the star.
Sym_tensor hij_auto
Deviation of the inverse conformal metric from the inverse flat metric generated principally by the ...
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
Base class for coordinate mappings.
Scalar logn_comp
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principally by...
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
const Map & get_mp() const
Returns the mapping.
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tensor field of valence 1.
Scalar logn_auto
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principally by...
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
Sym_tensor hij
Total deviation of the inverse conformal metric from the inverse flat metric.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the star) ...
Scalar lnq_auto
Scalar field generated principally by the star.
Scalar lnq_comp
Scalar field generated principally by the companion star.
int get_nzone() const
Returns the number of domains.
void extrinsic_curvature(double omega)
Computes tkij_auto and akcar_auto from beta_auto, nn and Q.
virtual void del_deriv() const
Deletes all the derived quantities.
Scalar logn
Logarithm of the lapse N .
Class for stars in binary system.
virtual const Scalar & determinant() const
Returns the determinant.
Metric gtilde
Conformal metric .
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
void update_metric(const Star_bin &comp, double omega)
Computes metric coefficients from known potentials, when the companion is another star...
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Scalar nn
Lapse function N .
Vector beta_comp
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
bool conf_flat
true if the 3-metric is conformally flat, false for a more general metric.
Vector beta_auto
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
virtual void change_triad(const Base_vect &new_triad)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Scalar & set(int)
Read/write access to a component.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Scalar psi4
Conformal factor .
Class intended to describe valence-2 symmetric tensors.