We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
diag
Diagonal
kron
1 parent e62193f commit 0fa47f6Copy full SHA for 0fa47f6
src/diagonal.jl
@@ -743,16 +743,16 @@ end
743
kron(A::Diagonal, B::Diagonal) = Diagonal(kron(A.diag, B.diag))
744
745
function kron(A::Diagonal, B::SymTridiagonal)
746
- kdv = kron(diag(A), B.dv)
+ kdv = kron(A.diag, B.dv)
747
# We don't need to drop the last element
748
- kev = kron(diag(A), _pushzero(_evview(B)))
+ kev = kron(A.diag, _pushzero(_evview(B)))
749
SymTridiagonal(kdv, kev)
750
end
751
function kron(A::Diagonal, B::Tridiagonal)
752
# `_droplast!` is only guaranteed to work with `Vector`
753
- kd = convert(Vector, kron(diag(A), B.d))
754
- kdl = _droplast!(convert(Vector, kron(diag(A), _pushzero(B.dl))))
755
- kdu = _droplast!(convert(Vector, kron(diag(A), _pushzero(B.du))))
+ kd = convert(Vector, kron(A.diag, B.d))
+ kdl = _droplast!(convert(Vector, kron(A.diag, _pushzero(B.dl))))
+ kdu = _droplast!(convert(Vector, kron(A.diag, _pushzero(B.du))))
756
Tridiagonal(kdl, kd, kdu)
757
758
0 commit comments