LAPACK中DGESV与DGBSV求解器输出不一致问题排查

我已经被这个问题困扰好几个小时了，希望有人能帮忙找出问题所在。我先用DGESV求解Ax=B，对这个结果有99%的把握是正确的；之后把矩阵A转换成带状形式，用DGBSV求解却得到了不一样的结果。

输入输出情况 #

• 使用DGESV时：A为9×9矩阵，B为-r（1×9向量），求解得到的x为sk（1×9向量），对应输入输出截图（标注：Inputs and Ouputs using DGESV）
• 使用DGBSV时：A=JBtemp，B=-r，求解得到的x为SK，对应输入输出截图（标注：Inputs and Ouputs using DGBSV）

两次求解得到的x值完全不匹配，以下是我调用这两个函数的代码：

call dgesv(n,1,J0,n,ipiv,r,n,info)

call dgbsv(n,KL,KU,1,JBtemp,KU+1+KL+KL,ipiv,r,n,info)

内容的提问来源于stack exchange，提问作者William Dennis