LORENE
|
Spacelike time slice of a 3+1 spacetime with conformal decomposition in the maximal slicing and Dirac gauge. More...
#include <time_slice.h>
Public Member Functions | |
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. More... | |
Tslice_dirac_max (const Map &mp, const Base_vect &triad, const Metric_flat &ff_in, int depth_in=3) | |
Constructor as standard time slice of flat spacetime (Minkowski). More... | |
Tslice_dirac_max (const Map &mp, const Base_vect &triad, const Metric_flat &ff_in, FILE *fich, bool partial_read, int depth_in=3) | |
Constructor from binary file. More... | |
Tslice_dirac_max (const Star_rot_Dirac &star, double pdt, int depth_in=3) | |
Construnction of a stationary slice from a rotating star. More... | |
Tslice_dirac_max (const Tslice_dirac_max &) | |
Copy constructor. More... | |
virtual | ~Tslice_dirac_max () |
Destructor. More... | |
void | operator= (const Tslice_dirac_max &) |
Assignment to another Tslice_dirac_max. More... | |
virtual void | set_hh (const Sym_tensor &hh_in) |
Sets the deviation of the conformal metric from the flat metric : . More... | |
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 approach. More... | |
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). More... | |
virtual void | set_AB_hh (const Scalar &A_in, const Scalar &B_in) |
Sets the potentials A and of the TT part of (see the documentation of Sym_tensor for details). More... | |
virtual void | set_trh (const Scalar &trh_in) |
Sets the trace, with respect to the flat metric ff , of . More... | |
virtual Scalar | solve_psi (const Scalar *ener_dens=0x0) const |
Solves the elliptic equation for the conformal factor $$ (Hamiltonian constraint). More... | |
virtual Scalar | solve_npsi (const Scalar *ener_dens=0x0, const Scalar *trace_stress=0x0) const |
Solves the elliptic equation for (maximal slicing condition + Hamiltonian constraint) More... | |
virtual Vector | solve_beta (int method=6) const |
Solves the elliptic equation for the shift vector from (Eq. More... | |
void | evolve (double pdt, int nb_time_steps, int niter_elliptic, double relax_elliptic, int check_mod, int save_mod, int method_poisson_vect=6, int nopause=1, const char *graph_device=0x0, bool verbose=true, const Scalar *ener_euler=0x0, const Vector *mom_euler=0x0, const Scalar *s_euler=0x0, const Sym_tensor *strain_euler=0x0) |
Time evolution by resolution of Einstein equations. More... | |
virtual double | adm_mass () const |
Returns the ADM mass at (geometrical units) the current step. More... | |
virtual const Sym_tensor & | hh (Param *par_bc=0x0, Param *par_mat=0x0) const |
Deviation of the conformal metric from the flat metric : . More... | |
virtual const Scalar & | trk () const |
Trace K of the extrinsic curvature at the current time step (jtime ). More... | |
virtual const Vector & | hdirac () const |
Vector which vanishes in Dirac gauge. More... | |
virtual const Scalar & | A_hh () const |
Returns the potential A of . More... | |
virtual const Scalar & | B_hh () const |
Returns the potential of . More... | |
virtual const Scalar & | trh () const |
Computes the trace h , with respect to the flat metric ff , of . More... | |
virtual void | set_psi_del_npsi (const Scalar &psi_in) |
Sets the conformal factor relating the physical metric to the conformal one: . More... | |
virtual void | set_psi_del_n (const Scalar &psi_in) |
Sets the conformal factor relating the physical metric to the conformal one: . More... | |
virtual void | set_npsi_del_psi (const Scalar &npsi_in) |
Sets the factor at the current time step (jtime ) and deletes the value of . More... | |
virtual void | set_npsi_del_n (const Scalar &npsi_in) |
Sets the factor at the current time step (jtime ) and deletes the value of N. More... | |
virtual void | set_hata (const Sym_tensor &hata_in) |
Sets the conformal representation of the traceless part of the extrinsic curvature: . More... | |
virtual void | set_hata_TT (const Sym_tensor_tt &hata_tt) |
Sets the TT part of (see member hata_evol ). More... | |
virtual void | set_hata_from_XAB (Param *par_bc=0x0, Param *par_mat=0x0) |
Sets the conformal representation of the traceless part of the extrinsic curvature from its potentials A, and . More... | |
virtual const Scalar & | nn () const |
Lapse function N at the current time step (jtime ) More... | |
virtual const Sym_tensor & | gam_dd () const |
Induced metric (covariant components ) at the current time step (jtime ) More... | |
virtual const Sym_tensor & | gam_uu () const |
Induced metric (contravariant components ) at the current time step (jtime ) More... | |
virtual const Sym_tensor & | k_dd () const |
Extrinsic curvature tensor (covariant components ) at the current time step (jtime ) More... | |
virtual const Sym_tensor & | k_uu () const |
Extrinsic curvature tensor (contravariant components ) at the current time step (jtime ) More... | |
virtual const Scalar & | A_hata () const |
Returns the potential A of . More... | |
virtual const Scalar & | B_hata () const |
Returns the potential of . More... | |
virtual const Scalar & | psi () const |
Conformal factor relating the physical metric to the conformal one: . More... | |
const Scalar & | psi4 () const |
Factor at the current time step (jtime ). More... | |
const Scalar & | ln_psi () const |
Logarithm of at the current time step (jtime ). More... | |
virtual const Scalar & | npsi () const |
Factor at the current time step (jtime ). More... | |
virtual const Metric & | tgam () const |
Conformal metric Returns the value at the current time step (jtime ). More... | |
virtual const Sym_tensor & | hata () const |
Conformal representation of the traceless part of the extrinsic curvature: . More... | |
virtual Sym_tensor | aa () const |
Conformal representation of the traceless part of the extrinsic curvature: . More... | |
virtual const Vector & | vec_X (int method_poisson=6) const |
Vector representing the longitudinal part of . More... | |
void | compute_X_from_momentum_constraint (const Vector &hat_S, const Sym_tensor_tt &hata_tt, int iter_max=200, double precis=1.e-12, double relax=0.8, int methode_poisson=6) |
Computes the vector from the conformally-rescaled momentum , using the momentum constraint. More... | |
virtual void | set_AB_hata (const Scalar &A_in, const Scalar &B_in) |
Sets the potentials A and of the TT part (see the documentation of Sym_tensor for details). More... | |
void | check_psi_dot (Tbl &tlnpsi_dot, Tbl &tdiff, Tbl &tdiff_rel) const |
Checks the relation. More... | |
void | set_scheme_order (int ord) |
Sets the order of the finite-differences scheme. More... | |
int | get_scheme_order () const |
Gets the order of the finite-differences scheme. More... | |
int | get_latest_j () const |
Gets the latest value of time step index. More... | |
const Evolution_std< double > & | get_time () const |
Gets the time coordinate t at successive time steps. More... | |
virtual const Vector & | beta () const |
shift vector at the current time step (jtime ) More... | |
const Metric & | gam () const |
Induced metric at the current time step (jtime ) More... | |
Tbl | check_hamiltonian_constraint (const Scalar *energy_density=0x0, ostream &ost=cout, bool verb=true) const |
Checks the level at which the hamiltonian constraint is verified. More... | |
Tbl | check_momentum_constraint (const Vector *momentum_density=0x0, ostream &ost=cout, bool verb=true) const |
Checks the level at which the momentum constraints are verified. More... | |
Tbl | check_dynamical_equations (const Sym_tensor *strain_tensor=0x0, const Scalar *energy_density=0x0, ostream &ost=cout, bool verb=true) const |
Checks the level at which the dynamical equations are verified. More... | |
void | save (const char *rootname) const |
Saves in a binary file. More... | |
Protected Member Functions | |
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 equation for and . More... | |
void | initialize_sources_copy () const |
Copy the sources source_A_XXX_evol and source_B_XXX_evol to all time-steps. More... | |
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 . More... | |
void | hh_det_one (const Sym_tensor_tt &hijtt, Param *par_mat=0x0) const |
Computes from the TT part using the condition , which fixes the trace of . More... | |
virtual ostream & | operator>> (ostream &) const |
Operator >> (virtual function called by the operator<<). More... | |
virtual void | sauve (FILE *fich, bool partial_save) const |
Total or partial saves in a binary file. More... | |
virtual void | del_deriv () const |
Deletes all the derived quantities. More... | |
void | set_der_0x0 () const |
Sets to 0x0 all the pointers on derived quantities. More... | |
Protected Attributes | |
Evolution_std< Scalar > | A_hh_evol |
The A potential of . More... | |
Evolution_std< Scalar > | B_hh_evol |
The potential of . More... | |
Evolution_std< Scalar > | source_A_hh_evol |
The A potential of the source of equation for . More... | |
Evolution_std< Scalar > | source_B_hh_evol |
The potential of the source of equation for . More... | |
Evolution_std< Scalar > | source_A_hata_evol |
The potential A of the source of equation for . More... | |
Evolution_std< Scalar > | source_B_hata_evol |
The potential of the source of equation for . More... | |
Evolution_std< Scalar > | trh_evol |
The trace, with respect to the flat metric ff , of . More... | |
const Metric_flat & | ff |
Pointer on the flat metric with respect to which the conformal decomposition is performed. More... | |
Evolution_std< Scalar > | psi_evol |
Values at successive time steps of the conformal factor relating the physical metric to the conformal one: . More... | |
Evolution_std< Scalar > | npsi_evol |
Values at successive time steps of the factor . More... | |
Evolution_std< Sym_tensor > | hh_evol |
Values at successive time steps of the components . More... | |
Evolution_std< Sym_tensor > | hata_evol |
Values at successive time steps of the components . More... | |
Evolution_std< Scalar > | A_hata_evol |
Potential A associated with the symmetric tensor . More... | |
Evolution_std< Scalar > | B_hata_evol |
Potential associated with the symmetric tensor . More... | |
Metric * | p_tgamma |
Pointer on the conformal metric at the current time step (jtime ) More... | |
Scalar * | p_psi4 |
Pointer on the factor at the current time step (jtime ) More... | |
Scalar * | p_ln_psi |
Pointer on the logarithm of at the current time step (jtime ) More... | |
Vector * | p_hdirac |
Pointer on the vector (which vanishes in Dirac gauge), at the current time step (jtime ). More... | |
Vector * | p_vec_X |
Pointer on the vector representing the longitudinal part of . More... | |
int | depth |
Number of stored time slices. More... | |
int | scheme_order |
Order of the finite-differences scheme for the computation of time derivatives. More... | |
int | jtime |
Time step index of the latest slice. More... | |
Evolution_std< double > | the_time |
Time label of each slice. More... | |
Evolution_std< Sym_tensor > | gam_dd_evol |
Values at successive time steps of the covariant components of the induced metric . More... | |
Evolution_std< Sym_tensor > | gam_uu_evol |
Values at successive time steps of the contravariant components of the induced metric . More... | |
Evolution_std< Sym_tensor > | k_dd_evol |
Values at successive time steps of the covariant components of the extrinsic curvature tensor . More... | |
Evolution_std< Sym_tensor > | k_uu_evol |
Values at successive time steps of the contravariant components of the extrinsic curvature tensor . More... | |
Evolution_std< Scalar > | n_evol |
Values at successive time steps of the lapse function N. More... | |
Evolution_std< Vector > | beta_evol |
Values at successive time steps of the shift vector . More... | |
Evolution_std< Scalar > | trk_evol |
Values at successive time steps of the trace K of the extrinsic curvature. More... | |
Evolution_full< Tbl > | adm_mass_evol |
ADM mass at each time step, since the creation of the slice. More... | |
Metric * | p_gamma |
Pointer on the induced metric at the current time step (jtime ) More... | |
Spacelike time slice of a 3+1 spacetime with conformal decomposition in the maximal slicing and Dirac gauge.
()
Definition at line 971 of file time_slice.h.
Lorene::Tslice_dirac_max::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.
lapse_in | lapse function N |
shift_in | shift vector |
ff_in | reference flat metric with respect to which the conformal decomposition is performed |
psi_in | conformal factor relating the physical metric to the conformal one: |
hh_in | deviation of the conformal metric from the flat metric : . must be such that . |
hata_in | conformal representation of the traceless part of the extrinsic curvature: , with K = 0 in the present case |
depth_in | number of stored time slices; this parameter is used to set the scheme_order member with scheme_order = depth_in - 1. scheme_order can be changed afterwards by the method set_scheme_order(int) . |
Definition at line 152 of file tslice_dirac_max.C.
Lorene::Tslice_dirac_max::Tslice_dirac_max | ( | const Map & | mp, |
const Base_vect & | triad, | ||
const Metric_flat & | ff_in, | ||
int | depth_in = 3 |
||
) |
Constructor as standard time slice of flat spacetime (Minkowski).
mp | Mapping on which the various Lorene fields will be constructed |
triad | vector basis with respect to which the various tensor components will be defined |
ff_in | reference flat metric with respect to which the conformal decomposition is performed |
depth_in | number of stored time slices; this parameter is used to set the scheme_order member with scheme_order = depth_in - 1. scheme_order can be changed afterwards by the method set_scheme_order(int) . |
Definition at line 166 of file tslice_dirac_max.C.
References A_hh_evol, B_hh_evol, Lorene::Time_slice::jtime, Lorene::Scalar::set_etat_zero(), source_A_hata_evol, source_A_hh_evol, source_B_hata_evol, source_B_hh_evol, Lorene::Time_slice::the_time, and trh_evol.
Lorene::Tslice_dirac_max::Tslice_dirac_max | ( | const Map & | mp, |
const Base_vect & | triad, | ||
const Metric_flat & | ff_in, | ||
FILE * | fich, | ||
bool | partial_read, | ||
int | depth_in = 3 |
||
) |
Constructor from binary file.
The binary file must have been created by method save
.
mp | Mapping on which the various Lorene fields will be constructed |
triad | vector basis with respect to which the various tensor components will be defined |
ff_in | reference flat metric with respect to which the conformal decomposition is performed |
fich | file containing the saved Tslice_dirac_max |
partial_read | indicates whether the full object must be read in file or whether the final construction is devoted to a constructor of a derived class |
depth_in | number of stored time slices; the given must coincide with that stored in the file. |
Definition at line 196 of file tslice_dirac_max.C.
References A_hh_evol, B_hh_evol, Lorene::Time_slice::depth, Lorene::fread_be(), Lorene::Map::get_mg(), Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, and Lorene::Time_slice::the_time.
Lorene::Tslice_dirac_max::Tslice_dirac_max | ( | const Star_rot_Dirac & | star, |
double | pdt, | ||
int | depth_in = 3 |
||
) |
Construnction of a stationary slice from a rotating star.
Definition at line 250 of file tslice_dirac_max.C.
References Lorene::Time_slice_conf::A_hata(), Lorene::Time_slice_conf::A_hata_evol, A_hh(), A_hh_evol, Lorene::Time_slice_conf::B_hata(), Lorene::Time_slice_conf::B_hata_evol, B_hh(), B_hh_evol, Lorene::Time_slice::beta_evol, compute_sources(), Lorene::Time_slice::depth, Lorene::exp(), Lorene::Star_rot_Dirac::get_aa(), Lorene::Star_rot_Dirac::get_ln_psi(), Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice_conf::hh_evol, initialize_sources_copy(), Lorene::Time_slice::jtime, Lorene::Time_slice::k_dd_evol, Lorene::Time_slice::n_evol, Lorene::Time_slice_conf::psi(), Lorene::Time_slice_conf::psi4(), Lorene::Time_slice_conf::psi_evol, Lorene::Scalar::std_spectral_base(), Lorene::Time_slice::the_time, trh_evol, Lorene::Time_slice::trk_evol, and Lorene::Evolution_std< TyT >::update().
Lorene::Tslice_dirac_max::Tslice_dirac_max | ( | const Tslice_dirac_max & | tin | ) |
Copy constructor.
Definition at line 302 of file tslice_dirac_max.C.
|
virtual |
Destructor.
Definition at line 317 of file tslice_dirac_max.C.
|
virtualinherited |
Returns the potential A of .
See the documentation of Sym_tensor
for details. Returns the value at the current time step (jtime
).
Definition at line 667 of file time_slice_conf.C.
References Lorene::Time_slice_conf::A_hata_evol, Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice::jtime, and Lorene::Time_slice::the_time.
|
virtual |
Returns the potential A of .
See the documentation of Sym_tensor
for details. Returns the value at the current time step (jtime
).
Definition at line 528 of file tslice_dirac_max.C.
References A_hh_evol, Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, and Lorene::Time_slice::the_time.
|
virtualinherited |
Conformal representation of the traceless part of the extrinsic curvature: .
Returns the value at the current time step (jtime
).
Definition at line 768 of file time_slice_conf.C.
References Lorene::Time_slice_conf::hata(), Lorene::Time_slice_conf::psi(), and Lorene::Time_slice_conf::psi4().
|
virtual |
Returns the ADM mass at (geometrical units) the current step.
Moreover this method updates adm_mass_evol
if necessary.
Reimplemented from Lorene::Time_slice_conf.
Definition at line 146 of file tslice_adm_mass.C.
References Lorene::Time_slice::adm_mass_evol, Lorene::Scalar::derive_con(), Lorene::Time_slice_conf::ff, Lorene::Vector::flux(), Lorene::Map::get_mg(), Lorene::Tensor::get_mp(), Lorene::Mg3d::get_nzone(), Lorene::Tbl::get_taille(), hh(), Lorene::Time_slice::jtime, Lorene::Time_slice_conf::psi(), Lorene::Tbl::set(), Lorene::Tbl::set_etat_qcq(), Lorene::Time_slice::the_time, Lorene::Tensor::trace(), and Lorene::Map::val_r().
|
virtualinherited |
Returns the potential of .
See the documentation of Sym_tensor_tt
for details. Returns the value at the current time step (jtime
).
Definition at line 681 of file time_slice_conf.C.
References Lorene::Time_slice_conf::A_hata_evol, Lorene::Time_slice_conf::B_hata_evol, Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice::jtime, and Lorene::Time_slice::the_time.
|
virtual |
Returns the potential of .
See the documentation of Sym_tensor_tt
for details. Returns the value at the current time step (jtime
).
Definition at line 539 of file tslice_dirac_max.C.
References B_hh_evol, Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, and Lorene::Time_slice::the_time.
|
virtualinherited |
shift vector at the current time step (jtime
)
Definition at line 90 of file time_slice_access.C.
References Lorene::Time_slice::beta_evol, and Lorene::Time_slice::jtime.
|
inherited |
Checks the level at which the dynamical equations are verified.
strain_tensor | : a pointer on the strain_tensor measured by the Eulerian observer of 4-velocity ; if this is the null pointer, it is assumed that = 0 (vacuum). |
energy_density | : a pointer on the energy density E (see check_hamiltonian_constraint ) |
ost | : output stream for a formatted output of the result |
Definition at line 142 of file tslice_check_einstein.C.
References Lorene::Tensor::annule_domain(), Lorene::Tbl::annule_hard(), Lorene::Time_slice::beta(), Lorene::Metric::con(), Lorene::contract(), Lorene::Tensor::derive_con(), Lorene::Scalar::derive_con(), Lorene::Sym_tensor::derive_lie(), Lorene::Time_slice::gam(), Lorene::Tbl::get_dim(), Lorene::Map::get_mg(), Lorene::Tensor::get_mp(), Lorene::Mg3d::get_nzone(), Lorene::Time_slice::jtime, Lorene::Time_slice::k_dd(), Lorene::Time_slice::k_dd_evol, Lorene::Time_slice::k_uu(), Lorene::maxabs(), Lorene::Time_slice::nn(), Lorene::Tensor_sym::position(), Lorene::Metric::ricci(), Lorene::Time_slice::scheme_order, Lorene::Itbl::set(), Lorene::Tensor::trace(), Lorene::Time_slice::trk(), and Lorene::Tensor::up_down().
|
inherited |
Checks the level at which the hamiltonian constraint is verified.
energy_density | : a pointer on the energy density E measured by the Eulerian observer of 4-velocity ; if this is the null pointer, it is assumed that E = 0 (vacuum). |
ost | : output stream for a formatted output of the result |
Definition at line 82 of file tslice_check_einstein.C.
References Lorene::contract(), Lorene::Scalar::dec_dzpuis(), Lorene::Time_slice::gam(), Lorene::Time_slice::k_dd(), Lorene::Time_slice::k_uu(), Lorene::maxabs(), Lorene::Metric::ricci_scal(), and Lorene::Time_slice::trk().
|
inherited |
Checks the level at which the momentum constraints are verified.
momentum_density | : a pointer on the momentum density measured by the Eulerian observer of 4-velocity ; if this is the null pointer, it is assumed that = 0 (vacuum). |
ost | : output stream for a formatted output of the result |
Definition at line 112 of file tslice_check_einstein.C.
References Lorene::Scalar::derive_con(), Lorene::Sym_tensor::divergence(), Lorene::Time_slice::gam(), Lorene::Tensor::get_index_type(), Lorene::Time_slice::k_uu(), Lorene::maxabs(), and Lorene::Time_slice::trk().
|
inherited |
Checks the relation.
tlnpsi_dot | [output] maximun value in each domain of |
tdiff | [output] maximum value in each domain of |
tdiff_rel | [output] relative error on the above relation in each domain. |
Definition at line 907 of file time_slice_conf.C.
References Lorene::abs(), Lorene::Time_slice::beta(), Lorene::contract(), Lorene::Scalar::dec_dzpuis(), Lorene::Vector::divergence(), Lorene::Time_slice_conf::ff, Lorene::Time_slice::jtime, Lorene::Time_slice_conf::ln_psi(), Lorene::max(), Lorene::Time_slice::n_evol, Lorene::Time_slice_conf::nn(), Lorene::Time_slice_conf::npsi(), Lorene::Time_slice_conf::npsi_evol, Lorene::Time_slice_conf::psi(), Lorene::Time_slice_conf::psi_evol, Lorene::Time_slice::scheme_order, and Lorene::Time_slice_conf::trk().
|
protected |
Computes the sources source_A_XXX_evol
and source_B_XXX_evol
, for the solution of the evolution equation for and .
strain_tensor | [input] : a pointer on the strain_tensor measured by the Eulerian observer of 4-velocity ; if this is the null pointer (default value), it is assumed that = 0 (vacuum). |
Definition at line 222 of file tslice_dirac_max_setAB.C.
References Lorene::Time_slice_conf::A_hata_evol, A_hh_evol, Lorene::Time_slice_conf::aa(), Lorene::Tensor::annule_domain(), Lorene::Time_slice_conf::B_hata_evol, B_hh_evol, Lorene::Time_slice::beta(), Lorene::Sym_tensor::compute_A(), Lorene::Sym_tensor::compute_tilde_B_tt(), Lorene::Metric::con(), Lorene::contract(), Lorene::Metric::cov(), Lorene::Scalar::derive_con(), Lorene::Tensor::derive_cov(), Lorene::Scalar::derive_cov(), Lorene::Tensor_sym::derive_cov(), Lorene::Sym_tensor::derive_lie(), Lorene::Vector::divergence(), Lorene::Time_slice_conf::ff, Lorene::Map::get_mg(), Lorene::Tensor::get_mp(), Lorene::Mg3d::get_nzone(), Lorene::Tensor::get_triad(), Lorene::Time_slice_conf::hata(), hh(), Lorene::Tensor::inc_dzpuis(), Lorene::Time_slice::jtime, Lorene::Time_slice_conf::ln_psi(), Lorene::log(), Lorene::maxabs(), Lorene::Time_slice_conf::nn(), Lorene::Vector::ope_killing_conf(), Lorene::Time_slice_conf::psi(), Lorene::Time_slice_conf::psi4(), Lorene::Tensor::set(), Lorene::Sym_tensor_tt::set_A_tildeB(), Lorene::Scalar::set_dzpuis(), Lorene::Tensor::set_etat_zero(), Lorene::Scalar::set_spectral_va(), source_A_hata_evol, source_A_hh_evol, source_B_hata_evol, source_B_hh_evol, Lorene::Scalar::std_spectral_base(), Lorene::Time_slice_conf::tgam(), Lorene::Time_slice::the_time, Lorene::Tensor::trace(), and Lorene::Valeur::ylm().
|
inherited |
Computes the vector from the conformally-rescaled momentum , using the momentum constraint.
Definition at line 840 of file time_slice_conf.C.
References Lorene::abs(), Lorene::contract(), Lorene::Tensor::get_index_type(), Lorene::max(), and Lorene::Vector::poisson().
|
protectedvirtualinherited |
Deletes all the derived quantities.
Reimplemented from Lorene::Time_slice.
Definition at line 350 of file time_slice_conf.C.
References Lorene::Time_slice::del_deriv(), Lorene::Time_slice_conf::p_hdirac, Lorene::Time_slice_conf::p_ln_psi, Lorene::Time_slice_conf::p_psi4, Lorene::Time_slice_conf::p_tgamma, Lorene::Time_slice_conf::p_vec_X, and Lorene::Time_slice_conf::set_der_0x0().
void Lorene::Tslice_dirac_max::evolve | ( | double | pdt, |
int | nb_time_steps, | ||
int | niter_elliptic, | ||
double | relax_elliptic, | ||
int | check_mod, | ||
int | save_mod, | ||
int | method_poisson_vect = 6 , |
||
int | nopause = 1 , |
||
const char * | graph_device = 0x0 , |
||
bool | verbose = true , |
||
const Scalar * | ener_euler = 0x0 , |
||
const Vector * | mom_euler = 0x0 , |
||
const Scalar * | s_euler = 0x0 , |
||
const Sym_tensor * | strain_euler = 0x0 |
||
) |
Time evolution by resolution of Einstein equations.
pdt | time step dt. |
nb_time_steps | number of time steps for the evolution |
niter_elliptic | number of iterations if the resolution of elliptic equations |
relax_elliptic | relaxation factor for the elliptic equations |
check_mod | determines the frequency of check of the constraint equations: they are checked every check_mod time step |
save_mod | determines the frequency of writing to file the monotoring quantities: they are written to file every save_mod time step |
method | method_poisson_vect to be used for solving vector Poisson equation (for the shift), see Vector::poisson(double, const Metric_flat&, int) const . |
nopause | = 1 if no pause between each time step, 0 otherwise |
graph_device | name of type of graphical device: 0x0 (default value) will result in interactive choice; "/xwin" in X-Window display and "/n" in no output. |
Definition at line 123 of file tslice_dirac_max_evolve.C.
References Lorene::Time_slice_conf::A_hata_evol, A_hh_evol, Lorene::Time_slice_conf::B_hata_evol, B_hh_evol, Lorene::Time_slice::beta(), Lorene::Time_slice::depth, Lorene::Time_slice_conf::ff, Lorene::Map::get_bvect_spher(), Lorene::Map::get_mg(), Lorene::Tensor::get_mp(), Lorene::Mg3d::get_np(), Lorene::Mg3d::get_nt(), Lorene::Mg3d::get_nzone(), Lorene::Scalar::get_spectral_base(), Lorene::Tensor::get_triad(), Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, Lorene::Base_val::mult_cost(), Lorene::Base_val::mult_x(), Lorene::Time_slice_conf::nn(), Lorene::Sym_tensor_tt::set_A_tildeB(), and Lorene::Scalar::std_spectral_base().
|
inherited |
Induced metric at the current time step (jtime
)
Definition at line 98 of file time_slice_access.C.
References Lorene::Time_slice::gam_dd(), and Lorene::Time_slice::p_gamma.
|
virtualinherited |
Induced metric (covariant components ) at the current time step (jtime
)
Reimplemented from Lorene::Time_slice.
Definition at line 611 of file time_slice_conf.C.
References Lorene::Time_slice::gam_dd_evol, Lorene::Time_slice::jtime, Lorene::Time_slice_conf::psi4(), Lorene::Time_slice_conf::tgam(), and Lorene::Time_slice::the_time.
|
virtualinherited |
Induced metric (contravariant components ) at the current time step (jtime
)
Reimplemented from Lorene::Time_slice.
Definition at line 622 of file time_slice_conf.C.
References Lorene::Time_slice::gam_uu_evol, Lorene::Time_slice::jtime, Lorene::Time_slice_conf::psi4(), Lorene::Time_slice_conf::tgam(), and Lorene::Time_slice::the_time.
|
inlineinherited |
Gets the latest value of time step index.
Definition at line 346 of file time_slice.h.
References Lorene::Time_slice::jtime.
|
inlineinherited |
Gets the order of the finite-differences scheme.
Definition at line 343 of file time_slice.h.
References Lorene::Time_slice::scheme_order.
|
inlineinherited |
Gets the time coordinate t at successive time steps.
Definition at line 349 of file time_slice.h.
References Lorene::Time_slice::the_time.
|
virtualinherited |
Conformal representation of the traceless part of the extrinsic curvature: .
Returns the value at the current time step (jtime
).
Definition at line 775 of file time_slice_conf.C.
References Lorene::Time_slice::beta(), Lorene::Sym_tensor::derive_lie(), Lorene::Vector::divergence(), Lorene::Time_slice_conf::ff, Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice_conf::hh(), Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, Lorene::Time_slice_conf::nn(), Lorene::Vector::ope_killing_conf(), Lorene::Time_slice_conf::psi(), Lorene::Time_slice_conf::psi4(), Lorene::Time_slice::scheme_order, and Lorene::Time_slice::the_time.
|
virtual |
Vector which vanishes in Dirac gauge.
It is null in the present case...
Reimplemented from Lorene::Time_slice_conf.
Definition at line 510 of file tslice_dirac_max.C.
References Lorene::Time_slice_conf::ff, Lorene::Metric::get_mp(), Lorene::Metric_flat::get_triad(), Lorene::Time_slice_conf::p_hdirac, and Lorene::Tensor::set_etat_zero().
|
virtual |
Deviation of the conformal metric from the flat metric : .
Returns the value at the current time step (jtime
).
Reimplemented from Lorene::Time_slice_conf.
Definition at line 477 of file tslice_dirac_max.C.
References A_hh_evol, B_hh_evol, hh_det_one(), Lorene::Time_slice_conf::hh_evol, and Lorene::Time_slice::jtime.
|
protected |
Computes from the values of A and and using the condition , which fixes the trace of .
j | time step at which the computation of is required. |
Definition at line 112 of file tslice_dirac_max_setAB.C.
References A_hh_evol, Lorene::Time_slice::adm_mass_evol, B_hh_evol, Lorene::Metric_flat::con(), Lorene::Time_slice_conf::ff, Lorene::Time_slice::gam_dd_evol, Lorene::Time_slice::gam_uu_evol, Lorene::Metric_flat::get_triad(), Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, Lorene::maxabs(), Lorene::Time_slice::p_gamma, Lorene::Time_slice_conf::p_hdirac, Lorene::Time_slice_conf::p_tgamma, Lorene::Sym_tensor_trans::set_AtBtt_det_one(), Lorene::Tensor::set_etat_zero(), Lorene::Sym_tensor::set_longit_trans(), Lorene::Time_slice_conf::tgam(), Lorene::Time_slice::the_time, Lorene::Sym_tensor_trans::the_trace(), and trh_evol.
|
protected |
Computes from the TT part using the condition , which fixes the trace of .
hijtt | : the TT part. |
Definition at line 171 of file tslice_dirac_max_setAB.C.
References Lorene::Time_slice::adm_mass_evol, Lorene::Time_slice_conf::ff, Lorene::Time_slice::gam_dd_evol, Lorene::Time_slice::gam_uu_evol, Lorene::Tensor::get_mp(), Lorene::Metric_flat::get_triad(), Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, Lorene::maxabs(), Lorene::Time_slice::p_gamma, Lorene::Time_slice_conf::p_hdirac, Lorene::Time_slice_conf::p_tgamma, Lorene::Tensor::set_etat_zero(), Lorene::Sym_tensor::set_longit_trans(), Lorene::Sym_tensor_trans::set_tt_part_det_one(), Lorene::Time_slice_conf::tgam(), Lorene::Time_slice::the_time, Lorene::Sym_tensor_trans::the_trace(), and trh_evol.
|
virtual |
Computes valid initial data by solving the constraint equations in the conformal thin-sandwich approach.
uu | value of (freely specifiable data). This quantity must be trace-free with respect to the conformal metric , reflecting the unimodular character of . |
trk_in | value of (freely specifiable data) |
trk_point | value of (freely specifiable data) |
pdt | time step, to be used in order to fill depth slices |
precis | convergence threshold required to stop the iteration |
method_poisson_vect | method to be used for solving vector Poisson equation (for the shift), see Vector::poisson(double, const Metric_flat&, int) const . |
graph_device | name of type of graphical device: 0x0 (default value) will result in interactive choice; "/xwin" in X-Window display and "/n" in no output. |
ener_dens | matter energy density E as measured by the Eulerian observer; this quantity is passed as a pointer, the null value of which (default) meaning E=0. |
mom_dens | matter momentum density J as measured by the Eulerian observer; this quantity is passed as a pointer, the null value of which (default) meaning J=0. |
trace_stress | trace of the matter stress S as measured by the Eulerian observer; this quantity is passed as a pointer, the null value of which (default) meaning S=0. |
Reimplemented from Lorene::Time_slice_conf.
Definition at line 355 of file tslice_dirac_max.C.
References Lorene::Time_slice::depth, Lorene::Scalar::get_etat(), Lorene::Map::get_mg(), Lorene::Tensor::get_mp(), Lorene::Mg3d::get_nzone(), Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice_conf::initial_data_cts(), and Lorene::Time_slice::jtime.
|
protected |
Copy the sources source_A_XXX_evol
and source_B_XXX_evol
to all time-steps.
Definition at line 414 of file tslice_dirac_max_setAB.C.
References Lorene::Time_slice::depth, Lorene::Time_slice::jtime, source_A_hata_evol, source_A_hh_evol, source_B_hata_evol, source_B_hh_evol, and Lorene::Time_slice::the_time.
|
virtualinherited |
Extrinsic curvature tensor (covariant components ) at the current time step (jtime
)
Reimplemented from Lorene::Time_slice.
Definition at line 633 of file time_slice_conf.C.
References Lorene::Time_slice::gam(), Lorene::Time_slice::jtime, Lorene::Time_slice::k_dd_evol, Lorene::Time_slice_conf::k_uu(), and Lorene::Time_slice::the_time.
|
virtualinherited |
Extrinsic curvature tensor (contravariant components ) at the current time step (jtime
)
Reimplemented from Lorene::Time_slice.
Definition at line 646 of file time_slice_conf.C.
References Lorene::Time_slice::gam(), Lorene::Time_slice_conf::hata(), Lorene::Time_slice::jtime, Lorene::Time_slice::k_uu_evol, Lorene::Time_slice_conf::psi(), Lorene::Time_slice_conf::psi4(), Lorene::Time_slice::the_time, and Lorene::Time_slice_conf::trk().
|
inherited |
Logarithm of at the current time step (jtime
).
Definition at line 722 of file time_slice_conf.C.
References Lorene::log(), Lorene::Time_slice_conf::p_ln_psi, Lorene::Time_slice_conf::psi(), and Lorene::Scalar::std_spectral_base().
|
virtualinherited |
Lapse function N at the current time step (jtime
)
Reimplemented from Lorene::Time_slice.
Definition at line 594 of file time_slice_conf.C.
References Lorene::Time_slice::jtime, Lorene::Time_slice::n_evol, Lorene::Time_slice_conf::npsi_evol, Lorene::Time_slice_conf::psi_evol, and Lorene::Time_slice::the_time.
|
virtualinherited |
Factor at the current time step (jtime
).
Definition at line 735 of file time_slice_conf.C.
References Lorene::Time_slice::jtime, Lorene::Time_slice::n_evol, Lorene::Time_slice_conf::npsi_evol, Lorene::Time_slice_conf::psi_evol, and Lorene::Time_slice::the_time.
void Lorene::Tslice_dirac_max::operator= | ( | const Tslice_dirac_max & | tin | ) |
Assignment to another Tslice_dirac_max.
Definition at line 324 of file tslice_dirac_max.C.
References A_hh_evol, B_hh_evol, Lorene::Time_slice_conf::operator=(), source_A_hata_evol, source_A_hh_evol, source_B_hata_evol, source_B_hh_evol, and trh_evol.
|
protectedvirtual |
Operator >> (virtual function called by the operator<<).
Reimplemented from Lorene::Time_slice_conf.
Definition at line 569 of file tslice_dirac_max.C.
References A_hh_evol, B_hh_evol, Lorene::Time_slice::jtime, Lorene::maxabs(), Lorene::Time_slice_conf::operator>>(), and trh_evol.
|
virtualinherited |
Conformal factor relating the physical metric to the conformal one: .
is defined by
Returns the value at the current time step (jtime
).
Definition at line 696 of file time_slice_conf.C.
References Lorene::Time_slice::jtime, Lorene::Time_slice::n_evol, Lorene::Time_slice_conf::npsi_evol, Lorene::Time_slice_conf::psi_evol, and Lorene::Time_slice::the_time.
|
inherited |
Factor at the current time step (jtime
).
Definition at line 710 of file time_slice_conf.C.
References Lorene::Time_slice_conf::p_psi4, Lorene::pow(), Lorene::Time_slice_conf::psi(), and Lorene::Scalar::std_spectral_base().
|
protectedvirtual |
Total or partial saves in a binary file.
This protected method is to be called either from public method save
or from method sauve
of a derived class.
fich | binary file |
partial_save | indicates whether the whole object must be saved. |
Reimplemented from Lorene::Time_slice_conf.
Definition at line 590 of file tslice_dirac_max.C.
References A_hh(), A_hh_evol, B_hh(), B_hh_evol, Lorene::Time_slice::depth, Lorene::fwrite_be(), Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, and Lorene::Time_slice_conf::sauve().
|
inherited |
Saves in a binary file.
The saved data is sufficient to restore the whole time slice via the constructor from file.
rootname | root for the file name; the current time step index will be appended to it. |
Definition at line 464 of file time_slice.C.
References Lorene::Time_slice::beta(), Lorene::Time_slice::depth, Lorene::fwrite_be(), Lorene::Map::get_mg(), Lorene::Tensor::get_mp(), Lorene::Tensor::get_triad(), Lorene::Time_slice::jtime, Lorene::Time_slice::nn(), Lorene::Base_vect::sauve(), Lorene::Time_slice::sauve(), Lorene::Mg3d::sauve(), and Lorene::Map::sauve().
|
virtualinherited |
Sets the potentials A and of the TT part (see the documentation of Sym_tensor
for details).
Sets the value at the current time step (jtime
).
Definition at line 894 of file time_slice_conf.C.
References Lorene::Time_slice_conf::A_hata_evol, Lorene::Time_slice_conf::B_hata_evol, Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice::jtime, Lorene::Time_slice::k_dd_evol, Lorene::Time_slice::k_uu_evol, and Lorene::Time_slice::the_time.
Sets the potentials A and of the TT part of (see the documentation of Sym_tensor
for details).
is not modified. Sets the value at the current time step (jtime
).
Definition at line 87 of file tslice_dirac_max_setAB.C.
References A_hh_evol, Lorene::Time_slice::adm_mass_evol, B_hh_evol, Lorene::Time_slice::gam_dd_evol, Lorene::Time_slice::gam_uu_evol, Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, Lorene::Time_slice::p_gamma, Lorene::Time_slice_conf::p_hdirac, Lorene::Time_slice_conf::p_tgamma, Lorene::Time_slice::the_time, and trh_evol.
|
protectedinherited |
Sets to 0x0
all the pointers on derived quantities.
Definition at line 364 of file time_slice_conf.C.
References Lorene::Time_slice_conf::p_hdirac, Lorene::Time_slice_conf::p_ln_psi, Lorene::Time_slice_conf::p_psi4, Lorene::Time_slice_conf::p_tgamma, and Lorene::Time_slice_conf::p_vec_X.
|
virtualinherited |
Sets the conformal representation of the traceless part of the extrinsic curvature: .
Sets the value at the current time step (jtime
), and updates the potentials A_hata_evol
, B_hata_evol
and p_vec_X
accordingly.
Definition at line 537 of file time_slice_conf.C.
References Lorene::Time_slice_conf::A_hata_evol, Lorene::Time_slice_conf::B_hata_evol, Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice::jtime, Lorene::Time_slice::k_dd_evol, Lorene::Time_slice::k_uu_evol, and Lorene::Time_slice::the_time.
|
virtualinherited |
Sets the conformal representation of the traceless part of the extrinsic curvature from its potentials A, and .
These potentials must be up-to-date. It sets the value at the current time step (jtime
).
Definition at line 569 of file time_slice_conf.C.
References Lorene::Time_slice_conf::A_hata_evol, Lorene::Time_slice_conf::B_hata_evol, Lorene::Time_slice_conf::ff, Lorene::Tensor::get_mp(), Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice::jtime, Lorene::Time_slice::k_dd_evol, Lorene::Time_slice::k_uu_evol, Lorene::Vector::ope_killing_conf(), Lorene::Time_slice_conf::p_vec_X, Lorene::Sym_tensor_tt::set_A_tildeB(), and Lorene::Time_slice::the_time.
|
virtualinherited |
Sets the TT part of (see member hata_evol
).
Sets the value at current time-step (jtime
) and updates the potentials A and .
Definition at line 549 of file time_slice_conf.C.
References Lorene::Time_slice_conf::A_hata_evol, Lorene::Time_slice_conf::B_hata_evol, Lorene::Sym_tensor::compute_A(), Lorene::Sym_tensor::compute_tilde_B_tt(), Lorene::Scalar::dec_dzpuis(), Lorene::Scalar::get_dzpuis(), Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice::jtime, Lorene::Time_slice::k_dd_evol, Lorene::Time_slice::k_uu_evol, and Lorene::Time_slice::the_time.
|
virtual |
Sets the deviation of the conformal metric from the flat metric : .
must be such that . Sets the value at the current time step (jtime
).
Reimplemented from Lorene::Time_slice_conf.
Definition at line 339 of file tslice_dirac_max.C.
References A_hh_evol, B_hh_evol, Lorene::Time_slice::jtime, Lorene::Time_slice_conf::set_hh(), source_A_hata_evol, source_A_hh_evol, source_B_hata_evol, source_B_hh_evol, and trh_evol.
Sets the potentials and of the TT part of (see the documentation of Sym_tensor_tt
for details).
The value of is then deduced from the unimodularity condition on the conformal metric. Sets the value at the current time step (jtime
).
Definition at line 411 of file tslice_dirac_max.C.
References Lorene::Scalar::dec_dzpuis(), Lorene::Time_slice_conf::ff, Lorene::Map::flat_met_spher(), Lorene::Map::get_bvect_spher(), Lorene::Tensor::get_mp(), Lorene::Metric_flat::get_triad(), Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, Lorene::Tensor::set_etat_zero(), Lorene::Sym_tensor_tt::set_khi_mu(), Lorene::Sym_tensor::set_longit_trans(), Lorene::Time_slice::the_time, Lorene::Sym_tensor_trans::trace_from_det_one(), and trh_evol.
|
virtualinherited |
Sets the factor at the current time step (jtime
) and deletes the value of N.
Definition at line 481 of file time_slice_conf.C.
References Lorene::Time_slice::adm_mass_evol, Lorene::Time_slice::jtime, Lorene::Time_slice::n_evol, Lorene::Time_slice_conf::npsi_evol, and Lorene::Time_slice::the_time.
|
virtualinherited |
Sets the factor at the current time step (jtime
) and deletes the value of .
Definition at line 456 of file time_slice_conf.C.
References Lorene::Time_slice::adm_mass_evol, Lorene::Time_slice::gam_dd_evol, Lorene::Time_slice::gam_uu_evol, Lorene::Time_slice::jtime, Lorene::Time_slice_conf::npsi_evol, Lorene::Time_slice::p_gamma, Lorene::Time_slice_conf::p_ln_psi, Lorene::Time_slice_conf::p_psi4, Lorene::Time_slice_conf::psi_evol, and Lorene::Time_slice::the_time.
|
virtualinherited |
Sets the conformal factor relating the physical metric to the conformal one: .
is defined by
Sets the value at the current time step (jtime
) and deletes the value of N.
Definition at line 431 of file time_slice_conf.C.
References Lorene::Time_slice::adm_mass_evol, Lorene::Time_slice::gam_dd_evol, Lorene::Time_slice::gam_uu_evol, Lorene::Time_slice::jtime, Lorene::Time_slice::n_evol, Lorene::Time_slice::p_gamma, Lorene::Time_slice_conf::p_ln_psi, Lorene::Time_slice_conf::p_psi4, Lorene::Time_slice_conf::psi_evol, and Lorene::Time_slice::the_time.
|
virtualinherited |
Sets the conformal factor relating the physical metric to the conformal one: .
is defined by
Sets the value at the current time step (jtime
) and deletes the value of .
Definition at line 407 of file time_slice_conf.C.
References Lorene::Time_slice::adm_mass_evol, Lorene::Time_slice::gam_dd_evol, Lorene::Time_slice::gam_uu_evol, Lorene::Time_slice::jtime, Lorene::Time_slice_conf::npsi_evol, Lorene::Time_slice::p_gamma, Lorene::Time_slice_conf::p_ln_psi, Lorene::Time_slice_conf::p_psi4, Lorene::Time_slice_conf::psi_evol, and Lorene::Time_slice::the_time.
|
inlineinherited |
Sets the order of the finite-differences scheme.
Definition at line 334 of file time_slice.h.
References Lorene::Time_slice::scheme_order.
|
virtual |
Sets the trace, with respect to the flat metric ff
, of .
Sets the value at the current time step (jtime
). Note that this method does not ensure that the conformal metric is unimodular.
Definition at line 440 of file tslice_dirac_max.C.
References Lorene::Time_slice::adm_mass_evol, Lorene::Time_slice::gam_dd_evol, Lorene::Time_slice::gam_uu_evol, Lorene::Time_slice_conf::hh_evol, Lorene::Time_slice::jtime, Lorene::Time_slice::p_gamma, Lorene::Time_slice_conf::p_hdirac, Lorene::Time_slice_conf::p_tgamma, source_A_hata_evol, source_A_hh_evol, source_B_hata_evol, source_B_hh_evol, Lorene::Time_slice::the_time, and trh_evol.
|
virtual |
Solves the elliptic equation for the shift vector from (Eq.
(73) of Bonazzola et al. 2004).
method | method to be used for solving vector Poisson equation (for the shift), see Vector::poisson(double, const Metric_flat&, int) const . |
Definition at line 227 of file tslice_dirac_max_solve.C.
References Lorene::Time_slice_conf::aa(), Lorene::Time_slice::beta(), Lorene::contract(), Lorene::Tensor::derive_con(), Lorene::Sym_tensor::divergence(), Lorene::Vector::divergence(), Lorene::Time_slice_conf::ff, hh(), Lorene::Tensor::inc_dzpuis(), Lorene::maxabs(), Lorene::Time_slice_conf::nn(), and Lorene::Vector::poisson().
|
virtual |
Solves the elliptic equation for (maximal slicing condition + Hamiltonian constraint)
ener_dens | conformal matter energy density , where E is measured by the Eulerian observer; this quantity is passed as a pointer, the null value of which (default) meaning E=0. |
trace_stress | trace of the conformal matter stress , where S is measured by the Eulerian observer; this quantity is passed as a pointer, the null value of which (default) meaning S=0. |
Definition at line 121 of file tslice_dirac_max_solve.C.
References Lorene::Time_slice_conf::aa(), Lorene::contract(), Lorene::Metric::cov(), Lorene::Tensor::derive_cov(), Lorene::Scalar::derive_cov(), Lorene::Tensor_sym::derive_cov(), Lorene::Time_slice_conf::ff, Lorene::Scalar::get_etat(), Lorene::Tensor::get_mp(), hh(), Lorene::Scalar::inc_dzpuis(), Lorene::Time_slice_conf::npsi(), Lorene::Scalar::poisson(), Lorene::Time_slice_conf::psi(), Lorene::Time_slice_conf::psi4(), Lorene::Scalar::set_etat_zero(), Lorene::Time_slice_conf::tgam(), and Lorene::Tensor::up_down().
Solves the elliptic equation for the conformal factor $$ (Hamiltonian constraint).
ener_dens | conformal matter energy density , where E is measured by the Eulerian observer; this quantity is passed as a pointer, the null value of which (default) meaning E=0. |
Definition at line 175 of file tslice_dirac_max_solve.C.
References Lorene::contract(), Lorene::Metric::cov(), Lorene::Tensor::derive_cov(), Lorene::Scalar::derive_cov(), Lorene::Tensor_sym::derive_cov(), Lorene::Time_slice_conf::ff, Lorene::Scalar::get_etat(), Lorene::Tensor::get_mp(), Lorene::Time_slice_conf::hata(), hh(), Lorene::Scalar::inc_dzpuis(), Lorene::Scalar::poisson(), Lorene::pow(), Lorene::Time_slice_conf::psi(), Lorene::Scalar::set_etat_zero(), Lorene::Time_slice_conf::tgam(), and Lorene::Tensor::up_down().
|
virtualinherited |
Conformal metric Returns the value at the current time step (jtime
).
Reimplemented in Lorene::Isol_hor.
Definition at line 750 of file time_slice_conf.C.
References Lorene::Metric_flat::con(), Lorene::Time_slice_conf::ff, Lorene::Time_slice_conf::hh(), and Lorene::Time_slice_conf::p_tgamma.
|
virtual |
Computes the trace h
, with respect to the flat metric ff
, of .
Returns the value at the current time step (jtime
).
Definition at line 550 of file tslice_dirac_max.C.
References hh_det_one(), Lorene::Time_slice::jtime, and trh_evol.
|
virtual |
Trace K of the extrinsic curvature at the current time step (jtime
).
It is null in the present case (maximal slicing)
Reimplemented from Lorene::Time_slice_conf.
Definition at line 494 of file tslice_dirac_max.C.
References Lorene::Time_slice_conf::ff, Lorene::Metric::get_mp(), Lorene::Time_slice::jtime, Lorene::Scalar::set_etat_zero(), Lorene::Time_slice::the_time, and Lorene::Time_slice::trk_evol.
|
virtualinherited |
Vector representing the longitudinal part of .
(see the documentation of hata_evol
)
Definition at line 828 of file time_slice_conf.C.
References Lorene::Time_slice_conf::ff, Lorene::Time_slice_conf::hata_evol, Lorene::Time_slice::jtime, and Lorene::Time_slice_conf::p_vec_X.
|
mutableprotectedinherited |
Potential A associated with the symmetric tensor .
(see documentation of Sym_tensor::p_aaa
).
Definition at line 550 of file time_slice.h.
|
mutableprotected |
The A potential of .
(see the documentation of Sym_tensor::p_aaa
for details).
Definition at line 980 of file time_slice.h.
|
mutableprotectedinherited |
ADM mass at each time step, since the creation of the slice.
At a given time step j
, adm_mass_evol
[j] is a 1-D Tbl
of size the number nz
of domains, containing the "ADM mass" evaluated at the outer boundary of each domain. The true ADM mass is thus the last value, i.e. adm_mass_evol
[j](nz-1).
Definition at line 236 of file time_slice.h.
|
mutableprotectedinherited |
Potential associated with the symmetric tensor .
(see documentation of Sym_tensor::p_tilde_b
).
Definition at line 555 of file time_slice.h.
|
mutableprotected |
The potential of .
(see the documentation of Sym_tensor::p_tilde_b
for details).
Definition at line 986 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the shift vector .
Definition at line 222 of file time_slice.h.
|
protectedinherited |
Number of stored time slices.
Definition at line 182 of file time_slice.h.
|
protectedinherited |
Pointer on the flat metric with respect to which the conformal decomposition is performed.
Definition at line 510 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the covariant components of the induced metric .
Definition at line 201 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the contravariant components of the induced metric .
Definition at line 206 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the components .
It is the conformal representation of the traceless part of the extrinsic curvature: . One can uniquely (up to the boundary conditions) define the decomposition: , where represents the longitudinal part and is the transverse-traceless part.
Definition at line 545 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the components .
It is the deviation of the conformal metric from the flat metric : .
Definition at line 533 of file time_slice.h.
|
protectedinherited |
Time step index of the latest slice.
Definition at line 193 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the covariant components of the extrinsic curvature tensor .
Definition at line 211 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the contravariant components of the extrinsic curvature tensor .
Definition at line 216 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the lapse function N.
Definition at line 219 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the factor .
Definition at line 525 of file time_slice.h.
|
mutableprotectedinherited |
Pointer on the induced metric at the current time step (jtime
)
Definition at line 242 of file time_slice.h.
|
mutableprotectedinherited |
Pointer on the vector (which vanishes in Dirac gauge), at the current time step (jtime
).
Definition at line 574 of file time_slice.h.
|
mutableprotectedinherited |
Pointer on the logarithm of at the current time step (jtime
)
Definition at line 569 of file time_slice.h.
|
mutableprotectedinherited |
Pointer on the factor at the current time step (jtime
)
Definition at line 566 of file time_slice.h.
|
mutableprotectedinherited |
Pointer on the conformal metric at the current time step (jtime
)
Definition at line 563 of file time_slice.h.
|
mutableprotectedinherited |
Pointer on the vector representing the longitudinal part of .
(see the documentation of hata_evol
)
Definition at line 580 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the conformal factor relating the physical metric to the conformal one: .
is defined by
Definition at line 520 of file time_slice.h.
|
protectedinherited |
Order of the finite-differences scheme for the computation of time derivatives.
This order is not constant and can be adjusted via set_scheme_order()
.
Definition at line 190 of file time_slice.h.
|
mutableprotected |
The potential A of the source of equation for .
(see the documentation of Sym_tensor::p_aaa
for details).
Definition at line 1004 of file time_slice.h.
|
mutableprotected |
The A potential of the source of equation for .
(see the documentation of Sym_tensor::p_aaa
for details).
Definition at line 992 of file time_slice.h.
|
mutableprotected |
The potential of the source of equation for .
(see the documentation of Sym_tensor::p_tilde_b
for details).
Definition at line 1010 of file time_slice.h.
|
mutableprotected |
The potential of the source of equation for .
(see the documentation of Sym_tensor::p_tilde_b
for details).
Definition at line 998 of file time_slice.h.
|
protectedinherited |
Time label of each slice.
Definition at line 196 of file time_slice.h.
|
mutableprotected |
The trace, with respect to the flat metric ff
, of .
Definition at line 1013 of file time_slice.h.
|
mutableprotectedinherited |
Values at successive time steps of the trace K of the extrinsic curvature.
Definition at line 227 of file time_slice.h.