标签： algorithm fortran sparse-matrix fortran90 matrix-inverse
我最近访问了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）都被访问，这与我所描述的尺寸相矛盾。
这里发生了什么？这个算法错了吗？我假设的矩阵大小是错误的吗？这个算法的输入矩阵是什么样的？
非常感谢任何帮助！
答案 0 :(得分：1)
cp(n)的计算方法是避免使用尴尬的if语句 - 它从未用于确定x。因此c(n)并未真正使用。
cp(n)
x
c(n)
a(1)也未按预期使用。
a(1)