- 在线时间
- 1957 小时
- 最后登录
- 2024-6-29
- 注册时间
- 2004-4-26
- 听众数
- 49
- 收听数
- 0
- 能力
- 60 分
- 体力
- 40960 点
- 威望
- 6 点
- 阅读权限
- 255
- 积分
- 23863
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 20501
- 主题
- 18182
- 精华
- 5
- 分享
- 0
- 好友
- 140
TA的每日心情 | 奋斗 2024-6-23 05:14 |
|---|
签到天数: 1043 天 [LV.10]以坛为家III
 群组: 万里江山 群组: sas讨论小组 群组: 长盛证券理财有限公司 群组: C 语言讨论组 群组: Matlab讨论组 |
<TABLE>6 M- R& R. ^9 e' S$ j6 l( J, A
2 u! |+ P K: R1 i1 y; y
<TR>0 b* h4 x0 h( N. d5 H/ T
<TD vAlign=center>
5 _ Y% L- D# q5 o7 s2 n<H1>Fortran 2003:完美还是虚幻?</H1></TD></TR></TABLE>
5 r' ^. Z$ \* v<DIV class="frame note">
3 Z- F" C( W9 C<DIV class=label>Note</DIV>
4 `% P3 F8 S. W<DIV class=content>本文发表于2004年10月《程序员》</DIV></DIV>) w6 W* P0 J/ E! ^6 U8 K: c5 m: z2 |
< >“Fortran 2003!面向对象的Fortran!自由、开放的Fortran!拥抱Fortran 2003吧!”我手舞足蹈、兴奋异常。</P>0 X& O% Z3 E, D. {3 N7 M
< >“Fortran 2003?你发什么颠?”地球物理系一位戴着宽边眼镜的博士生拦住了我,鄙夷的目光从厚重的眼镜片里直射过来,“我昨天才用Fortran 77做完了课题。你说的‘面向对象’、‘自由’、‘开放’什么的,都是你们程序员的鬼把戏吧?它们能让我的SCI论文数量翻番吗?能让我提前评上副教吗?”</P>& j# r( `" t" b C9 T! b% j
< >“……”</P>
# f+ {9 t( \; }0 q% x- D. Y4 O; R<TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>6 F3 f' d' P, ?1 _
9 P/ ]; l8 u% K }<TR>+ M# [! i: F$ ], G) a2 P6 j
<TD width=9 height=10></TD>( E7 l8 \! K" B2 q
<TD>
$ g' W% d# u6 a9 y<H3>1. 不断变革的Fortran</H3></TD>
" \' S" {% H2 d<TD vAlign=bottom align=right>5 I& S' J6 _9 Y9 B. Y9 U
<DIV ><a href="http://www.contextfree.net/wangyg/b/tech/fortran2003.html#top" target="_blank" > </A></DIV></TD></TR>3 g' t( K8 ~9 {: [ O( y6 S8 G
<TR>2 o8 ]% M1 T; l# c5 K5 l
<TD bgColor=#80c4ff colSpan=3><IMG src="http://www.contextfree.net/wangyg/skin/images/spacer.gif"></TD></TR></TABLE>- \# O! N, x% n4 D: w% O
<DIV class=section>& [4 x" T& Y' F
< >2004年5月,Fortran 2003,这个新世纪里诞生的Fortran语言新标准,终于在ISO、IEC的联合工作组JTC1/SC22/WG5以及美国Fortran委员会NCITS/J3的共同努力下,走完了编撰、修订的最后一步。不出意外的话,再过一两年,支持Fortran 2003的开发工具就会充斥Fortran世界;那些仍在坚守Fortran 77阵地的老用户们也会自然而然地把质疑的目光从Fortran 90/95投向更富有革命精神的Fortran 2003。</P>& O. @) ~2 J4 h5 {# Q6 I, u$ F
< >要了解Fortran 2003,我们有必要先简单回顾一下Fortran语言的历史。</P>1 w8 K: V/ K/ `4 j$ [, |' i) s
< >1954到1957年间,世界上第一种高级程序设计语言——Fortran诞生于IBM公司。Fortran这个名字的原意是IBM Mathematical FORmula TRANslation,其设计目的在于为科研人员提供一种符合数学思维习惯的高级语言,以满足科学计算的需要。1960年代,在美苏太空竞赛的大背景下,在国防、教育和科研领域对高性能计算工具的迫切需求下,Fortran语言蓬勃发展,成为了当时统治计算机世界的高级语言之王。</P>
1 a6 v$ x7 c9 i4 w! v' T) h< >1962年,为了统一不同公司、不同硬件平台上的Fortran语言,人们开始了Fortran语言标准化的尝试,这也是程序设计语言发展史上的第一次标准化历程。1972年,Fortran 66标准(标准编号来自标准草案的制定时间)正式发布。但因为标准文档过于简单,约束力不强,Fortran 66标准发布后,Fortran语言的统一问题并没有得到彻底解决。</P>. U9 y! ~9 m4 U/ @5 I! z9 F/ p
< >1978年,Fortran语言标准的第一个修订版本正式发布,这就是我们熟悉的Fortran 77了。Fortran 77细致描述了Fortran语言的各种特征,让Fortran成了一种真正规范、高效和强大的结构化程序设计语言。此后,无数性能优异的Fortran 77编译器和开发工具的问世更是让Fortran 77成为了几乎所有理工科学生的必修课。</P>
a" d) k1 c0 A! s+ c p6 X) e( J< >尽管Fortran 77的影响力一直延续到了今天,但Fortran语言不断变革的历程却从未停止过。为了改变Fortran 77那种老旧的、从打孔卡片遗传来的语言风格,为了给Fortran注入更多的现代特征,人们于1991年发布了崭新的Fortran 90标准。除了自由的代码风格外,Fortran 90还为Fortran语言引入了模块、接口、自定义(派生)数据类型和运算符、可动态分配和参与复杂运算的数组、泛型过程、指针、递归等重要的语法特征。这不但使结构化的Fortran语言更趋完善,也使其具备了少量的面向对象特性。</P>
- _8 r* M6 J5 E< >1997年发布的Fortran 95标准在Fortran 90的基础上,吸收了HPF语言(High Performance Fortran,Fortran语言在并行环境下的一个变种)的优点,提高了Fortran语言在并行任务中的表达和计算能力,并进一步完善了派生类型、指针、数组等要素的相关语法。</P>
0 \* V R S U< >以Fortran 66为基准,我们可以把后续的Fortran 77/90/95以及本文介绍的Fortran 2003均视为对Fortran语言标准的修订。在历次修订中,Fortran 77和Fortran 95是修订幅度相对较小的版本,而Fortran 90和Fortran 2003则是锐意变革的“大修”版本。如果说Fortran 90是在“形式”上让Fortran脱胎换骨的话,那么,Fortran 2003就近乎彻底地解决了Fortran语言现代化的问题:完整的面向对象机制、灵活的语法特征、统一的接口标准……这些难道还不足以让Fortran迷们跃跃欲试吗?</P></DIV>
( ]* Q! J6 F" n<TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>6 Q4 m* J$ F2 H U5 P! k
) P8 u% l7 I) L# ^ U# w9 T' Y/ i<TR>* K& Y" e. T3 k o! C+ D# S
<TD width=9 height=10></TD>+ a6 N7 L" k, `# \
<TD>7 H2 x& d6 G. t' v* t7 x
<H3>2. 面向对象的Fortran</H3></TD>( j3 |: O: a$ w6 Q4 v" }
<TD vAlign=bottom align=right>( r4 @0 e+ z# X
<DIV ><a href="http://www.contextfree.net/wangyg/b/tech/fortran2003.html#top" target="_blank" > </A></DIV></TD></TR>7 N" i8 F+ k# K) d
<TR># W' r% v3 j9 W, w6 G D4 J/ t' H
<TD bgColor=#80c4ff colSpan=3><IMG src="http://www.contextfree.net/wangyg/skin/images/spacer.gif"></TD></TR></TABLE>
5 z j1 ~! _0 k& Q3 Q<DIV class=section>
/ K8 i- ]9 P2 O( ]$ N8 H( O5 L" q< >从结构化程序设计到面向对象,这是大多数高级语言的发展、演变规律,Fortran也未能免俗。在Fortran 77中,由于缺乏对数据抽象和数据封装的基本支持,人们很难引入面向对象的理念和方法。在Fortran 90/95中,尽管我们可以用派生类型(Derived Type)和模块(Module)模拟一部分面向对象的特性,但却无法实现真正意义上的封装和继承。</P>" L( P4 A* U$ h7 v4 r
< >对数据和相关操作的封装是面向对象的基础。传统的结构化程序设计语言大多拥有简单的数据封装能力(如Pascal语言中的record),但缺乏对操作封装的支持。当我们为结构化程序设计语言引入面向对象机制时,在语法层面一般有两种思路:一种是引入一个全新的语法容器,允许用户借助该容器定义出封装了数据和操作的统一数据类型,例如,C++在C语言的基础上扩展出了class的概念,并以class为中心支持面向对象编程;另一种是通过某种机制,将原有的数据封装容器和表述操作的过程或函数捏合在一起,例如,Ada 95利用Ada 83中已有的概念,将private type和相关操作封装在package中,创建了完整的面向对象机制。Fortran 2003对Fortran 90/95的改造正是遵循着后一种思路进行的。</P>
3 n5 U2 p r, D& a( J< >首先,Fortran 2003第一次引入了过程指针(Procedure Pointer)的概念。Fortran 90/95只支持普通指针。为了实现回调机制,Fortran 90/95特意引入了用过程名充当形式参数(Dummy Procedure)的语法,但这种语法还无法解决在数据类型中内嵌相关操作的问题。在Fortran 2003中,我们可以使用过程指针将数据和操作绑定在一起:</P>< RE>TYPE point6 N( F* z" m2 {$ [
REAL :: x, y
1 _) X. y) \" V& | PROCEDURE (testif), POINTER :: test
/ }+ ^) u0 v/ }+ bEND TYPE
- g7 h6 q, Q6 B. \</PRE>' F7 i; Z+ y& ^$ {
< >这种绑定是一种运行时的动态绑定,我们可以根据需要为过程指针赋值,这类似于我们在C语言的struct中嵌入函数指针。</P>; T: e' ~6 r9 T9 x
< >其次,如果与类型相关的操作在编译时已经预知,我们就可以直接使用Fortran 2003提供的类型绑定过程(Type-bound Procedure),将数据和操作关联在一起:</P>< RE>TYPE point
, X( j' M9 v; W7 Z8 e8 O REAL :: x, y: N' R5 L/ |% u5 ]
CONTAINS& A1 s; }6 l; U/ ?+ ^1 \/ g
PROCEDURE :: test => my_test
# d$ {- l# o4 e& M2 E( R+ C! ~END TYPE6 i2 v' ]+ p6 D+ s" I, z2 V
</PRE>- W, n0 d' J# y; z
< >这样,在同一个模块中,类型和相关操作相互关联,用户可以直接使用CALL p1%test或类似的代码来访问特定操作,这已经非常近似于Ada 95的封装机制了。</P>; M; L% [' Q" C7 f1 U
< >和C++或Java语言相比,Fortran 2003的这种过程绑定方式较为靠近语言设计的底层。例如,对于绑定在point类型中的my_test过程,我们必须显示指明,该过程的第一个参数为point类型:</P>< RE>SUBROUTINE my_test(a)
: i" G7 `% T" ?1 L3 k& O' \/ l TYPE(point), INTENT(IN) :: a
! Y* ]7 o3 M' f1 i, _END SUBROUTINE! ?) I6 v! ]# h( a/ T8 N
</PRE>
% P9 s3 e4 ]+ r0 d1 K; K% N< >这里的参数a相当于C++或Java方法中隐含的this参数,其目的是将调用该操作的对象实例传入过程。我把这种显示实现this参数的做法称为“赤裸裸的面向对象”。尽管写起来稍嫌繁琐,但这种“赤裸裸”的做法也有不少好处。例如,我们可以在Fortran 2003中使用NOPASS属性禁止传递对象实例,或通过PASS属性改变对象实例参数在过程参数表中的位置。</P>
/ G! b0 I& \' ^" S* p< >Fortran 90/95已经允许用户自定义运算符过程。现在,Fortran 2003则进一步允许用户将运算符过程绑定到特定的类型,其方法和绑定普通过程类似:</P>< RE>TYPE point
% t% R: q6 p7 B' N, Q REAL :: x, y
8 [! ^& _6 j6 _; J, OCONTAINS$ N; H3 u% x% @1 k9 B0 F
GENERIC :: OPERATOR(+) => plus& N) F/ W1 h! w
GENERIC :: ASSIGNMENT(=) => assign
( v G7 R- H3 V$ U( U7 VEND TYPE
( R3 T4 _5 B+ o( M5 c, A* o</PRE>1 l% j' \: c ]3 {, Y4 G% T
< >有了完善的数据和操作封装机制,Fortran 90/95中简单的访问权限控制属性——PRIVATE和PUBLIC就不能满足要求了。Fortran 2003扩展了PRIVATE和PUBLIC的使用范围,允许PRIVATE类型中出现PUBLIC成员,同时又增加了PROTECTED属性,以限定某个成员的取值只能由自己所属的模块修改。</P>: \) S! ^* y$ V8 n
< >Fortran 2003增强了派生类型的初始化和构造语法。尽管我们在Fortran 2003中还不能像在C++里那样随心所欲地定义各种构造函数,但仍可以较精确地控制对象的初始化过程。此外,Fortran 2003允许我们为类型定义终结子例程(Final Subroutine),以完成对象销毁时的清理工作。终结子例程可以定义多个,像这样:</P>< RE>TYPE point
/ ]( `6 r6 c l$ B: Z REAL :: x, y0 G9 ~0 L& `& J
CONTAINS0 l, e' g; @8 J: Y
FINAL :: finish1, finish2- J, X; X" r- f, Q- U4 P7 h l; H- v9 R
END TYPE
5 c1 N: O: G+ v5 A2 x1 o9 t</PRE>
* _0 G: W" S2 H< >在Fortran 2003中,继承机制是通过类型扩展(Type Extension)实现的。相关的语法和大多数面向对象语言的继承机制没有什么不同:</P>< RE>TYPE, EXTENDS(point) :: funny_point+ K- _- v z: X/ n% `1 D
LOGICAL :: funny = .TRUE.5 H( ~- n8 A9 O# @& G3 I
END TYPE* y$ F. v3 r) g: j* e0 \/ l
</PRE>" i# q& S1 b0 {( P/ h9 ]
< >在扩展类型中,可以直接通过基类型的名称来引用基类型的成员,这类似于C++中引用基类成员的语法。</P>
5 r0 F3 \: R; {9 ?$ R<P>单从封装和继承这两个层面看来,Fortran 2003对Fortran 90/95的扩展就是革命性的。当年,Fortran 90刚问世的时候,许多人迫不及待地宣称Fortran已经是面向对象的语言了。殊不知,Fortran 90/95中的面向对象机制还不如Visual Basic 6.0中的类模块(Class Module)来得成熟,更不要说与真正面向对象的Fortran 2003相提并论了。</P></DIV>+ t: S# W8 ?4 g
<TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>. p% p$ X% d+ [) f. F
5 Z: B$ v& G2 J3 P, t4 s% U4 a! i<TR>
( U& E# ~0 N& d<TD width=9 height=10></TD>4 e4 Y$ b E" W/ `" j6 R
<TD>
( d0 b2 e7 z8 h$ h. V# g9 z; j<H3>3. 自由的Fortran</H3></TD>* t/ f# o* l0 o2 L9 e# a* ?
<TD vAlign=bottom align=right>
5 y) F: W- N& j<DIV ><a href="http://www.contextfree.net/wangyg/b/tech/fortran2003.html#top" target="_blank" > </A></DIV></TD></TR>
. c* P- b6 `8 F. {<TR>
& V6 y7 m& }* [5 e' A<TD bgColor=#80c4ff colSpan=3><IMG src="http://www.contextfree.net/wangyg/skin/images/spacer.gif"></TD></TR></TABLE>* X A1 Z6 C( c" e2 X
<DIV class=section>
3 z3 @8 w5 r. M0 b' V9 o. j<P>灵活性一直是Fortran语言的一个弱项。早期的Fortran甚至不支持内存的动态分配。Fortran 90引入了指针的概念,全面支持数据的动态存取,同时也引入了泛型过程(Generic Procedure)的机制,以支持与C++的函数重载类似的功能。在此基础上,Fortran 2003通过参数化派生类型、覆盖、多态型等机制,进一步增强了Fortran语言的灵活性。例如:</P><PRE>TYPE point(k, n)
* u* \2 G+ a* I& o' O8 k! \+ {7 D INTEGER, KIND :: k = KIND(0.0D0)
7 U& H$ U5 h* l INTEGER, LEN :: n: w J- S" E( X+ d8 }
REAL(KIND=k) :: x, y- x: {( j0 I) D
CHARACTER(LEN=n) :: name, Z" ?) i4 E' |5 S# d% J9 m) @
END TYPE
7 L* F/ q6 G. b* ]0 c/ {' B7 n. u</PRE>
( \. v+ j. ? U: {<P>基于上面这样的参数化派生类型(Parameterized Derived Type),我们可以在定义对象时,根据需要指定对象的具体类型参数,也可以在程序运行期间,动态指定可分配内存对象的类型参数(Deferred Type Parameter)。</P>
6 w* z$ ~) J, ~<P>除了Fortran 90已经实现的重载(Overloading)机制以外,Fortran 2003还允许扩展类型覆盖(Overriding)基类型中的类型绑定过程。更重要的是,Fortran 2003通过抽象类型(Abstract Type)、延迟绑定(Deferred Binding)、多态实体(Polymorphic Entity)、类型选择结构(SELECT TYPE)等要素完整地实现了面向对象语言必备的多态机制。举例来说,我们可以声明下面这样的抽象类型:</P><PRE>TYPE, ABSTRACT :: point 3 O. q3 V5 ^% m$ k- Q7 O/ {
REAL :: x, y& }& K2 u2 C4 J3 \; b+ x: U
CONTAINS
9 M1 n2 h2 G* T3 k" K PROCEDURE(testif), DEFERRED :: test: B7 g8 C f( e" O. b$ |$ ]
END TYPE
+ S- E! g4 h0 m$ o( T( XABSTRACT INTERFACE
0 U4 Q1 G2 |5 t% ~& U# S$ U2 b SUBROUTINE testif(p)2 Y1 U, F: [+ f0 u4 B x
CLASS(point), INTENT(in) :: p) W. z+ Y; L) g5 N* P9 x2 Z& R
END SUBROUTINE
4 j% q0 J$ Y& b5 E: D, Y) x5 sEND INTERFACE* o0 b6 [: d5 A# `
</PRE>
R3 d% `' A' {3 N4 h5 x0 [<P>上述代码相当于在C++或Java语言中声明了一个抽象类,其中的test过程相当于C++语言中的一个纯虚函数,只有接口声明,而没有实现代码。在此基础上,我们可以从抽象类型point扩展出具体的子类型。</P>
- d N5 d5 r( Q<P>当某个对象的具体类型在运行时才能确定时,我们可以用CLASS关键字(这里的CLASS指的是一组有继承关系的类型,与C++中的class并不等同)定义一个具有多态功能的指针或可分配对象:</P><PRE>CLASS(point), POINTER :: p$ `$ f, _ K" n) v
</PRE>8 ~) ~2 d; B- o' c7 O( D. o
<P>现在,指针p可以指向任何由point扩展的具体类型的对象,这种使用方式已经和我们在C++中用基类指针引用派生类对象的做法一模一样了。使用这种具有多态性特征的实体时,我们还可以用Fortran 2003提供的内部过程SAME_TYPE_AS和EXTENDS_TYPE_OF来判定对象的类型,这是一种简单的运行时类型识别(RTTI)机制,其功能和C++语言中的typeid关键字相仿。</P> H- {; B5 j& b; H) ~
<P>更有趣的是,Fortran 2003允许代码根据多态实体的实际类型,在运行时执行特定的程序流程:</P><PRE>CLASS(point), POINTER :: p* {$ R* s6 Q% L7 p6 {' Q
p => a
+ A$ v. r: E: h0 I( pSELECT TYPE ( pp => p ): z9 q4 ?) \1 g6 f1 z3 r
TYPE IS (point_3d)/ L' A0 }1 S) {, M4 h5 B
PRINT *, pp%z$ u* P2 C2 U! m
TYPE IS (point_2d) V/ ~- k7 K7 H; a: H5 g0 ]
PRINT *, pp%x
9 t6 ~: B" I8 T% [( ~" KEND SELECT- A( {& R+ P! l1 q0 p; A
</PRE>' ]$ n8 Y1 q' q0 E; x( s# I' X3 k
<P>这种语法在程序设计语言中比较少见。在Java或C#等具备较强的RTTI和Reflection功能的语言中,我们可以用switch语句结合类型标识实现类似的功能,如下面的C#代码:</P><PRE>void foo(Point p) {) o+ L/ D9 l( C. d: P, G
switch(p.GetType().FullName) {
( b- g- M; K- d& h( Q case "Science.Point2D":
% T; b S4 {# P2 } MessageBox.Show("2D point");; @2 N3 Y' i T+ T8 _, |; m
break;# x, I2 O! g) j
case "Science.Point3D":
+ q9 M S% l0 r8 Y! y$ S3 r MessageBox.Show("3D point");& ?7 |! a! f9 ^
break;
* L1 E {+ \" ~" I }1 v) S/ M% ?& }) e: I
}
& w& P- \; m( A</PRE>8 w# q! B( s" O1 [/ U
<P>但在Fortran这样不依赖于虚拟机和中间代码的纯编译型语言中,能够拥有如此方便的SELECT TYPE语句,实在是Fortran用户的一大幸事。</P># Z& T7 U' `; T0 ^3 c$ Q, O) U! a. b
<P>其他一些较小的语法功能扩充也能反映出Fortran 2003制定者对语言灵活性的追求。比方说,ASSOCIATE语句可以让Fortran用户体验到比C语言的宏定义更为便捷的复杂表达式替换功能:</P><PRE>ASSOCIATE ( Z => EXP(-(X**2+Y**2)) * COS(THETA) )
5 e5 z" p7 y: S( {% s' K) G& y PRINT *, A+Z, A-Z
: ^3 _' Q" c+ [END ASSOCIATE5 d) I O% M- \: S) S4 y
</PRE>4 {- h+ H% E }6 ~# M i% |
<P>再比方说,在Fortran 2003中,可以拥有ALLOCATABLE属性的变量已不再限于数组对象,参数化派生类型的对象也可以在定义时省略参数,而在动态分配时指明:</P><PRE>TYPE(point(KIND(0.0D0), n=20)) :: a
. v& O/ M3 Q4 _3 fTYPE(point(KIND(0.0D0), n= ),ALLOCATABLE :: b: {( z$ k/ T! I
ALLOCATE(b,SOURCE=a)
, s* V& q7 V9 p/ ^6 K/ k$ j+ o( V) J</PRE>" h0 s/ S1 i H# C$ d- R
<P>对于科学计算语言中最重要的数组(矩阵)操作,Fortran 2003也提供了更多的新功能。例如,为可动态分配的数组变量赋值时,目标数组的形状可以根据源数组的形状自动调整,并自动完成内存空间的重新分配,用户也可以使用新增的MOVE_ALLOC内部过程强制改变某个数组的形状:</P><PRE>REAL, ALLOCATABLE :: a( , temp( , ^! S. k! q. k4 t. k% S& G! ~
ALLOCATE(a(-n:n), H- d; g# K& c% m
ALLOCATE(temp(-2*n:2*n))
; W3 i5 _: s) o7 w/ t4 ttemp(::2) = a# `' k8 i/ C, n" v$ H9 `- T; S
CALL MOVE_ALLOC(TO=a, FROM=temp)
. _9 e5 A: s& B3 R- B</PRE># f( I) [% \! u# G; W9 ]9 `
<P>Fortran 2003允许我们为指针类型的参数增加INTENT属性,这为那些与动态数据结构(如链表)相关的算法程序提供了更多的选择空间。Fortran 2003引入了IMPORT语句,允许接口声明中引用接口所在模块内的相关定义。Fortran 2003还允许我们在USE语句中,改变被引用模块内自定义运算符的名称。</P>
! @ q! A9 D9 c# q' W9 V) J<P>总之,借助Fortran 2003提供的更加灵活的语法特性,我们能更容易地编写出可以适应不同数据类型、不同应用环境的可复用代码。对于Fortran语言拓展应用领域、发掘潜在用户的目标而言,这些灵活性都是至关重要和不可或缺的。</P></DIV>6 w7 F1 @: I) u) @9 l: K$ e$ j' v
<TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>; H! y' j$ V/ n( P/ ?" d- o
( T$ A4 z0 `# [+ u& w [! o: d
<TR>, i! h; A& q8 r
<TD width=9 height=10></TD>- Z* s+ \( y3 }, P/ e
<TD>3 H( P0 W4 S' F! D$ _; t
<H3>4. 高性能的Fortran</H3></TD>
4 a) y+ n3 H# h8 n0 v4 o; S<TD vAlign=bottom align=right>
! A1 d: c3 M5 K# \6 v6 S4 E<DIV ><a href="http://www.contextfree.net/wangyg/b/tech/fortran2003.html#top" target="_blank" > </A></DIV></TD></TR>$ Z4 Q, e) _* @; H. H$ `/ L7 u
<TR>$ Y) ^ R9 F9 s3 v
<TD bgColor=#80c4ff colSpan=3><IMG src="http://www.contextfree.net/wangyg/skin/images/spacer.gif"></TD></TR></TABLE>
1 p% A! G) @4 C7 r6 I2 u% E<DIV class=section>
1 @* d) k8 b) B$ j<P>与其他科学计算语言或工具相比,Fortran语言在性能上历来是出类拔萃的。因为语言本身专门针对数值计算、矩阵处理等功能进行了优化,大多数Fortran编译器产生的可执行代码在效率上甚至会超过以高效著称的C或C++语言。为了在高性能的并行处理系统(如IBM的“深蓝”和我国的“曙光”)上获得更出色的执行效率,Fortran 95还特意吸收了HPF语言的优点,为Fortran语言增添了若干支持并行计算的语法特征(比如著名的FORALL语句和PURE过程)。</P># z' d L* C# A
<P>继续改进Fortran语言的计算性能当然也是Fortran 2003的任务之一。一个最明显的例子是,Fortran 2003引入了VOLATILE属性。这个属性类似于C语言或Java语言中的volatile关键字。编写过并发或实时应用的程序员都知道这个含义为“易变”的关键字的价值:在并发系统中,如果没有这个关键字的帮助,我们就必须时刻警惕共享数据的取值是否已被正确刷新。</P>6 Z( E. T1 |% r5 ] P
<P>不过,相对而言,Fortran语言本身的并行计算机制仍不算十分健全。编写并发或实时程序时,Ada语言中的任务(Task)管理和同步(Synchronization)机制,或是Java语言中的多线程同步特性都可以为程序员提供更有力的支持。</P>1 q- K+ n L1 I& ^! v3 ?5 R! p
<P>实际上,Fortran 95和Fortran 2003陆续引入并行语法的目的之一是消除标准Fortran语言与已经存在并得到广泛应用的各种Fortran语言变种之间的隔阂,允许同一份Fortran代码在不同的语言环境间移植。至少到Fortran 2003为止,标准Fortran语言还没有能力完全替代以HPF为代表的“高性能”Fortran语言变种。仅就适应高性能并行计算环境的能力而言, HPF等Fortran变种也仍有足够的理由继续存在和发展,直到未来某一个大而全的Fortran标准把它们真正统一起来为止。</P>
# E$ T0 ~! b; E0 A- d+ O* l- ~% ^<P>为了适应现代CPU的时钟精度,Fortran 2003扩展了内部过程SYSTEM_CLOCK的功能,允许其COUNT_RATE参数为INTEGER或REAL类型。</P>! n% m! a: @# F% H# ]2 g
<P>对IEEE浮点数标准(IEEE 754)的完全支持是Fortran 2003在提高计算精确度和规范性方面的又一个努力。在Fortran 2003中,IEEE标准中描述的下溢(Underflow)可以被正确识别和处理,内部模块IEEE_ARITHMETIC中也为此增加了IEEE_SUPPORT_UNDERFLOW_CONTROL、IEEE_GET_UNDERFLOW_MODE、IEEE_SET_UNDERFLOW_MODE等几个相关的内部过程。</P>
8 t9 z) R1 V: B4 o<P>另一个可以大幅提升程序执行效率的改进是Fortran 2003新增加的异步I/O机制。也就是说,执行比较耗时的I/O操作时,程序不必阻塞在I/O语句上,可以继续执行后续的指令。这种机制类似于我们用WIN32函数操作磁盘文件时,可以借助FILE_FLAG_OVERLAPPED标记指明I/O操作应异步进行。</P>9 p8 d* ?' B1 h" p% X. d
<P>基于Fortran 2003,当我们在外部文件的OPEN语句中指明ASYNCHRONOUS=’YES’时,后续拥有ASYNCHRONOUS=’YES’属性的READ或WRITE语句即以异步方式执行。此后,我们可以继续执行其他语句,或使用WAIT语句等待异步I/O执行完毕。</P>
5 Z% q7 n5 s. T7 e& B+ _<P>异步I/O是高性能应用软件必备的特性之一。在Fortran 90/95的时代里,为了实现异步I/O功能,不同的开发商总会在自己的Fortran编译器或程序库中以各自不同的方式,增加对异步I/O的支持。Fortran 2003统一异步I/O语法的做法显然有助于Fortran应用的移植和推广。</P></DIV>
0 M1 r2 k- g4 ]9 t/ g, x<TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>3 h1 [6 m! Y4 `6 s9 X
9 q. h$ k( l# B, t' B C<TR>1 G8 j1 D) t& q; [2 { q
<TD width=9 height=10></TD>" Q1 X$ r0 I4 {3 K$ ?- x' Z; q
<TD>; j3 M6 Q" Y w5 C4 F6 j
<H3>5. 开放的Fortran</H3></TD>7 j3 |2 }/ }% W+ V$ i
<TD vAlign=bottom align=right>
% h) n+ R8 _+ E- D1 d<DIV ><a href="http://www.contextfree.net/wangyg/b/tech/fortran2003.html#top" target="_blank" > </A></DIV></TD></TR>
% c' v" l ?( y$ d3 h0 d/ p0 i% ]7 ?5 D<TR>
4 a' m9 a2 @9 D+ L5 Z& o<TD bgColor=#80c4ff colSpan=3><IMG src="http://www.contextfree.net/wangyg/skin/images/spacer.gif"></TD></TR></TABLE>
1 f* ^+ D6 U/ z, }- j<DIV class=section>
5 ~! t" ` V u9 s3 M<P>传统意义上的Fortran语言看上去并不那么开放,这恐怕和Fortran语言向来只注重科学计算,而忽视语言通用性的习惯有关。中国的Fortran用户最常抱怨的两件事情是:不同的Fortran编译器连接外部程序(如C语言程序或Matlab程序)的方法不尽相同,不同的Fortran编译器对中文的支持能力也强弱有别。人们迫切希望制定Fortran标准的专家们能体谅到最终用户的苦衷,尽早将Fortran纳入开放、规范的发展轨道。</P>0 K e9 j, J9 M) {( y0 v) p
<P>在开放性方面,Fortran 2003为我们带来的第一个福音是,新标准终于规范了Fortran语言与C语言的连接方式。Fortran 2003提供了一个名为ISO_C_BINDING的内部模块,该模块定义了Fortran与C语言连接时必需的类型常量。引用了ISO_C_BINDING模块后,我们就可以将Fortran变量定义成与C语言数据结构兼容的数据类型:</P><PRE>REAL(C_FLOAT), DIMENSION(100) :: ALPHA) T7 g$ ]# ?3 ]1 l) w
</PRE>, b H. ~1 J& R4 N8 x2 l b2 T; T6 M
<P>或者利用C_PTR、C_FUNPTR类型与C语言中的指针或函数指针交互操作:</P><PRE>TYPE(C_PTR) :: BETA1 G0 d5 G4 B3 ^5 Q
IF (.NOT. C_ASSOCIATED(BETA)) THEN5 K6 _- g a8 u' f7 q- P
BETA = C_LOC(ALPHA)7 A# {8 C Y$ ^% p7 E3 _
END IF
3 M5 F; M- y7 O. ?: P5 S. p</PRE>
% C; w' a2 G6 z! W<P>对于下面这样的C语言函数:</P><PRE>int foo(void* buf, int count, int *ret)$ _/ t- q' a" z$ ]
</PRE>
5 k# W/ T7 ~" n% t<P>它对应的Fortran函数接口是:</P><PRE>INTERFACE
- s. p" W) j9 ]4 I INTEGER (C_INT) FUNCTION foo &; L" ~6 A }; [, {
(buf, count, ret), BIND(C, NAME='foo')
2 s8 T- j) D9 [( L0 g USE ISO_C_BINDING
7 F7 V2 M1 J# n- y TYPE (C_PTR), VALUE :: buf6 I# x% U! X4 s2 a
INTEGER (C_INT), VALUE :: count
7 d! U: h" F& d3 }3 M D4 H TYPE (C_PTR), VALUE :: ret
$ {- Y$ z1 x1 V6 Z END FUNCTION
; d6 M! j5 _ k. X" @2 |END INTERFACE
' c$ n1 h1 j/ |</PRE>$ L& _! K$ _/ m9 o( i2 V8 z
<P>这样,C语言和Fortran语言就可以通过上述接口相互调用。当Fortran语言和C语言之间拥有了统一的连接方式后,Fortran语言与其他语言之间的连接也就不再是困难重重的事了。借助ISO_C_BINDING模块,Fortran语言可以直接(在二进制层面连接)或间接(以C语言为媒介)地与Delphi、C++、Ada、Java、C#等我们常见的通用编程语言交流、沟通。</P>
' h% L: s5 z4 Z0 f4 h+ `8 S<P>在国际化方面,Fortran 90已经为Fortran引入了不少国际化支持功能。Fortran 2003则全面支持ISO 10646标准,可以正确处理双字节或四字节的国际字符集。Fortran 2003程序可以通过SELECTED_CHAR_KIND内部过程指明字符数据的编码方式;打开文件时,可以在OPEN语句内利用ENCODING=’UTF-8’或类似的方式指明文件内字符的编码;可以在读写文件时完成ASCII、ISO 10646、UTF-8之间的转换;Fortran 2003代码本身的字符集,以及标识符和字符串常量的长度范围也有所扩充。</P>4 H: N5 Q% q1 e& c+ p) S
<P>另一个和开放性相关的新特性是Fortran 2003对流I/O(Stream Access Input/Output)的支持。作为一种以科学计算为目的的语言,传统的Fortran只提供了读写以记录为单位的外部文件的语法。在单纯的科学计算任务里,这种简单的I/O功能也许可以满足文件读写的需要,但当我们希望Fortran程序和外部应用交换数据的时候,有时就不得不面对如何读取没有固定记录结构的字节流数据的问题了。在Fortran 2003发布以前,我们只能通过编译器提供的扩展功能解决问题。现在,我们可以用标准的方式,在Fortran语言中随机访问外部文件的字节数据。应当说,直到Fortran 2003,Fortran语言才拥有了与C语言的I/O库大致相仿的I/O功能,才在I/O领域具备了通用语言的基本特征。</P>6 n0 H8 V5 R6 W6 H
<P>当我们在C语言中通过main函数的argc和argv参数获取命令行参数的时候,也许我们不会想到,在传统的Fortran语言中,要完成类似的任务,Fortran用户也必须求助于编译器提供的扩展功能。Fortran 2003显著增强了Fortran语言与外部环境的沟通能力。新引入的内部模块ISO_FORTRAN_ENV拥有INPUT_UNIT、OUTPUT_UNIT、ERROR_UNIT这一组对应于标准输入输出的常量(相当于C语言里的stdin、stdout和stderr),拥有GET_COMMAND、GET_COMMAND_ARGUMENT等获取命令行指令,以及GET_ENVIRONMENT_VARIABLE等访问环境变量或系统参数的内部过程。这些都是Fortran语言走向开放的重要标志。</P></DIV>
2 V4 ?3 W4 j. x<TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>0 _$ h6 s6 _. \% J
- S& @9 |' X M<TR>
/ ?, Q2 c* F% f* L<TD width=9 height=10></TD>; {3 c" v2 r+ {) n
<TD>" i7 ]4 w$ ~ L Z, [
<H3>6. 象牙塔里的Fortran</H3></TD>
* i/ e7 ~8 w# v3 O3 X<TD vAlign=bottom align=right>
9 q- T0 w* z' Z, t9 r<DIV ><a href="http://www.contextfree.net/wangyg/b/tech/fortran2003.html#top" target="_blank" > </A></DIV></TD></TR>
1 e, Y- a$ K: Z! e# R<TR>
! ^ W# m) {7 y s3 E) F<TD bgColor=#80c4ff colSpan=3><IMG src="http://www.contextfree.net/wangyg/skin/images/spacer.gif"></TD></TR></TABLE>
, W0 ]6 R' Q4 f! _4 D S<DIV class=section>' Q; A* B$ a, \6 _8 h7 _$ N
<P>无论Fortran 2003多么强调面向对象、灵活和开放,Fortran语言在根本上仍然是一种面向科学计算的高级程序设计语言。</P>
9 r0 Q, O2 g! W* o3 j<P>环顾一下四周,我们很容易发现,Fortran语言最常出现的地方也是科技创新最活跃、知识密集度最大的地方:在中科院大气物理研究所里,研究人员们正使用Fortran语言编写大气数据分析软件;在国家高性能计算中心里,Fortran语言正在新安装的曙光并行计算机上执行着各种科学计算任务;在一个又一个分子生物学、高能物理学、应用数学的国家重点实验室里,Fortran是研究者们最为倚重的工具之一;在相当数量的理工科学生的课程表中,Fortran是学生进入大学课堂后接触的第一门高级语言……</P>6 m! ?: m8 m, a# E( M- O
<P>但正如本文开头那位博士生所说的那样,国内的Fortran语言应用多半还局限在Fortran 77的层面,大多数研究者仍在使用打孔卡片式的代码风格编写Fortran程序。这件事背后的潜台词是:一方面,Fortran 77已经可以满足大多数科研人员的需要,另一方面,科研人员不是专业的程序员,不可能有充足的时间学习新语言。在这种情况下,我们有什么理由让科研人员学习更复杂的面向对象语法并转向更强大的Fortran 95或Fortran 2003呢?</P>) ~1 P! ?0 |2 a5 H. N
<P>我个人觉得,Fortran 2003为所有Fortran用户带来的不仅是新鲜的语法特性,也是一种新的选择和机遇。</P>! R' p2 A: r+ s+ ]& m# T) u( y
<P>对于利用Fortran语言开发商业软件包的公司(国内这类公司较少,但国外有相当数量的软件公司在开发和销售用Fortran语言开发的工具或程序库)而言,Fortran 2003提供的面向对象机制和开放性可以帮助他们优化大型软件的开发过程,提高软件的复用率,降低开发成本。</P>+ j' ?& {7 F3 W3 e6 ?5 ?
<P>对于将Fortran视为工具的科研人员来说,Fortran 2003既为我们提供了一种更自由的编程环境,也为我们提供了足够的向前兼容的保证——除了极少数“过时”或“危险”的语法特性被删除以外,Fortran 2003保持了与Fortran 77/90/95等语言版本的全面兼容。也就是说,你完全可以在Fortran 2003的环境中,一边编写你熟悉的Fortran 77代码,一边体验新的功能特性,直到你在不知不觉间喜欢上Fortran 2003为止。</P>
+ s& c, W$ o2 F1 ~$ F! n<P>对于理工科的学生来说,Fortran 2003也没有理由被拒之门外。与老旧、呆板的Fortran 77相比,Fortran 2003的语法更现代,更容易掌握;与C++、Java、Pascal等通用语言相比,Fortran 2003中强大的数据处理和矩阵运算功能显然可以帮我们节省许多开发时间;与Mathematica、Matlab、Macsyma,Mathcad等更“高级”的数学工具相比,Fortran 2003尽管在便捷性和可视化程度上略显不足,但Fortran语言编译后的执行效率却是这些工具软件难以企及的。</P>8 z1 [4 l g: k4 |( ^
<P>客观地说,Fortran 2003是一种在各方面都趋于“完美”的科学计算语言,是Fortran语言走向“现代化”的重要标志;不过,这种现代化改造大都是以Ada、C++、Java等现代主流程序设计语言为参照进行的,它们是否能与Fortran语言重视科学计算的初衷相适应,是否能赢得大多数科研人员的认可,还需要经过时间的检验——至少,我们应当等到HP(Compaq)、IBM、Salford、Intel、Fujitsu等主要的Fortran工具提供商推出了支持Fortran 2003的新版编译器,等到用Fortran 2003实现的程序库和应用系统越来越多的时候,再对Fortran 2003的表现做一个中肯的评价。</P></DIV>! y( b, T' l7 ^3 u T8 v
<TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
' g ?; a3 o+ @! F/ G* J/ H
/ O0 N" [7 L$ I: m3 {) b<TR>
! k! z! L$ O3 \4 `<TD width=9 height=10></TD>: c; I2 M1 G* S* q
<TD>
# A: z$ h6 F- L( n+ c1 A' R: f9 s<H3>7. 明天的Fortran</H3></TD>- i' T Q+ e0 Y; }! ?
<TD vAlign=bottom align=right>
# F5 v- M3 g4 k( U% D/ l- L<DIV ><a href="http://www.contextfree.net/wangyg/b/tech/fortran2003.html#top" target="_blank" > </A></DIV></TD></TR>0 L( i0 w$ L5 n: _/ B
<TR>
! _' A; C- _8 v4 C<TD bgColor=#80c4ff colSpan=3><IMG src="http://www.contextfree.net/wangyg/skin/images/spacer.gif"></TD></TR></TABLE>: Q' ?. F/ V% H* o! _4 y9 Y
<DIV class=section>
, S9 W$ o$ K P0 F* u- x) S<P>根据JTC1/SC22/WG5工作组和NCITS/J3委员会的规划,2004年5月到2009年8月是下一代Fortran标准的酝酿、编撰和定稿时间。到了2010年时,Fortran语言会变成什么样呢?</P>
+ ~! T4 o9 I7 N0 y3 P<P>在高级程序设计语言的发展史中,我们可以很容易地找出两种截然不同的语言风格。</P>, t5 r* A0 V0 f! E* G7 _; X
<P>第一种是LISP语言开创的“实验风格”,这种风格强调高级语言的语法逻辑应尽量服从计算机表达数据和操作的基本方式,以充分发挥计算机的“智慧”。例如,LISP语言中用“广义表”来统一描述代码和数据,以及Smalltalk语言用“对象-消息”这样的基本语法结构来实现面向对象机制的做法都是这种风格的完美体现。不幸的是,“实验风格”的语言尽管在表达能力和逻辑性上表现出色,但因为更接近电脑而不是更接近人,它们在与普通程序员的沟通上总显得力不从心,无法得到大多数程序员的认可。</P>
7 E$ E3 s9 w0 M<P>第二种是Fortran语言开创的“实用风格”,这种风格强调高级语言的语法应尽可能接近人类的自然语言或特定领域(如科学计算)的思维习惯,以简化程序员的编码或维护工作,降低程序员的学习难度。直接受Fortran影响的BASIC语言,以及间接受Fortran影响的C、Pascal、Ada、C++、Java、C#语言都是这种风格的杰出代表。在现代软件产业中,“实用风格”的语言占据了绝对的优势,但随着软件规模的增大,它们也逐渐暴露出了灵活度不足、表达能力不强、与计算机交流不顺畅的弱点。</P>
8 i: J( U5 I3 q5 m<P>现在看来,“实验风格”和“实用风格”在一定程度上的相互融合、相互促进是未来程序设计语言发展的总体趋势。这一点可以从Python、Ruby等新生代脚本语言的蓬勃发展中见到一些端倪。</P>
, L c) H7 T& [9 `+ T! f' L1 R+ e<P>根据这样的判断,Fortran语言在未来的几年里,是不是也会向更灵活、更富有表达能力的“实验派”靠拢呢?也许是,也许不是。但无论如何,Fortran语言注重科学计算、注重执行效率的基本特征都不会改变,语言设计者们仍要在Fortran语言天生的“科学属性”与“现代化”改造的矛盾中谨慎前行——没办法,谁让Fortran语言生来就活脱脱是一个爱科学、讲科学、用科学的“动脑筋爷爷”呢?</P></DIV>7 q4 m- ]! q" M" Y% w* ?( u
<P>[王咏刚,2004年8月]</P> |
zan
|