1 #include "rheolef/quadrature.h"
22 #include "rheolef/gauss_jacobi.h"
48 quadrature_option::family_type
f = opt.get_family();
52 if (
f == quadrature_option::equispaced) {
61 wx (
x(2*
T(
int(i))/r-1,2*
T(
int(j))/r-1), w);
74 check_macro (
f == quadrature_option::gauss,
75 "unsupported quadrature family \"" << opt.get_family_name() <<
"\"");
81 vector<T> zeta(
n), omega(
n);
85 wx (
x(zeta[i], zeta[j]), omega[i]*omega[j]);
90 #define _RHEOLEF_instanciation(T) \
91 template void quadrature_on_geo<T>::init_square (quadrature_option);
see the Float page for the full documentation
base::size_type size_type
void init_square(quadrature_option opt)
integrate_option quadrature_option
This file is part of Rheolef.
_RHEOLEF_instanciation(Float, sequential, std::allocator< Float >) _RHEOLEF_instanciation(Float
void gauss_jacobi(Size R, typename std::iterator_traits< OutputIterator1 >::value_type alpha, typename std::iterator_traits< OutputIterator1 >::value_type beta, OutputIterator1 zeta, OutputIterator2 omega)