24char 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)) ;
85 shift_tot = shift_tot + shift_auto ;
87 double indic = (orientation == 0) ? 1 : -1 ;
91 for (
int i=1 ; i<=3 ; i++) {
96 tbi.
set(1) = *shift_tot(1).get_spectral_va().c - indic *om * shift_tot.
get_mp().
ya ;
97 tbi.
set(2) = *shift_tot(2).get_spectral_va().c + indic *om * shift_tot.
get_mp().
xa ;
104 for (
int i=1 ; i<=3 ; i++)
105 derive_r.
set(i) = tbi(i).dsdr() ;
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) ;
131 enleve = fonction_radiale * val_hor ;
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))
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++) {
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) ;
209 enleve = fonction_radiale*val_hor ;
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) {
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.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Coord y
y coordinate centered on the grid
Coord ya
Absolute y coordinate.
Coord r
r coordinate centered on the grid
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
Coord x
x coordinate centered on the grid
double get_rot_phi() const
Returns the angle between the x –axis and X –axis.
Coord xa
Absolute x coordinate.
virtual double val_r(int l, double xi, double theta, double pphi) const =0
Returns the value of the radial coordinate r for a given in a given domain.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Tensor field of valence 0 (or component of a tensorial field).
Valeur & set_spectral_va()
Returns va (read/write version)
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
Vector beta_auto
Shift function .
Map_af & mp
Affine mapping.
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,...
double omega
Angular velocity in LORENE's units.
double regularise_one()
Corrects the shift in the innermost shell, so that it remains and that equals zero on the horizon.
Vector beta
Shift function .
Values and coefficients of a (real-value) function.
void set_etat_c_qcq()
Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl c ).
void set_base(const Base_val &)
Sets the bases for spectral expansions (member base )
void annule(int l)
Sets the Valeur to zero in a given domain.
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).
void coef_i() const
Computes the physical value of *this.
Base_val base
Bases on which the spectral expansion is performed.
void annule_hard()
Sets the Valeur to zero in a hard way.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Scalar & set(int)
Read/write access to a component.
Tbl norme(const Cmp &)
Sums of the absolute values of all the values of the Cmp in each domain.
Cmp pow(const Cmp &, int)
Power .
Tbl diffrelmax(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (max version).
const Map & get_mp() const
Returns the mapping.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
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).