Julia v1.9 (#123)

* Julia v1.9

* Update Project.toml

* Update Project.toml

* Update Project.toml

* Update Project.toml

* BLAS.axpy! -> axpy!

* Revert "BLAS.axpy! -> axpy!"

This reverts commit 99fab34.

* fixes to clenshaw

Author: dlfivefifty

.github/workflows/ci.yml

@@ -12,6 +12,7 @@ jobs:
         version:
           - '1.7'
           - '1'
+          - '^1.9.0-0'
         os:
           - ubuntu-latest
           - macOS-latest

Project.toml

@@ -1,7 +1,7 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
 authors = ["Sheehan Olver <[email protected]>"]
-version = "0.7.1"
+version = "0.7.2"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -26,22 +26,22 @@ SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 
 [compat]
 ArrayLayouts = "0.8"
-BandedMatrices = "0.17"
+BandedMatrices = "0.17.17"
 BlockArrays = "0.16.9"
 BlockBandedMatrices = "0.11.6"
-ContinuumArrays = "0.12.3"
+ContinuumArrays = "0.12.4"
 DomainSets = "0.5.6, 0.6"
 FFTW = "1.1"
 FastGaussQuadrature = "0.4.3, 0.5"
 FastTransforms = "0.14.9"
 FillArrays = "0.13"
 HypergeometricFunctions = "0.3.4"
-InfiniteArrays = "0.12.3"
+InfiniteArrays = "0.12.11"
 InfiniteLinearAlgebra = "0.6.7"
 IntervalSets = "0.5, 0.6, 0.7"
-LazyArrays = "0.22.16"
+LazyArrays = "0.22.18"
 LazyBandedMatrices = "0.8.5"
-QuasiArrays = "0.9"
+QuasiArrays = "0.9.6"
 SpecialFunctions = "1.0, 2"
 julia = "1.7"
 

src/ClassicalOrthogonalPolynomials.jl

@@ -10,7 +10,7 @@ using ContinuumArrays, QuasiArrays, LazyArrays, FillArrays, BandedMatrices, Bloc
 import Base: @_inline_meta, axes, getindex, unsafe_getindex, convert, prod, *, /, \, +, -,
                 IndexStyle, IndexLinear, ==, OneTo, tail, similar, copyto!, copy, setindex,
                 first, last, Slice, size, length, axes, IdentityUnitRange, sum, _sum, cumsum,
-                to_indices, _maybetail, tail, getproperty, inv, show, isapprox, summary,
+                to_indices, tail, getproperty, inv, show, isapprox, summary,
                 findall, searchsortedfirst
 import Base.Broadcast: materialize, BroadcastStyle, broadcasted, Broadcasted
 import LazyArrays: MemoryLayout, Applied, ApplyStyle, flatten, _flatten, adjointlayout,

src/clenshaw.jl

@@ -130,7 +130,7 @@ end
 Base.@propagate_inbounds function _clenshaw_next!(n, A::AbstractVector, B::AbstractVector, C::AbstractVector, x::AbstractMatrix, c, bn1::AbstractMatrix{T}, bn2::AbstractMatrix{T}) where T
     # bn2 .= B[n] .* bn1 .- C[n+1] .* bn2
     lmul!(-C[n+1], bn2)
-    BLAS.axpy!(B[n], bn1, bn2)
+    LinearAlgebra.axpy!(B[n], bn1, bn2)
     muladd!(A[n], x, bn1, one(T), bn2)
     view(bn2,band(0)) .+= c[n]
     bn2
@@ -164,7 +164,7 @@ end
 
 Base.@propagate_inbounds function _clenshaw_first!(A, B, C, X, c, bn1, bn2) 
     lmul!(-C[2], bn2)
-    BLAS.axpy!(B[1], bn1, bn2)
+    LinearAlgebra.axpy!(B[1], bn1, bn2)
     muladd!(A[1], X, bn1, one(eltype(bn2)), bn2)
     view(bn2,band(0)) .+= c[1]
     bn2
@@ -316,7 +316,7 @@ function materialize!(M::MatMulVecAdd{<:ClenshawLayout,<:PaddedLayout,<:PaddedLa
     _fill_lmul!(β,y)
     resizedata!(y, last(jkr))
     v = view(paddeddata(y),jkr)
-    BLAS.axpy!(α, Ax, v)
+    LinearAlgebra.axpy!(α, Ax, v)
     y
 end