Fortran Coder

标题: 函数传入和返回值均为包含动态数组的type [打印本页]

作者: Neo123    时间: 2022-6-1 16:59
标题: 函数传入和返回值均为包含动态数组的type
因为程序很大,简化之后的程序如下所示。原程序是采用ifortran采用参数化结构体编写,因为想对程序实现GPU并行,所以采用gfortran12.x编写。哪位大佬帮忙解答一下哪错了
[Fortran] 纯文本查看 复制代码
module ModParticle
  implicit none
  integer::x=4
  type ParticleType
    real(4) , allocatable :: MassGas(:)
  contains
    Procedure :: CopyParticle => CopyParticle

  end type
contains ! Routines of this module

  Subroutine CopyParticle( this , ID )
    integer(4), intent(in) ::  ID
    class(ParticleType) :: this
    this%MassGas = 0.1
  End Subroutine CopyParticle

  subroutine CopyType(TempParticlesx,TempParticles)
    type(ParticleType)::TempParticles
    type(ParticleType)::TempParticlesx
    allocate(TempParticles%MassGas(x))
      TempParticlesx%MassGas = TempParticles%MassGas
    deallocate(TempParticles%MassGas)
  end subroutine

end module

[Fortran] 纯文本查看 复制代码
PROGRAM Main
  use ModParticle
  implicit none

  integer::i
  type(ParticleType)::TempParticles
  type(ParticleType)::TempParticlesx
  allocate(TempParticles%MassGas(x))
  do i = 1, size(TempParticles%MassGas)
    call TempParticles%CopyParticle(i)
  end do
  write(*,*) TempParticles%MassGas

call CopyType(TempParticlesx,TempParticles)

write(*,*) TempParticlesx%MassGas
deallocate(TempParticles%MassGas)

END PROGRAM Main


作者: 风平老涡    时间: 2022-6-1 19:36
本帖最后由 风平老涡 于 2022-6-1 20:02 编辑

allocate(tempparticles%massgas(x)), 在 MAIN 和 SUBROUTINE COPYTYPE 中 重复 allocation。去掉 MAIN中的 allocate 语句。

作者: Neo123    时间: 2022-6-1 20:03
风平老涡 发表于 2022-6-1 19:36
allocate(tempparticles%massgas(x)), 在 MAIN 和 SUBROUTINE COPYTYPE 中 重复 allocation。

这不是问题所在,因为x在module ModParticle中已经声明了,main中引用了use ModParticle
作者: kyra    时间: 2022-6-2 10:31
Neo123 发表于 2022-6-1 20:03
这不是问题所在,因为x在module ModParticle中已经声明了,main中引用了use ModParticle ...

这就是问题所在。
重复分配了。一个可分配数组只能分配一次,要第二次分配必须先 Deallocate 释放。

与它在哪里声明无关。
作者: Neo123    时间: 2022-6-2 13:38
kyra 发表于 2022-6-2 10:31
这就是问题所在。
重复分配了。一个可分配数组只能分配一次,要第二次分配必须先 Deallocate 释放。

你看到的“风平老涡”的帖子是后来修改的,之前的内容是说我没定义x。
将ModParticle中的allocate删除后,可以在ifortran和gpfortran中正常编译,但是我在gfortran中还是报错undefined reference to ‘CopyType_’

打这些字的时候感觉是gfortran编译器的问题,刚才连之前用gfortran测试过的最简单的程序也都莫名奇妙报错,但吃个饭以后再用gfortran编译程序都可以正常编译了,不太明白是电脑问题还是软件问题
作者: zjk0112    时间: 2022-6-2 18:10
Neo123 发表于 2022-6-2 13:38
你看到的“风平老涡”的帖子是后来修改的,之前的内容是说我没定义x。
将ModParticle中的allocate删除后 ...

你要不把你测试报错的最简单程序发出来,我可以在我电脑上帮你测试一下是不是gfortran,或者gfortran版本的问题。
作者: Neo123    时间: 2022-6-2 22:13
zjk0112 发表于 2022-6-2 18:10
你要不把你测试报错的最简单程序发出来,我可以在我电脑上帮你测试一下是不是gfortran,或者gfortran版本 ...

我自己测试发现同一个程序codeblocks所带的旧版本的gfortran版本经常提示报错,有时会正常编译,但最新的gfortran 12.1是没有问题的
作者: zjk0112    时间: 2022-6-3 23:03
Neo123 发表于 2022-6-2 22:13
我自己测试发现同一个程序codeblocks所带的旧版本的gfortran版本经常提示报错,有时会正常编译,但最新的 ...

有时报错,有时正常?依我的经验,是不是初始化的问题。
作者: Neo123    时间: 2022-6-4 16:28
zjk0112 发表于 2022-6-3 23:03
有时报错,有时正常?依我的经验,是不是初始化的问题。

好的,谢谢你的建议,我根据你的建议再看看。确实是有时报错,有时正常。
作者: accding    时间: 2022-6-17 12:22
module里面,为何用两次contains?
作者: Neo123    时间: 2022-6-25 15:27
accding 发表于 2022-6-17 12:22
module里面,为何用两次contains?

其中一个是在派生type中的
作者: phost    时间: 2022-6-26 23:35
两次contain,一个是Fortran的03的新功能,派生类型的函数,类似c++,c#的class了。
不过感觉Fortran搞面向对象还是走入歧途,怎么搞也不可能搞过c#和java
还不如就加强数值计算功能得了
作者: 楚香饭    时间: 2022-6-27 08:07
phost 发表于 2022-6-26 23:35
两次contain,一个是Fortran的03的新功能,派生类型的函数,类似c++,c#的class了。
不过感觉Fortran搞面向 ...

话可不是这么说的。
能不能搞得过,是一回事。有没有,是另一回事。

毕竟,你要想借用Fortran的数值计算,又想要type过程捆绑。你不能简单方便的把C++/C的class特性拿过来用。是吧?
(不要提混编,混编的细粒度还不够)

我一直在用F2003的新特性,type捆绑过程,甚至抽象过程等等。且说和其他语言比较,还稍显稚嫩。但的确是对Fortran语法的重大补充。

恰如一台手机,他可以专注于摄像、专注于性能、专注于音乐。那么,它的屏幕可以差一些,但没有屏幕,那就是另外一回事了,对吧?




欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2