VC++++和Fortran混合編程借助於Fortran生成的DLL進行
(采用C默認的傳址方式進行函數參數傳遞)
1.Fortran 生成DLL
新建Fortran DLL程序test1.f
添加如下代碼:
! test1.f90
!
! FUNCTIONS/SUBROUTINES exported from test1.dll:
! test1 - subroutine
!示例沒有返回值的子例程
subroutine test1(a,b)
! Expose subroutine test1 to users of this DLL
!
!DEC$ ATTRIBUTES C,DLLEXPORT::test1
! Variables
! Body of test1
integer a,b
integer sum
sum=a+b
return
end subroutine test1
!示例有返回值的整數四則運算
!兩數相加
function add(a,b)
implicit none
!DEC$ ATTRIBUTES C,DLLEXPORT::add
integer a,b,add
add=a+b
return
end
!兩數相減
function abstract(a,b)
implicit none
!DEC$ ATTRIBUTES C,DLLEXPORT::abstract
integer a,b,abstract
abstract=a-b
return
end
!兩數相乘
function multiply(a,b)
implicit none
!DEC$ ATTRIBUTES C,DLLEXPORT::multiply
integer a,b,multiply
multiply=a*b
return
end
!兩數相除 (需要添加考慮被除數是否為0以及能否整除的判斷)
function divided(a,b)
implicit none
!DEC$ ATTRIBUTES C,DLLEXPORT::divided
integer a,b,divided
divided=a/b
return
end
編譯後生成test1.dll,test1.obj等文件。其中這兩個文件是我們在VC中調用所需要的。
查看test1.dll中生成的函數如下圖所示。
注意:使用偽注釋語句!DEC$ ATTRIBUTES C,DLLEXPORT::functionName後,生成的函數名與Fortran中定義的函數名一致,並沒有按照Fortran編譯器的默認屬性將 函數名都轉變為大些,如下圖所示。在後續的VC中調用的時候需要保持調用的函數名稱一致,否則會出現找不到函數的錯誤提示。
2.VC控制台調用Fortra生成的DLL
新建VC控制台應用程序,新建一個checktest1的工程。
注意:需要在工程的project菜單下的add to project子菜單的file對話框中添加上一步生成的test1.dll,test1.obj兩個文件,否則編譯能通過,鏈接的時候失敗。還需將以上 兩個文件拷貝到checktest1工程的debug目錄下,否則運行的時候出現找不到文件的錯誤提示。自己測試了一下,以上兩步是必須的。
添加如下代碼:(注意紅色的部分)
#include “stdafx.h”
#include “iostream.h”
//extern “C”{_stdcall TEST1(int* x,int* y);}
//extern “C”{_stdcall ADD(int* x,int* y);}
//extern “C”{_stdcall ABSTRACT(int* x,int* y);}
//extern “C”{_stdcall MULTIPLY(int* x,int* y);}
//extern “C”{_stdcall DIVIDED(int* x,int* y);}
//注意此處函數名稱要與DLL生成時保持一致(如下中的藍色部分),否則會出現找不到函數的錯誤提示。並且一定要記得去掉參數中的指針符號*。
extern “C”{_cdecl test1(int x,int y);}
extern “C”{_cdecl add(int x,int y);}
//采用C的傳值方式,則需要將_stdcall修改為_cdecl
//相應的Fortran DLL處要添加C的調用方式,即將!DEC$ ATTRIBUTES DLLEXPORT::add修改為:!DEC$ ATTRIBUTES C,DLLEXPORT::add
//適應偽注釋!DEC$ ATTRIBUTES C,DLLEXPORT::add後生成的DLL函數中只存在函數名為add的函數,ADD 和_ADD@8 均不存在,參見上圖中的DLL函數名稱
extern “C”{_cdecl abstract(int x,int y);}
extern “C”{_cdecl multiply(int x,int y);}
extern “C”{_cdecl divided(int x,int y);}
int main(int argc, char* argv[])
{
int a=35,b=5;
int sum=0;
int abs=0;
int mul=0;
int div=0;
//TEST1(&a,&b);
//sum=ADD(&a,&b);
//abs=ABSTRACT(&a,&b);
//mul=MULTIPLY(&a,&b);
//div=DIVIDED(&a,&b);
test1(a,b);
sum=add(a,b);
abs=abstract(a,b);
mul=multiply(a,b);
div=divided(a,b);
printf(“a+b= %dn”,sum);
printf(“a-b= %dn”,abs);
printf(“a*b= %dn”,mul);
printf(“a/b= %dn”,div);
printf(“Hello World!n”);
return 0;
}
然後編譯運行可以得出正確的結果: