29 char star_binupmetr_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/star_bin_upmetr.C,v 1.14 2014/10/13 08:53:38 j_novak Exp $" ;
85 #include "graphique.h"
86 #include "utilitaires.h"
105 if ( (comp.
logn_auto).get_etat() == ETATZERO ) {
130 if ( (comp.
lnq_auto).get_etat() == ETATZERO ) {
147 for(
int i=1; i<=3; i++)
148 for(
int j=i; j<=3; j++) {
183 for(
int i=1; i<=3; i++)
184 for(
int j=i; j<=3; j++) {
213 double* max_det =
new double[nz] ;
214 double* min_det =
new double[nz] ;
215 double* moy_det =
new double[nz] ;
217 for (
int i=0; i<nz; i++){
223 for (
int l=0; l<nz; l++)
224 for (
int k=0; k<np; k++)
225 for (
int j=0; j<nt; j++)
226 for (
int i=0; i<nr; i++){
237 cout <<
"average determinant of gtilde in each zone : " << endl ;
238 for (
int l=0; l<nz; l++){
239 cout << moy_det[l]/(nr*nt*np) <<
" " ;
244 cout <<
"maximum of the determinant of gtilde in each zone : " << endl ;
245 for (
int l=0; l<nz; l++){
246 cout << max_det[l] <<
" " ;
250 cout <<
"minimum of the determinant of gtilde in each zone : " << endl ;
251 for (
int l=0; l<nz; l++){
252 cout << min_det[l] <<
" " ;
278 double relax,
double om) {
291 if ( (comp.
logn_auto).get_etat() == ETATZERO ) {
317 if ( (comp.
lnq_auto).get_etat() == ETATZERO ) {
335 for(
int i=1; i<=3; i++)
336 for(
int j=i; j<=3; j++) {
346 double relaxjm1 = 1. - relax ;
356 for(
int i=1; i<=3; i++)
357 for(
int j=i; j<=3; j++) {
360 + relaxjm1 * (star_jm1.
hij_comp)(i,j) ;
392 for(
int i=1; i<=3; i++)
393 for(
int j=i; j<=3; j++) {
420 double* max_det =
new double[nz] ;
421 double* min_det =
new double[nz] ;
422 double* moy_det =
new double[nz] ;
424 for (
int i=0; i<nz; i++){
430 for (
int l=0; l<nz; l++)
431 for (
int k=0; k<np; k++)
432 for (
int j=0; j<nt; j++)
433 for (
int i=0; i<nr; i++){
444 cout <<
"average determinant of gtilde in each zone : " << endl ;
445 for (
int l=0; l<nz; l++){
446 cout << moy_det[l]/(nr*nt*np) <<
" " ;
451 cout <<
"maximum of the determinant of gtilde in each zone : " << endl ;
452 for (
int l=0; l<nz; l++){
453 cout << max_det[l] <<
" " ;
457 cout <<
"minimum of the determinant of gtilde in each zone : " << endl ;
458 for (
int l=0; l<nz; l++){
459 cout << min_det[l] <<
" " ;
Scalar logn
Logarithm of the lapse N .
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
Scalar psi4
Conformal factor .
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
Vector beta_auto
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
Class intended to describe valence-2 symmetric tensors.
void update_metric(const Star_bin &comp, double omega)
Computes metric coefficients from known potentials, when the companion is another star.
virtual const Scalar & determinant() const
Returns the determinant.
Cmp exp(const Cmp &)
Exponential.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
bool conf_flat
true if the 3-metric is conformally flat, false for a more general metric.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Class for stars in binary system.
Vector beta_comp
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
Metric gtilde
Conformal metric .
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
Tensor field of valence 0 (or component of a tensorial field).
Scalar nn
Lapse function N .
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Scalar logn_auto
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principally by...
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the star) .
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Sym_tensor hij
Total deviation of the inverse conformal metric from the inverse flat metric.
Scalar lnq_auto
Scalar field generated principally by the star.
void extrinsic_curvature(double omega)
Computes tkij_auto and akcar_auto from beta_auto, nn and Q.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Sym_tensor hij_comp
Deviation of the inverse conformal metric from the inverse flat metric generated principally by the ...
virtual void del_deriv() const
Deletes all the derived quantities.
Map & mp
Mapping associated with the star.
Scalar & set(int)
Read/write access to a component.
virtual void change_triad(const Base_vect &new_triad)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Scalar lnq_comp
Scalar field generated principally by the companion star.
int get_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Tensor field of valence 1.
Scalar logn_comp
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principally by...
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
Sym_tensor hij_auto
Deviation of the inverse conformal metric from the inverse flat metric generated principally by the ...
const Map & get_mp() const
Returns the mapping.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Base class for coordinate mappings.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.