从python调用LAPACK驱动程序/函数

我需要在python中使用lapack函数“DGGLSE”来解决带约束的最小二乘问题

我发现lapack的一些低级函数可以在scipy中找到，如下所示 但我感兴趣的那个并没有包括在内，而且似乎lapack的大多数驱动程序代码都没有包括在内

我想知道是否有办法做到这一点。谢谢

---

-Yan

我只有一点编写Fortran的经验，但据我所知，编写SciPy中尚未包含的LAPACK函数的标准方法是使用NumPy附带的工具f2py。有关基本说明，请访问

你也可以用Cython包装它。 通过这种方式，您可以编译Fortran文件并对其进行链接，如中所示 这里描述了如何使用Cython或ctypes包装Fortran子例程。 该特定示例可在以下网址获得：

---

请确保注意数组在内存中的排列方式（C连续与Fortran连续）。

？SciPy 1.0版中包含GGLSE

例程系列。您可以先调用

？gglse_lwork

以获得硬件上的最佳块大小，然后将其用作

？gglse

的参数

您可以通过以下方式在ipython控制台上检查签名：

scipy.linalg.lapack.dgglse?

或者使用“帮助”等等。

看一看，其中解释了如何使用Cython包装DGEMM，您可能可以使用类似的方法来使用DDGLSE。。。