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

Question

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

Answer 1

确实不可能设置 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] 到相同的值。