79 (
double*,
const int,
const int,
const int,
const double,
double*) ;
81 (
double*,
const int,
const int,
const double,
double*) ;
83 (
double*,
const int,
const double,
double*) ;
84 static int premier_appel = 1 ;
88 if (premier_appel == 1) {
94 som_phi[i/2] = som_phi_pas_prevu ;
96 som_tet[i] = som_tet_pas_prevu ;
97 som_r[i] = som_r_pas_prevu ;
121 assert (
etat != ETATNONDEF) ;
126 if (
etat == ETATZERO ) {
147 assert(
etat == ETATQCQ) ;
150 if (tbcf->
get_etat() == ETATZERO ) {
156 assert(tbcf->
get_etat() == ETATQCQ) ;
158 double* cf = tbcf->
t ;
161 double* trp =
new double [np+2] ;
162 double* trtp =
new double [(np+2)*nt] ;
169 som_tet[base_t](cf, nt, np, theta, trp) ;
176 som_r[base_r](cf, nr, nt, np, x, trtp) ;
177 som_tet[base_t](trtp, nt, np, theta, trp) ;
187 som_phi[base_p](trp, np, phi, &resu) ;
210 (
double*,
const int,
const int,
const int,
const double,
double*) ;
211 static int premier_appel = 1 ;
215 if (premier_appel == 1) {
220 som_r[i] = som_r_pas_prevu ;
232 assert (
etat != ETATNONDEF) ;
237 if (
etat == ETATZERO ) {
265 assert(
etat == ETATQCQ) ;
268 if (tbcf->
get_etat() == ETATZERO ) {
274 assert(tbcf->
get_etat() == ETATQCQ) ;
276 double* cf = tbcf->
t ;
279 double* coef_tp =
new double [(np+2)*nt] ;
285 som_r[base_r](cf, nr, nt, np, x, coef_tp) ;
290 double* pi = coef_tp ;
307 for (
int k=3 ; k<np+1 ; k+=2) {
308 int m_par = (k/2)%2 ;
310 for (
int j=0 ; j<nt ; j++) {
311 somt += (*pi) * tab_theta(m_par, j, j0) ;
314 resu += somt * tab_phi(k, k0) ;
323 cout <<
"Mtbl_cf::val_point_jk_asymy: unknown theta basis ! " #define MAX_BASE_2
Smaller maximum bases used for phi (and higher dimensions for now)
#define P_COSSIN
dev. standart
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
#define TRA_P
Translation en Phi, used for a bitwise shift (in hex)
#define TRA_T
Translation en Theta, used for a bitwise shift (in hex)
#define MSQ_P
Extraction de l'info sur Phi.
const Mg3d * mg
Pointer on the multi-grid Mgd3 on which this is defined.
#define T_COS
dev. cos seulement
int get_etat() const
Gives the logical state.
#define T_SIN
dev. sin seulement
#define TRA_R
Translation en R, used for a bitwise shift (in hex)
#define R_CHEBI
base de Cheb. impaire (rare) seulement
int etat
Logical state (ETATNONDEF , ETATQCQ or ETATZERO ).
#define R_CHEBP
base de Cheb. paire (rare) seulement
#define T_COS_P
dev. cos seulement, harmoniques paires
double val_point_jk_asymy(int l, double x, int j, int k) const
Computes the value of the field represented by *this at an arbitrary point in , but collocation point...
#define MSQ_T
Extraction de l'info sur Theta.
double * t
The array of double.
int * b
Array (size: nzone ) of the spectral basis in each domain.
#define MSQ_R
Extraction de l'info sur R.
#define T_SIN_P
dev. sin seulement, harmoniques paires
#define R_CHEBPIM_I
Cheb. pair-impair suivant m, impair pour m=0.
#define R_CHEBPIM_P
Cheb. pair-impair suivant m, pair pour m=0.
const Tbl & theta_functions(int l, int nt) const
Values of the theta basis functions at the theta collocation points.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
#define T_COSSIN_CI
cos impair-sin pair alternes, cos pour m=0
const Tbl & phi_functions(int l, int np) const
Values of the phi basis functions at the phi collocation points.
Base_val base
Bases of the spectral expansions.
#define P_COSSIN_I
dev. sur Phi = 2*phi, freq. impaires
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
#define R_CHEBU
base de Chebychev ordinaire (fin), dev. en 1/r
#define P_COSSIN_P
dev. sur Phi = 2*phi, freq. paires
#define T_COSSIN_CP
cos pair-sin impair alternes, cos pour m=0
#define MAX_BASE
Nombre max. de bases differentes.
double val_point_asymy(int l, double x, double theta, double phi) const
Computes the value of the field represented by *this at an arbitrary point, by means of the spectral ...
Tbl ** t
Array (size nzone ) of pointers on the Tbl 's which contain the spectral coefficients in each domain...
#define R_CHEB
base de Chebychev ordinaire (fin)