23char sxdsdx_1d_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/sxdsdx_1d.C,v 1.5 2014/10/13 08:53:27 j_novak Exp $" ;
60#include "type_parite.h"
82void _sxdsdx_1d_pas_prevu(
int nr,
double* tb,
double *res) {
83 cout <<
"sxdsdx pas prevu..." << endl ;
84 cout <<
" valeurs: " << tb <<
" " << res << endl ;
85 cout <<
"nr : " << nr << endl ;
96void _dsdx_1d_r_cheb(
int nr,
double* tb,
double *xo)
102 som = 2*(nr-1) * tb[nr-1] ;
104 for (
int i = nr-4 ; i >= 0 ; i -= 2 ) {
105 som += 2*(i+1) * tb[i+1] ;
108 som = 2*(nr-2) * tb[nr-2] ;
110 for (
int i = nr-5 ; i >= 0 ; i -= 2 ) {
111 som += 2*(i+1) * tb[i+1] ;
122void sxm1_1d_cheb (
int,
double*) ;
124void _sxmundsdx_1d_r_chebu(
int nr,
double* tb,
double *xo)
130 som = 2*(nr-1) * tb[nr-1] ;
132 for (
int i = nr-4 ; i >= 0 ; i -= 2 ) {
133 som += 2*(i+1) * tb[i+1] ;
136 som = 2*(nr-2) * tb[nr-2] ;
138 for (
int i = nr-5 ; i >= 0 ; i -= 2 ) {
139 som += 2*(i+1) * tb[i+1] ;
144 sxm1_1d_cheb (nr, xo) ;
151void _sxdsdx_1d_r_chebp(
int nr,
double* tb,
double *xo)
157 som = 8 * (nr-1) * tb[nr-1] ;
159 for (
int i = nr-4 ; i >= 0 ; i -= 2 ) {
160 som += 8 * (i+1) * tb[i+1] ;
164 som = 8 * (nr-2) * tb[nr-2] ;
166 for (
int i = nr-5 ; i >= 0 ; i -= 2 ) {
167 som += 8 * (i+1) * tb[i+1] ;
179void _sxdsdx_1d_r_chebi(
int nr,
double* tb,
double *xo)
186 som = 4 * (2*(nr-1)+1) * tb[nr-1] ;
188 for (
int i = nr-4 ; i >= 0 ; i -= 2 ) {
189 som += 4 * (2*(i+1)+1) * tb[i+1] ;
193 som = 4 * (2*(nr-2)+1) * tb[nr-2] ;
195 for (
int i = nr-5 ; i >= 0 ; i -= 2 ) {
196 som += 4 * (2*(i+1)+1) * tb[i+1] ;
208void sxdsdx_1d(
int nr,
double **tb,
int base_r)
212static void (*sxdsdx_1d[
MAX_BASE])(int,
double *,
double *) ;
219 sxdsdx_1d[i] = _sxdsdx_1d_pas_prevu ;
228 double *result =
new double[nr] ;
229 sxdsdx_1d[base_r](nr, *tb, result) ;
#define MAX_BASE
Nombre max. de bases differentes.
#define R_CHEBU
base de Chebychev ordinaire (fin), dev. en 1/r
#define R_CHEBI
base de Cheb. impaire (rare) seulement
#define TRA_R
Translation en R, used for a bitwise shift (in hex)
#define R_CHEB
base de Chebychev ordinaire (fin)
#define R_CHEBP
base de Cheb. paire (rare) seulement