翻转矩阵fortran

Question

我想将矩阵上下颠倒。 使得T（1,1）= C（2,1）

我已经编写了这个程序，我发现了一个应该在线执行的代码，该代码为C=T(2:1:-1, :)但是当尝试获取值C（1,1）时应该为3我得到1.3533635457363350E-306。 您如何翻转矩阵以使上升变为下降？

program main


implicit none
  integer iMax, jMax
  double precision, dimension(:,:), allocatable :: T,C

double precision x, dx,f,L2old,L2norm

integer i, j,n


 allocate(T(0:2, 0:2))
 allocate(C(0:2, 0:2))


T(1,1)=1
T(1,2)=2
T(2,1)=3
T(2,2)=4

write(*,*) T(2,2)

C=T(2:1:-1, :)

Write(*,*) C(1,2)


end program main

Answer 1

如果您分配合适大小的矩阵，则所有内容均应按预期工作。

例如这个程序

program main
  implicit none

  double precision, dimension(:, :), allocatable :: t, c
  integer :: i

  allocate (t(1:2, 1:2))
  allocate (c(1:2, 1:2))

  t = reshape([1, 3, 2, 4], shape(t))
  do i = 1, 2
    write (*, *) t(i, :)
  end do
  write (*, *) ""

  c = t(2:1:-1, :)
  do i = 1, 2
    write (*, *) c(i, :)
  end do
end program main

产生以下输出

   1.0000000000000000        2.0000000000000000
   3.0000000000000000        4.0000000000000000

   3.0000000000000000        4.0000000000000000
   1.0000000000000000        2.0000000000000000

---

或者，如果您确实想使用3x3矩阵，则错误在于C=T(2:1:-1, :) 。 它应该是C=T(2:0:-1, :) 。

program main
  implicit none

  double precision, dimension(:, :), allocatable :: t, c
  integer :: i

  allocate (t(0:2, 0:2))
  allocate (c(0:2, 0:2))

  t = reshape([1, 4, 7, 2, 5, 8, 3, 6, 9], shape(t))
  do i = 0, 2
    write (*, *) t(i, :)
  end do
  write (*, *) ""

  c = t(2:0:-1, :)
  do i = 0, 2
    write (*, *) c(i, :)
  end do
end program main

输出：

   1.0000000000000000        2.0000000000000000        3.0000000000000000
   4.0000000000000000        5.0000000000000000        6.0000000000000000
   7.0000000000000000        8.0000000000000000        9.0000000000000000

   7.0000000000000000        8.0000000000000000        9.0000000000000000
   4.0000000000000000        5.0000000000000000        6.0000000000000000
   1.0000000000000000        2.0000000000000000        3.0000000000000000

---

计数数组的元素时要小心。 一对一错误可能很难调试，因此最好始终从0或始终从1开始计数。为了安全起见，请始终在lbound和ubound内部函数的帮助下遍历数组，而不是使用显式边界作为上面已经完成了：

  do i = lbound(t, dim=1), ubound(t, dim=1)
    write (*, *) t(i, :)
  end do

Answer 2

您正在分配3 x 3数组（索引范围为0-2，含2和3）。 您只将值分配给数组T的九个元素中的四个，然后将该数组的2 x 3切片分配给C 这是无效的，因为赋值运算符右侧的数组表达式的形状与左侧的不同。

此外，您选择要写入的数组元素掩盖了对所涉及数组尺寸的持续困惑。 如果将T的尺寸分配为0：2、0：2，然后垂直翻转，则与T(2,2)对应的结果元素将为C(0,2) 。

有几种可能的方法来更正代码。 其中更可能的是：

• 如果要继续使用现在使用的尺寸，则翻转数组的表达式为T(2:0:-1, :) 。

• 如果要使用索引范围为1-2（含2和2）的2 x 2数组C=T(2:1:-1, :)与C=T(2:1:-1, :) ，则适当地分配数组（例如， allocate(T(2, 2)) ）。