Use IntervalContractors.jl for reverse functions (#91)

* Replace reverse_mode.jl with using IntervalContractors package

* Make alias for unexported rev_ops

* Remove reverse_mode.jl

* Rename rev_ops -> reverse_operations

* Add IntervalContractors to REQUIRE

Author: dpsanders

REQUIRE

@@ -1,5 +1,6 @@
 julia 0.5
 IntervalArithmetic 0.9
 IntervalRootFinding 0.1
+IntervalContractors 0.1
 MacroTools 0.3
 

src/IntervalConstraintProgramming.jl

@@ -2,7 +2,9 @@ __precompile__()
 
 module IntervalConstraintProgramming
 
-using IntervalArithmetic, IntervalRootFinding
+using   IntervalArithmetic,
+        IntervalRootFinding,
+        IntervalContractors
 
 using MacroTools
 
@@ -20,8 +22,8 @@ export
     Vol,
     show_code
 
+const reverse_operations = IntervalContractors.reverse_operations
 
-include("reverse_mode.jl")
 include("ast.jl")
 include("code_generation.jl")
 include("contractor.jl")

src/ast.jl

@@ -309,7 +309,7 @@ function process_call!(flatAST::FlatAST, ex, new_var=nothing)
 
     #@show op
 
-    if op ∈ keys(rev_ops)  # standard operator
+    if op ∈ keys(reverse_operations)  # standard operator
         if new_var == nothing
             new_var = make_symbol()
         end

src/code_generation.jl

@@ -53,7 +53,7 @@ function emit_backward_code(a::Assignment)
     args = isa(a.args, Vector) ? a.args : [a.args]
 
     return_args = [a.lhs, args...]
-    rev_op = rev_ops[a.op]  # find reverse operation
+    rev_op = reverse_operations[a.op]  # find reverse operation
 
     if rev_op == :()   # empty
         args = make_tuple(args)

src/reverse_mode.jl

@@ -163,9 +163,3 @@ end
 
     @test C3(A, x) == IntervalBox(sqrt(A / 16))
 end
-
-@testset "power_rev for odd power" begin
-    x = -∞..∞
-    a = -8..27
-    power_rev(a, x, 3)[2] == Interval(-2.0000000000000004, 3.0000000000000004)
-end