remove literal power special case

Author: dlfivefifty

Project.toml

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

src/ClassicalOrthogonalPolynomials.jl

@@ -212,8 +212,29 @@ gives the singularity structure of an expansion, e.g.,
 """
 singularities(w) = singularities_layout(MemoryLayout(w), w)
 
+
+
 struct NoSingularities end
 
+## default is to just assume no singularities
+singularitiesbroadcast(_...) = NoSingularities()
+
+for op in (:+, :*)
+    @eval singularitiesbroadcast(::typeof($op), A, B, C, D...) = singularitiesbroadcast(*, singularitiesbroadcast(*, A, B), C, D...)
+    @eval singularitiesbroadcast(::typeof($op), ::NoSingularities, ::NoSingularities) = NoSingularities()
+    @eval singularitiesbroadcast(::typeof($op), ::NoSingularities, ::NoSingularities, ::NoSingularities) = NoSingularities()
+end
+
+singularitiesbroadcast(::typeof(*), ::NoSingularities, b) = b
+singularitiesbroadcast(::typeof(*), a, ::NoSingularities) = a
+
+
+# for singularitiesbroadcast(literal_pow), ^, ...)
+# singularitiesbroadcast(F::Function, G::Function, V::SubQuasiArray, K) = singularitiesbroadcast(F, G, parent(V), K)[parentindices(V)...]
+# singularitiesbroadcast(F, V::SubQuasiArray...) = singularitiesbroadcast(F, map(parent,V)...)[parentindices(V...)...]
+# singularitiesbroadcast(F, V::NoSingularities...) = NoSingularities() # default is to assume smooth
+
+
 
 
 basis_axes(ax::Inclusion{<:Any,<:AbstractInterval}, v) = convert(AbstractQuasiMatrix{eltype(v)}, basis_singularities(ax, singularities(v)))

src/classical/jacobi.jl

@@ -94,24 +94,6 @@ singularities(a::AbstractAffineQuasiVector) = singularities(a.x)
 
 singularities(w::AbstractJacobiWeight) = w
 
-## default is to just assume no singularities
-singularitiesbroadcast(_...) = NoSingularities()
-
-for op in (:+, :*)
-    @eval singularitiesbroadcast(::typeof($op), A, B, C, D...) = singularitiesbroadcast(*, singularitiesbroadcast(*, A, B), C, D...)
-    @eval singularitiesbroadcast(::typeof($op), ::NoSingularities, ::NoSingularities) = NoSingularities()
-    @eval singularitiesbroadcast(::typeof($op), ::NoSingularities, ::NoSingularities, ::NoSingularities) = NoSingularities()
-end
-
-singularitiesbroadcast(::typeof(*), ::NoSingularities, b) = b
-singularitiesbroadcast(::typeof(*), a, ::NoSingularities) = a
-
-
-# for singularitiesbroadcast(literal_pow), ^, ...)
-singularitiesbroadcast(F::Function, G::Function, V::SubQuasiArray, K) = singularitiesbroadcast(F, G, parent(V), K)[parentindices(V)...]
-singularitiesbroadcast(F, V::SubQuasiArray...) = singularitiesbroadcast(F, map(parent,V)...)[parentindices(V...)...]
-singularitiesbroadcast(F, V::NoSingularities...) = NoSingularities() # default is to assume smooth
-
 
 abstract type AbstractJacobi{T} <: OrthogonalPolynomial{T} end
 

src/classical/legendre.jl

@@ -74,6 +74,8 @@ singularitiesbroadcast(::typeof(^), L::LegendreWeight, ::NoSingularities) = L
 singularitiesbroadcast(::typeof(/), ::NoSingularities, L::LegendreWeight) = L # can't find roots
 
 basis_singularities(ax::Inclusion, ::NoSingularities) = legendre(ax)
+basis_singularities(ax, sing) = basis_singularities(sing) # fallback for back compatibility
+
 
 """
     Legendre{T=Float64}(a,b)