1

我有一个关于 Fortran 和正确分配可分配用户派生类型的问题。

这是我的代码：

module polynom_mod
 implicit none

 type monomial
  integer,dimension(2) :: exponent
 end type

type polynom
  real, allocatable, dimension(:) :: coeff
  type(monomial),allocatable, dimension(:)   :: monom
  logical :: allocated
 !recursive type
  type(polynom),pointer :: p_dx,p_dy
 contains
  procedure :: init
  procedure :: init_dx
end type

在这里，我想导出一个类型多项式，我可以在其中执行以下操作：

p%coeff(1)=1.0 
p%monom(1)%exponent(1)=2

和类似的东西：

p%p_dx%coeff(1)=1.0 
p%p_dx%monom(1)%exponent(1)=2

所以我写了一些初始化类型绑定程序，我可以在其中初始化和分配我的类型：

contains

function init(this,num) result(stat)
  implicit none
  integer, intent(in)      :: num
  class(polynom),intent(inout) :: this
  logical :: stat

  allocate(this%coeff(num))
  allocate(this%monom(num))

  this%allocated = .TRUE.
  stat = .TRUE.
end function

function init_dx(this,num) result(stat)
  implicit none

  integer, intent(in)      :: num
  class(polynom),intent(inout) :: this

  logical :: stat

  allocate(this%p_dx%coeff(num))
  allocate(this%p_dx%monom(num))

  this%p_dx%allocated = .TRUE.
  stat = .TRUE.
 end function   
end module

program testpolytype
 use polynom_mod

 type(polynom) :: p

 if(p%init(2)) then
  print *,"Polynom allocated!"
 end if

 if(p%p_dx%init_dx(2)) then
  print *,"Polynom_dx allocated!"
 end if

结束程序

这将与 gfortran 4.6.3 一起编译，但是当我运行它时，我遇到了分段错误！

有没有办法分配递归可分配类型？

4

2 回答 2

4

您的代码的表面问题是，当p%p_dx%init_dx(2)计算表达式时，指针组件p%p_dx未定义，并且引发了分段错误。请注意，指针是未定义的，而不仅仅是未关联的。

现在我正在努力想出一个快速解决方案。长期的解决办法是解决我认为您的方法中的严重缺陷；请注意，这是我的意见，而不是非黑即白的问题，因此请仅在您关心我的意见时继续阅读。

函数init和init_dx并非没有副作用，实际上可以说它们几乎都是副作用——它们返回一个逻辑值，并且作为副作用初始化一个polynom变量。该程序似乎无法在polynom不评估的情况下初始化 a ，init并且init如果不将其包装到语句中就无法评估，例如

if (p%init(2)) then
end if

我想您可以将这些初始化函数重写为子例程，可能带有诸如

call initialise_polynom(p,2)

这至少会从您的代码中去除不纯函数的污点。但更好的方法是编写一个函数，例如：

function new_poly(num)
  implicit none
  integer, intent(in) :: num
  type(polynom) :: new_poly
  allocate(new_poly%coeff(num))
  allocate(new_poly%monom(num))
  allocate(new_poly%p_dx)
end function new_poly

哪个

a) 返回一个新的polynom；和

b) 分配组件p_dx；和

c) 无副作用。

然后，您可以使用表达式创建一个新polynom的，例如

p = new_poly(3)

并使用表达式初始化组件，例如

p%p_dx = new_poly(3)
于 2013-05-07T13:24:29.620 回答
1

回答我自己的问题，我想出了另一个解决方案，女巫也可以在没有指针的情况下工作，但它不像马克的那样优雅。

定义另一种类型：

type p_dx
 real, allocatable, dimension(:) :: coeff
 type(monomial),allocatable, dimension(:)   :: monom
 logical :: allocated
end type

然后将其用于：

type polynom
 real, allocatable, dimension(:) :: coeff
 type(monomial),allocatable, dimension(:)   :: monom
 type(p_dx) :: dx
 logical :: allocated
contains
 procedure     :: init
end type

所以你可以做类似的事情：

type(polynom) :: p

p%init(2)
p%dx%init_dx(3)
于 2013-05-08T10:54:26.163 回答