24 char single_regul_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Isol_hor/single_regul.C,v 1.4 2014/10/13 08:53:01 j_novak Exp $" ;
60 const Vector& shift_comp_temp,
double om) {
62 Vector shift_auto(shift_auto_temp) ;
64 Vector shift_comp(shift_comp_temp) ;
66 Vector shift_old (shift_auto) ;
69 assert ((orientation==0) || (orientation == M_PI)) ;
71 assert ((orientation_autre==0) || (orientation_autre == M_PI)) ;
73 int alignes = (orientation == orientation_autre) ? 1 : -1 ;
81 shift_tot.
set(1).
import(alignes*shift_comp(1)) ;
82 shift_tot.set(2).import(alignes*shift_comp(2)) ;
83 shift_tot.set(3).import(shift_comp(3)) ;
85 shift_tot = shift_tot + shift_auto ;
87 double indic = (orientation == 0) ? 1 : -1 ;
91 for (
int i=1 ; i<=3 ; i++) {
103 Vector derive_r (shift_auto.get_mp(), CON, *shift_auto.get_triad()) ;
104 for (
int i=1 ; i<=3 ; i++)
105 derive_r.
set(i) = tbi(i).dsdr() ;
110 Valeur val_hor (shift_auto.get_mp().get_mg()) ;
111 Valeur fonction_radiale (shift_auto.get_mp().get_mg()) ;
112 Scalar enleve (shift_auto.get_mp()) ;
115 for (
int comp=1 ; comp<=3 ; comp++) {
117 for (
int k=0 ; k<np ; k++)
118 for (
int j=0 ; j<nt ; j++)
119 for (
int i=0 ; i<nr ; i++)
120 val_hor.set(1, k, j, i) = derive_r(comp).
121 val_grid_point(1, k, j, 0) ;
123 double r_0 = shift_auto.get_mp().val_r (1, -1, 0, 0) ;
124 double r_1 = shift_auto.get_mp().val_r (1, 1, 0, 0) ;
126 fonction_radiale =
pow(r_1-shift_auto.get_mp().r, 3.)*
127 (shift_auto.get_mp().r-r_0)/
pow(r_1-r_0, 3.) ;
128 fonction_radiale.
annule(0) ;
129 fonction_radiale.annule(2,
nz-1) ;
131 enleve = fonction_radiale * val_hor ;
132 enleve.set_spectral_va().set_base (shift_auto(comp).
133 get_spectral_va().get_base()) ;
135 if (
norme(enleve)(1) != 0)
136 shift_auto.set(comp) = shift_auto(comp) - enleve ;
137 if (
norme(shift_auto(comp))(1) > 1e-5) {
139 if (erreur < diff(1))
144 shift_auto.change_triad(shift_auto.get_mp().get_bvect_spher()) ;
161 for (
int i=1 ; i<=3 ; i++) {
166 for (
int i=1 ; i<=3 ; i++)
167 shift(i).get_spectral_va().coef_i() ;
169 tbi.
set(1) = *shift(1).get_spectral_va().c -
omega*
mp.
y ;
170 tbi.
set(2) = *shift(2).get_spectral_va().c +
omega*
mp.
x ;
171 if (shift(3).get_etat() != ETATZERO)
172 tbi.
set(3) = *shift(3).get_spectral_va().c ;
183 for (
int i=1 ; i<=3 ; i++)
184 derive_r.set(i) = tbi(i).dsdr() ;
195 double r_0 =
mp.
val_r(1, -1, 0, 0) ;
196 double r_1 =
mp.
val_r(1, 1, 0, 0) ;
198 for (
int comp=1 ; comp<=3 ; comp++) {
199 val_hor.annule_hard() ;
200 for (
int k=0 ; k<np ; k++)
201 for (
int j=0 ; j<nt ; j++)
202 for (
int i=0 ; i<nr ; i++)
203 val_hor.set(1, k, j, i) = derive_r(comp).val_grid_point(1, k, j, 0) ;
205 fonction_radiale =
pow(r_1-
mp.
r, 3.)* (
mp.
r-r_0)/
pow(r_1-r_0, 3.) ;
206 fonction_radiale.
annule(0) ;
207 fonction_radiale.annule(2,
nz-1) ;
209 enleve = fonction_radiale*val_hor ;
210 enleve.set_spectral_va().base = shift(comp).get_spectral_va().base ;
212 Scalar copie (shift(comp)) ;
213 shift.
set(comp) = shift(comp)-enleve ;
216 assert (shift(comp).check_dzpuis(0)) ;
220 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