Update src/SemiclassicalOrthogonalPolynomials.jl

Co-authored-by: Sheehan Olver <solver@mac.com>
JuliaApproximation · Jun 24, 2024 · 5279d29 · 5279d29
1 parent dc70b6b
commit 5279d29
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/SemiclassicalOrthogonalPolynomials.jl b/src/SemiclassicalOrthogonalPolynomials.jl
@@ -342,7 +342,8 @@ Returns conversion operator from SemiclassicalJacobi `P` to SemiclassicalJacobi
 """
 function semijacobi_ldiv(Q::SemiclassicalJacobi, P::SemiclassicalJacobi)
     @assert Q.t ≈ P.t
-    (Q.t ≈ P.t) && (Q.a ≈ P.a) && (Q.b ≈ P.b) && (Q.c ≈ P.c) && return SymTridiagonal(Ones(∞),Zeros(∞))
+    T = promote_type(eltype(Q), eltype(P))
+    (Q.t ≈ P.t) && (Q.a ≈ P.a) && (Q.b ≈ P.b) && (Q.c ≈ P.c) && return SquareEye{T}(∞)
     Δa = Q.a-P.a
     Δb = Q.b-P.b
     Δc = Q.c-P.c