From 57756e6485566e00f89078132bf7d3ac32d029c0 Mon Sep 17 00:00:00 2001
From: Stephan Hilb <stephan@ecshi.net>
Date: Sun, 22 Aug 2021 17:33:44 +0200
Subject: [PATCH] improve lagrange point generator
should maybe do an explicit enumeration instead of filtering
---
src/operator.jl | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/operator.jl b/src/operator.jl
index 6221626..b86afcd 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
--
GitLab