Fortran 90中如何无拷贝缩小大型数组尺寸？

好问题！处理大数组时，避免不必要的拷贝确实是提升效率的关键。在Fortran里，你完全可以不用复制数据，就能“重塑”出只包含有效元素的数组视图，核心是利用数组指针或者假定形状数组来直接指向原数组的有效区域。

方案1：使用数组指针（最直接的方式） #

数组指针可以直接关联到原数组的一部分内存，不需要任何数据复制，操作指针就等同于操作原数组的对应区域，开销几乎为0。具体代码如下：

! 你的原数组声明（注意：Fortran默认是1-based索引，如果要使用0-based，需要编译器支持，比如gfortran加-fzero选项）
real, dimension(100000, 5) :: a
integer :: n = 30000  ! 填充完成后确定的有效行数上限

! 声明一个可分配的数组指针
real, dimension(:,:), pointer :: a_trimmed

! 将指针关联到原数组的有效部分（这里对应你说的a(0:n, :)）
a_trimmed => a(0:n, :)

这样a_trimmed的形状就是(n+1, 5)（因为从0到30000共30001行），如果你的实际需求是(30000,5)，只需要调整切片范围为a(1:30000, :)即可。之后你所有对a_trimmed的读写操作，都会直接作用在原数组的对应内存上，完全没有拷贝开销。

方案2：在子程序中使用假定形状数组 #

如果你的后续操作是在子程序中进行的，可以直接把原数组的有效切片传递给子程序，子程序用假定形状数组接收，同样不会产生拷贝：

! 定义处理数据的子程序，用假定形状数组作为参数
subroutine process_valid_data(arr)
    real, dimension(:,:), intent(inout) :: arr
    ! 在这里arr的形状就是你传递进来的切片的形状，比如传入a(0:30000,:)，arr就是(30001,5)
    ! 直接操作arr即可，所有修改都会反映到原数组上
end subroutine process_valid_data

! 主程序中调用子程序
call process_valid_data(a(0:30000, :))

Fortran的假定形状数组是通过数组描述符传递的，本质上也是直接指向原数组的内存区域，所以同样没有数据复制的开销。

注意事项 #

• 如果你的代码依赖0-based索引，要确保编译器支持（比如GCC的-fzero选项，Intel Fortran的-zero选项），否则建议切换为Fortran默认的1-based索引，避免潜在的兼容性问题。
• 如果你后续不再需要原数组的剩余空间，静态数组（比如你声明的dimension(100000,5)）无法动态释放内存。如果想更灵活地管理内存，建议一开始就使用可分配数组（real, dimension(:,:), allocatable :: a），填充完成后可以重新分配并复制数据，但这种方式还是会有一次拷贝，所以如果追求极致效率，指针方案依然是最优选择。

内容的提问来源于stack exchange，提问作者fronthem