如何高效的将一维数组重复赋值给二维数组？

mmSCI · 发表于 2020-9-10 10:55:29

有一个一维数组U包含N个元素，和二维数组sU N*N个元素，需要把U赋值给sU的每一行，如何快速高效的赋值？

li913 · 发表于 2020-9-10 10:59:26

do concurrent (i=1:n)
su(:,i) = u
end do

风平老涡 · 发表于 2020-9-10 21:24:53

li913 发表于 2020-9-10 10:59
do concurrent (i=1:n)
su(:,i) = u
end do

如果一定强调“每一行”的话，最后可以加
su = transpose(su)

mmSCI · 发表于 2020-9-11 11:09:21

风平老涡发表于 2020-9-10 21:24
如果一定强调“每一行”的话，最后可以加
su = transpose(su)

如果不加会影响效率还是报错呢？

风平老涡 · 发表于 2020-9-11 21:05:31

mmSCI 发表于 2020-9-11 11:09
如果不加会影响效率还是报错呢？

不影响也不报错。只不过是每一列sU含N个U元素。

mmSCI · 发表于 2020-9-12 10:34:06

风平老涡发表于 2020-9-11 21:05
不影响也不报错。只不过是每一列sU含N个U元素。

懂了，这应该和初始赋值规则有关吧，如果是sU(:,i)=U就是sU是由U一列一列||构成的，如果是sU(i,:)=U就是sU是由U一行一行构成的=，我理解的没错吧

风平老涡 · 发表于 2020-9-12 11:29:24

mmSCI 发表于 2020-9-12 10:34
懂了，这应该和初始赋值规则有关吧，如果是sU(:,i)=U就是sU是由U一列一列||构成的，如果是sU(i,:)=U就是s ...

对的。不过如果N值很大的话，用sU(i,:)=U 不是很高效，因为Fortran数组在内存中是按列存放的。

[数值问题] 如何高效的将一维数组重复赋值给二维数组？