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;
}
然後編譯運行可以得出正確的結果: