Fix a precision bug (#5)

* Update PolynomialSpace.jl * Update SpacesTest.jl
JuliaApproximation · Jul 2, 2019 · 77904ee · 77904ee
1 parent 51f505e
commit 77904ee
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/src/Spaces/PolynomialSpace.jl b/src/Spaces/PolynomialSpace.jl
@@ -265,7 +265,7 @@ clenshaw(sp::PolynomialSpace,c::AbstractMatrix,x) = clenshaw(c,x,ClenshawPlan(pr
 clenshaw!(sp::PolynomialSpace,c::AbstractVector,x::AbstractVector)=clenshaw!(c,x,ClenshawPlan(promote_type(eltype(c),eltype(x)),sp,length(x)))
 
 function clenshaw(sp::PolynomialSpace,c::AbstractVector,x)
-    N,T = length(c),promote_type(eltype(c),typeof(x))
+    N,T = length(c),promote_type(prectype(sp),eltype(c),typeof(x))
     TT = eltype(T)
     if isempty(c)
         return zero(T)

diff --git a/test/SpacesTest.jl b/test/SpacesTest.jl
@@ -6,6 +6,9 @@ import ApproxFunBase: space, SpaceOperator,
                     testfunctional
 
 @testset "Spaces" begin
+    @testset "evaluate" begin
+        @test evaluate([1],Legendre(),0)==evaluate([1.0],Legendre(),0)
+    end
     @testset "ChebyshevDirichlet" begin
         testtransforms(ChebyshevDirichlet{1,1}())