Added special case of degree=0 for local_element_end in FemAssembly…

…Grid
pyccel · Oct 6, 2022 · e031c07 · e031c07 · spossann · Sep 12, 2024
1 parent 327afb3
commit e031c07
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/psydac/fem/grid.py b/psydac/fem/grid.py
@@ -162,7 +162,10 @@ def __init__( self, space, start, end, *, quad_order=None, nderiv=1, parent_star
             local_element_start = self._spans.searchsorted( degree + start )
             local_element_end   = self._spans.searchsorted( degree + end   )
         else:
-            if end+1 >= degree:
+            if degree == 0:
+                local_element_start = self._spans.searchsorted( degree if start == 0   else 1 + start)
+                local_element_end   = self._spans.searchsorted( end )
+            elif end+1 >= degree:
                 local_element_start = self._spans.searchsorted( degree if start == 0   else 1 + start)
                 local_element_end   = self._spans.searchsorted( end if end   == n-1 else 1 + end )
             else: