ContinuumArrays v0.16 (#164)

* ContinuumArrays v0.16

* add tests

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "MultivariateOrthogonalPolynomials"
 uuid = "4f6956fd-4f93-5457-9149-7bfc4b2ce06d"
-version = "0.6"
+version = "0.6.1"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -29,7 +29,7 @@ BandedMatrices = "0.17.30"
 BlockArrays = "0.16.14"
 BlockBandedMatrices = "0.12.5"
 ClassicalOrthogonalPolynomials = "0.11"
-ContinuumArrays = "0.15"
+ContinuumArrays = "0.15, 0.16"
 DomainSets = "0.6"
 FastTransforms = "0.15"
 FillArrays = "1.0"

src/rect.jl

@@ -16,6 +16,8 @@ function show(io::IO, P::KronPolynomial)
     end
 end
 
+==(A::KronPolynomial, B::KronPolynomial) = length(A.args) == length(B.args) && all(map(==, A.args, B.args))
+
 const RectPolynomial{T, PP} = KronPolynomial{2, T, PP}
 
 

src/rectdisk.jl

@@ -30,6 +30,8 @@ end
 
 DunklXuDiskWeight(β::T) where T = DunklXuDiskWeight{float(T),T}(β)
 
+==(a::DunklXuDiskWeight, b::DunklXuDiskWeight) = a.β == b.β
+
 axes(P::DunklXuDiskWeight{T}) where T = (Inclusion(UnitDisk{T}()),)
 
 show(io::IO, P::DunklXuDiskWeight) = summary(io, P)

test/test_rect.jl

@@ -74,6 +74,9 @@ import ClassicalOrthogonalPolynomials: expand
         W² = RectPolynomial(W, W)
         P² = RectPolynomial(P, P)
         Q² = RectPolynomial(Q, Q)
+
+        @test P² == RectPolynomial(Jacobi(0,0), Jacobi(0,0))
+
         𝐱 = axes(W²,1)
         D_x,D_y = PartialDerivative{1}(𝐱),PartialDerivative{2}(𝐱)
         Δ = Q²\(D_x^2 + D_y^2)*W²

test/test_rectdisk.jl

@@ -24,6 +24,9 @@ import MultivariateOrthogonalPolynomials: dunklxu_raising, dunklxu_lowering, Ang
         WP = WeightedDunklXuDisk(β)
         WQ = WeightedDunklXuDisk(β+1)
 
+        @test WP ≠ WQ
+        @test WP == WP
+
         x, y = first.(axes(P, 1)), last.(axes(P, 1))
 
         L = WP \ WQ