Issue with `LanczosPolynomial` After Julia Update

[JiajieYao1106]

I just updated Julia to the latest version (v"1.11.5") and attempted to run a script I had previously written for testing purposes. While it worked as expected in Julia v1.10, it now throws an error when using the \ command with LanczosPolynomial in Julia v1.11. I attempted to follow the suggested fix by adding a method definition for copy(::Ldiv{...}), as indicated in the error message, but the error still persists.

(@v1.11) pkg> st ClassicalOrthogonalPolynomials
Status `C:\Users\jjyao\.julia\environments\v1.11\Project.toml`
  [b30e2e7b] ClassicalOrthogonalPolynomials v0.15.2

julia> using ClassicalOrthogonalPolynomials, LinearAlgebra, SparseArrays, BandedMatrices, LazyBandedMatrices

julia> x = Inclusion(0..1)
Inclusion(0 .. 1)

julia> Q1 = LanczosPolynomial(@. (1-x) * x * sqrt(2+x))
LanczosPolynomial{Float64} with weight with singularities LegendreWeight{Float64}() affine mapped to 0 .. 1

julia> Q2 = LanczosPolynomial(@. (1-x)^2 * x^2 * (sqrt(2+x))^3)
LanczosPolynomial{Float64} with weight with singularities LegendreWeight{Float64}() affine mapped to 0 .. 1

julia> D = Q2 \  (Derivative(x) * Q1)
ERROR: MethodError: copy(::ArrayLayouts.Ldiv{…}) is ambiguous.

julia> show(err)
1-element ExceptionStack:
MethodError: copy(::ArrayLayouts.Ldiv{BandedMatrices.BandedColumns{LazyArrays.LazyLayout}, LazyArrays.ApplyLayout{typeof(*)}, BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}, Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}}}) is ambiguous.

