10218阅读
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; ?D`h[ai
;`-@L
} <B3$ODGJp
ca!DZ%y
return ret_val; /yO|Q{C}M8
)0p7d:%mV
} /* dot_ */ a) 5;Od
+JRPd.B"@
qB44;!(
yOXL19d@p_
/* Main program alias */ int test_ () { MAIN__ (); } *HXx;:
0X[uXf
!w=6>B^
^h\Y.
g|PRk9
VLs%;|`5D
Iji9N!Yx
M|uWSG
x4nmDEpa
*Z C$DW!-
&WAJ;7f
Mg76v<mv<
二、编译c文件为可执行文件 ,T\)%q
iP^o]4[c
a>XlkkX
M~t;&po
从源程序我们可以看出，如果要这个C能编译通过的话 .0f6b
3pk`&'
:Vl2\H=P
e:kd0)9
A.要让编译器能找到f2c.h qJPEq%'Q
4J6,_8`U
<\C/;
RYem(%jq
B.要让编译器识别fortran的一些常量 g6Nw].{
2P4$^G[
|>L|7>J{<d
-P7JaH/Q
以Dev-C为例 r&?i>.Kz8
>xJh!w<pB
=[T_`*s&
>,s.!vpK
A.将工具包中的f2c.h拷贝到Dev-C安装目录下的include是最好的办法 U|}Bk/0.
AEr8^6
-P'KpX:]hd
dE0p>4F
B.在Dev-C的Project->Project Option->Parameters->linker中用“Add Library or Object”将工具包中的libf77.lib、libi77.lib添加到linker参数中 [ bB
Ua=w;h
XDcA&cM}p
{Wp5Ane
以上两步设置好就正常编译转换过来的C源程序了，以上的程序跑起来会是这个样子 rLzN#Zoi
@#= ail
VQ((c:+!
l!^+Xeg~
D:\Mydoc\DevC\misc>ls *.exe ( 17=|s
iF1E 5{dH
+9d]([Lx
22gk1'~dO
HelloWorld.exe& .. I!OV+utF
=& .KKr
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册