Fortran中是否有高斯消除的标准解决方案？

Question

我要在Fortran做一个标准的高斯消除程序，并且想知道是否已经有一些准备好的解决方案了吗？优选地，具有功能的形状。我看了很多有关数字的书籍，但似乎没有任何好的解决方案。

Answer 1

您正在寻找LAPACK子例程DGETRF。这执行LU分解。 LU分解的“U”部分对应于高斯消元的解。

请注意，对于所有LAPACK / BLAS例程，第一个字母对应于数据类型...有关详细信息，请参阅wikipedia。