Fortran内置函数sum()在real*4溢出情况下的表现

最新推荐文章于 2023-09-01 13:55:03 发布
最新推荐文章于 2023-09-01 13:55:03 发布
阅读量1.1k 收藏 1
点赞数 1
分类专栏： Fortran 文章标签： c++ 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_40759055/article/details/125440758
版权

 变量说明

（1）变量前缀为‘fortran’表示在Fortran中进行求和操作，且当后缀不为‘DFC’时表示为使用Fortran读取文件，数据类型为real*4；

（2）变量后缀为‘DFC’（data from c）表示为数据来自于C程序，数据类型为real*4；

（3）变量前缀为‘c’表示为使用C读取文件，数据类型为float，在C中进行求和操作；

（4）文件相同，读取位置相同，读取数据大小相同，共计1440万个float类型数据；

（4）数据最小值为1649.254272，最大值为12704.22168。

fortran_sum_float使用sum()对数据进行求和，结果存入real*4变量中
fortran_sum_double使用sum()对数据进行求和，结果存入real*8变量中
fortran_sum_auto_float自己编写sum()对数据进行求和，以real*4变量进行累加
fortran_sum_auto_double自己编写sum()对数据进行求和，以real*8变量进行累加
fortran_sum_float_DFC使用sum()对数据进行求和，结果存入real*4变量中
fortran_sum_double_DFC使用sum()对数据进行求和，结果存入real*8变量中
fortran_sum_auto_float_DFC自己编写sum()对数据进行求和，以real*4变量进行累加
fortran_sum_auto_double_DFC自己编写sum()对数据进行求和，以real*8变量进行累加
c_sum_auto_float自己编写sum()对数据进行求和，以real*4变量进行累加
c_sum_auto_double自己编写sum()对数据进行求和，以real*8变量进行累加

数据验证

编译器：mpiifort/ifort编译器：mpiicc
VALUEVALUE
fortran_sum_float74177954000
fortran_sum_double74177953792
fortran_sum_auto_float74177954000c_sum_auto_float75781332992
fortran_sum_auto_double75781374529c_sum_auto_double75781374529
fortran_sum_float_DFC74177954000
fortran_sum_double_DFC74177953792
fortran_sum_auto_float_DFC74177954000
fortran_sum_auto_double_DFC75781374529

        从数据可以看出不管是通过Fortran还是C进行文件读取，以double类型的变量进行累加后的结果一致，程序也对两组数据进行过一一比较，结果数据一致。

更换编译器

编译器：gfortran编译器：g++
VALUEVALUE
fortran_sum_float69373526000
fortran_sum_double69373526016
fortran_sum_auto_float69373526000c_sum_auto_float69373526016
fortran_sum_auto_double75781374529c_sum_auto_double75781374529
fortran_sum_float_DFC69373526000
fortran_sum_double_DFC69373526016
fortran_sum_auto_float_DFC69373526000
fortran_sum_auto_double_DFC75781374529

        在更换编译器后，Fortran的sum求和以及自己编写的sum求和与之前对比都发生了变化，可见不同的编译器对float数据溢出存在不同的处理方式，由此将导致不同程度的精度损失。

        以上测试Fortran代码与C代码在编译时都使用了-O3进行优化。

 取消-O3代码优化

编译器：gfortran编译器：g++
VALUEVALUE
fortran_sum_float69373526000
fortran_sum_double69373526016
fortran_sum_auto_float69373526000c_sum_auto_float69373526016
fortran_sum_auto_double75781374529c_sum_auto_double75781374529
fortran_sum_float_DFC69373526000
fortran_sum_double_DFC69373526016
fortran_sum_auto_float_DFC69373526000
fortran_sum_auto_double_DFC75781374529

         以上数据没有变化。

编译器：mpiifort/ifort编译器：mpiicc
VALUEVALUE
fortran_sum_float69373526000
fortran_sum_double69373526016
fortran_sum_auto_float69373526000c_sum_auto_float75781332992
fortran_sum_auto_double75781374529c_sum_auto_double75781374529
fortran_sum_float_DFC69373526000
fortran_sum_double_DFC69373526016
fortran_sum_auto_float_DFC69373526000
fortran_sum_auto_double_DFC75781374529

        从数据可以看出C的编译器不会因为-O3而改变对float数据溢出的处理方式，但mpiifort会受-O3的影响。

        以最终结果存入real*8变量的情况进行对比，增加-O3后的Fortran的sum()求和后的相对误差由0.084557037下降为0.021158507。

