Fix missing show method

* Fix issue #103

* Fix other test

* Add interval into printing for SemiclassicalJacobi, update README

* typo

* bump

Author: DanielVandH

Project.toml

@@ -1,7 +1,7 @@
 name = "SemiclassicalOrthogonalPolynomials"
 uuid = "291c01f3-23f6-4eb6-aeb0-063a639b53f2"
 authors = ["Sheehan Olver <[email protected]>"]
-version = "0.5.6"
+version = "0.5.7"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

README.md

@@ -11,8 +11,9 @@ This builds on top of [ClassicalOrthogonalPolynomials.jl](https://github.com/Jul
 
 For example, the following gives a half-range Chebyshev polynomial:
 ```julia
+julia> using SemiclassicalOrthogonalPolynomials
 julia> T = SemiclassicalJacobi(2, -1/2, 0, -1/2)
-SemiclassicalJacobi with weight x^-0.5 * (1-x)^0.0 * (2.0-x)^-0.5
+SemiclassicalJacobi with weight x^-0.5 * (1-x)^0.0 * (2.0-x)^-0.5 on 0..1
 
 julia> T[0.1,1:10] # evaluate first 10 OPs at 0.1
 10-element Array{Float64,1}:
@@ -28,5 +29,5 @@ julia> T[0.1,1:10] # evaluate first 10 OPs at 0.1
  -1.1219700834800677
 
 julia> U = SemiclassicalJacobi(2, 1/2, 0, 1/2, T) # last argument reuses computation from T
-SemiclassicalJacobi with weight x^0.5 * (1-x)^0.0 * (2.0-x)^0.5
+SemiclassicalJacobi with weight x^0.5 * (1-x)^0.0 * (2.0-x)^0.5 on 0..1
 ```

src/SemiclassicalOrthogonalPolynomials.jl

@@ -3,7 +3,7 @@ using ClassicalOrthogonalPolynomials: WeightedOPLayout
 using ClassicalOrthogonalPolynomials, FillArrays, LazyArrays, ArrayLayouts, QuasiArrays, InfiniteArrays, ContinuumArrays, LinearAlgebra, BandedMatrices,
         SpecialFunctions, HypergeometricFunctions, InfiniteLinearAlgebra
 
-import Base: getindex, axes, size, \, /, *, +, -, summary, ==, copy, sum, unsafe_getindex, convert, OneTo, diff
+import Base: getindex, axes, size, \, /, *, +, -, summary, show, ==, copy, sum, unsafe_getindex, convert, OneTo, diff
 
 import ArrayLayouts: MemoryLayout, ldiv, diagonaldata, subdiagonaldata, supdiagonaldata
 import BandedMatrices: bandwidths, AbstractBandedMatrix, BandedLayout, _BandedMatrix
@@ -53,9 +53,10 @@ function sum(P::SemiclassicalJacobiWeight{T}) where T
     return abs(convert(T, t^c*exp(loggamma(a+1)+loggamma(b+1)-loggamma(a+b+2)) * pFq((a+1,-c),(a+b+2, ), 1/t)))
 end
 
+show(io::IO, P::SemiclassicalJacobiWeight) = summary(io, P)
 function summary(io::IO, P::SemiclassicalJacobiWeight)
     t,a,b,c = P.t,P.a,P.b,P.c
-    print(io, "x^$a * (1-x)^$b * ($t-x)^$c")
+    print(io, "x^$a * (1-x)^$b * ($t-x)^$c on 0..1")
 end
 
 function ==(A::SemiclassicalJacobiWeight, B::SemiclassicalJacobiWeight)
@@ -219,11 +220,13 @@ axes(P::SemiclassicalJacobi{T}) where T = (Inclusion(UnitInterval{T}()),OneToInf
 
 orthogonalityweight(P::SemiclassicalJacobi) = SemiclassicalJacobiWeight(P.t, P.a, P.b, P.c)
 
+show(io::IO, P::SemiclassicalJacobi) = summary(io, P)
 function summary(io::IO, P::SemiclassicalJacobi)
     t,a,b,c = P.t,P.a,P.b,P.c
-    print(io, "SemiclassicalJacobi with weight x^$a * (1-x)^$b * ($t-x)^$c")
+    print(io, "SemiclassicalJacobi with weight x^$a * (1-x)^$b * ($t-x)^$c on 0..1")
 end
 
+
 jacobimatrix(P::SemiclassicalJacobi) = P.X
 
 """

test/runtests.jl

@@ -86,7 +86,7 @@ end
         W = SemiclassicalJacobi(2,  1/2, 0, -1/2, T)
         w_T = orthogonalityweight(T)
         w_W = orthogonalityweight(W)
-        @test summary(SemiclassicalJacobiWeight(2,1,1,1)) == "x^1 * (1-x)^1 * (2-x)^1"
+        @test summary(SemiclassicalJacobiWeight(2,1,1,1)) == "x^1 * (1-x)^1 * (2-x)^1 on 0..1"
         @test WeightedSemiclassicalJacobi{Float64}(2, -1/2, 0, 1/2) == SemiclassicalJacobiWeight{Float64}(2., -1/2, 0., 1/2) .* SemiclassicalJacobi{Float64}(2., -1/2, 0., 1/2)
 
         @test copy(T) == T
@@ -610,5 +610,14 @@ end
     end
 end
 
+@testset "Issue #103" begin 
+    w = SemiclassicalJacobiWeight(2.0, 2.3, 5.3, 0.4)
+    P = SemiclassicalJacobi(2.0, 2.3, 5.3, 0.4)
+    W = Weighted(P)
+    @test sprint(show, w) == "x^2.3 * (1-x)^5.3 * (2.0-x)^0.4 on 0..1"
+    @test sprint(show, P) == "SemiclassicalJacobi with weight x^2.3 * (1-x)^5.3 * (2.0-x)^0.4 on 0..1"
+    @test sprint(show, W) == "Weighted(SemiclassicalJacobi with weight x^2.3 * (1-x)^5.3 * (2.0-x)^0.4 on 0..1)"
+end
+
 include("test_derivative.jl")
 include("test_neg1c.jl")