Add more ForwardDiff related tests

JuliaIntervals · Dec 2, 2023 · 8a854e2 · 8a854e2
1 parent 6a8f853
commit 8a854e2
Showing 1 changed file with 27 additions and 0 deletions.
diff --git a/test/interval_tests/autodiff.jl → test/interval_tests/forwarddiff.jl b/test/interval_tests/autodiff.jl → test/interval_tests/forwarddiff.jl
@@ -23,4 +23,31 @@ end
         @test all(ForwardDiff.gradient(  v -> g(v[1]), [interval(-1, 1)]) .=== [interval(convert(Interval{Float64}, -2), convert(Interval{Float64}, 2), trv)])
         @test_broken all(ForwardDiff.hessian(   v -> g(v[1]), [interval(  0  )]) .=== [interval(convert(Interval{Float64}, -2), convert(Interval{Float64}, 2), trv)])
     end
+
+    @testset "sin" begin
+        x, y, w = interval(2), interval(-0.5, 0.5)
+        ϕ(t)    =  sin(x + (1+t)*w)
+        ϕ′(t)   =  cos(x + (1+t)*w) * w
+        ϕ′′(t)  = -sin(x + (1+t)*w) * w * w
+        ϕ′′′(t) = -cos(x + (1+t)*w) * w * w * w
+        dϕ(t)   = ForwardDiff.derivative(ϕ, t)
+        ddϕ(t)  = ForwardDiff.derivative(dϕ, t)
+        dddϕ(t) = ForwardDiff.derivative(ddϕ, t)
+
+        @test ϕ′(0)   === dϕ(0)
+        @test ϕ′′(0)  === ddϕ(0)
+        @test ϕ′′′(0) === dddϕ(0)
+
+        y = interval(1)
+        ψ(t)    =  sin(x + (y+t)*w)
+        ψ′(t)   =  cos(x + (y+t)*w) * w
+        ψ′′(t)  = -sin(x + (y+t)*w) * w * w
+        ψ′′′(t) = -cos(x + (y+t)*w) * w * w * w
+        dψ(t)   = ForwardDiff.derivative(ψ, t)
+        ddψ(t)  = ForwardDiff.derivative(dψ, t)
+        dddψ(t) = ForwardDiff.derivative(ddψ, t)
+        @test ψ′(0) === dψ(0)
+        @test_broken ψ′′(0)  === ddψ(0)
+        @test_broken ψ′′′(0) === dddψ(0)
+    end
 end