10217阅读
8回复

[资料共享]将FORTRAN源程序转换为C

离线febi

UID ：2537

注册：2007-05-11
登录：2016-05-05
发帖：123
等级：退休版主

0楼发表于: 2007-08-03 13:55:56

将FORTRAN源程序转换为C

来自：http://www.51fy.cn/program/other/200705/33522.htm

来源/ChinaUnix 文/blau 更新时间：2007-5-31 5:50:44

现在的大部分人听到fortran都会冷笑一番，特别是计算机系的同学，更是嗤之以鼻,多古老的语言啊，现代语言该有的特性它一点都没有,早该进回收站了！如果这些同学写过要用到很多浮点数、复数域的计算估计就不会有这样的感慨。大学毕业论文用fortran做过水力计算，所以至今仍认为在科学计算领域fortran是最合适的。

今天逛母校论坛发现一个学妹要把fotran的程序改写为C的，搞不清现在的老师为啥要学生把fortran转到c，可能顺应潮流吧，呵呵。一时兴起上google看了一下，老外早就写好现成的转换程序，还是开源的！顺手也玩了一把，:-)

一、从Fortran源程序转C源程序

先来一段很普通的Fortran程序

C FORTRAN USES POS. 1 TO MAX. 72: POS. 1 FOR COMMENT SIGN,

C POS 1 to 5 FOR LABELS, POS. 6 IS TO BE LEFT FREE EXCEPT OF LINE

C CONTINUATION

C2345*789#123456789#123456789#123456789#123456789#123456789#123456789#12

PROGRAM TEST

INTEGER*4 I

REAL*4 Y(10),Z(10)

DO 100 I=1,10

Z(I)=.75*I

Y(I)=.37*I

100 CONTINUE

DO 110 I=1,10

PRINT*,DOT(I,Y,Z)

110 CONTINUE

END

C THIS FUNCTION IS AN EXAMPLE FROM FORTRAN.DOC

FUNCTION DOT(N,X,Y)

INTEGER N

REAL X(N),Y(N)

DOT = 0

DO 10 I = 1, N

10 DOT = DOT + X(I)*Y(I)

END

Fortran源码是对位置敏感的：

A. 1~5列只能用来做注释（c打头）和标签

B. 第6列最好是空着不用

C. 一行不能超过72个字符

之所以有这样古怪的定义是因为，很久以前是没有键盘，数据只能靠纸带输入 ^_^

如果不符以上要求，转换时就会报

Error test.f 5 illegal continuation card (starts " P")

这样的错误

如果源程序格式无误，用工具包中的命令行工具

D:\f2cexe>f2c test.f

test.f:

MAIN test:

dot:

D:\f2cexe>

就可以转成test.c了

转成的c程序会是这样子的

/* test.f -- translated by f2c (version of 23 April 1993 18:34:30).

You must link the resulting object file with the libraries:

-lf2c -lm (in that order)

*/

#include "f2c.h"

/* Table of constant values */

static integer c__4 = 4;

static integer c__1 = 1;

/* FORTRAN USES POS. 1 TO MAX. 72: POS. 1 FOR COMMENT SIGN, */

/* POS 1 to 5 FOR LABELS, POS. 6 IS TO BE LEFT FREE EXCEPT OF LINE */

/* CONTINUATION */

/* 2345*789#123456789#123456789#123456789#123456789#123456789#123456789#12 */

/* Main program */ MAIN__()

{

/* System generated locals */

real r__1;

/* Builtin functions */

integer s_wsle(), do_lio(), e_wsle();

/* Local variables */

static integer i;

static real y[10], z[10];

extern doublereal dot_();

/* Fortran I/O blocks */

static cilist io___4 = { 0, 6, 0, 0, 0 };

for (i = 1; i <= 10; ++i) {

z[i - 1] = i * (float).75;

y[i - 1] = i * (float).37;

/* L100: */

}

for (i = 1; i <= 10; ++i) {

s_wsle(&io___4);

r__1 = dot_(&i, y, z);

do_lio(&c__4, &c__1, (char *)&r__1, (ftnlen)sizeof(real));

e_wsle();

/* L110: */

}

} /* MAIN__ */

/* THIS FUNCTION IS AN EXAMPLE FROM FORTRAN.DOC */

doublereal dot_(n, x, y)

integer *n;

real *x, *y;

{

/* System generated locals */

integer i__1;

real ret_val;

/* Local variables */

static integer i;

/* Parameter adjustments */

--y;

--x;

/* Function Body */

ret_val = (float)0.;

i__1 = *n;

for (i = 1; i <= i__1; ++i) {

/* L10: */

ret_val += x * y; yfR0vp<&
$#W^JWN1
} Uspv^O9_
i.?rom
return ret_val; '"rm66
) .' + {
} /* dot_ */ /RxqFpu|.
b E40^e
A=e1uBGA
\_(0V"
/* Main program alias */ int test_ () { MAIN__ (); } P6ka'!z
6cbV[!BL
eWD!/yr|
xy$aFPH!-
v4uQ0~k~X
|UQ[pas
KZ^>_K&
:P+7ti@
f4NN?"W)
CtD<%v3`
V$Oj@vI
b'R]DS{8
二、编译c文件为可执行文件 yT7{,Z7t
oY,{9H37b
M\7F1\ X
Dvl\o;
从源程序我们可以看出，如果要这个C能编译通过的话 ]iFW>N*a
&G\C[L
;4/ n~
Ah5o>ZtcO
A.要让编译器能找到f2c.h HJpx,NU'
EQZ/v gho
%]tW2s"
REi"Aj=
B.要让编译器识别fortran的一些常量 a[sdYZ
-O/[c
$M4_"!
xUa9>=JU{
以Dev-C为例 !O'p{dj][
hjk]?MC
,+gtr.
l^}5PHLd
A.将工具包中的f2c.h拷贝到Dev-C安装目录下的include是最好的办法 xt`a":lru
EqIs&){
{BaPK&x,
FJ?]|S.?,
B.在Dev-C的Project->Project Option->Parameters->linker中用“Add Library or Object”将工具包中的libf77.lib、libi77.lib添加到linker参数中 ^(7l!
5bI4' ;
HxLuJ
\Ov~ t
以上两步设置好就正常编译转换过来的C源程序了，以上的程序跑起来会是这个样子 ~G6xk/+n-m
,lFzL3'_0x
Txpj#JD
cSG(kFQ
D:\Mydoc\DevC\misc>ls *.exe gf4Hq&Rf
) R\";{`M
^%g8OP
j/E(*Hv
HelloWorld.exe& .. $`t2SD
pjSM7PhQ
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册