70 #include "time_slice.h" 90 for (
int l=0; l<nz; l++) {
91 double radius = mp.
val_r(l, 1., 0., 0.) ;
92 tmass->
set(l) = ww.
flux(radius, ff) / (16.* M_PI) ;
117 Tbl* tmass =
new Tbl(nz) ;
123 for (
int l=0; l<nz; l++) {
124 double radius = mp.
val_r(l, 1., 0., 0.) ;
125 tmass->
set(l) = - ww.
flux(radius,
ff) / (2.* M_PI) ;
152 Tbl* tmass =
new Tbl(nz) ;
159 for (
int l=0; l<nz; l++) {
160 double radius = mp.
val_r(l, 1., 0., 0.) ;
161 tmass->
set(l) = - ww.
flux(radius,
ff) / (2.* M_PI) ;
virtual const Scalar & psi() const
Conformal factor relating the physical metric to the conformal one: .
double flux(double radius, const Metric &met) const
Computes the flux of the vector accross a sphere r = const.
virtual double adm_mass() const
Returns the ADM mass at (geometrical units) the current step.
virtual const Vector & hdirac() const
Vector which vanishes in Dirac gauge.
Evolution_full< Tbl > adm_mass_evol
ADM mass at each time step, since the creation of the slice.
Tensor up(int ind, const Metric &gam) const
Computes a new tensor by raising an index of *this.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
double & set(int i)
Read/write of a particular element (index i) (1D case)
Flat metric for tensor calculation.
Base class for coordinate mappings.
const Tensor_sym & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
int jtime
Time step index of the latest slice.
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
Tensor field of valence 1.
virtual const Sym_tensor & gam_dd() const
Induced metric (covariant components ) at the current time step (jtime )
virtual const Sym_tensor & hh(Param *=0x0, Param *=0x0) const
Deviation of the conformal metric from the flat metric : .
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
virtual double val_r(int l, double xi, double theta, double pphi) const =0
Returns the value of the radial coordinate r for a given in a given domain.
int get_nzone() const
Returns the number of domains.
virtual double adm_mass() const
Returns the ADM mass (geometrical units) at the current step.
Evolution_std< double > the_time
Time label of each slice.
const Tensor & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
const Metric_flat & ff
Pointer on the flat metric with respect to which the conformal decomposition is performed.
int get_taille() const
Gives the total size (ie dim.taille)
virtual double adm_mass() const
Returns the ADM mass (geometrical units) at the current step.
virtual const Sym_tensor & hh(Param *par_bc=0x0, Param *par_mat=0x0) const
Deviation of the conformal metric from the flat metric : .
const Map & get_mp() const
Returns the mapping.
Tensor trace(int ind1, int ind2) const
Trace on two different type indices.