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);
75 "unsupported quadrature family \"" << opt.get_family_name() <<
"\"");
80 size_type n = n_node_gauss(opt.get_order());
81 vector<T> zeta(n), omega(n);
85 wx (
x(zeta[i], zeta[j]), omega[i]*omega[j]);
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)