4 #include "utilitaires.h"
9 #include "param_elliptic.h"
12 #include "excision_surf.h"
45 for (
int k=0; k<np; k++)
46 for (
int j=0; j<nt; j++) {
62 Scalar ff = lapse2*(c_theta_lap*thetaplus.
lapang() + c_theta_fin*(thetaplus - expa_fin));
72 Scalar ff_int = lapse2*(c_theta_lap*theta_int.
lapang() + c_theta_fin*(theta_int - expa_fin));
79 Scalar bound_theta = thetaplus + k_2;
119 for (
int k=0; k<np; k++)
120 for (
int j=0; j<nt; j++) {
137 Scalar K_1 = beta0*d_thetaplus + alpha0*d_thetaplus.
lapang() + gamma0*thetaplus ;
144 Scalar K_2 = beta0*d_thetaplus + alpha0*d_thetaplus.
lapang() + gamma0*thetaplus;
151 Scalar K_3 = beta0*d_thetaplus + alpha0*d_thetaplus.
lapang() + gamma0*thetaplus;
158 Scalar K_4 = beta0*d_thetaplus + alpha0*d_thetaplus.
lapang() + gamma0*thetaplus;
162 thetaplus =
expa + (1./6.)*
delta_t*(k_1 + 2.*k_2 + 2.*k_3 + k_4);
163 d_thetaplus =
dt_expa + (1./6.)*
delta_t*(K_1 + 2.*K_2 + 2.*K_3 + K_4);
Scalar & set_expa()
Sets the expansion function on the surface at time t (considering to protect this function)
int get_etat() const
Returns the logical state.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
Tensor field of valence 0 (or component of a tensorial field).
const Map & get_mp() const
Returns the mapping.
const Valeur & get_spectral_va() const
Returns va (read only version)
Valeur & set_spectral_va()
Returns va (read/write version)
Scalar & set_dt_expa()
Sets the time derivative of the expansion function on the surface at time t (considering to protect t...
const Scalar & lapang() const
Returns the angular Laplacian of *this , where .
void set_expa_hyperb(double alph0, double beta0, double gamma0)
Sets a new value for expansion rescaled over lapse (and its derivative), obtained by hyperbolic evolu...
double delta_t
The time step for evolution in parabolic drivers.
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
void annule_hard()
Sets the Scalar to zero in a hard way.
void set_expa_parab(double c_theta_lap, double c_theta_fin, Scalar &expa_fin)
Sets a new value for expansion rescaled over lapse (and its derivative), obtained by parabolic evolut...
Scalar expa
The 2d expansion, directly evolved from the initial excision with Einstein Equations.
Scalar lapse
The lapse defined on the 3 slice.
void ylm()
Computes the coefficients of *this.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
Scalar dt_expa
The time derivative of the expansion, derived from Einstein equations and arbitrary evolution.