permutedims, singularitiesbroadcast(*, ...)

Author: dlfivefifty

src/clenshaw.jl

@@ -290,6 +290,8 @@ end
 getindex(M::Clenshaw, k::Int, j::Int) = M[k:k,j][1]
 
 transposelayout(M::ClenshawLayout) = ClenshawLayout()
+# TODO: generalise for layout, use Base.PermutedDimsArray
+Base.permutedims(M::Clenshaw{<:Number}) = transpose(M)
 
 
 function materialize!(M::MatMulVecAdd{<:ClenshawLayout,<:PaddedLayout,<:PaddedLayout})

src/jacobi.jl

@@ -61,7 +61,7 @@ singularities(a::AbstractAffineQuasiVector) = singularities(a.x)
 singularitiesbroadcast(_, L::LegendreWeight) = L # Assume we stay smooth
 singularitiesbroadcast(::typeof(exp), L::LegendreWeight) = L
 singularitiesbroadcast(::typeof(Base.literal_pow), ::typeof(^), L::LegendreWeight, ::Val) = L
-for op in (:+, :-)
+for op in (:+, :-, :*)
     @eval begin
         singularitiesbroadcast(::typeof($op), ::LegendreWeight{T}, ::LegendreWeight{V}) where {T,V} = LegendreWeight{promote_type(T,V)}()
         singularitiesbroadcast(::typeof($op), L::LegendreWeight, ::NoSingularities) = L