24char nullite_plm_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/nullite_plm.C,v 1.8 2014/10/13 08:53:14 j_novak Exp $" ;
93#include "type_parite.h"
106int nullite_plm_t_leg (
int j,
int nt,
int k,
int np) {
108 int m = (k%2 == 0) ? k/2 : (k-1)/2 ;
109 int borne_sup = nt-1 ;
111 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
112 return 0 ;
else return 1 ;
120int nullite_plm_nonsym (
int j,
int nt,
int k,
int np) {
122 int m = (k%2 == 0) ? k/2 : (k-1)/2 ;
123 int borne_sup = (m%2 == 0) ? nt-1 : nt-2 ;
124 int borne_inf = (m%2 == 0) ? m/2 : (m-1)/2 ;
125 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
126 return 0 ;
else return 1 ;
133int nullite_plm_nonsym_anti (
int j,
int nt,
int k,
int np) {
135 int m = (k%2 == 0) ? k/2 : (k-1)/2 ;
137 int borne_sup = nt-2 ;
139 int borne_inf = (m%2 == 0) ? m/2 : (m-1)/2 ;
140 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
141 return 0 ;
else return 1 ;
151int nullite_plm_sym (
int j,
int nt,
int k,
int np) {
153 int m = (k%2 == 0) ? k : k-1 ;
154 int borne_inf = m/2 ;
155 int borne_sup = nt-1 ;
156 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
157 return 0 ;
else return 1 ;
165int nullite_plm_t_leg_ip(
int j,
int nt,
int k,
int np) {
167 int m = (k%2 == 0) ? k : k-1 ;
168 int borne_sup = nt-2 ;
169 int borne_inf = m/2 ;
170 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
171 return 0 ;
else return 1 ;
179int nullite_plm_t_leg_pi(
int j,
int nt,
int k,
int np) {
186 m = (k%2 == 0) ? k-1 : k ;
189 int borne_sup = nt-2 ;
190 int borne_inf = (m-1)/2 ;
191 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
192 return 0 ;
else return 1 ;
199int nullite_plm_t_leg_ii(
int j,
int nt,
int k,
int np) {
206 m = (k%2 == 0) ? k-1 : k ;
209 int borne_sup = nt-2 ;
210 int borne_inf = (m+1)/2 ;
211 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
212 return 0 ;
else return 1 ;
219int nullite_plm_t_leg_mp (
int j,
int nt,
int k,
int np) {
221 int m = (k%2 == 0) ? k : k-1 ;
223 int borne_sup = nt-1 ;
224 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
225 return 0 ;
else return 1 ;
233int nullite_plm_t_leg_mi (
int j,
int nt,
int k,
int np) {
235 int m = 2*( (k-1) / 2) + 1 ;
237 int borne_sup = nt-1 ;
238 if ((j<borne_inf) || (j>borne_sup) || (k==1) || (k>np))
239 return 0 ;
else return 1 ;
249int nullite_plm (
int j,
int nt,
int k,
int np, Base_val base) {
254 int base_t = (base.b[0] &
MSQ_T) ;
255 int base_p = (base.b[0] &
MSQ_P) ;
263 result = nullite_plm_nonsym (j, nt, k, np) ;
268 result = nullite_plm_nonsym_anti (j, nt, k, np) ;
272 result = nullite_plm_t_leg (j, nt, k, np) ;
276 cout <<
"nullite_plm : cas inconnu ..." << endl ;
284 result = nullite_plm_sym (j, nt, k, np) ;
289 result = nullite_plm_t_leg_ip (j, nt, k, np) ;
293 result = nullite_plm_t_leg_mp (j, nt, k, np) ;
297 cout <<
"nullite_plm : cas inconnu ..." << endl ;
305 result = nullite_plm_t_leg_pi (j, nt, k, np) ;
309 result = nullite_plm_t_leg_ii (j, nt, k, np) ;
313 result = nullite_plm_t_leg_mi (j, nt, k, np) ;
317 cout <<
"nullite_plm : cas inconnu ..." << endl ;
323 cout <<
"nullite_plm : cas inconnu ..." << endl ;
#define T_LEG_MP
fct. de Legendre associees avec m pair
#define P_COSSIN_P
dev. sur Phi = 2*phi, freq. paires
#define T_LEG_PI
fct. de Legendre associees paires avec m impair
#define T_LEG
fct. de Legendre associees
#define T_LEG_P
fct. de Legendre associees paires
#define T_LEG_IP
fct. de Legendre associees impaires avec m pair
#define P_COSSIN
dev. standart
#define P_COSSIN_I
dev. sur Phi = 2*phi, freq. impaires
#define T_LEG_MI
fct. de Legendre associees avec m impair
#define MSQ_T
Extraction de l'info sur Theta.
#define T_LEG_II
fct. de Legendre associees impaires avec m impair
#define T_LEG_I
fct. de Legendre associees impaires
#define MSQ_P
Extraction de l'info sur Phi.
#define T_LEG_PP
fct. de Legendre associees paires avec m pair