Fortran - 有效地复制矩阵

时间:2014-11-08 16:41:13

标签: performance matrix copy fortran

我的代码中有几个地方可以复制稀疏矩阵,大约有100,000,000个元素。我遇到的问题是,我在一个循环中处理矩阵,这个循环被调用大约1000次,并且处理矩阵需要大约2秒(完成复制花费大约一个小时花费大部分时间)。我使用等号(即A = B)复制矩阵。我使用以下内容初始化矩阵

data a/ myzp * 0.0/

这需要几乎时间来执行,其中' myzp'是矩阵中的元素,' a'是矩阵。

我的问题是在Fortran中复制矩阵比使用等号更有效吗?

更新: 这是有问题的代码需要2秒,

call cpu_time(t1)
a = b
call cpu_time(t2)

其中a和b都是矩阵。

0 个答案:

没有答案