fix computation of multiplicity

psydac/api/tests/test_api_feec_2d.py

@@ -729,9 +729,9 @@ def test_maxwell_2d_multiplicity():
         mult = 2
     )
 
-    TOL = 1e-6
-    ref = dict(error_l2_Ex = 4.35005708e-04,
-               error_l2_Ey = 4.35005708e-04,
+    TOL = 1e-5
+    ref = dict(error_l2_Ex = 4.350041934920621e-04,
+               error_l2_Ey = 4.350041934920621e-04,
                error_l2_Bz = 3.76106860e-03)
 
     assert abs(namespace['error_l2_Ex'] - ref['error_l2_Ex']) / ref['error_l2_Ex'] <= TOL

psydac/feec/tests/test_differentiation_matrices.py

@@ -911,10 +911,11 @@ def eval_div(fx, fy, fz, *eta):
 
 #==============================================================================
 if __name__ == '__main__':
-
-    test_Derivative_1D(domain=[0, 1], ncells=12, degree=3, periodic=False, seed=1, multiplicity=1)
+        
+    test_Derivative_1D(domain=[0, 1], ncells=3, degree=3, periodic=False, seed=1, multiplicity=1)
     test_Derivative_1D(domain=[0, 1], ncells=12, degree=3, periodic=True, seed=1, multiplicity=1)
 
+
     test_Gradient_2D(
         domain   = ([0, 1], [0, 1]),
         ncells   = (10, 15),

psydac/fem/splines.py

@@ -91,15 +91,13 @@ def __init__( self, degree, knots=None, grid=None, multiplicity=None, parent_mul
         if grid is None:
             grid = breakpoints(knots, degree)
 
-        indices = np.where(np.diff(knots[degree+1:-degree-1])>1e-15)[0]
+        indices = np.where(np.diff(knots[degree:len(knots)-degree])>1e-15)[0]
 
+        if len(indices)>0:
+            multiplicity = np.diff(indices).max(initial=1)
+        else:
+            multiplicity = max(1,len(knots[degree+1:-degree-1]))
 
-        if multiplicity is None: #V.C 20/10/23 : why computing the multiplicity again?
-            if len(indices)>0:
-                multiplicity = np.diff(indices).max(initial=1)
-            else:
-                multiplicity = max(1,len(knots[degree+1:-degree-1]))
-
         if parent_multiplicity is None:
             parent_multiplicity = multiplicity