fix jacobi evaluation (#90)

Author: jishnub

Project.toml

@@ -1,6 +1,6 @@
 name = "ApproxFunOrthogonalPolynomials"
 uuid = "b70543e2-c0d9-56b8-a290-0d4d6d4de211"
-version = "0.4.15"
+version = "0.4.16"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"

src/Spaces/Jacobi/JacobiOperators.jl

@@ -90,7 +90,7 @@ end
 
 function getindex(op::ConcreteEvaluation{<:Jacobi,<:Number},kr::AbstractRange)
     @assert op.order == 0
-    jacobip(eltype(op),kr-1,op.space.a,op.space.b,tocanonical(domain(op),op.x))
+    jacobip(eltype(op),kr.-1,op.space.a,op.space.b,tocanonical(domain(op),op.x))
 end
 
 

test/JacobiTest.jl

@@ -28,6 +28,13 @@ import ApproxFunOrthogonalPolynomials: jacobip
         @test norm((Fun(f,sp2)-f2).coefficients)<100eps()
         @test norm((Fun(f,sp3)-f3).coefficients)<100eps()
         @test norm((Fun(f,sp4)-f4).coefficients)<200eps()
+
+        f = Fun(Jacobi(0,0), Float64[1,2,3])
+        for x in [-1, 0, 1]
+            g = Evaluation(x) * f
+            @test ncoefficients(g) == 1
+            @test coefficients(g)[1] == f(x)
+        end
     end
 
     @testset "Conversion" begin