29char 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 $" ;
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++){
244 cout <<
"maximum of the determinant of gtilde in each zone : " <<
endl ;
245 for (
int l=0;
l<nz;
l++){
250 cout <<
"minimum of the determinant of gtilde in each zone : " <<
endl ;
251 for (
int l=0;
l<nz;
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++) {
356 for(
int i=1;
i<=3;
i++)
357 for(
int j=
i;
j<=3;
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++){
451 cout <<
"maximum of the determinant of gtilde in each zone : " <<
endl ;
452 for (
int l=0;
l<nz;
l++){
457 cout <<
"minimum of the determinant of gtilde in each zone : " <<
endl ;
458 for (
int l=0;
l<nz;
l++){
Time evolution with partial storage (*** under development ***).
Base class for coordinate mappings.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
virtual const Scalar & determinant() const
Returns the determinant.
Tensor field of valence 0 (or component of a tensorial field).
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
Class for stars in binary system.
Scalar lnq_auto
Scalar field generated principally by the star.
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.
Vector beta_comp
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
void update_metric(const Star_bin &comp, double omega)
Computes metric coefficients from known potentials, when the companion is another star.
Scalar logn_auto
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principally by...
Scalar logn_comp
Part of the lapse logarithm (gravitational potential at the Newtonian limit) generated principally by...
void extrinsic_curvature(double omega)
Computes tkij_auto and akcar_auto from beta_auto, nn and Q.
Scalar lnq_comp
Scalar field generated principally by the companion star.
Vector beta_auto
Part of the shift vector generated principally by the star (Spherical components with respect to the ...
Sym_tensor hij
Total deviation of the inverse conformal metric from the inverse flat metric.
Sym_tensor hij_auto
Deviation of the inverse conformal metric from the inverse flat metric generated principally by the ...
Scalar psi4
Conformal factor .
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the star) .
Metric gtilde
Conformal metric .
bool conf_flat
true if the 3-metric is conformally flat, false for a more general metric.
Scalar logn
Logarithm of the lapse N .
Scalar nn
Lapse function N .
Map & mp
Mapping associated with the star.
Class intended to describe valence-2 symmetric tensors.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Scalar & set(int)
Read/write access to a component.
Cmp exp(const Cmp &)
Exponential.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.