fix for Kronecker product scaling

JuliaSmoothOptimizers / LinearOperators.jl

Linear Operators for Julia

Other

150 stars 32 forks source link

Closed fguevaravas closed 1 year ago

fguevaravas commented 1 year ago

The Kronecker product code currently does not deal with multiplication by scalars correctly. This is illustrated below:

using LinearOperators, LinearAlgebra
A = kron(opEye(2),I(1))
Matrix(2A)

This results in:

2×2 Matrix{Float64}:
 1.0  0.0
 0.0  1.0

when it should be 2 times identity.

github-actions[bot] commented 1 year ago