Fortran学习笔记6(函数、子程序)


程序代码中,常常会在不同的地方重复用到某一功能和重复某一代码,这个时候就要使用函数。函数包括内嵌函数,自定义函数,子程序


子程序Subroutine

子程序之所以叫子程序,是因为在没有经过调用之前,它不会被执行。(当然,如果把编译器对主函数的调用也看做同样的话则另当别论。)

包含子程序的程序的格式:

program main 
......
......
call sub_1(...)
......
call sub_2(...)
......
end program main

subroutine sub_1(...)
...
end subroutine sub_1

subroutine sub_2(...)
...
end subroutine sub_2

主程序并不一定要放在最开始,它可以放在程序的任何地方(但是最好放在最开始或者结尾,如果你放在了一大段程序的中间......何苦为难自己!)子程序之间也可以互相调用。甚至可以自己调用自己,这个算法就叫递归,但是Fortran77不支持。调用语句为call

有关子程序,一个很重要的概念是,子程序独立的拥有属于自己的变量申明。主程序和子程序,不同的子程序之间的变量名可以相同,他们互不干扰。除了变量独立,子程序还独立拥有自己的行代码

在子程序的调用中会涉及到参数传递。Fortran中的参数传递使用的是传址调用。也就是说,实参和形参公用一块内存,两者改变其中之一,另一者也会发生变化,(其实是形参暂时的借住在实参家里。)


自定义函数Function

自定义函数的运行和上面的子程序没有什么两样,也需要调用才能被执行,也可以独立声明变量,参数传递的方法和子程序相同。尽管如此,自定义函数和子程序还是有两点不同:

  1. 调用自定义函数前要声明。
  2. 自定义函数调用后会返回一个数值

下面是一个简单的程序:用牛顿法来解方程。

program main
    implicit real*8(a-z)
    external f,g
    integer::max
    integer::i,k=1
    real*8::x0,x1,x2,dx,tol,f,g
    open(unit=11,file='fin.txt')
    open(unit=12,file='fout.txt')
    read(11,*)x0,max,tol
    
    do 
        x1=x0-f(x0)/g(x0)
        dx=abs(x1-x0)
        if(dx<=tol) exit
                k=k+1
        if(k>=max) exit
                x0=x1
    end do
      
    write(12,*)x0
    end program main

function f(x)
      implicit none
      real*8::f,x
      f=x**3+2*x**2+10*x-20
end function

function g(x)
      implicit none
      real*8::g,x !,dx
      g=3*x**2+4*x+10
end function

下面是一个更简单的例子:(摘自彭国伦《Fortran95程序设计》一书)

program main
    implicit  none
    real::a=1,b
    real add
    add(a,b)=a+b             
! 简单的子程序直接写在主程序里面,并且不用call调用。并且这样的函数只能在本主函数中使用,其他地方不能调用。
    write(*,*)add(a,3.0) 
    stop
end program

写函数的时候请注意,尽量不要去改变传入的参数的数值,最好另开辟一个内存来存储要用的应变量,而不要随意去改变传入的自变量。


全局变量COMMON

不同的程序之间,也就是在不同的函数之间,除了可以通过传递参数的方法来共享内存,还可以通过全局变量来让不同程序中声明出来的变量使用相同的内存位置。

program main
    implicit none
    integer::a,b
    common a,b
    a=1
    b=2
    call sub()
end program

subroutine sub()
    implicit none
    integer::num1,num2
    common num1,num2
    write(*,*)num1,num2
    return
end subroutine sub

由于全局变量使用的“地址对应”的方法在程序中共享数据,所以在程序设计时常常会出现一些不方便,比如在主程序中声明了6个全域变量,而你在子程序中只需要使用第六个全域变量,但是你不得不把前五个都写出来。这样很麻烦,一个幼小的解决办法是将全域变量分区。

program main
    implicit none
    integer::a,b,c,d
    common/group1/ a,b
    common/group2/ c,d
    a=1
    b=2
    c=3
    d=4
    call sub_1()
    call sub_2()
end program

subroutine sub_1()
    implicit none
    integer::num1,num2
    common/group1/ num1,num2
    write(*,*)num1,num2
    return
end subroutine sub_1

subroutine sub_2()
    implicit none
    integer::num1,num2
    common/group2/ num1,num2
    write(*,*)num1,num2
    return
end subroutine sub_2

BLOCK DATA

关于COMMON设置初值,不能直接在子程序或者主程序中使用data来设置初值,要在block data程序模块中使用data命令来设置初值。BLOCK data很像一段子程序,它也是一段独立的程序模块,也拥有自己的变量声明,不过它不需要别人调用就可以自己执行。事实上,这段程序会在主程序执行前生效,不过它的功能只在于设置全局变量的初值,不能有其他执行命令出现。具体语法如下:

 block data name    !   name可以省略
    implicit none       !   最好不要省略
    integer...               !  声明变量
    real...
    common ...          ! 把变量放在公共空间
    common/group1/...         !   公共空间分区  
    data var1,var2,...,varn/value1,value2,...,valuen/     !   data语句赋初值
    ...
    ...
    end block data name              ! 可以只写end或end block data

还有一点要注意,全局变量不能声明成常量。所以block data中不能出现parameter。

实际的例子如下:

 program main
    implicit none
    integer::a,b
    common a,b
    integer::c,d
    common/group1/c,d
    integer::e,f
    common/group2/e,f

    write(*,*)a,b,c,d,e,f
    stop
end  program main

 block data name    !   name可以省略
    implicit none       !   最好不要省略
    integer::a,b
    common a,b
    data a,b/1,2/
    integer::c,d
    common/group1/ c,d
    data c,d/3,4/
    integer::e,f
    common/group2/e,f
    data e,f/5,6/
    end block data name              ! 可以只写end或end block data

小结

COMMON语法在数据共享式会出现很多限制和不足,建议不要使用而采用module语句来实现数据共享。


最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 202,980评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,178评论 2 380
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 149,868评论 0 336
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,498评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,492评论 5 364
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,521评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,910评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,569评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,793评论 1 296
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,559评论 2 319
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,639评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,342评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,931评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,904评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,144评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,833评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,350评论 2 342

推荐阅读更多精彩内容