D 的快速线性系统求解器？ [关闭]

Question

我在哪里可以获得用 D 编写的快速线性系统求解器？它应该能够取一个方阵 A 和一个向量 b 并求解方程 Ax = b 为b，理想情况下，还对 A 执行显式反转。我有一个我自己写的，但它很慢，可能是因为它完全是缓存幼稚的。但是，对于我的用例，我需要具有以下 绝对、不可协商 要求的东西，即如果它不满足这些要求，那么我不在乎否则有多好：

1. 必须获得公共领域许可、Boost 许可或一些类似的许可许可。理想情况下，它不应该要求在二进制文件中进行归属（即不是 BSD），尽管这一点是可以协商的。

2. 必须用纯 D 编写或易于翻译成纯 D。难以理解的 Fortran 代码（即 LAPACK）无论多快都不是一个好的答案。

3. 必须针对大型（即 n > 1000）系统进行优化。我不希望为游戏程序员设计的东西真的非常快地解决 4x4 矩阵。

4. 不得与我不需要的庞大资料库密不可分。

编辑：这些看似疯狂的要求的原因是，我需要此代码用于获得许可的开源库，我不想拥有任何第三方依赖项。

Answer 1

如果您不喜欢 Fortran 代码，Eigen 是一个相当快速的 C++ 密集矩阵库，它具有适度的多核支持、编写良好的代码和良好的用户界面。将其代码转换为 D（或从中获取一些算法）应该很简单。

现在我的“考虑您的要求”：“每个人”（Mathematica、Matlab、Maple、SciPy、GSL、R、...）使用 ATLAS / LAPACK、UMFPACK、PARDISO、CHOLMOD 等是有原因的。编写快速、多线程、内存效率高、可移植且数值稳定的矩阵求解器是一项艰巨的工作（相信我，我已经尝试过）。很多这些艰苦的工作已经进入了 ATLAS 和其他领域。

所以我的方法是根据您的矩阵类型为相关库编写绑定，并从 D 链接到 C 接口。也许multiarray 中的绑定就足够了（我没试过）。否则，我建议查看另一个 C++ 库，即 uBlas 和各自的 bindings 以获得想法。