29 char tslice_adm_mass_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Time_slice/tslice_adm_mass.C,v 1.7 2014/10/13 08:53:47 j_novak Exp $" ;
67 #include "time_slice.h"
87 for (
int l=0; l<nz; l++) {
88 double radius = mp.
val_r(l, 1., 0., 0.) ;
89 tmass->
set(l) = ww.
flux(radius, ff) / (16.* M_PI) ;
114 Tbl* tmass =
new Tbl(nz) ;
120 for (
int l=0; l<nz; l++) {
121 double radius = mp.
val_r(l, 1., 0., 0.) ;
122 tmass->
set(l) = - ww.
flux(radius,
ff) / (2.* M_PI) ;
135 return tadm(tadm.get_taille()-1) ;
149 Tbl* tmass =
new Tbl(nz) ;
153 - 0.125* (
hh().
trace(ff)).derive_con(ff) ;
156 for (
int l=0; l<nz; l++) {
157 double radius = mp.
val_r(l, 1., 0., 0.) ;
158 tmass->
set(l) = - ww.
flux(radius, ff) / (2.* M_PI) ;
Base class for coordinate mappings.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
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.
Flat metric for tensor calculation.
int get_nzone() const
Returns the number of domains.
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
double & set(int i)
Read/write of a particular element (index i) (1D case)
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
int get_taille() const
Gives the total size (ie dim.taille)
virtual const Sym_tensor & hh(Param *=0x0, Param *=0x0) const
Deviation of the conformal metric from the flat metric : .
virtual const Vector & hdirac() const
Vector which vanishes in Dirac gauge.
virtual const Scalar & psi() const
Conformal factor relating the physical metric to the conformal one: .
const Metric_flat & ff
Pointer on the flat metric with respect to which the conformal decomposition is performed.
int jtime
Time step index of the latest slice.
Evolution_full< Tbl > adm_mass_evol
ADM mass at each time step, since the creation of the slice.
virtual const Sym_tensor & gam_dd() const
Induced metric (covariant components ) at the current time step (jtime )
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.
virtual double adm_mass() const
Returns the ADM mass at (geometrical units) 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 : .
Tensor field of valence 1.
double flux(double radius, const Metric &met) const
Computes the flux of the vector accross a sphere r = const.
const Map & get_mp() const
Returns the mapping.
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 Tensor_sym & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of this with respect to some metric , by raising the last inde...
Tensor up(int ind, const Metric &gam) const
Computes a new tensor by raising an index of *this.
Tensor trace(int ind1, int ind2) const
Trace on two different type indices.