backport: isnegative in intpow (#591)

* backport: isnegative in intpow

* Bump version to v0.8.58

* Add UInt test

Author: jishnub

Project.toml

@@ -1,6 +1,6 @@
 name = "ApproxFunBase"
 uuid = "fbd15aa5-315a-5a7d-a8a4-24992e37be05"
-version = "0.8.57"
+version = "0.8.58"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"

src/Fun.jl

@@ -594,6 +594,10 @@ end
 
 \(c::Number, f::Fun) = Fun(f.space, c \ f.coefficients)
 
+# eliminate the type-unstable 1/t branch by using an unsigned integer exponent
+isnegative(x) = x < zero(x)
+isnegative(::Unsigned) = false
+
 @static if VERSION >= v"1.8"
     Base.@constprop :aggressive intpow(f::Fun, k::Integer) = _intpow(f, k)
 else
@@ -612,10 +616,10 @@ end
         f * f * f * f
     else
         t = foldl(*, fill(f, abs(k)-1), init=f)
-        if k > 0
-            return t
-        else
+        if isnegative(k)
             return 1/t
+        else
+            return t
         end
     end
 end

test/SpacesTest.jl

@@ -70,6 +70,7 @@ using LinearAlgebra
                 @test ApproxFunBase.intpow(f, n) == f^n == reduce(*, fill(f, n))
             end
             @test ApproxFunBase.intpow(f,-2) == f^-2 == 1/(f*f)
+            @test f^2 == f^UInt(2)
 
             if VERSION >= v"1.8"
                 @test (@inferred (x -> x^1)(sp)) == sp