63 const Vector& shift_comp_temp,
double om) {
65 Vector shift_auto(shift_auto_temp) ;
67 Vector shift_comp(shift_comp_temp) ;
69 Vector shift_old (shift_auto) ;
72 assert ((orientation==0) || (orientation == M_PI)) ;
74 assert ((orientation_autre==0) || (orientation_autre == M_PI)) ;
76 int alignes = (orientation == orientation_autre) ? 1 : -1 ;
84 shift_tot.
set(1).
import(alignes*shift_comp(1)) ;
85 shift_tot.set(2).import(alignes*shift_comp(2)) ;
86 shift_tot.set(3).import(shift_comp(3)) ;
88 shift_tot = shift_tot + shift_auto ;
90 double indic = (orientation == 0) ? 1 : -1 ;
94 for (
int i=1 ; i<=3 ; i++) {
106 Vector derive_r (shift_auto.get_mp(), CON, *shift_auto.get_triad()) ;
107 for (
int i=1 ; i<=3 ; i++)
108 derive_r.
set(i) = tbi(i).dsdr() ;
113 Valeur val_hor (shift_auto.get_mp().get_mg()) ;
114 Valeur fonction_radiale (shift_auto.get_mp().get_mg()) ;
115 Scalar enleve (shift_auto.get_mp()) ;
118 for (
int comp=1 ; comp<=3 ; comp++) {
120 for (
int k=0 ; k<np ; k++)
121 for (
int j=0 ; j<nt ; j++)
122 for (
int i=0 ; i<nr ; i++)
123 val_hor.set(1, k, j, i) = derive_r(comp).
124 val_grid_point(1, k, j, 0) ;
126 double r_0 = shift_auto.get_mp().val_r (1, -1, 0, 0) ;
127 double r_1 = shift_auto.get_mp().val_r (1, 1, 0, 0) ;
129 fonction_radiale =
pow(r_1-shift_auto.get_mp().r, 3.)*
130 (shift_auto.get_mp().r-r_0)/
pow(r_1-r_0, 3.) ;
131 fonction_radiale.
annule(0) ;
132 fonction_radiale.annule(2,
nz-1) ;
134 enleve = fonction_radiale * val_hor ;
135 enleve.set_spectral_va().set_base (shift_auto(comp).
136 get_spectral_va().get_base()) ;
138 if (
norme(enleve)(1) != 0)
139 shift_auto.set(comp) = shift_auto(comp) - enleve ;
140 if (
norme(shift_auto(comp))(1) > 1e-5) {
142 if (erreur < diff(1))
147 shift_auto.change_triad(shift_auto.get_mp().get_bvect_spher()) ;
164 for (
int i=1 ; i<=3 ; i++) {
169 for (
int i=1 ; i<=3 ; i++)
170 shift(i).get_spectral_va().coef_i() ;
172 tbi.
set(1) = *shift(1).get_spectral_va().c -
omega*
mp.
y ;
173 tbi.
set(2) = *shift(2).get_spectral_va().c +
omega*
mp.
x ;
174 if (shift(3).get_etat() != ETATZERO)
175 tbi.
set(3) = *shift(3).get_spectral_va().c ;
186 for (
int i=1 ; i<=3 ; i++)
187 derive_r.set(i) = tbi(i).dsdr() ;
198 double r_0 =
mp.
val_r(1, -1, 0, 0) ;
199 double r_1 =
mp.
val_r(1, 1, 0, 0) ;
201 for (
int comp=1 ; comp<=3 ; comp++) {
202 val_hor.annule_hard() ;
203 for (
int k=0 ; k<np ; k++)
204 for (
int j=0 ; j<nt ; j++)
205 for (
int i=0 ; i<nr ; i++)
206 val_hor.set(1, k, j, i) = derive_r(comp).val_grid_point(1, k, j, 0) ;
208 fonction_radiale =
pow(r_1-
mp.
r, 3.)* (
mp.
r-r_0)/
pow(r_1-r_0, 3.) ;
209 fonction_radiale.
annule(0) ;
210 fonction_radiale.annule(2,
nz-1) ;
212 enleve = fonction_radiale*val_hor ;
213 enleve.set_spectral_va().base = shift(comp).get_spectral_va().base ;
215 Scalar copie (shift(comp)) ;
216 shift.
set(comp) = shift(comp)-enleve ;
219 assert (shift(comp).check_dzpuis(0)) ;
223 if (norm(1) > 1e-5) {
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
void annule(int l)
Sets the Cmp to zero in a given domain.
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.
double regularisation(const Vector &shift_auto, const Vector &shift_comp, double ang_vel)
Corrects shift_auto in such a way that the total is equal to zero in the horizon, which should ensure the regularity of .
Tensor field of valence 0 (or component of a tensorial field).
void coef_i() const
Computes the physical value of *this.
double omega
Angular velocity in LORENE's units.
virtual double val_r(int l, double xi, double theta, double pphi) const
Returns the value of the radial coordinate r for a given in a given domain.
Values and coefficients of a (real-value) function.
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
double get_rot_phi() const
Returns the angle between the x –axis and X –axis.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
void annule_hard()
Sets the Scalar to zero in a hard way.
Tbl norme(const Cmp &)
Sums of the absolute values of all the values of the Cmp in each domain.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Mtbl * c
Values of the function at the points of the multi-grid.
Vector beta_auto
Shift function .
Vector beta
Shift function .
int get_nzone() const
Returns the number of domains.
Cmp pow(const Cmp &, int)
Power .
double regularise_one()
Corrects the shift in the innermost shell, so that it remains and that equals zero on the horizon...
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
Map_af & mp
Affine mapping.
Coord y
y coordinate centered on the grid
Coord x
x coordinate centered on the grid
void set_etat_c_qcq()
Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl c )...
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Valeur & set_spectral_va()
Returns va (read/write version)
Scalar & set(int)
Read/write access to a component.
const Map & get_mp() const
Returns the mapping.
Tbl diffrelmax(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (max version).
const Valeur & get_spectral_va() const
Returns va (read only version)
Coord r
r coordinate centered on the grid