以文本方式查看主题 - Fortran中文网 (http://fortran.cn/bbs/index.asp) -- Fortran语言使用经验交流 (http://fortran.cn/bbs/list.asp?boardid=2) ---- Fortran常见问题汇总。(转自编程爱好者Fortran论坛的FAQ) (http://fortran.cn/bbs/dispbbs.asp?boardid=2&id=376) |
-- 作者:Evan -- 发布时间:2006/9/12 20:37:50 -- Fortran常见问题汇总。(转自编程爱好者Fortran论坛的FAQ) 以下操作,如无特别说明,都是以CVF6.6为例。 1. 如何加大Stack size? 选Project => Settings => Link => Category: Output => Stack allocations Reserve: 这里填新值(默认为1M,若需要10M,则填10000000) 2. 如何用Fortran批量生成文件? 设要生成4000个文件,文件名为AA1-AA4000,如何写循环生成文件,而不用写4000次write命令呢? 用内部文件: character(len=80) :: filename,form integer :: i do i=1,4000 select case (i) case (1:9) write(form,\'(i1)\') i case (10:99) write(form,\'(i2)\') i case (100:999) write(form,\'(i3)\') i case (1000:9999) write(form,\'(i4)\') i end select write(filename,*) "AA",trim(form),".TXT" open(10,file=filename) write(10,*) i close(10) end do stop end 3. 如何用Fortran动态生成输出格式? 设有一个数组data(100),输出时,希望每行输出num个数,而num由用户输入,如何实现? 用内部文件: character(len=80) :: form real :: data(100) integer :: i,num data = (/ (i,i=1,100) /)/10.0 read(*,*) num write(form,*) "(",num,"f10.3)" write(*,form) data stop end 4. MS POWERSTATION4.0是不是很垃圾? POWERSTATION4.0是垃圾,其中Bug太多,多到不可用的地步! 在这个主题里,换了CVF后问题就没了的人已有相当的数目。 如果你用POWERSTATION4.0,遇到莫名其妙的错误,建议换CVF6.6,这是一个比较成熟的编译器。 5. 如何用F90/95生成随机数? 注意: 现在计算机产生的随机数都是伪随机数。 random_number(x) 产生一个0到1之间的随机数(x可以是向量),但是每次总是那几个数。 用了random_seed ()后,系统根据日期和时间随机地提供种子,使得随机数更随机了。 program random implicit none real :: x call random_seed () ! 系统根据日期和时间随机地提供种子 call random_number (x) ! 每次的随机数就都不一样了 write(*,*) x stop end program random 6. 函数/子程序超载的例子 设要编一个两个变量值互换的子程序swap(a,b),哑元a,b可能是实型数,整型数,数组,矩阵,字符串,派生类型等等。但是希望只用一个子程序接口swap(a,b)来实现。F90可以用类属接口来实现这种子程序超载: module Utilities implicit none private I_Swap,R_Swap,RVec_Swap,RMat_Swap,Type_Swap public :: Swap interface Swap module procedure I_Swap,R_Swap,RVec_Swap,RMat_Swap,Type_Swap end interface contains subroutine i_swap (a,b) ! 整型置换 integer (ikind),intent(in out) :: a,b integer (ikind) :: t 。。。 ! 略 end subroutine i_swap subroutine r_swap (a,b) ! 实型置换 real (rkind), intent(in out) :: a,b real (rkind) :: t t = a a = b b = t return end subroutine r_swap subroutine RVec_swap (a,b) ! 实型向量置换 real (rkind), intent(in out) :: a(:),b(:) integer (ikind) :: i do i=1, size(a) call R_Swap (a(i),b(i)) end do return end subroutine RVec_swap subroutine RMat_swap (a,b) ! 实型矩阵置换 。。。 ! 略 end subroutine RMat_swap subroutine Type_swap (a,b) ! 派生类型置换 。。。 ! 略 end subroutine Type_swap end module Utilities 7. 如何在CVF中为多行代码加注释? 见:http://www.programfan.com/club/showbbs.asp?id=124318 8. 推荐好的代码风格 根据F90子集语言ELF90和F的要求整理(部分)。 “强迫用”的语言特性: + F90的自由格式的源代码。 + implicit none。 + 子过程的哑元都要有intent属性。 + 函数子程序的哑元必须指定为intent(in)。 + 所有子程序和函数都放在模块(module)中,然后引用(use)该模块;或者放在program中。 + 数组哑元要求是假定形状的,或者有固定的维数和大小。字符哑元要求是假定长度的。 + 对于recursive function(递归函数)语句,必须有result子句。 + 在所有派生类型(type)的定义语句中,必须用双冒号分隔符(::)。 + 主程序要求有program语句。 + 在程序单元的end语句中要求后跟程序单元的类型和名称。 + 在end type语句中要求后跟类型的名称。 + end program前必须有stop语句以表示停止执行。 + 子过程中必须有return语句,以表示返回。 + subroutine s( )并且call s( ),即必须有括号。 “不得用”的语言特性: - allocatable、intent、pointer、save、dimension、parameter和target语句形式。(用属性形式代替。) - external语句形式。(用显式的接口代替。) - assign、赋值go to、交错return、continue、entry、和计算go to 语句。 - include文件。(用模块代替。) - data和block data。(在类型声明语句中进行初始化或赋值。) - common块。(将全局数据放在模块中,用模块代替。) - equivalence。(被认为是许多难以查找的编程错误的来源。) - double precision语句。(用real语句声明双精度的实型数。) - 语句函数。(用内部函数代替。) - 专用固有函数。(用类属函数代替。) - 假定大小数组。(用假定形状数组代替。) - do n (其中n为语句标号)。(用do和end do代替。) - 非整数do变量和表达式。 - 同一行上多条语句。 - 逻辑型case表达式。 - 从if块外面分支到end if。 - where语句形式。(用where结构形式。) - 在open和inquire语句中的blank= 说明符。 - 双字关键词之间要求有空格:in out,go to。不能写为inout,goto。 9. 将字符串改为大写的子程序 subroutine UpCase (str) !========================================= ! change to upper case !========================================= character(len=*),intent(in out) :: str integer(4) :: icha,LL,icval integer(4),parameter :: diff = ichar(\'a\') - ichar(\'A\') LL = len_trim(str) do icha=1,LL icval = ichar(str(icha:icha)) if (icval>=ichar(\'a\') .and. icval<=ichar(\'z\')) then str(icha:icha) = char(icval-diff) end if end do return end subroutine UpCase 10. CVF中源代码信息浏览 默认情况下,代码信息浏览是失效的(Go to Definition/Reference都不起作用),你可以将其启用: 1. Project/Settings/Fortran, 选中Generate Source Browse Information。 2. 选BrowseInfo卡页,选中Build Browse info file,点击OK。 3. Build程序。代码浏览器仅当成功Build程序后才可应用。 4. Tools/Source Browser。 5. 在Browse对话框底部,找到Case sensitive复选框。Fortran是大小写不敏感的语言,去除Case sensitive的选中。点击OK。 6. 现在可以用Browse对话框来浏览函数调用关系以及变量声明和引用关系。 7. 光标置于一个变量名或函数名中,鼠标右键:Go to Definition使光标跳到该变量声明或函数定义语句;Go to Reference使光标跳到该变量的引用或函数引用语句。 注:(Addison提供) IVF中现在不适用,因为此技术是属于MS的,不是INTEL的,所以,intel正在与MS就有关问题进行磋商,目前的进展比较顺利,可能要在9.1以后的版本中支持,目前还没有时间表。 [此贴子已经被作者于2006-9-12 20:41:24编辑过]
|
-- 作者:euler -- 发布时间:2006/10/16 21:27:29 -- |
-- 作者:dentain -- 发布时间:2006/10/18 15:17:40 -- 菜鸟问题: 我的程序运行中出现数组下标超出范围--array bounds exceeded,请问怎么处理? |
-- 作者:ruoti -- 发布时间:2007/1/4 15:01:28 -- 运行debug,按照箭头指示的地方你看看数组的定义对不对 |
-- 作者:liucunxi -- 发布时间:2007/1/7 10:11:28 -- 谢谢啊] |
-- 作者:qqz2001 -- 发布时间:2007/4/24 19:36:14 -- 这个帖子见过,呵呵。。。 不过的确是个好帖子 |
-- 作者:kittyxue80 -- 发布时间:2007/5/9 10:35:00 -- 我是FORTRAN新手,问题很多啊, |
-- 作者:ninegreatsun -- 发布时间:2007/5/15 7:53:27 -- 回答得太好了。!!!!! |
-- 作者:ywhwjg -- 发布时间:2007/7/20 16:22:56 -- 请问goto(1001,1002,1003,1004,1005,1006,1007),igeom 是什么意思,我怎么看不懂,一般goto直接跟转向语句标号,这个却语句标号放到括号里,后面还跟着igeom 请在线的各位高手解释以下 |
-- 作者:mickmade -- 发布时间:2007/8/6 9:04:43 -- 这个是goto的带条件的用法,类似于switch啦,后面那个变量igeom=1时,跳到1001行,等于2时,跳到1002,跳到3时,1003,如此类推。 |