122 (
double*,
const int,
const int,
const int,
const double,
double*) ;
124 (
double*,
const int,
const int,
const double,
double*) ;
126 (
double*,
const int,
const double,
double*) ;
127 static int premier_appel = 1 ;
131 if (premier_appel == 1) {
137 som_phi[i/2] = som_phi_pas_prevu ;
139 som_tet[i] = som_tet_pas_prevu ;
140 som_r[i] = som_r_pas_prevu ;
176 assert (
etat != ETATNONDEF) ;
181 if (
etat == ETATZERO ) {
202 assert(
etat == ETATQCQ) ;
205 if (tbcf->
get_etat() == ETATZERO ) {
211 assert(tbcf->
get_etat() == ETATQCQ) ;
213 double* cf = tbcf->
t ;
216 double* trp =
new double [np+2] ;
217 double* trtp =
new double [(np+2)*nt] ;
224 som_tet[base_t](cf, nt, np, theta, trp) ;
231 som_r[base_r](cf, nr, nt, np, x, trtp) ;
232 som_tet[base_t](trtp, nt, np, theta, trp) ;
242 som_phi[base_p](trp, np, phi, &resu) ;
265 (
double*,
const int,
const int,
const int,
const double,
double*) ;
266 static int premier_appel = 1 ;
270 if (premier_appel == 1) {
275 som_r[i] = som_r_pas_prevu ;
290 assert (
etat != ETATNONDEF) ;
295 if (
etat == ETATZERO ) {
323 assert(
etat == ETATQCQ) ;
326 if (tbcf->
get_etat() == ETATZERO ) {
332 assert(tbcf->
get_etat() == ETATQCQ) ;
334 double* cf = tbcf->
t ;
337 double* coef_tp =
new double [(np+2)*nt] ;
343 som_r[base_r](cf, nr, nt, np, x, coef_tp) ;
348 double* pi = coef_tp ;
353 for (
int j=0 ; j<nt ; j++) {
354 somt += (*pi) * tab_theta(0, j, j0) ;
357 resu = somt * tab_phi(0, k0) ;
377 for (
int k=2 ; k<np+1 ; k++) {
379 for (
int j=0 ; j<nt ; j++) {
380 somt += (*pi) * tab_theta(0, j, j0) ;
383 resu += somt * tab_phi(k, k0) ;
395 for (
int k=2 ; k<np+1 ; k++) {
396 int m_par = (k/2)%2 ;
398 for (
int j=0 ; j<nt ; j++) {
399 somt += (*pi) * tab_theta(m_par, j, j0) ;
402 resu += somt * tab_phi(k, k0) ;
408 cout <<
"Mtbl_cf::val_point_jk: unknown theta basis ! " << endl ;
443 for (
int i=0; i<nr; i++)
444 resu +=
operator()(l, k0, j0, i) ;
466 for (
int i=0; i<nr; i++) {
#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 R_CHEBPI_I
Cheb. pair-impair suivant l impair pour l=0.
const Tbl & operator()(int l) const
Read-only of the Tbl containing the coefficients in a given domain.
#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
#define T_COSSIN_SP
sin pair-cos impair alternes, sin pour m=0
int get_etat() const
Gives the logical state.
#define R_LEGP
base de Legendre paire (rare) seulement
#define R_LEGI
base de Legendre impaire (rare) seulement
#define R_JACO02
base de Jacobi(0,2) ordinaire (finjac)
#define T_SIN
dev. sin seulement
double val_out_bound_jk(int l, int j, int k) const
Computes the angular coefficient of index j,k of the field represented by *this at by means of the s...
#define TRA_R
Translation en R, used for a bitwise shift (in hex)
double val_in_bound_jk(int l, int j, int k) const
Computes the angular coefficient of index j,k of the field represented by *this at by means of the s...
#define T_COS_I
dev. cos seulement, harmoniques impaires
double val_point_jk(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 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
#define MSQ_T
Extraction de l'info sur Theta.
#define T_COSSIN_C
dev. cos-sin alternes, cos pour m=0
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 T_COSSIN_SI
sin impair-cos pair alternes, sin pour m=0
#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.
double val_point(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 ...
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 R_CHEBPI_P
Cheb. pair-impair suivant l pair pour l=0.
#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.
#define T_SIN_I
dev. sin seulement, harmoniques impaires
Tbl ** t
Array (size nzone ) of pointers on the Tbl 's which contain the spectral coefficients in each domain...
#define T_COSSIN_S
dev. cos-sin alternes, sin pour m=0
#define R_LEG
base de Legendre ordinaire (fin)
#define R_CHEB
base de Chebychev ordinaire (fin)