23char som_asymy_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/som_asymy.C,v 1.3 2014/10/13 08:53:26 j_novak Exp $" ;
81 const double x,
double* trtp) {
89 double* cheb =
new double [nr] ;
92 for (i=2; i<nr; i++) {
93 cheb[i] = 2*x* cheb[i-1] - cheb[i-2] ;
103 for (k=3 ; k<np+1 ; k+=2) {
104 for (j=0 ; j<nt ; j++) {
107 for (i=0 ; i<nr ; i++) {
108 *po += (*pi) * cheb[i] ;
129 const double x,
double* trtp) {
137 double* cheb =
new double [nr] ;
140 for (i=2; i<nr; i++) {
141 cheb[i] = 2*x* cheb[i-1] - cheb[i-2] ;
151 for (k=3 ; k<np+1 ; k+=2) {
152 for (j=0 ; j<nt ; j++) {
155 for (i=0 ; i<nr ; i++) {
156 *po += (*pi) * cheb[i] ;
178 const int np,
const double x,
double* trtp) {
186 double* cheb =
new double [2*nr] ;
189 for (i=2 ; i<2*nr ; i++) {
190 cheb[i] = 2*x* cheb[i-1] - cheb[i-2] ;
201 for (k=3 ; k<np+1 ; k+=2) {
203 for (j=0 ; j<nt ; j++) {
206 for (i=0 ; i<nr ; i++) {
207 *po += (*pi) * cheb[2*i + m] ;
228 const int np,
const double x,
double* trtp) {
236 double* cheb =
new double [2*nr] ;
239 for (i=2 ; i<2*nr ; i++) {
240 cheb[i] = 2*x* cheb[i-1] - cheb[i-2] ;
250 for (k=3 ; k<np+1 ; k+=2) {
252 for (j=0 ; j<nt ; j++) {
255 for (i=0 ; i<nr ; i++) {
256 *po += (*pi) * cheb[2*i + 1 - m] ;
282 const double tet,
double* to) {
290 double* cossin =
new double [2*nt] ;
291 for (j=0 ; j<2*nt ; j +=2) {
292 cossin[j] =
cos(j * tet) ;
293 cossin[j+1] =
sin((j+1) * tet) ;
302 for (k=3 ; k<np+1 ; k+=2) {
305 for (j=0 ; j<nt ; j++) {
306 *po += (*pi) * cossin[2*j + m] ;
326 const double tet,
double* to) {
334 double* cossin =
new double [2*nt] ;
335 for (j=0 ; j<2*nt ; j +=2) {
336 cossin[j] =
cos((j+1) * tet) ;
337 cossin[j+1] =
sin(j * tet) ;
346 for (k=3 ; k<np+1 ; k+=2) {
349 for (j=0 ; j<nt ; j++) {
350 *po += (*pi) * cossin[2*j + m] ;
370void som_phi_cossin_asymy(
double* ti,
const int np,
const double phi,
376 for (
int k=3 ; k<np+1 ; k +=2 ) {
378 *xo += ti[k] *
sin(m * phi) ;
Cmp sin(const Cmp &)
Sine.
Cmp cos(const Cmp &)
Cosine.
void som_r_cheb_asymy(double *, const int, const int, const int, const double, double *)
Cas R_CHEB ///.
void som_tet_cossin_cp_asymy(double *, const int, const int, const double, double *)
Cas T_COSSIN_CP ///.
void som_r_chebpim_p_asymy(double *, const int, const int, const int, const double, double *)
Cas R_CHEBPIM_P ///.
void som_r_chebpim_i_asymy(double *, const int, const int, const int, const double, double *)
Cas R_CHEBPIM_I ///.
void som_tet_cossin_ci_asymy(double *, const int, const int, const double, double *)
Cas T_COSSIN_CI ///.
void som_r_chebu_asymy(double *, const int, const int, const int, const double, double *)
Cas R_CHEBU ///.