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) ;
112 const Map& mp = psi().get_mp() ;
114 Tbl* tmass =
new Tbl(nz) ;
118 + 0.125* ( hdirac() - (hh().trace(ff)).derive_con(ff) ) ;
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) ;
147 const Map& mp = psi().get_mp() ;
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) ;
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.
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.
Tensor field of valence 1.
virtual const Sym_tensor & gam_dd() const
Induced metric (covariant components ) at the current time step (jtime )
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...
int get_taille() const
Gives the total size (ie dim.taille)
const Map & get_mp() const
Returns the mapping.
Tensor trace(int ind1, int ind2) const
Trace on two different type indices.