Fortran Coder

楼主: lookbook
打印 上一主题 下一主题

[派生类型] 参数化类型如何绑定

[复制链接]

178

帖子

15

主题

0

精华

大宗师

F 币
4973 元
贡献
1152 点
15#
发表于 2018-11-16 22:42:50 | 只看该作者
本帖最后由 liudy02 于 2018-11-16 22:44 编辑

顺便说句,我终于现实解决了之前说的想定义不同kind但是功能几乎一样的类型,
甚至包括不同的数据类型,例如有Real又有Complex的办法了,不过在LINUX下才好操作这个
此法既可以减小文件和module的大小,又可以保证所有定义各种派生类型的module的一致性
简单来说就是在一开始只写一个Module文件,在其开始时定义kind常数
Module当然采用Private语句,Type的数据成员和绑定进程基本也全Private,
只用Interface和operator暴露几个想要被引用的进程
然后写个sh脚本,自动copy出想要的几个文件名,
然后sh脚本里继续用sed命令逐个文件自动替换掉kind常数定义和定义的type类型名
写个含有program的主文件,啥都不用干,就是USE所有的module
最后写个makefile文件,将sh脚本的运行也加上去,让它能够用主文件和所有的module文件编译出个程序就OK
这样也能看出来哪里有错了,哈哈

178

帖子

15

主题

0

精华

大宗师

F 币
4973 元
贡献
1152 点
14#
发表于 2018-11-16 22:31:46 | 只看该作者
lookbook 发表于 2018-11-15 22:08
版主您好,我去试了下,我的编译器(gfortran8)不行,错误提示是:
procedure:: barsub => barsub4 , ba ...

哈哈,版主也犯错了啊
把后面两个进程里对哑元定义的关键字TYPE换成CLASS应该就好了

79

帖子

17

主题

0

精华

专家

齊天大聖

F 币
433 元
贡献
266 点
13#
 楼主| 发表于 2018-11-15 22:08:47 | 只看该作者
楚香饭 发表于 2018-11-9 00:55
kind的参数必须要编译时确定,不能假定
len的参数可以假定。
所以,多个kind,就要写多个函数。

版主您好,我去试了下,我的编译器(gfortran8)不行,错误提示是:
procedure:: barsub => barsub4 , barsub8
               1
Error: Non-polymorphic passed-object dummy argument of ‘barsub4’ at (1)

79

帖子

17

主题

0

精华

专家

齊天大聖

F 币
433 元
贡献
266 点
12#
 楼主| 发表于 2018-11-10 13:52:33 | 只看该作者
我的理解是,不同的编译器可能对kind的表示不同,所以才会有selected_real_kind(),我用的是gfotran,4就代表real,8就代表double。
我测试了一下:
module style
      implicit none
      integer,parameter :: nnn=8
end module

use style
implicit none
real(kind=nnn) :: a
     
a=1d0
print*,a,kind(a)

end
结果是:
   1.0000000000000000                8

所以是可行的。
不过我在编写程序的时候是严格的,是使用了Selected_Real_Kind()的,
希望它的移植性能更好。
不知道这么理解对不对,如有错误千万不要怪我呀

178

帖子

15

主题

0

精华

大宗师

F 币
4973 元
贡献
1152 点
11#
发表于 2018-11-9 16:20:57 | 只看该作者
lookbook 发表于 2018-11-9 15:58
其实有一个方法,只需要多一个module就行。就是在最底层定义一个kind常量,然后每个module都包含使用它即 ...

哈哈,这个果然是个好办法,根据自己实际用哪个来决定Kind常数
不过这里有一个问题是如果某个程序里同时需要用多种Kind的版本肿么办呢?
感觉上这样做的话似乎实际用的过程中还需要弄个中间产物
比如写个Make文件,可以自动生成各种kind的.o文件和.mod文件
每次定义Type的文件修改了之后,make一下就可以自动产生这三个文件
凡是需要用到这个Type的程序,每次copy最新的.o和.mod文件就可以了
嗯,可惜我最近要写的那个东西需要用到BLAS,要不然就能马上用你这一招了,非常感谢
另外,你这写的不完整吧,定义kind常数的不应该写成:
integer,parameter :: nnn=kind(0.0d0)
或者:
integer,parameter :: Selected_Real_Kind(14)吗?

