Add a test for issue JuliaMath#4

pabloferz · Aug 1, 2017 · 3005a12 · 3005a12
1 parent 21e075b
commit 3005a12
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/test/runtests.jl b/test/runtests.jl
@@ -44,3 +44,10 @@ end
       end
       @test HCubature.countevals(HCubature.Trivial()) == 1
 end
+
+@testset "axischoosing" begin
+    #Issue #4
+    let f = x -> 1.0 + (x[1] * x[3] * sin(x[2]))^2
+        @test hcubature(f, SVector(0,0,-0.2), SVector(0.2,2π,0.2))[1] ≈ (22502//140625)*π
+    end
+end