Fortran 03/08（gfortran 编译器）中使用无限多态类型进行数组操作

我想通过以下方式实现有用的数组操作（添加元素、删除元素、通过可分配/指针/二叉树结构实现不同的实现）class(*)特征（无限多态性）。我使用 gfortran 5.0 应该可以处理这样的功能。我需要它，以免为我使用的每种类型重复相同的代码。

这应该看起来像

function add_element(array,element)
  class(*),intent(in)::array(:)
  class(*),intent(in)::element
  class(*)::add_element(size(array)+1)
       add_element=[array,element]
end function

问题是，当我尝试将此函数与某些确定类型一起使用时，返回结果时出现错误。我无法分配class(*)到一些明确的类型变量没有select type，我当然不希望每次使用它时都选择类型结构。在子例程中，我不应该知道我想要使用的任何类型，因为我将创建其中的许多类型。

我尝试了一些变体move_alloc，来源，尝试使用子例程intent(out)争论等等。它不起作用。我认为它应该在参数属性中定义，与大小相同（带有源关键字？），但在标准中没有找到此类结构的示例或定义。当然，我会更多地研究这个标准（我不是专业程序员，而是物理学家，试图使我的程序可测试、可检查并且更容易更改），并且现在将简单地重复此代码以等待更好的解决方案，但也许有人知道在哪里在标准或某本书中搜索它？我认为这不仅与数组有关，还与使用class(*)无论如何，我认为应该有一些不知道类型的方法......

不知道我是否应该添加该子例程的其他无效形式的示例或它对错误的说明 - 否则问题将失去焦点。它可以被编译，但在所有情况下，在调用中分配给明确的类型都不起作用。为了论证intent(out) or (inout)它不能从虚拟论证转变为实际论证。从源重新分配会生成一个具有类型的对象（以及在我的示例中分配的结果），但该类型是隐藏的......并且我无法在子例程中使用选择类型，因为我不知道该类型。

另外，我不知道可以检查“相同类型”或在这种情况下的某些内容的构造......

这不是一个简单的问题您可以使用select type，但是 Fortran 没有类似的东西type is(type_of(x))。另一方面，还有SAME_TYPE_AS() and EXTENDS TYPE_OF()内在函数，但不能将它们用作类型保护。

有必要确保两者的动态类型array and element是相同的。

我认为这是标准的缺陷。

但是，您的方法仍然存在错误。您应该使函数结果可分配，以便能够将其分配给正确的动态类型：

class(*), allocatable ::add_element(:)

您可能会想到以下内容：（未经测试！使用 gfortran-4.9 ifort14 进行编译）

allocate(add_element(size(array)+1), mold=array)

但如何实际转移这些值我不知道，我担心如果不诉诸一些肮脏的伎俩，这可能是不可能的。

你甚至不能使用transfer这就是我看到的真正的缺陷。即使你可以用多态模具调用转移

transfer(element, add_element(1))

你无法将它分配给数组元素

add_element(1) = transfer(element, add_element(1))

我的观点是 Fortran 缺乏类型保护选项，只能确保两个变量具有相同的动态类型。

您可能会想到以下内容：（未经测试！使用 gfortran-4.9 ifort14 进行编译）

function add_element(array,element)
  use iso_c_binding
  implicit none
  class(*),intent(in)::array(:)
  class(*),intent(in)::element
  class(*), allocatable ::add_element(:)
  type(c_ptr) :: tmp

  interface
    function memcpy(dest, src, n) bind(c)
      use iso_c_binding
      integer(c_intptr_t),value :: dest, src
      integer(c_size_t) :: n
      type(c_ptr) :: memcpy
    end function
  end interface

  allocate(add_element(size(array)+1), mold=array)

  tmp = memcpy(loc(add_element(size(array)+1)), &
               loc(array), &
               size(array, kind=c_size_t) * storage_size(array, c_size_t)/8_c_size_t )
  tmp = memcpy(loc(add_element(size(array)+1)), &
               loc(array(1)), &
               storage_size(element, c_size_t)/8_c_size_t )

end function