Fix ScaledOperator concretize!

SciML · Jul 21, 2023 · 36d2ef2 · 36d2ef2
1 parent 958e00a
commit 36d2ef2
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/basic.jl b/src/basic.jl
@@ -260,7 +260,7 @@ has_ldiv(L::ScaledOperator) = has_ldiv(L.L) & !iszero(L.λ)
 has_ldiv!(L::ScaledOperator) = has_ldiv!(L.L) & !iszero(L.λ)
 
 function concretize!(A, L::ScaledOperator{T}, α, β) where {T}
-    concretize!(A, L.L, L.λ * α, β)
+    concretize!(A, L.L, convert(Number, L.λ) * α, β)
 end
 
 function cache_internals(L::ScaledOperator, u::AbstractVecOrMat)

diff --git a/test/basic.jl b/test/basic.jl
@@ -329,7 +329,7 @@ end
 
 @testset "ConcretizedOperator" begin
     A = rand(2, 2); B = rand(2, 2);
-    L = MatrixOperator(A; update_func=(u,p,t)->t * A) + MatrixOperator(B; update_func=(u,p,t)->t * B)
+    L = MatrixOperator(A; update_func=(u,p,t)->t * A) + 2 * MatrixOperator(B; update_func=(u,p,t)->t * B)
     C = ConcretizedOperator(L)
     v = rand(2)
     C * v ≈ L * v