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) {
37 wx (
x(1/
T(3),1/
T(3),0), 1);
44 wx (
x(
T(
int(i))/r,
T(
int(j))/r,2*
T(
int(k))/r-1), w);
54 check_macro (
f == quadrature_option::gauss,
55 "unsupported quadrature family \"" << opt.get_family_name() <<
"\"");
57 switch (opt.get_order()) {
72 size_t r = opt.get_order();
76 vector<T> zeta0(n0), omega0(n0);
77 vector<T> zeta1(n1), omega1(n1);
78 vector<T> zeta2(n2), omega2(n2);
87 T eta_0 = (1+zeta0[i])*(1-zeta1[j])/4;
88 T eta_1 = (1+zeta1[j])/2;
91 wx (
x(eta_0,eta_1,eta_2), J*omega0[i]*omega1[j]*omega2[k]);
101 #define _RHEOLEF_instanciation(T) \
102 template void quadrature_on_geo<T>::init_prism (quadrature_option);