Add tests

Author: eeshan9815

test/newton1d.jl

@@ -22,17 +22,24 @@ three_halves_pi = 3*big_pi/2
     f′(x) = 2*x*exp(x^2) + sin(x)
     f1(x) = x^4 - 10x^3 + 35x^2 - 50x + 24
     f1′(x) = 4x^3 - 30x^2 + 70x - 50
-
+    f2(x) = 4567x^2 - 9134x + 4567
+    f2′(x) = 9134x - 9134
+    f3(x) = (x^2 - 2)^2
+    f3′(x) = 4x * (x^2 - 2)
     for autodiff in (false, true)
         if autodiff
             rts1 = newton1d(sin, -5..5)
             rts2 = newton1d(f, -∞..∞)
             rts3 = newton1d(f1, -10..10)
+            rts4 = newton1d(f2, -10..11)
+            rts5 = newton1d(f3, -10..10)
 
         else
             rts1 = newton1d(sin, cos, -5..5)
             rts2 = newton1d(f, f′, -∞..∞)
             rts3 = newton1d(f1, f1′, -10..10)
+            rts4 = newton1d(f2, f2′, -10..11)
+            rts5 = newton1d(f3, f3′, -10..10)
         end
 
         @test length(rts1) == 3
@@ -42,12 +49,20 @@ three_halves_pi = 3*big_pi/2
         end
 
         @test length(rts2) == 1
-        @test (0..0) == rts2[1].interval && :unique == rts2[1].status
+        @test (0..0) == rts2[1].interval && :unknown == rts2[1].status
 
         @test length(rts3) == 4
         L = [1, 2, 3, 4]
         for i = 1:length(rts3)
             @test L[i] in rts3[i].interval && :unique == rts3[i].status
         end
+
+        @test length(rts4) == 1
+        @test 1 in rts4[1].interval && :unknown == rts4[1].status
+
+        L1 = [-sqrt(2), -sqrt(2), sqrt(2), sqrt(2)]
+        for i = 1:length(rts5)
+            @test L1[i] in rts5[i].interval && :unknown == rts5[i].status
+        end
     end
 end