fortran第7章2讲课讲稿.ppt

fortran第7章2; Fortran90提供了丰富的数组操作函数，但它们的函数值可以是数组，但不一定全是数组。 （1）函数值个数≥2，函数值一定为数组 （2）函数值个数=1，一般函数值不是数组， 但maxloc、minloc、shape函数值是数组 ;INTEGER,DIMENSION(1:2,1:3)::A ;例1：阅读下列程序，写出运行结果 PROGRAM main IMPLICIT NONE INTEGER,DIMENSION(2,3,2)::M INTEGER,DIMENSION(2,2)::S M=RESHAPE((/1,2,3,4,5,6,7,8,9,1, 2,3/),(/2,3,2/)) S=SUM(M,2,MASK=M3) PRINT*,S END ;练习1：阅读下列程序，写出运行结果 PROGRAM main IMPLICIT NONE INTEGER,DIMENSION(2,-1:1,2)::a INTEGER::i a= RESHAPE((/(i, i=1,12)/),(/2,3,2/)) print*, sum(a,3) end;例:INTEGER,DIMENSION(2,3)::A size(a) size(a,DIM=1) size(a,2) ;例: INTEGER,DIMENSION(2,3)::A MAXVAL(A) MAXVAL(A,2) MAXVAL(A,1,MASK=A2);例：INTEGER,DIMENSION(2,3)::A INTEGER,DIMENSION(2)::B A= RESHAPE((/1,2,6,3,4,5/),(/2,3/)) B=MAXLOC(A,A4) 例：INTEGER,DIMENSION(2,3)::A INTEGER,DIMENSION(3)::B A= RESHAPE(/1,2,6,3,4,5/,(/2,3/)) B=MAXLOC(A,1) ;INTEGER,DIMENSION(0:2)::A INTEGER,DIMENSION(1)::B A= (/1,2,6/) B=MAXLOC(A) PRINT*,B END ;7、 PRODUCT (数组名[,DIM=n][,MASK=数组逻辑表达式]) 求数组各元素之积 注意：没有满足MASK条件的元素时，结果为1 MASK=可省略不写 例如：PRODUCT(A, 2, A3) ;8、 shape (数组名) 求数组的形状 形状：各维维长的序列 判定维数，求各维维长d2-d1+1 函数值为数组 ;练习2：阅读下列程序，写出运行结果 IMPLICIT NONE INTEGER，DIMENSION（0:1，-1:1，2）：：M INTEGER，DIMENSION(3)::A M=RESHAPE（（/1，2，3，4，5，6，7，8，9， 1，2，3/），（/2，3，2/）） A= SHAPE（M） PRINT*，A END;例2：阅读下列程序，写出运行结果 PROGRAM MAIN IMPLICIT NONE INTEGER,DIMENSION(3)::M INTEGER,DIMENSION(1)::A M=(/1,2,3/) A= SHAPE(M) PRINT*,A END;练习3：阅读下列程序，写出运行结果 PROGRAM main INTERFACE SUBROUTINE S(b) INTEGER,DIMENSION(:,:)::b END SUBROUTINE ENDINTERFACE INTEGER,DIMENSION(2,3)::a a=RESHAPE((/1,2,3,4,5,6/),(/2,3/)) Call s(a) End SUBROUTINE S(b) INTEGER,DIMENSION(:,:)::b Print*,shape(b); End;9、 DOT_PRODUCT (数组名1, 数组名2) 求两个同形状的一维数组（向量）的点积 点积： REAL，DIMENSION(3)::a REAL，DIMENSION(0:2)::b DOT_PRODUCT (a,b): a(1)*b(0)+a(2)