这里将详细介绍一下在C++中如何调用用Fortran语言编写函数的问题,即Fortran与C++的混合编程问题。
通常情况下,C++与Fortran的混合编程问题是利用动态链接库的方式进行的,换句话说,如果在C/C++中调用Fortran函数,应将Fortran函数按照一定的协议编译为动态链接库,然后即可实现二者的混合编程问题。实现这一目的有两种方式:显示链接与隐式链接,下面通过两个非常简单的实例分别介绍之。
编译环境:Fortran:推荐使用Compaq Visual Fortran 6.0以上的版本,C++:Microsoft Visual C++6.0。
1. 显式链接,
(a) 找开CVF编译器,然后新建一个Fortran DLL工程(选择Fortran Dynamic Link Library),并指定工程名,如下图所示:
按确定按钮,然后新建一个Fortran 源文件,并输入以下的Fortran代码:
(这里仅为说明问题,实际的问题比下面的代码要复杂得多。)
SUBROUTINE OUTPUT(N)
!必须声明本函数为输出函数:DLLEXPORT
!DEC$ ATTRIBUTES DLLEXPORT::OUTPUT
IMPLICIT NONE
INTEGER N
N=N+10
WRITE(*,*) "N=",n
END SUBROUTINE OUTPUT
经编译,链接后,将在Debug目录下生成两个文件,即dll文件与lib文件。
(b) 启动Visual C++6.0,然后新建一个console工程(即Win32 Console Application,当然在MFC中也是完成可行的),如下图所示:
按确定按钮,新建一个空的工程后,再新建一个C++源文件,并输入以下的代码:
#include <iostream.h>
#include <windows.h>
int main()
{
typedef void (_stdcall * wndProc)(int& );
HINSTANCE hLibrary=LoadLibrary("pp.dll"); //加载动态库文件
if(hLibrary==NULL)
{
cout<<"can't find the dll file"<<endl;
return -1;
}
wndProc test=(wndProc)GetProcAddress(hLibrary,"OUTPUT"); //获得Fortran导出函数的地址
if(test==NULL)
{
cout<<"can't find the function file."<<endl;
return -2;
}
int n=3;
test(n); //调用fortran函数
FreeLibrary(hLibrary); //卸载动态库文件
return 0;
}
当编译通过后,将由fortran编译器生成的动态库文件(本算例为pp.dll)复制到C++的Debug文件夹中(即mm\debug),然后点击执行即可,程序的动行结果为:
2. 隐式链接
与上文所介绍的显示链接相比,隐式链接要相对容易一点。下面也通过一个算例进行说明。
(1) 建立一个Fortran动态库文件,其方法与上述完全相同,然后在CVF编译器中输入以下的代码:
SUBROUTINE OUTPUT(N)
!下面为对Fortran函数的声明
!ms$if .not. defined(LINKDIRECT)
!ms$attributes dllexport :: OUTPUT
!ms$endif
IMPLICIT NONE
INTEGER N
N=N+10
WRITE(*,*) "N=",n
END SUBROUTINE OUTPUT
编译通过后,将在Fortran工程文件夹的Debug文件夹中生成两个文件,即*.dll文件与*.lib文件。
(2) 新建一个C++ Console项目,并新建一个C++文件,然后输入以下的代码:
#include <iostream.h>
#include <windows.h>
//声明函数OUTPUT为extern型的,即是从外部调用的。
extern "C" void _stdcall OUTPUT(int& n);
int main()
{
int n=3;
OUTPUT(n);
return 0;
}
然后将Fortran编译器生成的两个文件(dll文件与lib文件)复制到C++的当前目录下,并将lib文件加入到当前的C++工程项目中,如下图所示:
点击菜单“工程-添加工程-Files”,然后选中pp.lib文件即可,如下图所示:
此时即可通过编译,执行此程序,其输出结果如下所示:
总结
本文通过两个简单的实例详细介绍了如何在C++中调用Fortran函数的两种方法,即显式链接与隐式链接。当然实际中我们所遇到的问题将会比本实例要复杂得多,本文仅作为抛砖引玉之用。显式链接与隐式链接两种方法均各自有其优缺点,但由于隐式链接要比显示链接容易得多,也易于理解,实际的大部分算例均采用这种方式。
在linux中操作
对于通常的用户,接促C语言的较多,但是,C语言的开发者,有时还想利用高效的Fortran 数据包 或者是Fortran语言的开发者想借用C提供的强大辅助功能,为此,为了方便大家对二语言相互调用的学习。这里给出了一些简单的入门性技术介绍。
1. C 调用Fortran
编辑Fortran 源文件 add.f95
subroutine add(a, b, c)
implicit none
real:: a, b, c
!
! the add routine
! c = a+b
!
c = a + b
write(*,*) a, ’+’, b, ’=’, c
return
end subroutine add
这里部介绍Fortran的语法格式,add.f95的功能是建立一个add routine 也就是函数。
编译 add.f95
f95 -o add.o -c add.f95
编辑C 主程序main.c
1 #include <stdio.h>
2
3 extern void add_(float *a, float *b, float *c);
4
5 int main(int argc, char *argvs[])
6 {
7 float a, b, c;
8 a = 5.0;
9 b = 7.0;
10 add_(&a, &b, &c);
11
12 return 0;
13 }
3: 声明你要调用的函数,这里注意,是add_,调用的时候用指针,即变量地址,你应该明白了Fortran的调用是改变参数值的。
10:调用fortran函数
编译与执行
#gfortran -o main main.c add.o
#./main
5.000000 + 7.000000 = 12.00000
NOTES; if you use the gcc tool, it will generate error
#gcc -o main main.c add.o
add.o: In function `add_’:
add.f95:(.text+0x4c): undefined reference to `_gfortran_st_write’
add.f95:(.text+0x69): undefined reference to `_gfortran_transfer_real’
add.f95:(.text+0x87): undefined reference to `_gfortran_transfer_character’
add.f95:(.text+0xa4): undefined reference to `_gfortran_transfer_real’
add.f95:(.text+0xc2): undefined reference to `_gfortran_transfer_character’
add.f95:(.text+0xdf): undefined reference to `_gfortran_transfer_real’
add.f95:(.text+0xed): undefined reference to `_gfortran_st_write_done’
collect2: ld returned 1 exit status
Some Wrong with your PATH for you gcc library
OK! you are clever and has known how call fortran routines in the C progaram. Next part I show you how to call C routines in the Fortran program!