[原创]在FORTRAN中调用matlab编译生成的库

taohe

（潜水太久，上来换口气）
拙作《深入浅出matlab7.x混合编程》中讲述了若干matlab和其他编程语言的混合编程问题，其中讲到了如何在c/c++中使用mcc编译生成的库以实现使用matlab计算、画图的功能。其实也可以在fortran中实现。这个帖子用一个简单的例子演示如何在fortran中使用这样的库。作为一个补充材料吧。

试验环境：windows server 2003 std，matlab 6.5.1, msvc 6.0, MinGW gcc 3.4.5.

这里仍然使用拙作《深入浅出matlab7.x混合编程》第七章中的mcctest01.m作为例子，这是一个简单的m文件，其功能是画一个简单的二维图。代码如下：

1. 
   
2. function mcctest01
   
3. % MCCTEST01 Create some data and draw a 2D x-y plot
   
4. %   mcctest01 is used to demonstrates how to convert M-functions into   
   
5. %   library files and how to use created library files with MFC
   
6. %   applications.
   
7. %
   
8. %   Inputs:
   
9. %       None
   
10. %   
    
11. %   Output:
    
12. %       None
    
13. 
    
14. %
    
15. %   Copyright (c) 2005 Dong Weiguo
    
16. 
    
17. % create data
    
18. x = 0:.01:20;
    
19. y = x.*sin(x);
    
20. 
    
21. % draw a plot
    
22. plot(x,y);
    
23. xlabel('x');
    
24. ylabel('y');
    

复制代码

假设你在matlab中已经设置过mbuild使用msvc6，那么可以在matlab中使用如下命令把mcctest01.m编译成一个库函数：

1. 
   
2. mcc -B sgl -t -W libhg:mcctestlib -T link:lib mcctest01 libmwsglm.mlib -h
   

复制代码

结果将产生mcctestlib.lib, mcctestlib.h, mcctestlib.dll等。把这三个文件拷贝到fortran源文件所在的目录。

然后在fortran源文件目录，首先产生一个c/c++ 的wrap文件：

1. 
   
2. // File: mcctest01wrap.cpp
   
3. // Description: A simple testing c++ souce file. A function
   
4. // is defined, which will be called in a FORTRAN 77 code.
   
5. //
   
6. // History: Created on 1 September, 2006
   
7. // Created by: Dong Weiguo
   
8. //////////////////////////////////////////////////////////////
   
9. 
   
10. #include "mcctestlib.h"
    
11. 
    
12. // By default, most of fortran compilers would add an
    
13. // undrescore at the end of function name
    
14. // So let's emulate this bahavior
    
15. extern "C" void mcctest01_()
    
16. {
    
17.      mcctestlibInitialize();
    
18.         mlfMcctest01();
    
19.      mcctestlibTerminate();
    
20. }
    

复制代码

然后创建一个fortran源文件f77mcctst.for：

1. 
   
2. C      A test program to demonstrate how to call c++
   
3. C      functions from a Fortran program.
   
4. C      Synopsis: f77mcctst
   
5. C      History: 1 September, 2006
   
6. C      Created by: Dong Weiguo
   
7. 
   
8.        program f77mcctst
   
9.       
   
10.        call mcctest01()
    
11.   
    
12.        stop
    
13.        end
    
14. 
    

复制代码

最后创建一个Makefile:

1. 
   
2. # makefile for f77mcctst on Windows by using MinGW
   
3. 
   
4. CC = g77
   
5. CFLAGS = -Wall -g -IH:\MATLAB6p5p1\extern\include
   
6. FFLAGS = #-fno-underscoring -fno-second-underscore
   
7. LDFLAGS = -LH:\MATLAB6p5p1\extern\lib\win32\microsoft\msvc60
   
8. 
   
9. MYEXE = f77mcctst.exe
   
10. MYOBJS = f77mcctst.o mcctest01wrap.o
    
11. 
    
12. #default target = all
    
13. all :  $(MYEXE)
    
14. 
    
15. $(MYEXE) : $(MYOBJS)
    
16.         g77 -o $(MYEXE) $(MYOBJS) mcctestlib.lib libmwsglm.lib libmmfile.lib libmatlb.lib libmx.lib sgl.lib libut.lib $(LDFLAGS) -lm -lstdc++
    
17.        
    
18. mcctest01wrap.o : mcctest01wrap.cpp
    
19.         g++ -c $(CFLAGS) mcctest01wrap.cpp
    
20. 
    
21. f77mcctst.o : f77mcctst.for
    
22.         g77 -c $(FFLAGS) f77mcctst.for
    
23. 
    
24. 
    
25. clean:
    
26.         del *.o
    
27.         del $(MYEXE)
    
28. 
    

复制代码

注：为了测试方便，我把matlab下面的msvc6所对应的lib文件拷贝到fortran源文件目录。

最后在fortran源文件目录，在dos命令行用如下命令编译：

1. 
   
2. mingw32-make
   

复制代码

结果生成f77mcctst.exe，运行后就可以看见matlab的图形窗口了。

附：由于时间仓促，少了很多说明文字。比如本来应该先说一下如何在fortran中调用c/c++等。以后有时间就会补上。欢迎大家讨论。

再附：这是偶的第1000贴。本来想写得更好一些，无奈时间不够用。只能这样了。抱歉。

bainhome

幸亏今天睡的晚！
第一千贴，说过要祝贺，早等得心焦了。
不过还没看呢^_^，鼓掌吆喝先！！

不知道taohe老大试过matfor没有，今天大概看了看本贴，从功能的实现上有类似之处。
以下是一个实例，从实现上的不同烦请taohe和scie老大给点儿评价
1.完成一系列库文件配置（我用的是CVF6.6）
2.在FORTRAN中使用控制台程序console application
3.F90文件写如下代码：

1. PROGRAM MAIN
   
2. USE FML
   
3. USE FGL
   
4. IMPLICIT NONE
   
5. TYPE(mfArray)::nx,ny,nz
   
6. TYPE(mfArray)::x,y,z,c,tet
   
7. nx=mfLinspace(-2,2.2d0,21)
   
8. ny=mfLinspace(-2,2.25d0,17)
   
9. nz=mfLinspace(-1.5d0,1.6d0,31)
   
10. call msMeshgrid(mfOut(y,x,z),ny,nx,nz)
    
11. c=2*mfCos(x**2)*mfExp(-(y**2)-(z**2))
    
12. tet=mfGetDelaunay3(x,y,z)
    
13. call msTetSliceXYZ(tet,x,y,z,c,&
    
14.                 mf((/1.0d0,1.0d0/)),mf(0),mf(-0.75d0))
    
15. call msViewPause()
    
16. call msFreeArgs(nx,ny,nz,x,y,z,c)
    
17. END PROGRAM MAIN

复制代码

4.编译后得到如下结果：

两位老大如果有兴趣可以尝试一下，关于速度，易用性等也请给一些评介，谢谢！

[ 本帖最后由 bainhome 于 2006-9-1 22:12 编辑 ]

taohe

完成了第一个1000贴，那就朝2000、10000或更多贴努力前进:)

感谢bainhome兄关于matfor的介绍和分享。我没有用过matfor，不过看起来好像挺不错得样子，如果有人得主要编程语言是fortran的话，倒是不妨考虑选用。我本人对fortran也不熟悉，用得顺手的还是C/C++以及matlab。

penghust

把这个老帖子翻出来，有没有人能把fortran调用c++的步骤给解释一下啊？