dvbbs
收藏本页
联系我们
论坛帮助
dvbbs

>> Fortran语言使用经验交流
搜一搜相关精彩主题 
Fortran中文网Fortran中文网—Fortran语言经验交流Fortran语言使用经验交流 → Fortran常见问题汇总。(转自编程爱好者Fortran论坛的FAQ)

您是本帖的第 31084 个阅读者
树形 打印
标题:
Fortran常见问题汇总。(转自编程爱好者Fortran论坛的FAQ)
Evan
帅哥哟,离线,有人找我吗?
等级:论坛游民
文章:82
积分:755
门派:无门无派
注册:2006年4月4日
楼主
 用支付宝给Evan付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给Evan

发贴心情
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编辑过]
ip地址已设置保密
2006/9/12 20:37:50
euler
帅哥哟,离线,有人找我吗?
等级:新手上路
文章:1
积分:201
门派:无门无派
注册:2006年10月16日
2
 用支付宝给euler付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给euler

发贴心情
ip地址已设置保密
2006/10/16 21:27:29
dentain
帅哥哟,离线,有人找我吗?
等级:新手上路
文章:2
积分:233
门派:无门无派
注册:2006年10月17日
3
 用支付宝给dentain付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给dentain

发贴心情

菜鸟问题:

    我的程序运行中出现数组下标超出范围--array bounds exceeded,请问怎么处理?

ip地址已设置保密
2006/10/18 15:17:40
ruoti
美女呀,离线,留言给我吧!
等级:论坛游民
文章:21
积分:519
门派:无门无派
注册:2006年12月30日
4
 用支付宝给ruoti付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给ruoti

发贴心情
运行debug,按照箭头指示的地方你看看数组的定义对不对
ip地址已设置保密
2007/1/4 15:01:28
liucunxi
帅哥哟,离线,有人找我吗?
等级:新手上路
文章:1
积分:211
门派:无门无派
注册:2007年1月7日
5
 用支付宝给liucunxi付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给liucunxi

发贴心情

谢谢啊]

ip地址已设置保密
2007/1/7 10:11:28
qqz2001
美女呀,离线,留言给我吧!
等级:新手上路
文章:9
积分:280
门派:无门无派
注册:2006年12月20日
6
 用支付宝给qqz2001付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给qqz2001

发贴心情

这个帖子见过,呵呵。。。

不过的确是个好帖子

ip地址已设置保密
2007/4/24 19:36:14
kittyxue80
美女呀,离线,留言给我吧!
等级:论坛游民
文章:16
积分:442
门派:无门无派
注册:2007年4月18日
7
 用支付宝给kittyxue80付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给kittyxue80

发贴心情
我是FORTRAN新手,问题很多啊,
ip地址已设置保密
2007/5/9 10:35:00
ninegreatsun
帅哥哟,离线,有人找我吗?
等级:论坛游民
文章:13
积分:197
门派:无门无派
注册:2007年5月13日
8
 用支付宝给ninegreatsun付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给ninegreatsun

发贴心情
回答得太好了。!!!!!
ip地址已设置保密
2007/5/15 7:53:27
ywhwjg
帅哥哟,离线,有人找我吗?
等级:新手上路
文章:1
积分:201
门派:无门无派
注册:2007年7月20日
9
 用支付宝给ywhwjg付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给ywhwjg

发贴心情

请问goto(1001,1002,1003,1004,1005,1006,1007),igeom

是什么意思,我怎么看不懂,一般goto直接跟转向语句标号,这个却语句标号放到括号里,后面还跟着igeom

请在线的各位高手解释以下

ip地址已设置保密
2007/7/20 16:22:56
mickmade
帅哥哟,离线,有人找我吗?
等级:新手上路
文章:1
积分:201
门派:无门无派
注册:2007年8月6日
10
 用支付宝给mickmade付款或购买其商品,支付宝交易免手续费、安全、快捷! 点击这里发送电子邮件给mickmade

发贴心情
这个是goto的带条件的用法,类似于switch啦,后面那个变量igeom=1时,跳到1001行,等于2时,跳到1002,跳到3时,1003,如此类推。
ip地址已设置保密
2007/8/6 9:04:43

 28   10   1/3页      1   2   3   尾页 
网上贸易 创造奇迹! 阿里巴巴 Alibaba
Powered By Dvbbs Version 7.1.0 Sp1
Copyright ©2005 - 2008 www.fortran.cn
页面执行时间 0.20313 秒, 5 次数据查询
京ICP备05056801号