Fortran语言的动态链接库及同其它语言的混合编译

FORTRAN语言动态链接库的使用及混合编译

转载自网上搜索到的一篇介绍FORTRAN动态链接库的资料——《Fortran语言的动态链接库的开发及应用》，同上海汽轮机期刊中《水和水蒸汽热力性质动态链接库的开发及在汽轮机中的应用》一文一致。仅供参考。

1 引言
由于FORTRAN在开发图形界面、可视计算、面对对象、事件驱动、数据库应用方面存在局限性，如果我们利用其它语言来重写，将费时费力，同时又损失了FORTRAN语言在科学计算上的优越性。经过努力，笔者寻求到一种折衷的方法：利用Microsoft Fortran PowerStation 4.0将FORTRAN源代码编译为动态链接库，然后用其它语言建立应用程序的项目文件，在适当的时候调用该动态链接库。这样就把DOS环境中的软件包括移植到Windows环境，使问题得以圆满解决。
以下以最常用的子例程PT为例来说明动态链接库的建立和在各种流行程序语言中的调用方法。

8 |0 T* y/ y$ a. y2 基于Fortran语言的动态链接库的建立
2.1 建立动态链接库的项目
在Fortran PowerStation开发环境中，建立动态链接库的项目文件(IFCSW.MDP)。
2.2 编写动态链接库的源代码
6 K/ T3 u9 b) L/ K  V' r编写动态链接库的源代码文件(IFCSW.F90)，并将源代码文件插入到动态链接库的项目文件中。
在动态链接库源代码(IFCSW.F90)中子例程PT的定义为
( ^, _! k( d4 NSubroutine PT(P,T,H,S,V,X,R)
!MS$Attributes Stdcall , Dllexport , Alias : ’PT’ :: PT
!MS$Attributes Value :: P,T
. e$ y3 n- u; r+ k!MS$Attributes Reference :: H,S,V,X,R
. d0 I" e* g) J( X* }" rReal(4) P,T,H,S,V,X
# @3 D5 c6 M5 S5 J! dInteger(2) R
…………
+ t- {. u$ U! s& SEnd Subroutine
- |% T" Y6 V0 w4 s* \
' \* z$ C7 @- U. u8 M+ v, C. `$ _! b其中：P为压力，T为温度，H为比焓，S为比熵，V为比容，X为干度，R为所处的区间；如果调用前已知所处的区间，可直接给R赋区间号后再调用子例程，可避免区间判断而加快调用速度。当所处的区间未知或不确定时，由子例程内部进行区间判断后，针对不同区间再调用不同的其它例程，并将所处区间号传回。
由于编译后动态链接库中命名方式的差异，为避免在VC++中隐式加载该DLL时产生连接错误，则专为VC++定义如下的子例程(仅供VC++使用)，即不要使用命名别名的Alias语句：
Subroutine PTVC(P,T,H,S,V,X,R)
!MS$Attributes Stdcall , Dllexport :: PTVC
/ _1 V2 J% z) J4 i& n& n!MS$Attributes Value :: P,T
!MS$Attributes Reference :: H,S,V,X,R
Real(4) P,T,H,S,V,X
$ F3 s/ m4 b, V& o7 I* C- ~& }# ZInteger(2) R
, B, t) u; |' V/ u: C2 T: V$ }( K…………
$ z) O  ^# O. t( f: @; C  E6 r6 ~2 TEnd Subroutine

' X! U5 t3 Z& Z& [5 m& }" i0 Q/ s2.3 创建动态链接库
1 e4 b8 R/ Q1 E0 s6 m) |$ E对动态链接库的项目文件进行编译连接，获得动态链接库(IFCSW.DLL)和接入库(IFCSW.LIB，该接入库通常用于隐性加载动态链接库，且只能用于Microsoft的开发语言，对于其它公司的产品，需利用专用程序生成对应的接入库)。
: }! k; h# m) [, Q- a1 G
3 动态链接库的使用
一旦生成了动态链接库，就可以在各种编程语言环境中调用该动态链接库中的各种子程序或子例程。
) d; N" p" n9 F1 ?- u3 G. W动态链接库有隐式和显式两种装载方式，各有优缺点。采用隐式比显式友好，需要较少的编码，装载快速、可靠，但缺乏控制，影响动态使用，资源使用不经济；采用显式则正好相反。
$ v2 @+ y5 M& _8 x2 {- }/ X采用隐式通常需要接入库和函数申明。
采用显式则应遵循以下步骤：
  j  }! y. u& p+ ]8 j* ^(1) 创建一个函数指针，其指针类型与要调用的DLL引出函数相吻合。
