[size=large]不知道是封装还是什么的问题,昨天提到的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
看看很其他矩阵库的速度对比如何,如果不错,就去修修他的封装...[/size]
相关资源：使用python/opengl进行二维重力坝显示
