Tests pass

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "MultivariateOrthogonalPolynomials"
 uuid = "4f6956fd-4f93-5457-9149-7bfc4b2ce06d"
-version = "0.2.3"
+version = "0.2.4"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -28,17 +28,17 @@ ArrayLayouts = "0.7"
 BandedMatrices = "0.16"
 BlockArrays = "0.16.7"
 BlockBandedMatrices = "0.11"
-ClassicalOrthogonalPolynomials = "0.4.2"
-ContinuumArrays = "0.8, 0.9"
+ClassicalOrthogonalPolynomials = "0.5"
+ContinuumArrays = "0.10"
 DomainSets = "0.5"
-FastTransforms = "0.12"
-FillArrays = "0.11, 0.12"
+FastTransforms = "0.13"
+FillArrays = "0.12"
 HarmonicOrthogonalPolynomials = "0.2"
 InfiniteArrays = "0.12"
 InfiniteLinearAlgebra = "0.6"
-LazyArrays = "0.21.8, 0.22"
+LazyArrays = "0.22"
 LazyBandedMatrices = "0.7"
-QuasiArrays = "0.8, 0.9"
+QuasiArrays = "0.9"
 SpecialFunctions = "1"
 StaticArrays = "1"
 julia = "1.6"

src/MultivariateOrthogonalPolynomials.jl

@@ -6,11 +6,11 @@ using ClassicalOrthogonalPolynomials, FastTransforms, BlockBandedMatrices, Block
       LazyArrays, SpecialFunctions, LinearAlgebra, BandedMatrices, LazyBandedMatrices, ArrayLayouts,
       HarmonicOrthogonalPolynomials
 
-import Base: axes, in, ==, *, ^, \, copy, OneTo, getindex, size, oneto, all
+import Base: axes, in, ==, *, ^, \, copy, OneTo, getindex, size, oneto, all, resize!
 import DomainSets: boundary
 
 import QuasiArrays: LazyQuasiMatrix, LazyQuasiArrayStyle
-import ContinuumArrays: @simplify, Weight, grid, plotgrid, TransformFactorization, Expansion
+import ContinuumArrays: @simplify, Weight, grid, plotgrid, TransformFactorization, ExpansionLayout
 
 import ArrayLayouts: MemoryLayout, sublayout, sub_materialize
 import BlockArrays: block, blockindex, BlockSlice, viewblock
@@ -20,7 +20,7 @@ import LazyArrays: arguments, paddeddata, LazyArrayStyle, LazyLayout
 import LazyBandedMatrices: LazyBandedBlockBandedLayout, AbstractBandedBlockBandedLayout, AbstractLazyBandedBlockBandedLayout
 import InfiniteArrays: InfiniteCardinal
 
-import ClassicalOrthogonalPolynomials: jacobimatrix, Weighted, orthogonalityweight, HalfWeighted
+import ClassicalOrthogonalPolynomials: jacobimatrix, Weighted, orthogonalityweight, HalfWeighted, WeightedBasis
 import HarmonicOrthogonalPolynomials: BivariateOrthogonalPolynomial, MultivariateOrthogonalPolynomial, Plan,
                                           PartialDerivative, AngularMomentum, BlockOneTo, BlockRange1, interlace
 

src/triangle.jl

@@ -569,7 +569,7 @@ getindex(P::JacobiTriangle, xy::SVector{2}, JR::BlockOneTo) =
 # Clenshaw
 ###
 
-function getindex(f::Expansion{T,<:JacobiTriangle}, xy::SVector{2}) where T
+function getindex(f::ApplyVector{T,typeof(*),<:Tuple{JacobiTriangle,AbstractVector}}, xy::SVector{2}) where T
     P,c∞ = arguments(f)
     c = paddeddata(c∞)
     N = blocksize(c,1)
@@ -605,7 +605,7 @@ function getindex(f::Expansion{T,<:JacobiTriangle}, xy::SVector{2}) where T
     γ1[1]
 end
 
-getindex(f::Expansion{T,<:JacobiTriangle}, xys::AbstractArray{<:SVector{2}}) where T =
+getindex(f::ApplyVector{T,typeof(*),<:Tuple{JacobiTriangle,AbstractVector}}, xys::AbstractArray{<:SVector{2}}) where T =
     [f[xy] for xy in xys]
 
 # getindex(f::Expansion{T,<:JacobiTriangle}, x::AbstractVector{<:Number}) where T =