“…Typically, the rule is required to be exact, that is, R m (f ) ≡ 0 for each element of a predefined linear function space L. Moreover, the rule is said to be Gaussian if m is the minimal number of nodes t i ∈ R n at which f has to be evaluated. There is an extensive number of various quadrature rules depending on n (f is univariate [15], bivariate [28,39], multivariate [16]), domain shape (disc, hypercube, simplex) [37], and the type of the linear space (polynomials [15], splines [4,6,29,33], rational functions [38], smooth non-polynomial [7,27]). For polynomial multivariate integration, the field is well studied and the reader is referred to [37].…”