1 #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);
58 check_macro (
f == quadrature_option::gauss,
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) \
76 template void quadrature_on_geo<T>::init_hexahedron (quadrature_option);