consistency with ClassicalOPs 0.12.4

Project.toml

@@ -1,7 +1,7 @@
 name = "SemiclassicalOrthogonalPolynomials"
 uuid = "291c01f3-23f6-4eb6-aeb0-063a639b53f2"
 authors = ["Sheehan Olver <solver@mac.com>"]
-version = "0.5.4"
+version = "0.5.5"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -20,7 +20,7 @@ SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 [compat]
 ArrayLayouts = "1"
 BandedMatrices = "0.17, 1"
-ClassicalOrthogonalPolynomials = "0.11, 0.12"
+ClassicalOrthogonalPolynomials = "0.12.4"
 ContinuumArrays = "0.15, 0.16, 0.17"
 FillArrays = "1"
 HypergeometricFunctions = "0.3.4"

src/SemiclassicalOrthogonalPolynomials.jl

@@ -142,9 +142,9 @@ function semiclassical_jacobimatrix(t, a, b, c)
     P = Normalized(jacobi(b, a, UnitInterval{T}()))
     iszero(c) && return jacobimatrix(P)
     if isone(c)
-        return cholesky_jacobimatrix(x->(t-x),P)
+        return cholesky_jacobimatrix(Symmetric(P \ ((t.-axes(P,1)).*P)), P)
     elseif isone(c/2)
-        return qr_jacobimatrix(x->(t-x),P)
+        return qr_jacobimatrix(Symmetric(P \ ((t.-axes(P,1)).*P)), P)
     elseif isinteger(c) && c ≥ 0 # reduce other integer c cases to hierarchy
         return SemiclassicalJacobi.(t, a, b, 0:Int(c))[end].X
     else # if c is not an integer, use Lanczos
@@ -153,7 +153,7 @@ function semiclassical_jacobimatrix(t, a, b, c)
     end
 end
 
-function semiclassical_jacobimatrix(Q::SemiclassicalJacobi, a, b, c)
+|function semiclassical_jacobimatrix(Q::SemiclassicalJacobi, a, b, c)
     Δa = a-Q.a
     Δb = b-Q.b
     Δc = c-Q.c