SVector space in Chebyshev Dirichlet (#143)

Author: jishnub

Project.toml

@@ -17,6 +17,7 @@ IntervalSets = "8197267c-284f-5f27-9208-e0e47529a953"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
+StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
 [compat]

src/ApproxFunOrthogonalPolynomials.jl

@@ -93,6 +93,8 @@ import SpecialFunctions: sinpi, cospi, airy, besselh,
                     abs, sign, log, expm1, tan, abs2, sqrt, angle, max, min, cbrt, log,
                     atan, acos, asin, erfc, inv
 
+using StaticArrays: SVector
+
 points(d::IntervalOrSegmentDomain{T},n::Integer) where {T} =
     fromcanonical.(Ref(d), chebyshevpoints(float(real(eltype(T))), n))  # eltype to handle point
 bary(v::AbstractVector{Float64},d::IntervalOrSegmentDomain,x::Float64) = bary(v,tocanonical(d,x))

src/Spaces/Chebyshev/ChebyshevOperators.jl

@@ -117,8 +117,14 @@ function getindex(op::ConcreteEvaluation{<:Chebyshev{DD,RR},typeof(rightendpoint
     scal!(cst,ret)
 end
 
+convert_vector_or_svector(v::AbstractVector) = convert(Vector, v)
+convert_vector_or_svector(t::Tuple) = SVector(t)
+
 @inline function _Dirichlet_Chebyshev(S, order)
-    order == 0 && return ConcreteDirichlet(S,ArraySpace([ConstantSpace.(Point.(endpoints(domain(S))))...]),0)
+    order == 0 && return ConcreteDirichlet(S,
+        ArraySpace(convert_vector_or_svector(
+            ConstantSpace.(Point.(endpoints(domain(S)))))),
+        0)
     default_Dirichlet(S,order)
 end
 @static if VERSION >= v"1.8"