如何在 julia 的 SymTridiagonal 矩阵中分配非对角线条目?

How to assign off-diagonal entries in SymTridiagonal matrix in julia?

Julia 中的 SymTridiagonal 数据类型不允许我将非对角线值分配给除零以外的任何值。我收到此错误:ArgumentError:无法设置非对角线条目 (2, 1)。 我需要分配非对角线值,因为我正在尝试实现 ImplicitSymmetricQRStep 算法,该算法需要在此过程中执行此操作。

确实不可能设置 SymTridiagonal 矩阵的非对角线值 - 我不能说为什么要做出这个决定。

我现在看到两个选择:

1) 在 Julia 中,结构的字段不会被隐藏,因此可以通过这种方式更改值。但这很危险,因为该矩阵的内部结构可能会在未来版本中发生变化而没有任何警告。这是您如何执行此操作的示例:

using LinearAlgebra: SymTridiagonal

a = SymTridiagonal([1 2 0; 2 1 2; 0 2 1)] # 1 on diagonal, 2 on off diagonals
a.ev[1] = 4 # a[1, 2] == 4 and a[2, 1] == 4

2) 你也可以使用 Tridiagonal 矩阵类型,它也在 LinearAlgebra 包中;这种类型允许设置关闭对角线条目。然后你只需要确保自己没有违反该矩阵的对称属性,即如果你设置 a[i, j] 那么你还必须将 a[j, i] 设置为相同的值。