(2) 通过Win32 API函数LoadLibrary()显式地调用DLL，此函数返回DLL的实例句柄。
/ I, _5 c8 {0 Z& D5 `(3) 通过Win32 API函数GetProAddress()获得要调用的DLL的函数地址，把结果赋予自定义函数指针类型的指针。
$ o5 Y, ~6 p8 @9 i8 d+ T% s(4) 使用函数指针来调用DLL中函数。
(5) 调用完成后，通过Win32 API函数FreeLibrary()释放DLL。
7 d4 L  p' K- S( _) n
实际使用中可将隐式和显式两者结合使用，采用混合模式。

3.1 在MS Fortran PowerStation 4.0中调用
① 隐式调用
将接入库(IFCSW.LIB)包含到项目文件中，另在源文件的开始部分加入以下的声明语句
Interface
. w/ L/ ]. d9 B2 _Subroutine PT(P,T,H,S,V,X,R)
!MS$Attributes Stdcall , Dllimport , Alias : ’PT’ :: PT
; v7 c$ a4 @2 N' p1 y' XReal(4) P,T,H,S,V,X
9 C2 k$ o, L+ {' ]9 w7 U: BInteger(2) R
& h$ P. F% C* n' l, W!MS$Attributes Value :: P,T
, \6 G8 u9 w* V1 B!MS$Attributes Reference :: H,S,V,X,R
4 s& f6 [& F  e  f. kEnd Subroutine
. N0 P: v& M6 I8 ^7 [End Interface
此后，在需要的地方可用以下语句调用
Call PT(P,T,H,S,V,X,R)

; o- i* L" \; v, h) o% v② 显式调用
+ y: z5 C( x7 U, V9 C动态装载首先要用LoadLibrary函数将动态链接库装入；然后用GetProcAddress函数获取动态链接库中要调用的函数或过程的地址；由于Fortran 90没有提供函数指针，不能直接用地址进行函数的调用，通过编写一个用函数名或过程名作参数的子程序CallProc来解决；将函数名或过程名作为参数传递时，实际传递的是四字节的函数或过程的入口地址，CallProc子程序接收该地址后，根据所声明的函数接口，以正确的方式对函数或过程进行调用；动态链接库使用完后，要及时用FreeLibrary函数释放动态链接库。
Subroutine CallProc(Pn,P,T,H,S,V,X,R)
- n+ i" t0 j5 I' z! O* NInterface
1 R. h. ^! A! ?5 ~  TSubroutine Pn(P,T,H,S,V,X,R)
2 K* U. @# ^2 A!MS$Attributes Stdcall , Dllimport :: Pn
8 q, G& q; |& g& S2 lReal(4) P,T,H,S,V,X
Integer(2) R
' d6 x6 j" y4 K! s9 T2 K!MS$Attributes Value :: P,T
  x4 x9 S! ?: A% ]8 N! X!MS$Attributes Reference :: H,S,V,X,R
End Interface
Real(4) P,T,H,S,V,X
# ?! |8 C: I/ [Integer(2) R
Call Pn(P,T,H,S,V,X,R)
End Subroutine
/ @- e9 N* V6 U2 J( f5 m
# P; A, o9 k6 l  ?9 ?& c2 sCallProc(PT, P,T,H,S,V,X,R)

3.2 在Visual Basic和Visual Basic Application中调用
( a# s: r6 ~! Y9 D2 ~VB和VBA中由于显式加载过于繁琐，通常不使用，以下仅介绍隐式加载。
. E4 \9 Q! `/ Y9 [) y0 cVB和VBA中隐式加载不需要接入库，只需在模块代码的声明部分增加Declare语句，声明需调用的函数在外部动态链接库中。
Decalre Sub PT Lib “IFCSW.DLL” (ByVal P as Double , ByVal T as Double , ByRef H as Double , ByRef S as Double , ByRef V as Double , ByRef X as Double , ByRef R as Integer)
3 }, E9 w) R& `而后，即可在任何需要的地方用以下语句调用子程序PT。
CALL PT(P,T,H,S,V,X,R)
7 C# c; C3 r0 n$ u# W( Y$ B! M
3.3 在Delphi(Visual Pascal)中调用
3 C( {, L" N# d) x! Y2 L