Add tests to handle Derivative when zero degree

src/compmec/nurbs/calculus.py

@@ -1,5 +1,5 @@
 from fractions import Fraction
-from typing import Any, Callable, Optional, Tuple
+from typing import Any, Callable, Optional
 
 import numpy as np
 
@@ -81,11 +81,6 @@ def nonrational_spline(curve: Curve) -> Curve:
         assert isinstance(curve, Curve)
         assert curve.weights is None
         assert curve.degree + 1 != curve.npts
-        if curve.degree == 0:
-            newknotvector = curve.knotvector.limits
-            anypoint = curve.ctrlpoints[0]
-            newctrlpoints = [0 * anypoint]
-            return Curve(newknotvector, newctrlpoints)
 
         knotvector = tuple(curve.knotvector)
         matrix = heavy.Calculus.derivate_nonrational_spline(knotvector)
@@ -216,8 +211,7 @@ def density(
         method: Optional[str] = None,
         nnodes: Optional[int] = None,
     ) -> float:
-        """Computes the integral :math:`I`
-        Computes the integral :math:`I`
+        """Computes the integral I
 
         The operation ``@`` is needed cause ``norm(curve(u)) = numpy.sqrt(curve(u) @ curve(u))``
 
@@ -290,7 +284,7 @@ def function(
         nnodes: Optional[int] = None,
     ) -> float:
         """Computes the integral I
-        
+
         .. math::
             I = \int_{a}^{b} g ( u ) \ du
 

src/compmec/nurbs/heavy.py

@@ -84,8 +84,6 @@ def number_type(number: Union[int, float, Fraction]):
     except TypeError:
         if isinstance(number, (int, np.integer)):
             return int
-        if isinstance(number, (float, np.floating)):
-            return float
         if isinstance(number, Fraction):
             return Fraction
         return float
@@ -1235,19 +1233,6 @@ def add_spline_curve(
         matrixb = Operations.matrix_transformation(knotvectorb, knotvectorc)
         return totuple(matrixa), totuple(matrixb)
 
-    @staticmethod
-    def sub_spline_curve(
-        knotvectora: Tuple[float], knotvectorb: Tuple[float]
-    ) -> Tuple["Matrix2D"]:
-        assert KnotVector.is_valid_vector(knotvectora)
-        assert KnotVector.is_valid_vector(knotvectorb)
-        assert knotvectora[0] == knotvectorb[0]
-        assert knotvectora[-1] == knotvectorb[-1]
-
-        matrixa, matrixb = MathOperations.add_spline_curve(knotvectora, knotvectorb)
-        matrixb = -np.array(matrixb, dtype="object")
-        return matrixa, totuple(matrixb)
-
     @staticmethod
     def mul_spline_curve(
         knotvectora: Tuple[float], knotvectorb: Tuple[float]

tests/test_calculus.py

@@ -123,7 +123,7 @@ def test_begin(self):
     def test_bezier(self):
         deltau = 1e-6
         usample = np.linspace(2 * deltau, 1 - 2 * deltau, 5)
-        for degree in range(1, 4):
+        for degree in range(0, 4):
             vector = GeneratorKnotVector.bezier(degree)
             curve = Curve(vector)
             points = np.random.uniform(-1, 1, curve.npts)
@@ -140,7 +140,7 @@ def test_bezier(self):
     )
     def test_spline(self):
         deltau = 1e-6
-        for degree in range(1, 4):
+        for degree in range(0, 4):
             for npts in range(degree + 3, degree + 9):
                 vector = GeneratorKnotVector.uniform(degree, npts)
                 knots = vector.knots
@@ -169,7 +169,7 @@ def test_spline(self):
     def test_rationalbezier(self):
         deltau = 1e-6
         usample = np.linspace(2 * deltau, 1 - 2 * deltau, 5)
-        for degree in range(1, 4):
+        for degree in range(0, 4):
             vector = GeneratorKnotVector.bezier(degree)
             curve = Curve(vector)
             points = np.random.uniform(-1, 1, curve.npts)
@@ -194,7 +194,7 @@ def test_rationalbezier(self):
     )
     def test_rationalspline(self):
         deltau = 1e-6
-        for degree in range(1, 4):
+        for degree in range(0, 4):
             for npts in range(degree + 3, degree + 9):
                 vector = GeneratorKnotVector.uniform(degree, npts)
                 knots = vector.knots
@@ -310,7 +310,9 @@ def test_lenght_integral(self):
         curve.ctrlpoints = points
 
         test = Integrate.lenght(curve)
-        good = sum(np.linalg.norm(points[i+1] - points[i]) for i in range(curve.npts-1))
+        good = sum(
+            np.linalg.norm(points[i + 1] - points[i]) for i in range(curve.npts - 1)
+        )
 
         assert abs(test - good) < 1e-9
 

tests/test_splinecurve.py

@@ -991,9 +991,26 @@ def test_fractions(self):
                 assert type(tval) is frac
                 assert tval == gval
 
+    @pytest.mark.order(5)
+    @pytest.mark.timeout(10)
+    @pytest.mark.dependency(depends=["TestOthers::test_begin"])
+    def test_fraction_function(self):
+        knotvector = GeneratorKnotVector.uniform(3, 10)
+        curve = Curve(knotvector)
+        curve.ctrlpoints = np.random.uniform(-1, 1, curve.npts)
+
+        num, den = curve.fraction()
+        newcurve = num / den
+        assert curve == newcurve
+
     @pytest.mark.order(5)
     @pytest.mark.dependency(
-        depends=["TestOthers::test_begin", "TestOthers::test_others"]
+        depends=[
+            "TestOthers::test_begin",
+            "TestOthers::test_others",
+            "TestOthers::test_fractions",
+            "TestOthers::test_fraction_function",
+        ]
     )
     def test_end(self):
         pass
@@ -1009,7 +1026,7 @@ def test_end(self):
         "TestSplitUnite::test_end",
         "TestDegreeOperations::test_end",
         "TestSumSubtract::test_end",
-        "TestOthers::test_others",
+        "TestOthers::test_end",
     ]
 )
 def test_end():