冰冰冰冰淇淋君
关注
专栏目录
Python 标准库:: 内置类型
Simon112119的博客
05-20 495
内置类型 以下部分描述了解释器中内置的标准类型。 主要内置类型有数字、序列、映射、类、实例和异常。 有些多项集类是可变的。 它们用于添加、移除或重排其成员的方法将原地执行，并不返回特定的项，绝对不会返回多项集实例自身而是返回 None。 有些操作受多种对象类型的支持；特别地，实际上所有对象都可以比较是否相等、检测逻辑值，以及转换为字符串（使用 repr() 函数或略有差异的 str() 函数）。 后一个函数是在对象由 print() 函数输出时被隐式地调用的。 逻辑值检测 任何对象都可以进行逻辑值的检测，以
numpy学习
li_xin_smile的博客
07-26 663
Python是一种计算机程序设计语言. 完成同一个任务，C语言要写1000行代码，Java只需要写100行，而Python可能只要20行. 代码少的代价是运行速度慢，C程序运行1秒钟，Java程序可能需要2秒，而Python程序可能就需要10秒. *********** 命令行交互模式 *********** 在Windows开始菜单选择“命令提示符”，就进入到命令行模式，它的提示符类似 C:\> 在命令行模式下敲命令python ，然后就进入到Python交互模式
描述常用统计变量计算的ＦＯＲＴＲＡＮ程序
05-06
调用计算最大值和最小值的冒泡程序，计算平均值，标准差，计偏态、峰态系数的函数等 可以计算的统计变量包括下列各项： sample size 样本大小 平均值 中位数 方差 样本标准差 standard deviation　 样本平均数标准误差：standard error of mean the maximum of the sample values the minumum of the smaple values 样本值总和 SUM 极差， 全距 变异系数 Coefficient of variation 离散系数 Coefficient of dispersin 偏态系数，偏度系数 kurtosis coefficient 峰度系数，峰态系数，峭度系数 百分位数 5，10，15，20，25，....99% 百分位数值 也包括计算频数分布表
Fortran大数输出*科学计数法*指数溢出问题
weixin_41671554的博客
06-18 2137
Fortran大数输出*科学计数法*指数溢出问题指数溢出问题解决方案 指数溢出问题 由于输出格式的限制，在输出数据时，格式长度不够，导致将科学计数法的E标志丢失，引起数据的错误 解决方案 修改输出格式，将 e24.16 改为 e24.15e3, 从而使得指数上的长度满足要求 图示为将原输出数据格式进行了修改：将指数上的长度扩展成三位；为保持数据长度，将小数部分缩减一位。 修改后数据输出为 不再溢出，满足数据计算要求。 ...
一个简单的Fortran程序
最新发布
wu5_liu6_qi7的博客
09-01 301
这个程序使用Fortran语言的基本语法，包括声明变量、为变量赋值、计算和输出结果等。Fortran是高阶语言，适合用于科学计算和工程领域。
Fortran 浮点型變量（REAL）的類型
LookIn的专栏
10-09 1万+
!單精度浮點數的聲明方法（使用4個字節） real(kind=4) x_float !Fortran 90 real(4) y_float !Fortran 77 real*4 z_float !Fortran 77 !雙精度浮點數的聲明方法（使用8個字節） real(kind
Fortran学习6：数组3
faltas的博客
10-02 1656
Fortran学习
Fortran 程序集
11-08
real t, sum i = 1 t = 1 sum = 0 do while (abs(t) >= 1e-6) sum = sum + t t = t / i i = i + 1 end do print *, 'e =', sum end program ``` 该程序通过循环累加每一项的值来计算自然对数底数e的...
Matlab基本语法（一）
既然选择远方，便只顾风雨兼程！
07-20 8475
MATLAB（矩阵实验室）是用于数字计算，可视化和编程的第四代高级编程语言和交互式环境。MATLAB由MathWorks开发。https它允许矩阵操作；功能和数据绘制；算法的实现；创建用户界面；与使用其他语言（包括C，C++，Java和FORTRAN）编写的程序进行交互；分析数据；开发算法；并创建模型和应用程序。它具有许多内置命令和数学函数，可帮助您进行数学计算，生成图和执行数值方法。MATLAB环境的行为就像一个超级复杂的计算器。您可以在>>命令提示符下输入命令。MATLAB是一个解释型环境。......
Lu脚本教程
Lu的专栏
06-28 2102
欢迎访问Lu程序设计 Lu脚本教程 目 录 [页首] 1概述 2安装OpenLu 3Lu基础知识 3.1第一个程序 3.2简单的例子 3.3标识符3.4表达式及函数简介3.5源代码格式简介3.6o函数 3.7赋值语句 3.8常用算术运算符 3.9逗号（冒号、分号）运算符和括号运算符 4Lu常量 5Lu变量 5.1自变量 5.2动态变量 5.3静态变量 5.4模块变量 5....
识别FORTRAN中的double类型数据
01-10
本资源是编译原理课程设计中关于double数据识别的源码
2-编程练习
weixin_43515797的博客
10-08 181
编写一个求和函数sum() function sum() { var result=0; for(var i=0;i<arguments.length;i++){ result+=parseFloat(arguments[i])||0;//去除非数字参数 } return result.toFixed(3)/1;//去掉溢出小数 } var a=sum(1,2,3,4,5); var b=sum(0.1,0.2,false,true,'a','F',6)
sql中常用函数的用法
Hebe
12-24 6274
1、sql中sum和count的区别                        ----sql聚合函数       （1）首先，sum是对一个字段求和，hive中字段的类型一般是string或者是int，如果是int当然没有问题，如果是string类型但是全部是数字也没有问题，但是如果string类型中包含一个字母，sum的结果将会是0       例如如下的表：        user
fortran 学习笔记
plkolili的博客
10-11 3051
20180907运行 CRTL+F5 编译器:一般编译器本身都是控制台程序(console)，没有界面。 命令行 IDE :集成开发环境，用于提供程序开发环境的应用程序，一般包括代码编辑器、编译器、调试器和图形用户界面等工具。如Microsoft Visual Studio和Code::Blocks等 编辑器 : UEStudio Intel Visual Fortran 配合...
浮点数的存储(IEEE754)
秦黎明的博客
06-21 1108
VS2008 MSDN 中有 浮点数格式的详细介绍，本文只是一个总结性的东西，欲 详细研究请参考 （1）MSDN中输入：IEEE，找到浮点数相关部分 （2）可以参考网上一些文章          http://www.cnblogs.com/jillzhang/archive/2007/06/24/793901.html          http://blog.csdn.net/gan
Fortran入门 - 基本语法一篇通
热门推荐
sodaxyh的博客
08-19 4万+
0. 前些日子为了准备汇报学了Fortran这个语言，最近又要用到发现有些些忘记了，所以写一篇博客整理记录学了些啥就当复习~ 主要内容还是比较基本的语法，一些小细节等以后有机会来补上吧（下次一定.jpg） 1.固定格式&自由格式 （图源自网络 侵删） 简而言之： 固定格式——比较老，老古董代码一般都是这种格式，需要读的时候能读懂就行，不用会写。 自由格式——目前的使用主流。 2.输入输出 入门一个语言从Helloworld开始！顺便熟悉一下Fortran的程序结构。 ..
Fortran简明自学手册
蒸发杰作的博客
10-27 6430
Fortran简明自学手册 Fortran即Formula Translation（公式翻译）。是世界上诞生最早的高级语言之一，其易读，好写，很适合科学计算，一诞生便很快统治了数值计算界。因此也不难理解为何Fortran中保留了一些特殊的基础类型，如虚数。过去学习Fortran也是一种风尚。很多科学界的前辈，往往就是从这门语言开始了对编程的学习。 如其字面意，本手册是一本简明扼要学习Fortran...
Fortran：数组函数详解
chd_lkl的博客
08-29 1万+
声明：本博文翻译自：https://www.tutorialspoint.com/fortran/location.htm等相关网页 上一篇博文讲了数组的基本用法。这一篇博文着重讲解关于数组内置函数的用法。 1. 向量与矩阵乘法 向量的乘法其实就是数学里的内积运算；而矩阵乘法matmul是指数学中的矩阵乘法运算。而矩阵(向量)A*B则是两个大小相同的矩阵对应位置元素相乘。 dot_produc...
Matlab中sum函数使用
SeisTang的博客
11-18 2222
今天使用matlab中sum函数，记录一下: sum(A)函数主要功能为针对A的元素求和。 1 如果A为矩阵(matrix)： 例如： A=[1 2 3; ... 4 9 6;] ; （1）sum(A,1) 为沿着最快维fast dimension，第一维，进行求和，针对Matlab的数据存储方式，其与Fortran类似，相当于是以矩阵的列为最快维进行存储，此时 sum(A,...
Fortran内部函数错误处理教程
2. Fortran内部函数概念：Fortran语言中的内部函数指的是那些内置在Fortran编译器中的函数，可以直接使用而无需用户自定义。这些函数覆盖了数学运算、字符处理、数组操作等多个方面。了解这些内部函数对于编写高效、...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值