TL: solving numerical issue

Parallel-in-Time · Jun 21, 2024 · 1c4aecd · 1c4aecd
1 parent 239f3b5
commit 1c4aecd
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/qmat/qcoeff/collocation.py b/qmat/qcoeff/collocation.py
@@ -30,8 +30,9 @@ def __init__(self, nNodes, nodeType, quadType, tLeft=0, tRight=1):
         nodes *= (tRight-tLeft)
         nodes += tLeft
         self.tLeft, self.tRight = tLeft, tRight
-        # Rounding for safety
-        np.round(nodes, 14, out=nodes) # TODO : check if necessary ...
+        # Safety when bound should be included ...
+        if np.allclose(tLeft, nodes[0]): nodes[0] = tLeft
+        if np.allclose(tRight, nodes[-1]): nodes[-1] = tRight
         self._nodes = nodes
 
         # Lagrange approximation based on nodes