c++的矩阵库getfem with python

不知道是封装还是什么的问题,昨天提到的pyulbas显得很慢.

我又找来一个c++的矩阵库,看起来不错

http://home.gna.org/getfem/

不过他的python封装不是很好,首先是只封装了一个运算用的稀疏矩阵,这种矩阵赋值的很慢的,不过好在还可以从文件中读取数据,而且比较快.

然后就是csc_val(self)
Return the array of values of all non-zero entries of M.
这个函数居然会报一个错,不过好在也可以保持到文件,囧

给出一点测试数据:

AMD Opteron(tm) Processor 2350 ,32G内存

边长100000的方阵,5000000的数据

然后跑了一次运算的速度如下

11.8575739861 loading cost

a+a 4.88795804977
a.T 11.9869470596
a*a 198.991373062

看看很其他矩阵库的速度对比如何,如果不错,就去修修他的封装...