MetMan's Blog

标签搜索

搜索到 13 篇与的结果 ———

Fortran数组转置

Fortran数组转置

数组转置（transpose）是一个常见的操作。二维数组转置Fortran标准有一个内置函数transpose可以对Rank-2数组进行转置操作。program main implicit none integer :: a(3,4) = reshape([1,2,3,4,5,6,7,8,9,10,11,12], [3,4]) integer,allocatable :: b(:,:) b = transpose(a) print*, 'shape(b) = ', shape(b) print*, 'b = ', b end program更高维数组转置但对于高维数组，如何转置？显式循环使用显式嵌套循环进行数组转置，下面使用do concurrent结构对一个4维数组转职。dims = shape(x) allocate(res(dims(4), dims(3), dims(2), dims(1))) do concurrent(i = 1:dims(1), j = 1:dims(2), k = 1:dims(3), l = 1:dims(4)) ...

2025年03月22日

75 阅读

0 评论

Fortran Associate结构学习

Fortran Associate结构学习

只要代码写得多，总能Get新知识。Fortran 2003标准中associate结构应该不常用，但associate结构用好可以简化复杂语句，增加代码的可读性。可以将associate关联名字看成是要指向的表达式或者变量的别名（引用），供后续代码使用。语法[name:] ASSOCIATE (assoc-entity[, assoc-entity]...) block END ASSOCIATE [name]其中name是可选的，associate结构名称，如果指定，后面的end assoicate也必须加上相同的名称。assoc-entity语法为assoicate-name => selector，其中selector是一个表达式或者变量，而associate-name是关联selector的别名。block是零条或者多条Fortran语句/结构。示例：program associateTest implicit none integer :: a=1,b=1 associate( x => a*b ) ! 表达式的引用 ...

2025年03月22日

71 阅读

0 评论

ODC库编程实践

ODC库编程实践

前面提到过odb-2数据格式看成是二维表格型数据，但表格中每列的数据类型可以不同。开发人员可能考虑到通用性以及读写性能，在实现时将表格数据统一处理成双精度的二维数组（也可以看成更一般的一维数组）。real(8), target :: data(ncols, nrows) rc = encoder%set_data(data, column_major)此时需要解决两个问题：二维数组在内存中的布局，即以列主元还是行主元的内存布局。如果column的类型不是双精度，那么就存在数据类型转换，对于数值类型还算简单，但字符类型如何处理？第一个问题，我们从ODC API看起。C API odc_encoder_set_data_array函数用于设置编码的输入数组。函数第二个参数是指定的输入数组，这里使用const void*类型符，通用类型指针。第三个参数是指数组（表格）的宽度，以字节为单位。第四个参数是数组的行数。最后一个参数是指列主元布局每列宽度，以字节为单位，一般为8字节，如果是行主元，传入0值。对应的Fortran API，则是类对象encoder的成员函数set_data，查看...

2024年12月02日

101 阅读

0 评论

如何查阅Fortran语言标准

如何查阅Fortran语言标准

如果你对Fortran的语法细节感兴趣，希望了解标准如何定义的，可以阅读官方Fortran标准说明文档。各大编译器厂商基本上按照Fortran标准实现自家的编译器，再加上一些各自的扩展功能。由于ISO组织拥有Fortran最终标准的版权，无法免费获取正式标准的电子文档，需要花钱购买它们的纸质或电子文档。但有一个网上公开合法获取Fortran标准说明文档途径是美国Fortran编程语言标准技术委员会（J3）提供的内部工作版本。Fortran标准由"Working Group 5 of Sub-Committee 22 of the Joint Technical Committee 1 of the International Organization for Standardization and the International Electrotechnical Commission (IEC)"（ISO/IEC JTC1/SC22/WG5）负责开发，这个小组简称为WG5。官方网站：https://wg5-fortran.org而前面提到的J3是ISO/IEC JTC1/SC...

2024年11月03日

166 阅读

0 评论

Fortran标准库系列：查询变量占用内存大小

Fortran标准库系列：查询变量占用内存大小

