Use broadcast

Author: dlfivefifty

src/ChebyshevJacobiPlan.jl

@@ -140,7 +140,7 @@ function ForwardChebyshevJacobiPlan{T}(c_jac::AbstractVector{T},α::T,β::T,M::I
     @inbounds for i=1:N+1 tempcosβsinα[i] = tempcos[i]^(β+1/2)*tempsin[i]^(α+1/2) end
 
     # Initialize normalizing constant
-    cnαβ = Cnαβ(0:N,α,β)
+    cnαβ = Cnαβ.(0:N,α,β)
     cnmαβ = zero(cnαβ)
 
     # Get indices
@@ -182,11 +182,11 @@ function BackwardChebyshevJacobiPlan{T}(c_cheb::AbstractVector{T},α::T,β::T,M:
     @inbounds for i=1:2N+1 tempcosβsinα[i] = tempcos[i]^(β+1/2)*tempsin[i]^(α+1/2) end
 
     # Initialize normalizing constant
-    cnαβ = Cnαβ(0:2N,α,β)
+    cnαβ = Cnαβ.(0:2N,α,β)
     cnmαβ = zero(cnαβ)
 
     # Initialize orthonormality constants
-    anαβ = Anαβ(0:N,α,β)
+    anαβ = Anαβ.(0:N,α,β)
 
     # Get indices
     CJI = ChebyshevJacobiIndices(α,β,CJC,tempmindices,cfs,tempcos,tempsin,tempcosβsinα)

src/ChebyshevUltrasphericalPlan.jl

@@ -139,7 +139,7 @@ function ForwardChebyshevUltrasphericalPlan{T}(c_ultra::AbstractVector{T},λ::T,
     tempsinλm = similar(tempsinλ)
 
     # Initialize normalizing constant
-    cnλ = Cnλ(0:N,λ)
+    cnλ = Cnλ.(0:N,λ)  # broadcast is overloaded
     cnmλ = similar(cnλ)
 
     # Get indices
@@ -180,7 +180,7 @@ function BackwardChebyshevUltrasphericalPlan{T}(c_ultra::AbstractVector{T},λ::T
     tempsinλm = similar(tempsinλ)
 
     # Initialize normalizing constant
-    cnλ = Cnλ(0:2N,λ)
+    cnλ = Cnλ.(0:2N,λ) # broadcast is overloaded
     cnmλ = similar(cnλ)
 
     # Initialize orthonormality constants

test/runtests.jl

@@ -8,7 +8,7 @@ n = 0:1000_000
 @time FastTransforms.Cnλ.(n,λ);
 
 x = linspace(0,20,81);
-@test norm((FastTransforms.Λ.(x)-FastTransforms.Λ.(big.(x)))./FastTransforms.Λ.(big.(x)),Inf) < 2eps()
+@test norm((FastTransforms.Λ.(x).-FastTransforms.Λ.(big.(x)))./FastTransforms.Λ.(big.(x)),Inf) < 2eps()
 
 x = 0:0.5:10_000
 λ₁,λ₂ = 0.125,0.875
@@ -20,11 +20,11 @@ n = 0:1000
 α = 0.125
 β = 0.375
 
-@time FastTransforms.Cnαβ(n,α,β);
-@test norm(FastTransforms.Cnαβ(n,α,β)./FastTransforms.Cnαβ(n,big(α),big(β))-1,Inf) < 3eps()
+@time FastTransforms.Cnαβ.(n,α,β);
+@test norm(FastTransforms.Cnαβ.(n,α,β)./FastTransforms.Cnαβ.(n,big(α),big(β))-1,Inf) < 3eps()
 
-@time FastTransforms.Anαβ(n,α,β);
-@test norm(FastTransforms.Anαβ(n,α,β)./FastTransforms.Anαβ(n,big(α),big(β))-1,Inf) < 4eps()
+@time FastTransforms.Anαβ.(n,α,β);
+@test norm(FastTransforms.Anαβ.(n,α,β)./FastTransforms.Anαβ.(n,big(α),big(β))-1,Inf) < 4eps()
 
 
 N = 20