fortran问题

imsl7.0中用use linear_operators这句话是会出错的，当时intel论坛上也有人问，后来intel 给出了个X64 imsl的补丁，但是32位的没有。

要使用包含在linear_operators这个库中的函数时，要用use+原函数

例如：上面那个例子把use linear_operators改成use operation_xt即可

另外，imsl7.0引用函数和以前版本不一样的

补充一点：

linear_operators这个文件是这样子的（一看就知道原因了）：

modulelinear_operators

usecond_int

usedet_int

usediag_int

usediagonals_int

usefft_int

useifft_int

useeye_int

uselin_eig_self_int

uselin_sol_self_int

usenorm_int

useoperation_i

useoperation_ix

useoperation_t

useoperation_h

useoperation_tx

useoperation_hx

useoperation_x

useoperation_xi

useoperation_xt

useoperation_xh

useorth_int

userand_int

userank_int

usesvd_int

useunit_int

useeig_int

usechol_int

useisnan_int

end module

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

doi=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

callrandom_seed () ! 系统根据日期和时间随机地提供种子

callrandom_number (x) ! 每次的随机数就都不一样了

write(*,*) x

stop

end program random

6. 函数/子程序超载的例子

设要编一个两个变量值互换的子程序swap(a,b)，哑元a,b可能是实型数，整型数，数组，矩阵，字符串，派生类型等等。但是希望只用一个子程序接口swap(a,b)来实现。F90可以

用类属接口来实现这种子程序超载：

module Utilities

implicit none

privateI_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

subroutinei_swap (a,b) ! 整型置换

integer (ikind),intent(in out) :: a,b

integer (ikind) :: t

。。。! 略

end subroutine i_swap

subroutiner_swap (a,b) ! 实型置换

real (rkind), intent(in out) :: a,b

real (rkind) :: t

t = a

a = b

b = t

return

end subroutine r_swap

subroutineRVec_swap (a,b) ! 实型向量置换

real (rkind), intent(in out) :: a(:),b(:)