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);
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]);
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
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)