152 double pi = 3.141592653589793;
157 std::cerr <<
"CLENSHAW_CURTIS_POINTS - Fatal error!\n";
158 std::cerr <<
" Order is less than 1.\n";
167 for ( index = 1; index <= n; index++ )
169 x[index-1] = std::cos ( (
double ) ( n - index ) * pi
170 / (
double ) ( n - 1 ) );
173 if ( ( n % 2 ) == 1 )
215 double pi = 3.141592653589793;
221 std::cerr <<
"CLENSHAW_CURTIS_WEIGHTS - Fatal error!\n";
222 std::cerr <<
" Order is less than 1.\n";
231 for ( i = 1; i <= n; i++ )
233 theta = ( double ) ( i - 1 ) * pi / ( double ) ( n - 1 );
237 for ( j = 1; j <= ( n - 1 ) / 2; j++ )
239 if ( 2 * j == ( n - 1 ) )
248 w[i-1] = w[i-1] - b * std::cos ( 2.0 * (
double ) ( j ) * theta )
249 / ( double ) ( 4 * j * j - 1 );
253 w[0] = w[0] / ( double ) ( n - 1 );
254 for ( i = 1; i < n - 1; i++ )
256 w[i] = 2.0 * w[i] / ( double ) ( n - 1 );
258 w[n-1] = w[n-1] / ( double ) ( n - 1 );
300 double pi = 3.141592653589793;
305 std::cerr <<
"FEJER2_POINTS - Fatal error!\n";
306 std::cerr <<
" Order is less than 1.\n";
315 for ( index = 1; index <= n; index++ )
317 x[index-1] = std::cos ( (
double ) ( n + 1 - index ) * pi
318 / (
double ) ( n + 1 ) );
320 if ( ( n % 2 ) == 1 )
361 double pi = 3.141592653589793;
367 std::cerr <<
"FEJER2_WEIGHTS - Fatal error!\n";
368 std::cerr <<
" N < 1.\n";
382 for ( i = 1; i <= n; i++ )
384 theta = ( double ) ( n + 1 - i ) * pi
385 / ( double ) ( n + 1 );
389 for ( j = 1; j <= ( ( n - 1 ) / 2 ); j++ )
391 w[i-1] = w[i-1] - 2.0 * std::cos ( 2.0 * (
double ) ( j ) * theta )
392 / ( double ) ( 4 * j * j - 1 );
394 p = 2.0 * ( double ) ( ( ( n + 1 ) / 2 ) ) - 1.0;
395 w[i-1] = w[i-1] - std::cos ( ( p + 1.0 ) * theta ) / p;
397 for ( i = 0; i < n; i++ )
399 w[i] = 2.0 * w[i] / ( double ) ( n + 1 );
void hce_weights(int n, int dim, double w[])
void nco_compute_weights(int n, double w[])
void patterson_lookup_weights(int n, double w[])
void fejer2_points(int n, int dim, double x[])
void laguerre_compute_points(int order, double x[])
void nco_weights(int n, int dim, double w[])
void gen_laguerre_compute_points(int order, double alpha, double x[])
void hce_points(int n, int dim, double x[])
void gen_hermite_compute_weights(int order, double alpha, double w[])
void hcc_compute_points(int n, double x[])
void patterson_points(int n, int dim, double x[])
void nco_points(int n, int dim, double x[])
void jacobi_compute_points(int order, double alpha, double beta, double x[])
void ncc_weights(int n, int dim, double w[])
void ccn_points(int n, int dim, double x[])
void laguerre_weights(int n, int dim, double w[])
void hce_compute_points(int n, double x[])
void jacobi_compute_weights(int order, double alpha, double beta, double w[])
void hermite_compute_points(int order, double x[])
void jacobi_weights(int n, int dim, double w[])
double parameter(int dim, int offset)
void hermite_points(int n, int dim, double x[])
void hcc_weights(int n, int dim, double w[])
void legendre_compute_weights(int order, double w[])
void laguerre_points(int n, int dim, double x[])
void clenshaw_curtis_weights(int n, int dim, double w[])
void ccn_weights(int n, int dim, double w[])
void patterson_lookup_points(int n, double x[])
void laguerre_compute_weights(int order, double w[])
void gen_laguerre_compute_weights(int order, double alpha, double w[])
void gen_hermite_compute_points(int order, double alpha, double x[])
void clenshaw_curtis_points(int n, int dim, double x[])
void hermite_weights(int n, int dim, double w[])
void ccn_compute_points(int n, double x[])
void hcc_points(int n, int dim, double x[])
void hce_compute_weights(int n, double w[])
void gen_laguerre_points(int n, int dim, double x[])
void ncc_points(int n, int dim, double x[])
void hcc_compute_weights(int n, double w[])
void gen_laguerre_weights(int n, int dim, double w[])
void nco_compute_points(int n, double x[])
void hermite_genz_keister_lookup_points(int n, double x[])
void ncc_compute_weights(int n, double w[])
void patterson_weights(int n, int dim, double w[])
void gen_hermite_points(int n, int dim, double x[])
void gen_hermite_weights(int n, int dim, double w[])
void hermite_genz_keister_lookup_weights(int n, double w[])
void legendre_points(int n, int dim, double x[])
void hermite_genz_keister_points(int n, int dim, double x[])
void hermite_genz_keister_weights(int n, int dim, double w[])
void legendre_compute_points(int order, double x[])
void hermite_compute_weights(int order, double w[])
void ncc_compute_points(int n, double x[])
void jacobi_points(int n, int dim, double x[])
void ccn_compute_weights(int n, double w[])
void legendre_weights(int n, int dim, double w[])
void fejer2_weights(int n, int dim, double w[])