28char tslice_dirax_max_setAB_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Time_slice/tslice_dirac_max_setAB.C,v 1.11 2014/10/13 08:53:48 j_novak Exp $" ;
77#include "time_slice.h"
128 wzero.set_etat_zero() ;
158 "Max. of absolute value of deviation from det tgam = 1") ;
163 "Max. of absolute value of deviation from det tgam = 1") ;
185 wzero.set_etat_zero() ;
212 "Max. of absolute value of deviation from det tgam = 1") ;
223 const Map& map =
hh().get_mp() ;
225 int nz = map.
get_mg()->get_nzone() ;
234 aij.annule_domain(nz-1) ;
243 log_N.std_spectral_base() ;
250 a_hat.annule_domain(nz-1) ;
260 -
beta().ope_killing_conf(ff) + 0.6666666666666667*
div_beta*
hh() ;
262 for (
int i=1;
i<=3;
i++)
263 for (
int j=
i;
j<=3;
j++)
267 tmp.set_spectral_va().ylm() ;
268 tmp.annule_domain(nz-1) ;
273 tmp.set_spectral_va().ylm() ;
274 tmp.annule_domain(nz-1) ;
294 for (
int i=1;
i<=3;
i++) {
295 for (
int j=1;
j<=
i;
j++) {
297 for (
int k=1;
k<=3;
k++) {
298 for (
int l=1;
l<=3;
l++) {
306 for (
int i=1;
i<=3;
i++) {
307 for (
int j=1;
j<=
i;
j++) {
309 for (
int k=1;
k<=3;
k++) {
310 for (
int l=1;
l<=3;
l++) {
311 for (
int m=1;
m<=3;
m++) {
312 for (
int n=1;
n<=3;
n++) {
355 for (
int i=1;
i<=3;
i++) {
356 for (
int j=1;
j<=
i;
j++) {
358 for (
int k=1;
k<=3;
k++) {
359 for (
int l=1;
l<=3;
l++) {
371 for (
int i=1;
i<=3;
i++) {
372 for (
int j=1;
j<=
i;
j++) {
374 for (
int k=1;
k<=3;
k++) {
375 for (
int l=1;
l<=3;
l++) {
391 for (
int i=1;
i<=3;
i++)
392 for (
int j=
i;
j<=3;
j++)
399 tmp.annule_domain(nz-1) ;
403 tmp.annule_domain(nz-1) ;
411void Tslice_dirac_max::initialize_sources_copy()
const {
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 Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Metric for tensor calculation.
Tensor field of valence 0 (or component of a tensorial field).
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
Transverse symmetric tensors of rank 2.
void set_AtBtt_det_one(const Scalar &a_in, const Scalar &tbtt_in, const Scalar *h_prev=0x0, Param *par_bc=0x0, Param *par_mat=0x0, double precis=1.e-14, int it_max=100)
Assigns the derived member A and computes from its TT-part (see Sym_tensor::compute_tilde_B_tt() ).
const Scalar & the_trace() const
Returns the trace of the tensor with respect to metric *met_div.
void set_tt_part_det_one(const Sym_tensor_tt &hijtt, const Scalar *h_prev=0x0, Param *par_mat=0x0, double precis=1.e-14, int it_max=100)
Assignes the TT-part of the tensor.
Transverse and traceless symmetric tensors of rank 2.
Class intended to describe valence-2 symmetric tensors.
Symmetric tensors (with respect to two of their arguments).
Evolution_std< Sym_tensor > hh_evol
Values at successive time steps of the components .
virtual const Scalar & nn() const
Lapse function N at the current time step (jtime )
virtual const Sym_tensor & hata() const
Conformal representation of the traceless part of the extrinsic curvature: .
const Scalar & ln_psi() const
Logarithm of at the current time step (jtime ).
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 Sym_tensor aa() const
Conformal representation of the traceless part of the extrinsic curvature: .
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 ).
virtual const Metric & tgam() const
Conformal metric Returns the value at the current time step (jtime ).
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 .
int jtime
Time step index of the latest slice.
Metric * p_gamma
Pointer on the induced metric at the current time step (jtime)
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 .
virtual const Vector & beta() const
shift vector at the current time step (jtime )
int depth
Number of stored time slices.
Evolution_std< double > the_time
Time label of each slice.
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 .
Evolution_std< Scalar > A_hh_evol
The A potential of .
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 .
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_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).
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 log(const Cmp &)
Neperian logarithm.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
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.
Standard units of space, time and mass.