UInt exponents in JacobiWeight defaultConversion (#90)

Author: jishnub

src/JacobiWeightOperators.jl

@@ -306,23 +306,26 @@ conversion_rule(A::JacobiWeight, B::Space{<:IntervalOrSegmentDomain}) =
 function defaultConversion(A::JacobiWeight{<:Any,<:IntervalOrSegmentDomain},
         B::JacobiWeight{<:Any,<:IntervalOrSegmentDomain})
 
-    @assert isapproxinteger(A.β-B.β) && isapproxinteger(A.α-B.α)
+    Aβ, Aα = A.β, A.α
+    Bβ, Bα = B.β, B.α
 
-    if isapprox(A.β,B.β) && isapprox(A.α,B.α)
+    @assert isapproxinteger(Aβ-Bβ) && isapproxinteger(Aα-Bα)
+
+    if isapprox(Aβ,Bβ) && isapprox(Aα,Bα)
         ConversionWrapper(SpaceOperator(Conversion(A.space,B.space),A,B))
     else
-        @assert A.β≥B.β && A.α≥B.α
+        @assert Aβ≥Bβ && Aα≥Bα
         # first check if a multiplication by JacobiWeight times B.space is overloaded
         # this is currently designed for Jacobi multiplied by (1-x), etc.
-        βdif,αdif=round(Int,A.β-B.β),round(Int,A.α-B.α)
+        βdif,αdif=round(Int,Aβ-Bβ),round(Int,Aα-Bα)
         d=domain(A)
         M=Multiplication(jacobiweight(βdif,αdif,d), A.space)
 
         if rangespace(M)==JacobiWeight(βdif,αdif,A.space)
             # M is the default, so we should use multiplication by polynomials instead
             x=Fun(identity,d)
             y=mobius(d,x)   # we use mobius instead of tocanonical so that it works for Funs
-            m=(1+y)^βdif*(1-y)^αdif
+            m = (1+y)^UInt(βdif) * (1-y)^UInt(αdif)
             MC=promoterangespace(Multiplication(m,A.space), B.space)
 
             ConversionWrapper(SpaceOperator(MC,A,B))# Wrap the operator with the correct spaces