关于fortran的程序问题，看不懂了，求各位大神指导

DDOT=DTEMP

// DOUBLE PRECISION相当于 C语言的 double

相当于 C语言的

double result=DTEMP；

。。。。。。。

return result；//函数返回时才执行，平常DDOT这样的函数名，代表函数返回值而已。

并不表示函数返回。

只有递归调用，和函数调用，和用函数作为参数时，函数声明和定义时才是一个函数！

其他时候在函数内部，就直接相当于函数返旦含回值！！！

在调用函数的表达式中整个函数调用才相当于函数返回值！！！

有点绕！！！

DOUBLE PRECISION FUNCTION DDOT(N,DX,INCX,DY,INCY)

INTEGER INCX,INCY,N

DOUBLE PRECISION DX(*),DY(*)

DOUBLE PRECISION DTEMP

DTEMP = 0.0d0

!-------------------------------------------------------------------------------------------

DDOT=DTEMP 函数结束，函数返回值不变了,此时和C语言的 return DTEMP是一样的。

end function

!---------------------------------------------

如果在赋值处，函数结束，和C语言的 return DTEMP是一样的。

如又重新薯虚赋值，

DOUBLE PRECISION FUNCTION DDOT(N,DX,INCX,DY,INCY)

INTEGER INCX,INCY,N

DOUBLE PRECISION DX(*),DY(*)

DOUBLE PRECISION DTEMP

DTEMP = 0.0d0

DDOT=DTEMP ！函数返回值暂时是DTEMP，也就是 0.0d0双精度0

。。。。。

DDOT=DTEMP+200 ！现在是数迟燃函数返回值200了！

end function ！ 函数结束，函数返回值就是200了这两条语句联手，

!相当于return DTEMP+200 ；

double randomz (int ia, int ib)

{

double x//返回值

static int initial[15]={0}

static double iz,iy[15]={0.0}

//使用static类型，为下次调用保留值，不然每次都要从头开始

switch(ia)

{

case 1: //ia参数为1，从键盘输入种子初始化随机数生成器

iz=100000001.0

read(5,'(5i8)') initial 这句话直接翻成C很难，

//意思就是从键盘(5号设备)连续读入5个整数（难道不是15？

// 我怀疑源码写错了，应该是15i8），存入整型数组 initial ；

//而且Fortran的格式描述符i8强制每个整数都是8位（千万位），

//如果输入不足8位就切换到下一个数进行输入，则Fortran会

//自动在不足8位的数的右侧补上足够的0，以放大到千万位

iy=dble(initial) 这句话直接翻成C很难，

//本句用到Fortran90/95特色的数组整体 *** 作，C/C++要用循环；

//是把整型数组initial逐项复制给double数组iy，转换成double型

x = iy[ib]*1.0E-8

case 2: //旦缓租ia参数为2，继续使用已经初始化好了的生成器

iy[ib] = 329.0 * iy[ib] % iz

x = iy[ib]*1.0E-8

case 3: //ia参数为3，重新初始化，但还使用原先的种子

iy=dble(initial) 这句话直接翻成C很难，用循环完成。

}

//switch结束

return x

}

=================================

算法的主要思想就是“线性同余法”，

linear-congruential method

其基本迭代公式为

X[n+1] = ( A * X[n] + B ）% C

X的初始值随便取

在你给的源代码里：

A= 329.0

B= 0

C= 100000001.0

源代码中模兆的关键是 iy[ib] = 329.0 * iy[ib] % iz

另外 x = iy[ib]*1.0E-8 是为了将结果归一化到0~1之间再返回

你可以自己手工算几个数，就能看出这个算法的奥妙了。

另外需要指出的是——源代哪伍码里用static就是为了每次case 2时候的调用，都是在对上一次的结果进行迭代。而ib参数的用处就是保持有几组不同的独立迭代序列可用，防止不同用途的几处“生成伪随机数”调用互相干扰。

---