29char 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"
80 int nz = mp.
get_mg()->get_nzone() ;
82 tmass->set_etat_qcq() ;
85 -
gam_dd().trace(ff).derive_con(ff) ;
87 for (
int l=0;
l<nz;
l++) {
88 double radius = mp.
val_r(
l, 1., 0., 0.) ;
112 const Map& mp =
psi().get_mp() ;
113 int nz = mp.
get_mg()->get_nzone() ;
115 tmass->set_etat_qcq() ;
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.) ;
135 return tadm(tadm.get_taille()-1) ;
147 const Map& mp =
psi().get_mp() ;
148 int nz = mp.
get_mg()->get_nzone() ;
150 tmass->set_etat_qcq() ;
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.) ;
Time evolution with partial storage (*** under development ***).
virtual void update(const TyT &new_value, int j, double time_j)
Sets a new value at a given time step.
Base class for coordinate mappings.
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.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Flat metric for tensor calculation.
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.