10216阅读
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; 9u?Eb~#$
?0uOR*y'
} wePI*."]
+\srZ<67
return ret_val; X0!Bs-WFp
G"{4'LlA
} /* dot_ */ e) ]RA?bF
D/cg7
C_-%*]*,j
aj?2jU~Pq
/* Main program alias */ int test_ () { MAIN__ (); } 7MoR9,(
5 >'66gZ
gU+ss
.(JE-upJ"
ygMd$0:MN
b]"2VN
ma1(EJ/
%n4@[fG%K
p%tE v
$.``OxJk%
b'&LBT7
\zBd<H4S:
二、编译c文件为可执行文件 ^u3*hl}YKy
J0Jr BXCh
c5<kbe
xUQdVrFU
从源程序我们可以看出，如果要这个C能编译通过的话 /9P^{OZ;y
)sRN!~
^)Smv\Md
7,f:Qi@g
A.要让编译器能找到f2c.h !;TR2Zcn
`(P "u
)ni"qv~J
:-(U%`a[
B.要让编译器识别fortran的一些常量 E1#H{)G
33*^($bE&
#]k0Z~Bl
dJloH)uJZ>
以Dev-C为例 GI}4,!^N
+I*k0"gj6
GahaZ F
"pOqd8>]
A.将工具包中的f2c.h拷贝到Dev-C安装目录下的include是最好的办法 3T"2S[gT
J0&zb'1
3(MoXA*
@8QFP3\1
B.在Dev-C的Project->Project Option->Parameters->linker中用“Add Library or Object”将工具包中的libf77.lib、libi77.lib添加到linker参数中 j[_t6Z
dUI3erO
pi|\0lH6W
>9 q]>fJ
以上两步设置好就正常编译转换过来的C源程序了，以上的程序跑起来会是这个样子 p 7 ,f6kG
lB=(8.
C\y[&egww
ThjUiuWe
D:\Mydoc\DevC\misc>ls *.exe ,m2A p\l
o^8*aH)I>Y
Jw2B&)k/
k&WUv0
HelloWorld.exe& .. 5P-K *C&
c]GQU
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册