23char multxpun_1d_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/multxpun_1d.C,v 1.3 2014/10/13 08:53:24 j_novak Exp $" ;
57#include "type_parite.h"
78void _multxpun_1d_pas_prevu(
int nr,
double* tb,
double *res) {
79 cout <<
"multxpun pas prevu..." << endl ;
80 cout <<
" valeurs: " << tb <<
" " << res << endl ;
81 cout <<
"nr : " << nr << endl ;
90void _multxpun_1d_r_jaco02 (
int nr,
double* tb,
double *xo)
94 for (
int i = 1 ; i < nr-1 ; i++) {
95 xo[i] = i*(i+2)/
double((i+1)*(2*i+1))*tb[i-1] + (i*i+3*i+3)/
double((i+1)*(i+2))*tb[i] + (i+1)*(i+3)/double((i+2)*(2*i+5))*tb[i+1] ;
97 xo[0] = 1.5*tb[0] + 0.3*tb[1] ;
98 xo[nr-1] = (nr*nr-1)/
double((nr)*(2*nr-1))*tb[nr-2] + (1+1/double((nr)*(nr+1)))*tb[nr-1] ;
106void multxpun_1d(
int nr,
double **tb,
int base_r)
110static void (*multxpun_1d[
MAX_BASE])(int,
double *,
double *) ;
117 multxpun_1d[i] = _multxpun_1d_pas_prevu ;
123 double *result =
new double[nr] ;
124 multxpun_1d[base_r](nr, *tb, result) ;
#define MAX_BASE
Nombre max. de bases differentes.
#define R_JACO02
base de Jacobi(0,2) ordinaire (finjac)
#define TRA_R
Translation en R, used for a bitwise shift (in hex)