Overload _getindex instead of getindex for orthogonal polynomials (#251)

* Overload `_getindex` instead of `getindex` for orthogonal polynomials

* Update clenshaw.jl

* Update clenshaw.jl

* v0.15.9

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
-version = "0.15.8"
+version = "0.15.9"
 authors = ["Sheehan Olver <[email protected]>"]
 
 [deps]

src/clenshaw.jl

@@ -51,10 +51,7 @@ function copyto!(dest::AbstractVector, V::SubArray{<:Any,1,<:OrthogonalPolynomia
     dest
 end
 
-getindex(P::OrthogonalPolynomial, x::Number, n::AbstractVector) = layout_getindex(P, x, n)
-getindex(P::OrthogonalPolynomial, x::AbstractVector, n::AbstractVector) = layout_getindex(P, x, n)
-getindex(P::SubArray{<:Any,1,<:OrthogonalPolynomial}, x::AbstractVector) = layout_getindex(P, x)
-Base.@propagate_inbounds function getindex(P::OrthogonalPolynomial, x::Number, n::Number)
+Base.@propagate_inbounds function _getindex(::Type{IND}, P::OrthogonalPolynomial, (x,n)::IND) where IND
     @boundscheck checkbounds(P, x, n)
     Base.unsafe_getindex(P, x, n)
 end