C语言提供sizeof()运算符计算数据类型对象占用内存字节数，常用于malloc分配对象空间时使用。那么Fortran语言中是否有类似函数计算对象大小？Fortran标准实际上提供了STORAGE_SIZE和C_SIZEOF内置函数可以实现该功能，此外GNU编译器提供了额外的扩展SIZEOF支持该功能。STORAGE_SIZEFortran 2008标准内置函数，函数返回参数A的存储大小，单位是bits。RESULT = STORAGE_SIZE(A [, KIND])参数A是任意类型的标量或者数组。可选参数KIND是标量整型常量表达式，用于指定返回值的整型KIND。需要注意的是，如果传入的参数A是数组，返回的是数组一个元素占有的位数，而不是整个数组存储空间大小。SIZEOFGNU扩展，不是Fortran标准内置函数。计算参数X占用的存储字节数，单位是字节（bytes）。函数假设数组连续。N = SIZEOF(X)返回值类型是INTEGER(C_SIZE_T)，值是参数占用的字节数。如果参数是POINTER属性，返回指向的存储区域字节数。如果参数是包含POINTER或ALLOC...

2024年11月03日

132 阅读

0 评论

Fortran Sequential二进制文件读入问题诊断

Fortran Sequential二进制文件读入问题诊断

问题现有一个二进制文件rain_clmn_202306.0000.dat及对应的GrADS CTL数据描述文件，使用GrADS绘图正常显示，但通过Fortran程序读入数据会报错。CTL数据描述文件如下所示：DSET ^rain_clmn_202306.0000.dat OPTIONS sequential big_endian UNDEF 9999.000 TITLE Result for CTRL,MWR, XDEF 45 linear 5.0000 9.0000 YDEF 8 levels 3.0000 6.0000 9.0000 12.0000 15.0000 18.0000 21.0000 24.0000 ZDEF 1 levels 0.000 TDEF 30 linear 00Z01JUN2023 1dy VARS 12 CTRL1 1 0 CTRL (TS) CTRL2 1 0 CTRL (MISSED) CTRL3...

2024年09月22日

185 阅读

0 评论

Fortran标准库系列：MOD vs MODULO

Fortran标准库系列：MOD vs MODULO

取余运算 vs 取模运算取余运算（remainder）在计算商时，会将商向0方向舍入。这意味着取余运算尽可能让余数的绝对值最小。取模运算（modulus）在计算商时，会将商向负无穷方向舍入。取模运算尽可能让商的值最小。当两个数都是正数时，取余运算和取模运算的结果是一样的，因为它们的商都是向0方向舍入的。然而，当两个数符号不同时，取余运算和取模运算的结果会不同。Fortran标准提供取余函数MOD和取模函数MODULO。MOD函数计算A除以P的余数。R=MOD(A,P)参数A是INTEGER或者REAL类型的标量。参数P必须与A相同类型和精度，不能为0。函数返回结果是$A-(INT(A/P)* P)$ ，返回值R符号与参数A相同，且$|R|<|P|$。MODULO计算A对P取模结果。R=MODULO(A,P)参数A是INTEGER或REAL类型的标量。参数P必须与A相同类型和精度,。返回结果类型和精度与参数相同，且返回值R符号与P相同,且$|R|<|P|$。。如果A和P是INTEGER类型，返回结果R满足$A=Q*P+R$ ，其中Q是一个整数，$R\in [0,P)$ ...

2024年08月10日

335 阅读

0 评论

Fortran面向对象（三）

Fortran面向对象（三）

继续Fortran面向对象的学习。Unlimited polymorphic typesclass(*) 类型称为无限制多态类型，必须拥有pointer或allocatable属性。括号中用*代替具体的类名表示可以是任意类。（联想C语言中void *通用指针概念）一般不直接使用无限制多态类型，常用作过程的哑元，使用select type结构转换为指定类型使用。CLASS(*), POINTER :: unlimited_ptr unlimited_ptr => base_target SELECT type(unlimited_ptr) TYPE IS (base) PRINT *, "base type: component value: ", unlimited_ptr%i TYPE IS (child) PRINT *, "child type: component values: ", unlimited_ptr%i,unlimited_ptr%j END SELECTpass ...

2024年08月10日

390 阅读

1 评论

查看更多

MetMan博主神仙

访问 66773 ℃

145 文章数

5 评论量

互动读者

标签云

最新回复

tqymnonccc打酱油
2024-09-27

怎么收藏这篇文章？
toibdpojay打酱油
2024-09-22

博主真是太厉害了！！！
yvctxyevvw打酱油
2024-09-22

博主真是太厉害了！！！
frezhwzwuq打酱油
2024-09-22

博主真是太厉害了！！！

繁

登录

X