fortran95教程9.docx
《fortran95教程9.docx》由会员分享,可在线阅读,更多相关《fortran95教程9.docx(25页珍藏版)》请在冰豆网上搜索。
fortran95教程9
第9章 基本计算
(一)赋值语句
从计算机的实际状态来看,给出表达式并不意味着计算的开始,真正能够驱动计算的是可执行语句,而最直接驱动表达式的计算过程的就是赋值语句,完成一个赋值步骤,就意味着机器的状态的局部或全局发生了一个根本的变化。
从语法的角度来看,一个表达式还只是一个比较复杂的复合词汇,还不能构成一个完整的语句,而赋值语句则是一个在语法意义上的完整的语句,表达式在赋值语句中扮演关键的语法作用。
从计算问题的角度来看,赋值可以说是基本的计算步骤,考虑一下我们人工计算的过程就可以发现,任何一个完整的计算步骤,都可以说就是一个赋值步骤,特别是公式演算和数值计算,完成一个局部运算过程的标志,往往就是求出某个中间变量的数值或表达式。
FORTRAN语言作为一种以公式翻译为初衷的高级语言,它的赋值语句的一般形式就是一个数学等式,当然不是那种公式恒等变换得到的恒等式,而是要把运算的结果赋予一个变量,因此它的一般句法形式为:
variable=expression
variable=>expression
可以看到一个赋值语句作为一个完整的语句,分为三个部分:
●被赋值的变量
●赋值符号
●计算(表达式)
其中被赋值的变量,通过赋值,就拥有了明确的取值形式,如果该变量含有下标,片断下标或子串范围,则赋值的前提是它们都已经获得具体的取值。
基于表达式结果的不同种类,赋值语句分为:
● 固有赋值:
固有赋值是把任意类型的值赋予一个非指针变量,或把一个指针变量赋予一个与之相关联的目标。
● 自定义赋值:
自定义赋值得以构成的前提是存在一个可访问的子例行程序,它包含一个具有ASSIGNMENT形式的赋值界面,其属性与自定义赋值语句里的变量和表达式的属性保持一致。
● 指针赋值:
指针赋值把一个指针变量关联到它的目标对象,或者说把一个目标对象赋予到一个指针变量。
在形式上它使用符号=>。
● 过滤数组赋值:
这个赋值过程是对满足一定条件的数组元素进行赋值,而不是对数组整体进行赋值。
● 并行指标数组赋值。
这种赋值形式给出了一种有效的并行机制,能够大规模地对多重指标变量进行赋值。
下面我们将分节讨论这5类赋值形式,首先给出一些例子如下:
【例9-1】
X=X-2
实型的固有赋值
CHAR(5:
8)=“M_30”
逻辑型的固有赋值
SAMPLE=NOTE_2
结构的固有赋值
STRING=”MYSAMPLE”
可变串结构的自定义赋值
WHERE(X/=0.0)
A=B/X
ENDWHERE
过滤数组赋值
FORALL(I=0:
N,J=1:
N)
A(I,J)=3.0/(I-J+2)
B(I,J)%PTR=>C(I:
N,J:
N)
ENDFORALL
PTR=>X
并行指标结构
并行指标赋值
并行指标指针赋值
指针赋值
9.1 固有赋值语句
所谓固有赋值的主要特征就是被赋值的变量不能是指针变量,赋值语句本身就是对变量的定义或重定义,变量的取值就是直接计算赋值语句右边的表达式所得。
不需要额外的子程序来加以说明。
下面说明固有赋值语句的用法:
● 固有赋值中变量与表达式所能够具有的类型,以及它们的一致性要求见表9-1。
表9-1固有赋值中变量与表达式的类型
变量的类型
表达式的类型
INTEGER
INTEGER,REAL,COMPLEX
REAL
INTEGER,REAL,COMPLEX
COMPLEX
INTEGER,REAL,COMPLEX
CHARACTER
变量为具有相同种别参数的CHARACTER
LOGICAL
LOGICAL
派生类型
变量为相同的派生类型
● 如果变量是标量,那么表达式也必须是标量。
● 如果变量是数组,那么表达式或者是标量,或者是相同形状的数组。
● 如果变量是显形数组,那么变量的形状可以通过说明语句说明。
● 如果变量是待定形数组,那么它的形状由ALLOCATE语句,或指针赋值语句确定。
● 如果变量是哑形数组,那么它的形状由变量里的片断下标,或实元确定。
变量不能是哑尺度数组,除非存在含有最后一个维度上的上界的片断下标,或数组的下标向量,或者标量下标。
● 表达式的形状由算元的形状,表达式里的运算,以及其中的函数引用决定。
● 如果变量是一个指针,那么它必定已经关联到一个目标,赋值语句把表达式的值赋予指针的目标。
指针所关联的目标可以是一个数组,指针决定数组的秩,而每个维度的宽度则由目标决定。
● 在赋值过程本身被执行之前,赋值语句右边表达式以及表达式和变量里面包含的下标与下标片断表达式,都必须预先求值完毕。
● 如果变量的类型以及种别参数和表达式的不一致,那么在执行赋值语句之前,还必须完成必要的针对表达式的类型转换,固有转换函数见下表9-2:
表9-2 固有转换函数
变量的类型
被赋的值
INTEGER
INT(expression,KIND(variable))
REAL
REAL(expression,KIND(variable))
COMPLEX
CMPLX(expression,KIND(variable))
LOGICAL
LOGICAL(expression,KIND(variable))
● 表达式里面可以使用赋值语句左边变量的部分值。
【例9-2】
DATE(5:
10)=DATE(2:
6)
注意这是FORTAN现代版本所增加的功能,在FORTRAN77里面是不允许的。
● 如果变量和表达式都是字符型的,那么它们必须具有相同的种别参数值。
● 如果变量和表达式是具有不同长度属性的字符型对象,那么赋值过程遵循以下规则:
● 如果变量的长度比表达式的短,那么把表达式的右边长出来的部分截除;
● 如果变量的长度比表达式的长,那么在表达式的右边添加空格,直到跟变量一样长。
● 如果在赋值语句左边的变量里面含有表达式,例如下标表达式,那么这个表达式的计算与取值完全与赋值语句右边的表达式无关,因为左边变量的任何表达式都必须是预先完成计算的。
● 如果把一个标量赋值给一个数组,为了保持形状的一致性,需要把该标量扩充为一个与左边数组同样形状的数组,该数组的每个元素都是那个标量。
【例9-3】
INTEGERX(100)
X=5.0
经过赋值后X的取值为含有100个元素的数组,每个元素都是5.0。
● 数组的赋值过程是按照元素位置一一对应进行的。
【例9-4】设有两个相同形状的数组X和Y,有赋值语句:
X=Y
那么具体的赋值过程就是Y的元素赋予X的相同位置的元素,至于不同位置的元素是按照什么顺序来进行赋值的,并不需要考虑,可以认为所有元素的赋值是同时进行的。
● 对于派生类型的对象的固有赋值,变量与表达式必须是属于同一个派生类型,而赋值方法类似于数组,在相同位置的成员之间进行赋值,如果某个成员是指针,则该成员的赋值过程按照指针赋值的规则进行。
9.2 自定义赋值语句
如果需要进行赋值的变量和表达式是不符合表9-1,表9-2所规定的一致性要求的固有类型或派生类型,那么就需要使用自定义赋值语句,它提供赋值界面与相应子例行程序以供访问,从而完成赋值。
自定义赋值是由具有赋值说明符ASSIGNMENT(=)的子例行程序来完成赋值操作的。
用户可以在子例行程序当中通过定义新的规则,来扩充可以进行赋值的类型。
下面介绍自定义赋值运算的用法:
● 自定义赋值运算由具有两个哑元的子例行程序声明。
(子例行程序可以是外部或模块子例行程序里的对象。
)
● 子例行程序的哑元按照其出现的前后顺序,分别表示自定义赋值的变量与表达式。
其中第一个哑元必须具有INTENT的OUT或INOUT的属性,而第二个哑元必须具有INTENT的IN的属性。
● 具有ASSIGNMENT(=)形式的通用说明符的子例行程序必须带有界面块。
● 变量与表达式的类型以及种别参数必须与哑元保持一致。
● 对于一个非基本子例行程序,变量与表达式的秩必须与哑元的秩相同。
● 对于一个基本子例行程序,变量必须是数组,而表达式与之保持一致;或者两者都是标量。
● 如果变量和表达式同时与一个基本子例行程序和一个非基本子例行程序的界面匹配,那么赋值运算由非基本子例行程序定义。
● 两个哑元或者有一个是派生类型,或者都是固有类型,但不满足固有赋值的一致性条件。
● 对变量的自定义赋值的结果由所引用的子例行程序决定。
【例9-5】
INTERFACEASSIGNMENT(=)
ELEMENTALSUBROUTINGRATIONAL_TO_REAL(L,R)
USERATIONAL_MODULE
TYPE(RATIONAL),INTENT(IN):
:
R
REAL,INTENT(OUT) :
:
L
ENDSUBROUTINGRATIONAL_TO_REAL
ELEMENTALSUBROUTINGREAL_TO_RATIONAL(L,R)
USERATIONAL_MODULE
REALINTENT(IN) :
:
R
TYPE(RATIONAL),INTENT(OUT):
:
L
ENDSUBROUTINGREAL_TO_RATIONAL
ENDINTERFACE
上面的界面块通过两个外部子例行程序给出了自定义赋值,一个把派生类型RATIONAL的对象赋值给实型对象,另一个把实型对象赋值给派生类型RATIONAL的对象。
有了上面的界面块,就可以进行下面的自定义赋值:
REALR_VALUE,R_ARRAY(20)
TYPE(RATIONAL)RAT_VALUE,RAT_ARRAY(20)
R_VALUE=RATIONAL(5,8)
RAT_VALUE=6.2
9.3 指针赋值
所谓指针赋值,实际上就是把指针变量关联到一个具有TARGET属性的对象上,即使得指针变量成为目标对象的一个“别名”。
如果指针所关联的目标的状态发生改变,例如去关联,或去定义,那么指针的状态也发生相应的改变。
一旦指针赋值语句被执行之后,指针的关联状态就不会改变了,除非执行另外一个指针赋值语句来改变其状态,或执行ALLOCATE,DEALLOCATE,NULLIFY语句对指针重定义。
指针赋值的一般形式(R736)为:
pointer-object=>target
其中指针对象(pointer-object)的一般形式(R630)为:
variable-name
structure-component
而目标(target)的一般形式(R737)为:
variable
expression
下面说明指针赋值的用法:
● 如果指针对象为变量名称或结构成员,那么它们都必须具有POINTER属性。
● 如果目标是命名变量,那么它必须具有TARGET或POINTER属性。
● 如果目标是子对象指示符,那么它的父结构必须具有TARGET或POINTER属性,或者它必须是一个在最右边成员具有POINTER属性的结构引用。
● 目标与指针对象的类型,种别参数,长度参数,秩都必须一样。
● 如果=>右边的变量具有TARGET属性,那么=>左边的指针对象就被关联到该目标。
● 如果=>右边的变量具有POINTER属性,并且处于关联状态,那么=>左边的指针对象在指针赋值语句执行之后,也被关联到目标指针所关联的同一个数据对象。
● 如果=>右边的变量具有POINTER属性,并且处于去关联状态,或者右边的表达被引用到固有函数NULL,那么=>左边的指针对象也成为去关联状态。
● 如果=>右边的变量具有POINTER属性,并且处于去定义关联状态,那么=>左边的指针对象也成为去定义关联状态。
● 指针赋值语句终止了指针对象此前所具有的任何关联状态,从而产生一个新的关联状态。
● 如果指针对象是待定形数组,那么指针赋值语句就能够确立指针数组每个维度的宽度,即其目标的相应维度的宽度,除非目标本身是去关联指针或去定义指针。
【例9-6】执行以下语句:
INTEGER,TARGET:
:
T(11:
20)
INTEGER,POINTER:
:
P1(:
),P2(:
)
P1=>T
P2=>T(:
)
得到P1的宽度是11和20,而P2的宽度是1和10。
● 目标不能是属于哑尺度数组的命名变量。
● 如果目标是一个哑尺度数组的数组片断,那么它必须具有下标,或者具有在最后一维给出了上界的下标三元组。
● 如果目标为数组片断,那么它不能具有下标向量。
● 如果目标是一个表达式,那么它必须给出一个指针结果,这意味着该表达式引用了固有函数NULL,或者是引用了能够得到指针结果的用户自定义函数或自定义运算。
因为在固有函数或固有运算里,只有NULL得到指针结果。
● 如果一个指针的目标没有被引用或定义,那么该指针也不能被引用或定义。
● 如果一个结构的某些成员具有POINTER属性,而该结构被固有派生类型赋值语句进行了赋值,那么对于每个具有POINTER属性的成员,是运用指针赋值语句进行赋值的。
自定义赋值语句可以用来在结构的成员之间进行指针赋值。
【例9-7】
MONTH=>DAYS(1:
30)
PTR=>X(:
5)
NUMBER=>INDEXONE%INDEXTWO
FIRST_RESLUT=>NULL()
【例9-8】下面的目标也是指针:
REAL,POINTER:
:
PTR,P
REAL,TARGET:
:
X
REALY
X=5.7
P=>X
PTR=>P
Y=PTR-8.1
【例9-9】下面的目标为表达式:
INTERFACE
FUNCTIONPOINTER_FCN(X)
REALX
REAL,POINTER:
:
POINTER_FCN
ENDFUNCTION
REAL,POINTER:
:
P
REALA
P=>POINTER_FCN(A)
注意以下情形:
●在指针与目标之间建立关联,除了指针赋值语句之外,使用ALLOCATE语句也能做到,而使得指针去关联则也可使用DEALLOCATE或NULLIFY语句。
●当在表达式或过程里引用指针时,会得到与该指针关联的目标的值,这点与在指针赋值语句里不同。
●当指针是作为一个与具有POINTER属性的哑元相应的实元出现的时候,那么该引用是到达指针本身,而不是它的值。
9.4 过滤数组赋值
对于数组的赋值,除了全数组的一一对应的赋值方式之外,很多时候,还需要做到在数组里挑选符合一定条件的某些元素,只对这些特定元素进行赋值。
这时已有的赋值语句都无法做到,因此FORTRAN专门提供了过滤数组赋值方式,用来完成这类任务。
所谓过滤数组赋值,实际上就是建立WHERE块,或ELSEWHERE块,或者是WHERE语句,然后在其中使用固有赋值语句,它的变量为数组时,就可以对其元素进行选择性赋值。
在WHERE语句或WHERE结构里对数组元素进行过滤,是根据一个逻辑数组表达式的值来进行的,该逻辑数组表达式是从数组取值而得到逻辑型值的条件表达式,作用于:
● 基本固有运算;
● 基本固有函数引用;
● 基本自定义运算;
● 基本自定义函数引用;
● 固有赋值;
● 基本自定义赋值。
注意所谓WHERE语句或WHERE结构非常类似于选择控制结构IF结构或IF语句,但它们有一个主要的差别,就是在WHERE结构里,所有的赋值语句都是要执行的,即使是不满足条件的元素,它的赋值同样执行了,只是不满足条件的元素的赋值对随后的程序不其任何作用而已。
但是在IF结构里,不满足条件的情形是不被执行的。
9.4.1 WHERE语句
WHERE语句的句法形式(R738)为:
WHERE(logical-expression)array-assignment-statement
其中逻辑表达式(logical-expression)对所有数组元素起作用,就相当于一个过滤网,使得表达式结果为TURE的数组元素被赋值给变量,而表达式结果为FALSE的数组元素不被赋值给变量。
【例9-10】
WHERE(CURRENT>1.00)PASS=CURRENT
WHERE(TEMP<0.40)USEFULL=TEMP
9.4.2 WHERE结构
WHERE结构的句法形式(R739)为:
[where-construct-name:
]WHERE(logical-expression)
[where-body-construct]…
[ELSEWHERE(logical-expression)[where-construct-name]
[where-body-construct]…]…
[ELSEWHERE[where-construct-name]
[where-body-construct]…]
ENDWHERE[where-construct-name]
其中where体结构的句法形式(R741)为:
assignment-statement
where-construct
where-statement
在初始WHERE语句当中出现的逻辑表达式构成了一个过滤网,可以控制在WHERE体结构里给出的数组赋值语句的表达式取值和赋值。
如果逻辑表达式在ELSEWHERE语句当中也出现,那么该语句就成为了过滤ELSEWHERE语句。
【例9-11】
INTEGER:
:
N(9)=(/1,2,3,4,5,6,7,8,9/)
WHERE(MOD(N,2)==0)
N=2 !
现在N=(/1,2,3,2,5,2,7,2,9/)
ELSEWHERE(MOD(N,3)==0)
N=3 !
现在N=(/1,2,3,2,5,2,7,2,3/)
ELSEWHERE(MOD(N,5)==0)
N=5 !
现在N=(/1,2,3,2,5,2,7,2,3/)
ELSEWHERE
N=0 !
现在N=(/0,2,3,2,5,2,0,2,3/)
ENDWHERE
下面说明WHERE结构的用法:
● 一个结构名称不能单独出现在初始WHERE语句,或ENDWHERE语句里,必须成对出现。
如果一个结构名称出现在ELSEWHERE语句里,或者是过滤ELSEWHERE语句里,那么它也必定出现在WHERE结构语句里。
● 在WHERE结构的每个赋值语句里,被定义的变量必须和过滤网具有相同的形状,如果在WHERE结构的其他位置也出现了过滤网,它们的形状必须都相同。
● WHERE结构里的语句按照其出现的顺序执行。
● 每个过滤网的有效取值只有一次,随后的逻辑表达式里的取值变化都对控制过滤网的值没有影响。
● WHERE结构里的自定义赋值必须由基本子例行程序定义。
● 只有在控制过滤网取真值时,赋值语句里的基本运算或函数才得到计算。
【例9-12】
REALX(1,100)
…
WHERE(X>0.0)
SQRTX=SQRT(X)
ENDWHERE
● 如果在结构的逻辑表达式或赋值语句当中出现了数组构造器,那么它将不经过过滤而得到完全计算。
● 结构里的表达式或变量或逻辑表达式里的非基本函数引用将得到完全计算,尽管有可能结果数组的所有元素都用不上。
【例9-13】
REALA(2,3),B(3,10),C(2,10),D(2,10)
IMTRINSICMATMUL
…
WHERE(D<0.0)
C=MATMUL(A,B)
ENDWHERE
其中矩阵乘积的所有元素都得到了,但是只有小于0的元素才得到赋值。
● 在WHERE结构里,只有WHERE语句才能成为分支目标。
9.4.3 WHERE结构的运行
当WHERE结构开始运行的时候,其中的控制过滤网和候选过滤网就都已经建立了,而控制过滤网主宰了随后的语句块的执行。
● 如果WHERE结构不是嵌套WHERE结构,那么控制过滤mask取得逻辑表达式的值,而候选过滤网的值为.NOT.mask。
【例9-14】
WHERE(A1) !
语句1
CALLB1 !
块1
ELSEWHERE(A2) !
语句2
CALLB2 !
块2
ELSEWHERE !
语句3
CALLB3 !
块3
ENDWHERE
在执行语句1之后,控制过滤网的值为A1,而候选过滤网的值为.NOT.A1;
在执行语句2之后,控制过滤网的值为(.NOT.A1).AND.A2,而候选过滤网的值为(.NOT.A1).AND(.NOT..A2);
在执行语句3之后,控制过滤网的值为(.NOT.A1).AND(.NOT..A2)。
● 如果WHERE结构是嵌套WHERE结构,还是以上面的例子来说明运行过程:
如果最外层主宰着块内的语句运行的初始控制过滤网是outer-mask,那么块1的控制过滤网取值为outer-mask.AND.A1,从这个初始状态开始,就可以得到向内嵌套过程的每个控制过滤网取值,即总是向内取逻辑与运算。
到运行内部ENDWHERE语句时,控制过滤网就回复到outer-mask。
9.4.4 WHERE结构与控制结构的差别
WHERE结构与控制结构的主要差别在前面9.4已经讨论过,除此之外,还有个差别,就是在WHERE结构里,不存在从WHERE块或ELSEWHERE块里转移出来的机