intel MKL pardiso求解大型稀疏矩阵,是不是很浪费内存?
我在使用intel fortran中的pardiso的时候,求解非对称矩阵,发现在call pardiso的时候,内存使用量特别大,是所用数组的几倍,这限制了矩阵规模!这个问题可以解决吗?
------解决方案--------------------稀疏矩阵在分解的时候,会产生更多的非零元素,需要增加内存的使用。
程序是不是在使用最新版本的MKL. MKL 10.2中对pardiso 内存使用,有所优化。
如果矩阵实在太大, 可以使用 paridso 的 out of core 的模式 (将中间计算结果存储于一个文件中),克服内存不足的问题。