将数组变矩阵

630767695 · 发表于 2022-9-11 22:27:49

strain（nblock*6）这个数组的大小是这个，我想讲其变为矩阵STR(nblock,6)，不知道这样的操作是否正确

[Fortran] 纯文本查看 复制代码

call vgetvrm('PE',strain,jData,cData,jStatus)？
          e11=strain(k)
          e22=strain(nblock+k)
          e33=strain(nblock*2+k)
          e12=strain(nblock*3+k)
          e23=strain(nblock*4+k)
          e13=strain(nblock*5+k)
          stateNew(k,2)=e11
          stateNew(k,3)=e22
          stateNew(k,4)=e33
          stateNew(k,5)=e12
          stateNew(k,6)=e23
          stateNew(k,7)=e13
          do r=1,6
            STR(k,r)=strain(k+(r-1)*nblock)
          end do

li913 · 发表于 2022-9-12 17:42:58

do r=1,6
STR(k,r)=strain(k+(r-1)*nblock)
end do
没问题，也可以更简单：
str = reshape(strain,[nblck,6])

630767695 · 发表于 2022-9-13 10:09:50

li913 发表于 2022-9-12 17:42
do r=1,6
STR(k,r)=strain(k+(r-1)*nblock)
end do

谢谢

[求助] 将数组变矩阵