【讨论】请教如何选取高性能的开放源代码数学库？

lxuejin

最近在琢磨一个问题，需要自己编写程序进行科学计算，比如进行几十万自由度的有限元计算
有没有比较好的开放源代码库可以使用呢？
通过搜索，找到了mtl和lapack++等等面向c++的库，种类很多，不知道该如何选取

fireghg

呵呵，这个没做过！
认真研究一下ansys，看看它们用什么做的～

refeihc

MTL是基于C++的开源代码，曾有人对它评价不错，于是我花了不少时间想拿下它。后来还是放弃了，主要原因有3个，一是它的类结构太繁琐，主要采用宏来生成类，而不是派生、抽象的方式，因此不方便调试跟踪。二是从2.0版以后，它就没有进展了，虽然号称4.0版正在开发中，但是开发人员彻底换了一批，我担心新的这批人消化不了老版本的代码。三是我曾用2.0试算了一些简单例子，但是有错。

Lapack系列是目前最有名的开放源代码，Lapack++它是在Lapack基础上发展起来的，其实就是Fortran内核加一个C++壳。这一系列代码可以免费得到，而且协议友好（没有GPL的传染性）。主要优点是速度快，运算稳定，得到全世界的普遍认同。不过它最大的缺点是没有稀疏求解器，一般到了十万级以上自由度的有限元求解就有些吃力了，建议还补充一个UMFPack。