diff --git a/src/operator.jl b/src/operator.jl index 62216268db3cb3ed4861dfd01b3046bf9f476185..b86afcd8fdc3b184dbde9e50308e48e819b10287 100644 --- a/src/operator.jl +++ b/src/operator.jl @@ -127,8 +127,8 @@ function quadrature_composite_lagrange_midpoint(p) points = Matrix{Float64}(undef, 2, n) k = 0 - for I in Iterators.product(ntuple(_ -> 0:p, d_ + 1)...) - I[1] + I[2] + I[3] != p && continue + for I in Iterators.product(ntuple(_ -> 0:p, d_)...) + sum(Tuple(I)) > p && continue k += 1 weights[k] = 1 / n points[1, k] = I[1] / p