28char tslice_dirac_max_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Time_slice/tslice_dirac_max.C,v 1.26 2014/10/13 08:53:48 j_novak Exp $" ;
135#include "time_slice.h"
136#include "utilitaires.h"
175 tmp.set_etat_zero() ;
204 "Constructor of Tslice_dirac_max from file: the case of partial reading\n"
205 <<
" is not ready yet !"
248 :
Time_slice_conf(star.get_nn(), star.get_beta(), star.get_mp().flat_met_spher(),
249 0.*star.get_nn(), star.get_hh(), 0.*star.get_aa(),
257 tmp.std_spectral_base() ;
293 initialize_sources_copy() ;
301 A_hh_evol(
tin.A_hh_evol),
302 B_hh_evol(
tin.B_hh_evol),
303 source_A_hh_evol(
tin.source_A_hh_evol),
304 source_B_hh_evol(
tin.source_B_hh_evol),
305 source_A_hata_evol(
tin.source_A_hata_evol),
306 source_B_hata_evol(
tin.source_B_hata_evol),
363 int nz =
trk_in.get_mp().get_mg()->get_nzone() ;
371 if (
tmp.get_etat() != ETATZERO) {
372 assert(
tmp.get_mp().get_mg()->get_type_r(nz-1) == UNSURR) ;
373 tmp.annule_domain(nz-1) ;
380 if (
tmp.get_etat() != ETATZERO) {
381 assert(
tmp.get_mp().get_mg()->get_type_r(nz-1) == UNSURR) ;
382 tmp.annule_domain(nz-1) ;
389 "Tslice_dirac_max::initial_data_cts : variation of A and tilde(B) for J = "
404 initialize_sources_copy() ;
426 wzero.set_etat_zero() ;
440 cout <<
"Tslice_dirac_max::set_trh : #### WARNING : \n"
441 <<
" this method does not check whether det(tilde gamma) = 1"
496 resu.set_etat_zero() ;
570 flux <<
"Dirac gauge and maximal slicing" <<
'\n' ;
578 if (trh_evol.is_known(
jtime)) {
591 "Tslice_dirac_max::sauve : the partial_save case is not ready yet !"
Vectorial bases (triads) with respect to which the tensorial components are defined.
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.
void downdate(int j)
Suppresses a stored value.
bool is_known(int j) const
Checks whether the value a given time step has been set.
Base class for coordinate mappings.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Metric_flat & flat_met_spher() const
Returns the flat metric associated with the spherical coordinates and with components expressed in th...
Flat metric for tensor calculation.
Tensor field of valence 0 (or component of a tensorial field).
Class for relativistic rotating stars in Dirac gauge and maximal slicing.
const Sym_tensor get_aa() const
Returns .
const Scalar & get_ln_psi() const
Returns .
Transverse symmetric tensors of rank 2.
Transverse and traceless symmetric tensors of rank 2.
Class intended to describe valence-2 symmetric tensors.
Spacelike time slice of a 3+1 spacetime with conformal decomposition.
Evolution_std< Sym_tensor > hh_evol
Values at successive time steps of the components .
Evolution_std< Sym_tensor > hata_evol
Values at successive time steps of the components .
virtual const Scalar & A_hata() const
Returns the potential A of .
void operator=(const Time_slice_conf &)
Assignment to another Time_slice_conf.
virtual const Scalar & psi() const
Conformal factor relating the physical metric to the conformal one: .
Metric * p_tgamma
Pointer on the conformal metric at the current time step (jtime)
virtual const Scalar & B_hata() const
Returns the potential of .
virtual void del_deriv() const
Deletes all the derived quantities.
Vector * p_hdirac
Pointer on the vector (which vanishes in Dirac gauge), at the current time step (jtime).
const Scalar & psi4() const
Factor at the current time step (jtime ).
Evolution_std< Scalar > psi_evol
Values at successive time steps of the conformal factor relating the physical metric to the conform...
Evolution_std< Scalar > A_hata_evol
Potential A associated with the symmetric tensor .
Evolution_std< Scalar > B_hata_evol
Potential associated with the symmetric tensor .
virtual void set_hh(const Sym_tensor &hh_in)
Sets the deviation of the conformal metric from the flat metric : .
int jtime
Time step index of the latest slice.
Evolution_std< Scalar > trk_evol
Values at successive time steps of the trace K of the extrinsic curvature.
Metric * p_gamma
Pointer on the induced metric at the current time step (jtime)
Evolution_std< Sym_tensor > k_dd_evol
Values at successive time steps of the covariant components of the extrinsic curvature tensor .
virtual void sauve(FILE *fich, bool partial_save) const
Total or partial saves in a binary file.
Evolution_full< Tbl > adm_mass_evol
ADM mass at each time step, since the creation of the slice.
Evolution_std< Sym_tensor > gam_dd_evol
Values at successive time steps of the covariant components of the induced metric .
Evolution_std< Sym_tensor > gam_uu_evol
Values at successive time steps of the contravariant components of the induced metric .
int depth
Number of stored time slices.
Evolution_std< Scalar > n_evol
Values at successive time steps of the lapse function N.
Evolution_std< double > the_time
Time label of each slice.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator<<).
Evolution_std< Vector > beta_evol
Values at successive time steps of the shift vector .
Spacelike time slice of a 3+1 spacetime with conformal decomposition in the maximal slicing and Dirac...
virtual const Scalar & trh() const
Computes the trace h, with respect to the flat metric ff , of .
virtual const Sym_tensor & hh(Param *par_bc=0x0, Param *par_mat=0x0) const
Deviation of the conformal metric from the flat metric : .
Evolution_std< Scalar > source_B_hh_evol
The potential of the source of equation for .
Evolution_std< Scalar > source_A_hh_evol
The A potential of the source of equation for .
virtual const Scalar & B_hh() const
Returns the potential of .
virtual void initial_data_cts(const Sym_tensor &uu, const Scalar &trk_in, const Scalar &trk_point, double pdt, double precis=1.e-12, int method_poisson_vect=6, const char *graph_device=0x0, const Scalar *ener_dens=0x0, const Vector *mom_dens=0x0, const Scalar *trace_stress=0x0)
Computes valid initial data by solving the constraint equations in the conformal thin-sandwich approa...
Evolution_std< Scalar > A_hh_evol
The A potential of .
virtual ~Tslice_dirac_max()
Destructor.
virtual const Scalar & trk() const
Trace K of the extrinsic curvature at the current time step (jtime ).
void hh_det_one(int j, Param *par_bc=0x0, Param *par_mat=0x0) const
Computes from the values of A and and using the condition , which fixes the trace of .
Evolution_std< Scalar > source_B_hata_evol
The potential of the source of equation for .
virtual void set_hh(const Sym_tensor &hh_in)
Sets the deviation of the conformal metric from the flat metric : .
void compute_sources(const Sym_tensor *strain_tensor=0x0) const
Computes the sources source_A_XXX_evol and source_B_XXX_evol , for the solution of the evolution equa...
virtual void set_khi_mu(const Scalar &khi_in, const Scalar &mu_in)
Sets the potentials and of the TT part of (see the documentation of Sym_tensor_tt for details).
virtual void sauve(FILE *fich, bool partial_save) const
Total or partial saves in a binary file.
virtual void set_trh(const Scalar &trh_in)
Sets the trace, with respect to the flat metric ff , of .
virtual const Scalar & A_hh() const
Returns the potential A of .
virtual const Vector & hdirac() const
Vector which vanishes in Dirac gauge.
Tslice_dirac_max(const Scalar &lapse_in, const Vector &shift_in, const Metric_flat &ff_in, const Scalar &psi_in, const Sym_tensor_trans &hh_in, const Sym_tensor &hata_in, int depth_in=3)
Constructor from conformal decomposition.
Evolution_std< Scalar > B_hh_evol
The potential of .
Evolution_std< Scalar > source_A_hata_evol
The potential A of the source of equation for .
Tensor field of valence 1.
Cmp exp(const Cmp &)
Exponential.
int fread_be(int *aa, int size, int nb, FILE *fich)
Reads integer(s) from a binary file according to the big endian convention.
int fwrite_be(const int *aa, int size, int nb, FILE *fich)
Writes integer(s) into a binary file according to the big endian convention.
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Tbl maxabs(const Tensor &aa, const char *comment=0x0, ostream &ost=cout, bool verb=true)
Maxima in each domain of the absolute values of the tensor components.