linux fortran 大数组,Fortran动态数组

这篇博客介绍了Fortran编程中的动态数组特性,包括如何使用allocatable声明和allocate/deallocate函数进行内存管理。还展示了如何使用data语句初始化数组,以及where语句在有条件地修改数组元素时的应用。示例代码详细解释了这些概念并展示了运行结果。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

动态数组是一种数组,其尺寸在编译时不知道,而是在执行时才已知/确定的。

动态数组的属性使用 allocatable 声明。 .

例如,

real,dimension(:,:),allocatable::darray

数组的秩,即尺寸,必须提到但是,分配内存以这样的阵列,可以使用allocate函数。

allocate(darray(s1,s2))

该阵列使用后,在该程序中,所创建的存储器应该使用 deallocate 函数解除

deallocate(darray)

示例

下面的例子演示了上面讨论的概念。

program dynamic_arrayimplicitnone!rankis2,but sizenotknown

real,dimension(:,:),allocatable::darray

integer::s1,s2

integer::i,jprint*,"Enter the size of the array:"read*,s1,s2!allocate memory

allocate(darray(s1,s2))doi=1,s1doj=1,s2

darray(i,j)=i*jprint*,"darray(",i,",",j,") = ",darray(i,j)enddoenddodeallocate(darray)endprogram dynamic_array

当上述代码被编译和执行时,它产生了以下结果:

Enter the size of the array: 3,4

darray( 1 , 1 ) = 1.00000000

darray( 1 , 2 ) = 2.00000000

darray( 1 , 3 ) = 3.00000000

darray( 1 , 4 ) = 4.00000000

darray( 2 , 1 ) = 2.00000000

darray( 2 , 2 ) = 4.00000000

darray( 2 , 3 ) = 6.00000000

darray( 2 , 4 ) = 8.00000000

darray( 3 , 1 ) = 3.00000000

darray( 3 , 2 ) = 6.00000000

darray( 3 , 3 ) = 9.00000000

darray( 3 , 4 ) = 12.0000000

使用data语句

data 语句可用于初始化多个阵列,或用于阵列部分的初始化。

data 语句的语法是:

data variable/list/...

示例

下面的例子演示了这一概念:

program dataStatementimplicitnone

integer::a(5),b(3,3),c(10),i,j

data a/7,8,9,10,11/data b(1,:)/1,1,1/data b(2,:)/2,2,2/data b(3,:)/3,3,3/data(c(i),i=1,10,2)/4,5,6,7,8/data(c(i),i=2,10,2)/5*2/Print*,'The A array:'doj=1,5print*,a(j)enddoPrint*,'The B array:'doi=lbound(b,1),ubound(b,1)write(*,*)(b(i,j),j=lbound(b,2),ubound(b,2))enddoPrint*,'The C array:'doj=1,10print*,c(j)enddoendprogram dataStatement

当上述代码被编译和执行时,它产生了以下结果:

The A array:

7

8

9

10

11

The B array:

1 1 1

2 2 2

3 3 3

The C array:

4

2

5

2

6

2

7

2

8

2

使用where语句

where语句可以使用数组中的某些元素在一个表达式,根据一些逻辑条件的结果。它允许表达的执行在一个元素上,如果给定的条件为真。

例子

下面的例子演示了这一概念:

program whereStatementimplicitnone

integer::a(3,5),i,jdoi=1,3doj=1,5a(i,j)=j-ienddoenddoPrint*,'The A array:'doi=lbound(a,1),ubound(a,1)write(*,*)(a(i,j),j=lbound(a,2),ubound(a,2))enddowhere(a<0)a=1elsewhere

a=5endwherePrint*,'The A array:'doi=lbound(a,1),ubound(a,1)write(*,*)(a(i,j),j=lbound(a,2),ubound(a,2))enddoendprogram whereStatement

当上述代码被编译和执行时,它产生了以下结果:

The A array:

0 1 2 3 4

-1 0 1 2 3

-2 -1 0 1 2

The A array:

5 5 5 5 5

1 5 5 5 5

1 1 5 5 5

¥ 我要打赏

纠错/补充

收藏

加QQ群啦,易百教程官方技术学习群

注意:建议每个人选自己的技术方向加群,同一个QQ最多限加 3 个群。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值