Update show

Author: dlfivefifty

Project.toml

@@ -29,7 +29,7 @@ ArrayLayouts = "1.0.1"
 BandedMatrices = "0.17.17"
 BlockArrays = "0.16.9"
 BlockBandedMatrices = "0.12"
-ContinuumArrays = "0.12.4"
+ContinuumArrays = "0.13"
 DomainSets = "0.5.6, 0.6"
 FFTW = "1.1"
 FastGaussQuadrature = "0.4.3, 0.5"
@@ -41,7 +41,7 @@ InfiniteLinearAlgebra = "0.6.16"
 IntervalSets = "0.5, 0.6, 0.7"
 LazyArrays = "1.0.1"
 LazyBandedMatrices = "0.8.5"
-QuasiArrays = "0.9.6"
+QuasiArrays = "0.10"
 SpecialFunctions = "1.0, 2"
 julia = "1.7"
 

src/ClassicalOrthogonalPolynomials.jl

@@ -161,7 +161,7 @@ gives the singularity structure of an expansion, e.g.,
 `JacobiWeight`.
 """
 singularities(::WeightLayout, w) = w
-singularities(lay::BroadcastLayout, a) = singularitiesbroadcast(call(a), map(singularities, arguments(lay, a))...)
+singularities(lay::BroadcastLayout, a::AbstractQuasiVector) = singularitiesbroadcast(call(a), map(singularities, arguments(lay, a))...)
 singularities(::WeightedBasisLayouts, a) = singularities(BroadcastLayout{typeof(*)}(), a)
 singularities(::WeightedOPLayout, a) = singularities(weight(a))
 singularities(w) = singularities(MemoryLayout(w), w)

src/classical/jacobi.jl

@@ -156,7 +156,8 @@ axes(::AbstractJacobi{T}) where T = (Inclusion{T}(ChebyshevInterval{real(T)}()),
 ==(A::Legendre, B::Weighted{<:Any,<:AbstractJacobi}) = A == B.P
 ==(A::Weighted{<:Any,<:AbstractJacobi}, B::Legendre) = A.P == B
 
-
+show(io::IO, w::AbstractJacobiWeight) = summary(io, w)
+show(io::IO, P::AbstractJacobi) = summary(io, P)
 summary(io::IO, P::Jacobi) = print(io, "Jacobi($(P.a), $(P.b))")
 
 ###

src/classical/legendre.jl

@@ -35,8 +35,8 @@ singularitiesbroadcast(::typeof(^), L::LegendreWeight, ::NoSingularities) = L
 singularitiesbroadcast(::typeof(/), ::NoSingularities, L::LegendreWeight) = L # can't find roots
 
 singularities(::AbstractJacobi{T}) where T = LegendreWeight{T}()
-singularities(::Inclusion{T,<:AbstractInterval}) where T = LegendreWeight{T}()
-singularities(d::Inclusion{T,<:Interval}) where T = LegendreWeight{T}()[affine(d,ChebyshevInterval{T}())]
+singularities(::Inclusion{T,<:ChebyshevInterval}) where T = LegendreWeight{T}()
+singularities(d::Inclusion{T,<:AbstractInterval}) where T = LegendreWeight{T}()[affine(d,ChebyshevInterval{T}())]
 singularities(::AbstractFillLayout, P) = LegendreWeight{eltype(P)}()
 
 _basis(::LegendreWeight{T}) where T = Legendre{T}()
@@ -166,22 +166,3 @@ function _sum(P::Legendre{T}, dims) where T
 end
 
 _sum(p::SubQuasiArray{T,1,Legendre{T},<:Tuple{Inclusion,Int}}, ::Colon) where T = parentindices(p)[2] == 1 ? convert(T, 2) : zero(T)
-
-
-###
-# dot
-###
-
-_dot(::Inclusion{<:Any,<:AbstractInterval}, a, b) = __dot(MemoryLayout(a), MemoryLayout(b), a, b)
-function __dot(::ExpansionLayout, ::ExpansionLayout, a, b)
-    P,c = basis(a),coefficients(a)
-    Q,d = basis(b),coefficients(b)
-    c' * (P'Q) * d
-end
-__dot(_, _, a, b) = dot(expand(a), expand(b))
-__dot(::ZerosLayout, _, a, b) = __dot(ZerosLayout(), ZerosLayout(), a, b)
-__dot(_, ::ZerosLayout, a, b) = __dot(ZerosLayout(), ZerosLayout(), a, b)
-function __dot(::ZerosLayout, ::ZerosLayout, a, b)
-    axes(a) == axes(b) || throw(DimensionMismatch("axes must match"))
-    zero(promote_type(eltype(a), eltype(b)))
-end

src/normalized.jl

@@ -298,7 +298,7 @@ function mul(A::Derivative, B::Weighted{<:Any,<:SubQuasiArray{<:Any,2,<:Abstract
 end
 mul(Ac::QuasiAdjoint{<:Any, Weighted{<:Any,<:SubQuasiArray{<:Any,2,<:AbstractQuasiMatrix,<:Tuple{<:AbstractAffineQuasiVector,<:Any}}}}, Bc::QuasiAdjoint{<:Any,<:Derivative}) = mul(Bc', Ac')'    
 
-summary(io::IO, Q::Weighted) = print(io, "Weighted($(Q.P))")
+show(io::IO, Q::Weighted) = print(io, "Weighted($(Q.P))")
 
 __sum(::AbstractNormalizedOPLayout, A, dims) = __sum(ApplyLayout{typeof(*)}(), A, dims)
 function __sum(::WeightedOPLayout, A, dims)

test/runtests.jl

@@ -60,12 +60,16 @@ end
 
 @testset "basis" begin
     for x in (Inclusion(ChebyshevInterval()), Inclusion(1 .. 2))
-        @test sum(x) == last(x)-first(x)
-        # sum(x .^ 2)
-        # sum(exp.(x))
+        a,b = first(x),last(x)
+        @test sum(x) == b-a
+        @test sum(x .^ 2) ≈ (b^3 - a^3)/3
+        @test sum(exp.(x)) ≈ exp(b) - exp(a)
         @test dot(x, x) ≈ sum(expand(x .^2))
         @test dot(x.^2, x.^2) ≈ sum(expand(x .^4))
         @test dot(exp.(x), x.^2) ≈ sum(expand(x .^2 .* exp.(x)))
         @test dot(x, exp.(x)) ≈ dot(exp.(x), x)
     end
+
+    # A = x .^ (0:2)'
+    # sum(A; dims=1)
 end