回 帖 发 新 帖 刷新版面

主题:intel MKL pardiso能求解多大规模的矩阵?

不知道有人用过pardiso解大型稀疏矩阵吗?我估算好内存,构造矩阵a,但是在call pardiso的时候,内存使用突然增加,应该是在矩阵分解时浪费了更多的内存!我用OOC也不管用!希望使用过的指教我一下!我的a数组大小为4*10^7个复数,所占内存也就是16*4*10^7字节,还不到1G呢,就算不了了,帮我看看,这么大的矩阵能计算吗?用pardiso,谢谢啊!
我使用的参数为:      iparm(1) = 1 
      iparm(2) = 2 
      iparm(4) = 0 
      iparm(5) = 0 
      iparm(6) = 0 
      iparm(7) = 0 
      iparm(8) = 2 
      iparm(9) = 0 
      iparm(10) = 13
      iparm(11) = 1 
      iparm(12) = 0 
      iparm(13) = 1 
      iparm(14) = 0 
      iparm(15) = 0 
      iparm(16) = 0 
      iparm(17) = 0 
      iparm(18) = -1 
      iparm(19) = -1 
      iparm(20) = 0 
     下面是错误信息: 
      Error!***Memory allocation failed for SetUpCoarseGraph:gtata.Requested size:
      449357668 bytes

回复列表 (共2个回复)

沙发

我在CSDN上给你说过了,最好贴出你的运行报告,比如出错提示之类的,或者文字或者截屏,看看你有无遗漏掉什么内容。

板凳


谢谢啊!你在CSDN上叫什么名字啊?

我来回复

您尚未登录,请登录后再回复。点此登录或注册