21#include "rheolef/quadrature.h"
22#include "rheolef/gauss_jacobi.h"
30 quadrature_option::family_type
f = opt.get_family();
34 if (
f == quadrature_option::equispaced) {
44 wx (
x(2*
T(
int(i))/r-1,2*
T(
int(j))/r-1,2*
T(
int(k))/r-1), w);
59 "unsupported quadrature family \"" << opt.get_family_name() <<
"\"");
65 vector<T> zeta(
n), omega(
n);
70 wx (
x(zeta[i], zeta[j], zeta[k]), omega[i]*omega[j]*omega[k]);
75#define _RHEOLEF_instanciation(T) \
76template void quadrature_on_geo<T>::init_hexahedron (quadrature_option);
see the Float page for the full documentation
base::size_type size_type
void init_hexahedron(quadrature_option opt)
integrate_option quadrature_option
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
This file is part of Rheolef.
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)
_RHEOLEF_instanciation(Float, sequential, std::allocator< Float >) _RHEOLEF_instanciation(Float