一个关于无类型指针的问题

liudy02

lookbook 发表于 2019-7-3 22:39
谢谢绝地，
但这个代码在gfortran下测试并不成功，c_f_pointer不接受多态的指针，pgi是可以的，我猜ivf也 ...

难道不是c_f_pointer这个子进程属于iso_c_binding这个内部模块？

用父类确实可以实现一定程度的多态，但我认为还是很不理想的
多态进程中涉及到的成员变量在某个子类里有变化可能就会要重载一堆方法
而且多数时候事情没那么复杂，可能就是想省码点包含基本运算的代码
结果却要层层设置类树，设置一堆运算符之类的，麻烦要死
所以我很欣赏有位大佬在这里推荐的在fortran里写模块的办法
只是那个办法也是建立在C预处理的基础上的
如果不是和C编译器继承在一起的fortran编译器，事情就麻烦的多了……

randomzoom

liudy02 发表于 2019-7-3 13:40
不懂，我现在渐渐觉得fortran确实没啥不可替代性了
似乎一般的数学库也都有C语言版了
甚至连python里的nu ...

其他语言，你会不知道怎么死的，特别是业务逻辑多的。
fortran的优点就是太傻了，不按照他的语法就是不过。

fcode

fortran语法设计得比较严谨，不确切的、可能因疏忽导致莫名其妙错误的，都被严格禁止了。

lookbook

liudy02 发表于 2019-7-3 22:53
难道不是c_f_pointer这个子进程属于iso_c_binding这个内部模块？

用父类确实可以实现一定程度的多态，但 ...

我的本意是为什么要用c_f_pointer这个函数

fcode

gfortran的 c_loc 不支持多态的话，可能就麻烦点，得用不规范的 loc

  Subroutine CnvrtEndian( byteIn , nByte )
    use , intrinsic :: ISO_C_Binding
    class(*) , target , Intent(INOUT) :: byteIn
    integer  , Intent(IN)             :: nByte
    character , pointer :: p(:)
    type(c_ptr) :: cp
    cp = transfer( loc(byteIn) , cp )
    call c_f_pointer( cp , p , [nByte] )
    p = p(nByte:1:-1)
  End Subroutine CnvrtEndian

liudy02

lookbook 发表于 2019-7-4 11:05
我的本意是为什么要用c_f_pointer这个函数

我没用过c和 fortran的交互
印象中，似乎fortran 的位操作都是假操作，还限定整数类型
所以我虽然没看懂版主是怎么做的，但是应该是把fortran的指针和c的指针相互转换
然后利用c来不分类型地对位进行倒置

fcode

差不多是这样。fortran语法太严苛，保证安全的同时，也有失灵活。
大多数时候，没必要这么做。
一部分过程，如果传入类型不重要，只是单纯的做别的事情，把传入参数直接传走，可以考虑用class(*)

weixing1531

楼主的设想还是可以实现的
只是class(*)多态指针p的断开引用必须在select type内部才行     
  但Fortran2018标准已经放松该限制  
  摘自电子书《Modern Fortran Explained  Incorporating Fortran 2018》
  下载网址：http://www.doc88.com/p-0601727992809.html

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

program main
    implicit none
    integer,target::a=3
    real,target::b=8.15
    class(*),pointer::p
    
    ALLOCATE(p,source=a) !f2018标准支持 p=a
    call write_p(p)
    p=>b
    call write_p(p)
    ALLOCATE(p,source=b) !f2018标准支持 p=b
    call write_p(p)
    p=>a
    call write_p(p)
contains
    subroutine write_p(p)
        class(*),pointer,intent(in)::p

        select type(p)
        type is(integer)
            write(*,"(I3)")p
        type is(real)
            write(*,"(F8.2)")p
        end select
    end subroutine write_p
end program

lookbook

weixing1531 发表于 2019-7-6 21:41
楼主的设想还是可以实现的
只是class(*)多态指针p的断开引用必须在select type内部才行     
  但Fortran20 ...

我没明白你的意思，看了下代码，这还是要在select type里才能进行操作(包括绑定的过程的调用)啊。