fortran代码（fortran代码转python）

今天给各位分享fortran代码的知识，其中也会对fortran代码转python进行解释，如果能碰巧解决你现在面临的问题，别忘了关注本站，现在开始吧！

本文目录一览：

• 1、fortran代码
• 2、如何运行fortran源代码
• 3、一个公式写成Fortran语言代码
• 4、fortran代码优化从那里入手？

fortran代码

double precision a,b,h,eps,x1,y1,x2,y2,x

f(x)=x**3-x**2-1

data a,b,h,eps/0,3,0.01,0.00000001/

x1=a

y1=f(x1)

20 x2=x1+h

y2=f(x2)

if (x1.lt.b) then

if (y1*y2.le.0) then

a1=x1

b1=x2

10 x=(a1+b1)/2.0

if (y1*f(x).gt.0)then

a1=x

else

b1=x

endif

if (b1-a1.gt.eps) goto 10

write(*,*)x

endif

x1=x2

y1=y2

goto 20

endif

end

如何运行fortran源代码

你已的编译器已经能够完全正常运行，而且程序也可以运行了呢，只是你不知道如何运行它而已！

在unix或linux下编译fortran，一般系统没有自带fortran的编译器，可以自行安装g77或gfortran编译器，它们是gcc一个系列的，也可以下载单独的g95编译器。

假定你的fortran源程序为hello.for

则输入如下的命令来编译它（g77、gfortran、g95）

g77

hello.for

这样就会在当前目录生成a.out的可执行文件，这样既可执行

./a.out

如果需要指定生成的文件名，可以加个参数-o

来制定文件名，像这样

g77

-o

hello

hello.for

./hello

明白了吗？unix类操作系统的执行档，并不需要.exe这样的扩展名的。

一个公式写成Fortran语言代码

program baidu

integer::I,J,N

real*8::Cr,T

real*8,dimension(:),allocatable ::P,X

N=3 !变量X的个数

Cr=5.0d0 !常量Cr，自己设定

T=4.0d0 !常量T，自己设定

allocate(P(N),X(N))

! =======读入变量X的值

do I=1,N

write(*,*)"请输入第 ",I," 个变量的值："

read(*,*)X(I)

enddo

! =======读入变量X的值

do I=1,N

P(I)=(-4.2d0/Cr**2*X(I)+2.9/Cr)*T

write(*,*)“第 ”,I," 个变量X对应结果：",P(I)

enddo

end

fortran代码优化从那里入手？

从数组开始和循环开始，因为程序的运行性能主要看怎么处理迭代和对大量数据的存储。

如：

do i=1,L

do j=1,N

do k=1,M

write(*,*) C(i,j,k)

end do

end do

end do

do i=1,L

do j=1,N

do k=1,M

write(*,*) C(k,j,i)

end do

end do

end do

显然第二个要比第一个优化。

因为fortran的数组存储顺序为：

a(1,1) a(2,1) a(3,1) a(1,2) a(2,2) a(3,2)

这个一个连续的内存空间，我们在读取它的时候连续读取肯定要比跳跃性读取好！

fortran代码的介绍就聊到这里吧，感谢你花时间阅读本站内容，更多关于fortran代码转python、fortran代码的信息别忘了在本站进行查找喔。