Fix primitive /

Author: tansongchen

benchmark/case1.jl

@@ -1,4 +1,4 @@
-using TaylorDiff, BenchmarkTools
+using TaylorSeries, TaylorDiff, BenchmarkTools
 
 function f(t, p, α, s)
     x = 1.0 / (1.0 - s*(t+1)/(t-1))
@@ -16,5 +16,7 @@ p, α, s = (rand(N), rand(N), rand())
 p ./= sum(p)
 
 # Test
-t = TaylorScalar{Float64, 20}(0.)
+t0 = Taylor1(eltype(p),m)
+t = TaylorScalar{Float64, 21}(0.)
+@btime f(t0, p, α, s);
 @btime f(t, p, α, s);

src/primitive.jl

@@ -138,7 +138,7 @@ end
         ex = quote
             $ex
             $(Symbol('v', i)) = (va[$i] - +($([
-                :($(binomial(i - 1, j)) * $(Symbol('v', j)) * vb[$i + 1 - $j])
+                :($(binomial(i - 1, j - 1)) * $(Symbol('v', j)) * vb[$i + 1 - $j])
             for j = 1:i-1]...))) / vb[1]
         end
     end