79

帖子

17

主题

0

精华

专家

齊天大聖

F 币
433 元
贡献
266 点
10#
 楼主| 发表于 2018-11-9 16:00:03 | 只看该作者
楚香饭 发表于 2018-11-9 00:55
kind的参数必须要编译时确定,不能假定
len的参数可以假定。
所以,多个kind,就要写多个函数。

大佬就是大佬,学习了,我去试试。

79

帖子

17

主题

0

精华

专家

齊天大聖

F 币
433 元
贡献
266 点
9#
 楼主| 发表于 2018-11-9 15:58:36 | 只看该作者
liudy02 发表于 2018-11-9 12:38
原来如此,版主的解释让我终于明了Kind参数和Len参数为什么在编程上有不同之处了
写成不同的Type寻找上应 ...

其实有一个方法,只需要多一个module就行。就是在最底层定义一个kind常量,然后每个module都包含使用它即可。

module style
integer,parameter :: nnn=8
end module

module another

type anything
real(kind=nnn) :: a,b
balabala.....
.......
这样的话,关于绑定的函数也只要写一遍就行,但记得更改nnn之后要clean一下,然后再make。
我是这样处理的。

178

帖子

15

主题

0

精华

大宗师

F 币
4973 元
贡献
1152 点
8#
发表于 2018-11-9 12:38:26 | 只看该作者
本帖最后由 liudy02 于 2018-11-9 12:40 编辑
楚香饭 发表于 2018-11-9 11:16
不同的kind,需要用不同的函数。这是肯定的,毕竟二进制代码(指令)完全不同。
如果你有很高的通用编程需 ...

原来如此,版主的解释让我终于明了Kind参数和Len参数为什么在编程上有不同之处了
写成不同的Type寻找上应该是要困难点,所以同一类型不同kind就要分别去找文件
另外如果对与这几种Kind的相关过程有一致性要求,确实也容易顾此失彼,不过反正自己一个人维护的话还好
但是写在不同文件里边,从编辑灵活性上还是好些啊
非常感谢版主的理论解释和经验分享,对我相当具有启发性

712

帖子

4

主题

0

精华

大师

农村外出务工人员

F 币
608 元
贡献
311 点

新人勋章爱心勋章水王勋章元老勋章热心勋章

7#
发表于 2018-11-9 11:16:55 | 只看该作者
不同的kind,需要用不同的函数。这是肯定的,毕竟二进制代码(指令)完全不同。
如果你有很高的通用编程需求,可以考虑用C++,模板会让通用编程更方便。
(事实上,模板也是产生了多个目标代码的函数的)

Fortran在OOP方面还刚刚起步,尚需时间和契机来成熟。
这主要还是取决于需求,至少,我还没有太多这方面的需求。

至于哪个合适,就要自己斟酌了。写成不同type,可能在维护性上会有欠缺。
PS:我在实际应用中,参数化type也很少使用。

178

帖子

15

主题

0

精华

大宗师

F 币
4973 元
贡献
1152 点
6#
发表于 2018-11-9 10:17:10 | 只看该作者
本帖最后由 liudy02 于 2018-11-9 10:18 编辑
楚香饭 发表于 2018-11-9 00:55
kind的参数必须要编译时确定,不能假定
len的参数可以假定。
所以,多个kind,就要写多个函数。

因为IVF和gfortran都没有支持对参数化派生类型中(Len=:)对象的Allocate语句
所以原来我已经觉得参数化派生类型废掉了,
版主却向我展示了它的另外一种用法,确实可以用一种类型涵盖事实上的多种数据
不过这里我有一点点疑惑:
不同的Kind在这里的处理方式实际上还是用的不同的函数
并没有让程序写起来变得更为简洁
从减小单个文件或单个Module的角度,怎么感觉还是写成不同的TYPE封装进不同的Module文件更合适呢
这种参数化派生类型的优势或者是巧妙应用是什么呢……
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-5-19 14:36

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表