Merge pull request #207 from araujoms/main

fix power with negative integers

Author: JeffreySarnoff

.gitignore

@@ -1,5 +1,6 @@
+Manifest.toml
 *.jl.cov
 *.jl.*.cov
 *.jl.mem
 deps/deps.jl
-.DS_Store
+.DS_Store

src/math/elementary/explog.jl

@@ -123,13 +123,13 @@ function Base.:(^)(r::DoubleFloat{T}, n::DoubleFloat{T}) where {T <: IEEEFloat}
     if isinteger(n)
         return r^Int64(Float64(n)) # convert n to Float64 first
     else
-       return exp(n * log(r))
+        return exp(n * log(r))
     end
 end
 
 function Base.:(^)(r::Int, n::DoubleFloat{T}) where {T<:IEEEFloat}
     if isinteger(n)
-        return r^Int64(Float64(n)) # convert n to Float64 first
+        return DoubleFloat{T}(r)^Int64(Float64(n)) # convert n to Float64 first
     else
         return exp(n * log(r))
     end

test/functions.jl

@@ -17,6 +17,13 @@ end
     @test Double64(10.0)^0 == Double64(1.0)
     @test Double64(10.0)^1 == Double64(10.0)
     @test Double64(0.0)^0 == Double64(1.0)
+    @test 2^Double64(-3) == 0.125
+    @test isa(2^Double16(1), Double16)
+    @test isa(2^Double16(-1), Double16)
+    @test isa(2^Double32(1), Double32)
+    @test isa(2^Double32(-1), Double32)
+    @test isa(2^Double64(1), Double64)
+    @test isa(2^Double64(-1), Double64)
 end
 
 f64a = sqrt(2.0)