Numbers as derivative orders (#39)

* Numbers as derivative orders

* Add order test

Project.toml

@@ -1,6 +1,6 @@
 name = "ApproxFunSingularities"
 uuid = "f8fcb915-6b99-5be2-b79a-d6dbef8e6e7e"
-version = "0.3.3"
+version = "0.3.4"
 
 [deps]
 ApproxFunBase = "fbd15aa5-315a-5a7d-a8a4-24992e37be05"

src/JacobiWeight.jl

@@ -45,7 +45,7 @@ WeightedJacobi(β,α) = JacobiWeight(β,α,Jacobi(β,α))
 Fun(::typeof(identity), S::JacobiWeight) =
     isapproxinteger(S.β) && isapproxinteger(S.α) ? Fun(x->x,S) : Fun(identity,domain(S))
 
-order(S::JacobiWeight{Ultraspherical{Int,D,R},D,R}) where {D,R} = order(S.space)
+order(S::JacobiWeight{<:Ultraspherical{<:Any,D,R},D,R}) where {D,R} = order(S.space)
 
 
 spacescompatible(A::JacobiWeight,B::JacobiWeight) =

src/JacobiWeightOperators.jl

@@ -186,7 +186,8 @@ end
 
 Derivative(S::JacobiWeight{SS,DDD}) where {SS,DDD<:IntervalOrSegment} = jacobiweightDerivative(S)
 
-function Derivative(S::JacobiWeight{SS,DD}, k::Integer) where {SS,DD<:IntervalOrSegment}
+function Derivative(S::JacobiWeight{SS,DD}, k::Number) where {SS,DD<:IntervalOrSegment}
+    @assert Integer(k) == k "order must be an integer"
     if k==1
         Derivative(S)
     else

test/runtests.jl

@@ -5,6 +5,7 @@ using ApproxFunBase: HeavisideSpace, PointSpace, ArraySpace, DiracSpace, Piecewi
 using ApproxFunBaseTest: testbandedoperator, testtransforms, testfunctional,
                         testbandedblockbandedoperator
 using ApproxFunOrthogonalPolynomials
+using ApproxFunOrthogonalPolynomials: order
 using ApproxFunSingularities
 using LinearAlgebra
 using IntervalSets
@@ -167,6 +168,12 @@ end
         @test norm(Fun(cos,Chebyshev)-Fun(cos,JacobiWeight(0,0,Jacobi(-0.5,-0.5))))<100eps()
         @test norm(Fun(cos,Jacobi(-0.5,-0.5))-Fun(cos,JacobiWeight(0,0,Jacobi(-0.5,-0.5))))<100eps()
     end
+
+    @testset "Ultraspherical order" begin
+        us = Ultraspherical(0.5)
+        s = JacobiWeight(1, 1, us)
+        @test order(s) == order(us)
+    end
 end
 
 @testset "Ray and Line" begin