Candidates:
  copy(L::ArrayLayouts.Ldiv{<:Union{LazyArrays.AbstractLazyBandedLayout, ArrayLayouts.BidiagonalLayout{<:LazyArrays.AbstractLazyLayout}, ArrayLayouts.HermitianLayout{<:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.HermitianLayout{<:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.HermitianLayout{<:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.SymTridiagonalLayout{<:LazyArrays.AbstractLazyLayout}, ArrayLayouts.SymmetricLayout{<:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.SymmetricLayout{<:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.SymmetricLayout{<:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'N', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'U', 'N', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'L', 'U', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'U', 'U', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'L', 'N', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'N', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'U', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'U', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'N', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'N', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'U', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'U', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TridiagonalLayout{<:LazyArrays.AbstractLazyLayout}, BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}, BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}})
    @ LazyArraysBandedMatricesExt C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\ext\LazyArraysBandedMatricesExt.jl:595
  copy(L::ArrayLayouts.Ldiv{<:Union{LazyArrays.AbstractLazyBandedLayout, ArrayLayouts.BidiagonalLayout{<:LazyArrays.AbstractLazyLayout}, ArrayLayouts.HermitianLayout{<:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.HermitianLayout{<:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.HermitianLayout{<:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.SymTridiagonalLayout{<:LazyArrays.AbstractLazyLayout}, ArrayLayouts.SymmetricLayout{<:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.SymmetricLayout{<:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.SymmetricLayout{<:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'N', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'U', 'N', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'L', 'U', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'U', 'U', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'L', 'N', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'N', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'U', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'U', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'N', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'N', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'U', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'U', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TridiagonalLayout{<:LazyArrays.AbstractLazyLayout}, BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}, BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, <:LazyArrays.AbstractLazyLayout})
    @ LazyArraysBandedMatricesExt C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\ext\LazyArraysBandedMatricesExt.jl:596
  copy(L::ArrayLayouts.Ldiv{<:Any, LazyArrays.ApplyLayout{typeof(*)}})
    @ LazyArrays C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\linalg\inv.jl:119
  copy(L::ArrayLayouts.Ldiv{<:Any, <:LazyArrays.AbstractLazyLayout})
    @ LazyArrays C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\linalg\inv.jl:129

Possible fix, define
  copy(::ArrayLayouts.Ldiv{<:Union{LazyArrays.AbstractLazyBandedLayout, ArrayLayouts.BidiagonalLayout{<:LazyArrays.AbstractLazyLayout}, ArrayLayouts.HermitianLayout{<:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.HermitianLayout{<:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.HermitianLayout{<:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.SymTridiagonalLayout{<:LazyArrays.AbstractLazyLayout}, ArrayLayouts.SymmetricLayout{<:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.SymmetricLayout{<:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.SymmetricLayout{<:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'N', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'U', 'N', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'L', 'U', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'U', 'U', <:LazyArrays.AbstractLazyBandedLayout}, ArrayLayouts.TriangularLayout{'L', 'N', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'N', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'U', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'U', <:BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'N', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'N', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'L', 'U', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TriangularLayout{'U', 'U', <:BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, ArrayLayouts.TridiagonalLayout{<:LazyArrays.AbstractLazyLayout}, BandedMatrices.BandedColumns{<:LazyArrays.AbstractLazyLayout}, BandedMatrices.BandedRows{<:LazyArrays.AbstractLazyLayout}}, LazyArrays.ApplyLayout{typeof(*)}})

Stacktrace:
  [1] materialize
    @ C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\ldiv.jl:22 [inlined]
  [2] ldiv
    @ C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\ldiv.jl:98 [inlined]
  [3] \(x::BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, A::LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}, Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}}; kwds::@Kwargs{})
    @ ArrayLayouts C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\ldiv.jl:183
  [4] copy(M::ArrayLayouts.Mul{LazyArrays.ApplyLayout{typeof(\)}, ArrayLayouts.UnknownLayout, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}}}, LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}}})
    @ LazyArrays C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\linalg\inv.jl:146
  [5] copy(L::ArrayLayouts.Mul{LazyArrays.ApplyLayout{typeof(\)}, LazyArrays.ApplyLayout{typeof(*)}, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}}}, LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}}})
    @ LazyArrays C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\linalg\inv.jl:148
  [6] materialize
    @ C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:137 [inlined]
  [7] mul
    @ C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:138 [inlined]
  [8] *(A::LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}}}, B::LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}})
    @ ArrayLayouts C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:223
  [9] _broadcast_mul_mul(::typeof(\), ::Tuple{LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}}}}, C::LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}})
    @ LazyArrays C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\lazybroadcasting.jl:412
 [10] copy
    @ C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\lazybroadcasting.jl:432 [inlined]
 [11] materialize(M::ArrayLayouts.Mul{LazyArrays.BroadcastLayout{typeof(\)}, LazyArrays.ApplyLayout{typeof(*)}, LazyArrays.BroadcastMatrix{Float64, typeof(\), Tuple{LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}}}}}, LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}}})
    @ ArrayLayouts C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:137
 [12] mul
    @ C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:138 [inlined]
 [13] *
    @ C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:223 [inlined]
 [14] copy(M::ArrayLayouts.Mul{LazyArrays.ApplyLayout{typeof(\)}, ArrayLayouts.UnknownLayout, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{ClassicalOrthogonalPolynomials.LanczosConversion{Float64}, LazyArrays.BroadcastMatrix{Float64, typeof(\), Tuple{LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}}}}}}}, LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}}})
    @ LazyArrays C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\linalg\inv.jl:146
 [15] copy
    @ C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\linalg\inv.jl:148 [inlined]
 [16] materialize(M::ArrayLayouts.Mul{LazyArrays.ApplyLayout{typeof(\)}, LazyArrays.ApplyLayout{typeof(*)}, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{ClassicalOrthogonalPolynomials.LanczosConversion{Float64}, LazyArrays.BroadcastMatrix{Float64, typeof(\), Tuple{LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}}}}}}}, LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}}})
    @ ArrayLayouts C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:137
 [17] mul
    @ C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:138 [inlined]
 [18] *(A::LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{ClassicalOrthogonalPolynomials.LanczosConversion{Float64}, LazyArrays.BroadcastMatrix{Float64, typeof(\), Tuple{LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}, LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}}}}}}}, B::LazyArrays.ApplyArray{Float64, 2, typeof(*), Tuple{Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}})
    @ ArrayLayouts C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\mul.jl:223
 [19] _quad_matmul(A::LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{ClassicalOrthogonalPolynomials.LanczosConversion{Float64}, LazyArrays.BroadcastMatrix{Float64, typeof(\), Tuple{LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}, LazyArrays.ApplyArray{Float64, 2, typeof(inv), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}}}}}}}, B::BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}, C::Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, D::ClassicalOrthogonalPolynomials.LanczosConversion{Float64})
    @ LinearAlgebra C:\Users\jjyao\.julia\juliaup\julia-1.11.5+0.x64.w64.mingw32\share\julia\stdlib\v1.11\LinearAlgebra\src\matmul.jl:1212
 [20] *(A::LazyArrays.ApplyMatrix{Float64, typeof(\), Tuple{ClassicalOrthogonalPolynomials.LanczosConversion{Float64}, LazyArrays.BroadcastMatrix{Float64, typeof(\), Tuple{LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}, LazyArrays.ApplyArray{Float64, 2, typeof(inv), Tuple{BandedMatrix{Float64, LazyArrays.ApplyArray{Float64, 2, typeof(vcat), Tuple{Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(-), Tuple{LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, Zeros{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, Adjoint{Float64, LazyArrays.BroadcastVector{Float64, typeof(/), Tuple{Float64, InfiniteArrays.InfStepRange{Float64, Float64}}}}}}, InfiniteArrays.OneToInf{Int64}}}}}}}}, B::BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}, C::Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, D::ClassicalOrthogonalPolynomials.LanczosConversion{Float64})
    @ LinearAlgebra C:\Users\jjyao\.julia\juliaup\julia-1.11.5+0.x64.w64.mingw32\share\julia\stdlib\v1.11\LinearAlgebra\src\matmul.jl:1201
 [21] _copy_ldiv_mul
    @ C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\linalg\inv.jl:115 [inlined]
 [22] copy
    @ C:\Users\jjyao\.julia\packages\LazyArrays\ltmzk\src\linalg\inv.jl:121 [inlined]
 [23] copy
    @ C:\Users\jjyao\.julia\packages\ClassicalOrthogonalPolynomials\INrHP\src\normalized.jl:134 [inlined]
 [24] materialize(M::ArrayLayouts.Ldiv{ClassicalOrthogonalPolynomials.LanczosLayout, LazyArrays.ApplyLayout{typeof(*)}, LanczosPolynomial{Float64, QuasiArrays.ApplyQuasiVector{Float64, typeof(*), Tuple{Weighted{Float64, QuasiArrays.SubQuasiArray{Float64, 2, Legendre{Float64}, Tuple{ContinuumArrays.AffineMap{Float64, Inclusion{Float64, IntervalSets.ClosedInterval{Int64}}, Inclusion{Float64, ChebyshevInterval{Int64}}}, Base.Slice{InfiniteArrays.OneToInf{Int64}}}, false}}, LazyArrays.ApplyArray{Float64, 1, typeof(vcat), Tuple{Vector{Float64}, Zeros{Float64, 1, Tuple{InfiniteArrays.OneToInf{Int64}}}}}}}, Normalized{Float64, QuasiArrays.SubQuasiArray{Float64, 2, Legendre{Float64}, Tuple{ContinuumArrays.AffineMap{Float64, Inclusion{Float64, IntervalSets.ClosedInterval{Int64}}, Inclusion{Float64, ChebyshevInterval{Int64}}}, Base.Slice{InfiniteArrays.OneToInf{Int64}}}, false}, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}}, QuasiArrays.ApplyQuasiMatrix{Float64, typeof(*), Tuple{QuasiArrays.SubQuasiArray{Float64, 2, Ultraspherical{Float64, Float64}, Tuple{ContinuumArrays.AffineMap{Float64, Inclusion{Float64, IntervalSets.ClosedInterval{Int64}}, Inclusion{Float64, ChebyshevInterval{Int64}}}, Base.Slice{InfiniteArrays.OneToInf{Int64}}}, false}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}, Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}}})
    @ ArrayLayouts C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\ldiv.jl:22
 [25] ldiv
    @ C:\Users\jjyao\.julia\packages\ArrayLayouts\B2wRU\src\ldiv.jl:98 [inlined]
 [26] \(A::LanczosPolynomial{Float64, QuasiArrays.ApplyQuasiVector{Float64, typeof(*), Tuple{Weighted{Float64, QuasiArrays.SubQuasiArray{Float64, 2, Legendre{Float64}, Tuple{ContinuumArrays.AffineMap{Float64, Inclusion{Float64, IntervalSets.ClosedInterval{Int64}}, Inclusion{Float64, ChebyshevInterval{Int64}}}, Base.Slice{InfiniteArrays.OneToInf{Int64}}}, false}}, LazyArrays.ApplyArray{Float64, 1, typeof(vcat), Tuple{Vector{Float64}, Zeros{Float64, 1, Tuple{InfiniteArrays.OneToInf{Int64}}}}}}}, Normalized{Float64, QuasiArrays.SubQuasiArray{Float64, 2, Legendre{Float64}, Tuple{ContinuumArrays.AffineMap{Float64, Inclusion{Float64, IntervalSets.ClosedInterval{Int64}}, Inclusion{Float64, ChebyshevInterval{Int64}}}, Base.Slice{InfiniteArrays.OneToInf{Int64}}}, false}, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}}, B::QuasiArrays.ApplyQuasiMatrix{Float64, typeof(*), Tuple{QuasiArrays.SubQuasiArray{Float64, 2, Ultraspherical{Float64, Float64}, Tuple{ContinuumArrays.AffineMap{Float64, Inclusion{Float64, IntervalSets.ClosedInterval{Int64}}, Inclusion{Float64, ChebyshevInterval{Int64}}}, Base.Slice{InfiniteArrays.OneToInf{Int64}}}, false}, BandedMatrix{Float64, Fill{Float64, 2, Tuple{Base.OneTo{Int64}, InfiniteArrays.OneToInf{Int64}}}, InfiniteArrays.OneToInf{Int64}}, Diagonal{Float64, LazyArrays.Accumulate{Float64, 1, typeof(*), Vector{Float64}, AbstractVector{Float64}}}, ClassicalOrthogonalPolynomials.LanczosConversion{Float64}}})
    @ QuasiArrays C:\Users\jjyao\.julia\packages\QuasiArrays\hFvLm\src\matmul.jl:34
 [27] top-level scope
    @ REPL[11]:1

[DanielVandH]

It probably wouldn't've been an update from Julia 1.10 to 1.11, but rather an update to one of the packages. Do you have access to the package versions you were using when it worked, and the versions you have now with it broken?

[JiajieYao1106]

It probably wouldn't've been an update from Julia 1.10 to 1.11, but rather an update to one of the packages. Do you have access to the package versions you were using when it worked, and the versions you have now with it broken?

@DanielVandH Thank you for your comment. I didn’t have the earlier package versions saved, but everything is up to date now.

[dlfivefifty]

I'm about to merge a fix.

Ambiguity errors are super easy to fix so a PR like this should be doable without me. In particular, you can click on the links to see what the various overloads are and deduce what's missing.

Next time you encounter an error like this let me help you walk you through the fix

[JiajieYao1106]

I'm about to merge a fix.

Ambiguity errors are super easy to fix so a PR like this should be doable without me. In particular, you can click on the links to see what the various overloads are and deduce what's missing.

Next time you encounter an error like this let me help you walk you through the fix

@dlfivefifty Really appreciate your help, thank you!