我已经使用 repa 实现了欧拉辛方法。这是一个非常简单的方法。可悲的是，我得到的表现并没有我希望的那么好。repa 的人非常乐于助人（尤其是 Ben Lippmeier），我遵循了他关于内联和添加严格性注释的建议。不幸的是，这些似乎没有帮助。我正在生产核心
ghc --make -O2 RepaOnly.hs -main-is RepaOnly -ddump-prep -dsuppress-all -dppr-case-as-let -dppr-cols200 > RepaOnly.hcr
可悲的是，这似乎很大，我目前不确定如何进行。我认为更广泛的 Haskell 社区可能对如何调试或指出我哪里出错有一些想法。
代码在这里：https ://gist.github.com/idontgetoutmuch/6209752
如果有帮助，我也可以将核心发布在要点中。