BLAS/LAPACK例程用于进行高斯消元

Question

我是BLAS/Lapack的新用户，我想知道是否有一种例程可以执行高斯消元或高斯-约旦消元？我已经在谷歌上搜索并查看了它们的文档，但仍然找不到。

非常感谢你的帮助！

- Kelvin Lee

---

它必须是高斯/高斯-约旦算法吗？还是可以使用任何能够找到线性方程的精确解的算法？ - Matt Ball

@Matt Ball 我认为我可能更喜欢高斯/高斯-约旦算法，因为生成的矩阵将被放置在梯形形式中。其他算法是否也会产生行梯形形式？ - Kelvin Lee

@KelvinLee 只是好奇，行阶梯形式还有哪些其他操作会有帮助？（我只想学得更多，而不是暗示这不是一个有效的用例。） - Justin Meiners

Answer 1

高斯消元法基本上与LU分解相同。例程xGETRF计算LU分解（例如，DGETRF 用于实数双精度矩阵）。U因子对应于经过高斯消元后的矩阵。U因子存储在矩阵A的上三角部分（包括对角线）中，在退出时返回。

LU分解/高斯消元常用于解线性方程组。一旦您计算出LU分解，就可以使用xGETRS例程来解决线性系统。

- Jitse Niesen