Fortran 90 数组按引用传递问题

发布于 2025-01-05 07:41:31 字数 334 浏览 5 评论 0原文

如果我这样分配 F90 指针会发生什么：

real, pointer :: abc(:)

allocate abc (nx*ny*nz)

我将 abc 传递给子例程，在子例程中我将其重新定义为

real arg1(nx,ny,xz)

这似乎工作正常。

但是如果我重新定义为二维数组，就会出现段错误。

real arg1(nx,ny)

使用上面重新排序的数组，它应该可以工作。为什么会失败？任何帮助将不胜感激。

谢谢。

原文

What happens if I allocate a F90 pointer thus:

real, pointer :: abc(:)

allocate abc (nx*ny*nz)

I pass abc to a subroutine, where i redefine it as

real arg1(nx,ny,xz)

This seems to work fine.

But if I redefine as 2D array, I get a segfault.

real arg1(nx,ny)

With reordered arrays as above, it should work. Why does it fail? Any help will be appreciated.

Thanks.

分享到QQ

分享到微博

如果你对这篇内容有疑问，欢迎到本站社区发帖提问参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。

发布评论

需要登录才能够评论，你可以免费注册一个本站的账号。

撩发小公举 2025-01-12 07:41:31

它失败是因为你在数组的大小上向编译器撒谎。

在 Fortran 77 中，人们总是使用这些技巧，因为别无选择。在我们当前开明的时代，我们永远不应该做这些伎俩——只有 Fortran77 和 C 程序员才不得不诉诸这种诡计。

根据更改 fortran 中的数组尺寸的答案，如果您想调整数组大小，或者只是使用不同的形状创建对它的引用，您可以使用 RESHAPE 内在函数，或使用数组指针。

回复收藏 0 原文

木有鱼丸 2025-01-12 07:41:31

我不太明白你的问题，你所做的应该有效，因为在这种情况下只传递了对数组的引用。

program test
    integer :: nx, ny
    real, pointer :: abc(:)
    nx = 2
    ny = 3
    allocate(abc(nx**3))
    call my_sub(abc, nx, ny)
    write(*,*) abc
end program

subroutine my_sub(abc, nx, ny)
    integer :: nx, ny
    real :: abc(nx,ny)
    abc = 1.0
end subroutine

要了解数组是如何传递的，您可以阅读此页面。在底部您可以找到包含可能情况的表格。这里适用左下角的情况。

I don't really understand your problem, what you do should work, since in this case only a reference to the array is passed.

program test
    integer :: nx, ny
    real, pointer :: abc(:)
    nx = 2
    ny = 3
    allocate(abc(nx**3))
    call my_sub(abc, nx, ny)
    write(*,*) abc
end program

subroutine my_sub(abc, nx, ny)
    integer :: nx, ny
    real :: abc(nx,ny)
    abc = 1.0
end subroutine

To know how the arrays are passed, you can read this page. At the bottom you can find a table with possible situations. Here the bottom left case applies.

回复收藏 0 原文

~没有更多了~