我的代码中有几个地方可以复制稀疏矩阵,大约有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都是矩阵。