wiki 上的 Fortran TDMA 算法

Fortran TDMA algorithm on wiki

我最近访问了 wiki 站点:

https://en.wikibooks.org/wiki/Algorithm_Implementation/Linear_Algebra/Tridiagonal_matrix_algorithm#Fortran_90

表示a,b,c分别是次对角线、对角线和上对角线。如果'n'是b的大小,那么'a'和'c'的大小不就是n-1吗?通过检查可以清楚地看到 c(1) 和 c(n) 都被访问了,这与我描述的大小相矛盾。

这是怎么回事?这个算法错了吗?我假设的矩阵大小是错误的吗?该算法的输入矩阵是什么样的?

非常感谢任何帮助!

cp(n) 的计算只是为了避免尴尬的 if 语句——它从不用于确定 x。所以c(n)并没有真正用到

正如预期的那样,

a(1) 也未使用。