Fortran Coder

查看: 8769|回复: 2
打印 上一主题 下一主题

[原创] Fortran、Python、C求余数问题

[复制链接]

3

帖子

1

主题

0

精华

入门

F 币
110 元
贡献
47 点
跳转到指定楼层
楼主
发表于 2020-8-3 10:27:21 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
本帖最后由 watermelon818 于 2020-8-3 10:53 编辑

昨天群里面有同学讨论求余数的计算问题,Fortran中有专门的求余数函数mod和modulo,其中,两个函数说明如下:(参考资料《Fortran95程序设计》 彭国伦 著)
MOD(a, b)
a, b = INTEGER / REAL
return = INTEGER / REAL
功能:计算a/b余数。当参数为浮点数时,返回(a-int(a/b)*b)的值
MODULO(a, b)
a, b = INTEGER / REAL
RETURN = INTEGER / REAL
功能:同样计算a/b的余数,使用和MOD不同的公式来计算。当参数为整数时,返回a-FLOOR(REAL(a)/REAL(b))*b, 参数为浮点数时,返回a - FLOOR(a/b)*b

小弟主要讨论一下对整数的求余方法的计算过程,抛砖引玉,希望有老师们多多批评指正!
之前max同学在群里发送过一个图片,显示了Fortran求余过程中的计算方法,如下图所示:

对于mod函数的计算方法,是非常的清晰明了了。
接下来看看Fortran和Python以及C语言中的取余数过程:
实验平台:win10 x64, Intel Visual Fortran2013学生版, Visual Studio 2013 Community, Python3.8(CPython)
Fortran实验代码:
[Fortran] 纯文本查看 复制代码
    program main
    implicit none
        write(*, *) mod(139, 10)    ! 139 - 10*(int(139 // 10)) => 139 - 10*int(13.9) => 139 - 130 = 9
        write(*, *) mod(139, -10)   ! 139 - (-10)*(int(139 // -10) => 139 + 10*int(-13.9) => 139 + (-130) = 9
        write(*, *) mod(-139, 10)   ! -139 - 10*(int(-139 // 10)) => -139 - 10*int(-13.9) => -139 + 130 = -9
        write(*, *) mod(-139, -10)  ! -139 - (-10)*(int(-139 // -10)) => -139 - (-10)*(int(13.9)) => -139 + 130 = -9
    end program

运行结果:
[Fortran] 纯文本查看 复制代码
           9
           9
          -9
          -9
请按任意键继续. . .
依据代码的计算过程,我们可以看出,Fortran在取余的计算过程中使用int()进行取整数时,实现的是正数负数均向0取整的计算方法。
接下来我们继续看一下python3.8中对于取余数的计算方法:
Python实验代码:
[Python] 纯文本查看 复制代码
print(139 % 10)   # 139 - 10*(int(139 // 10)) => 139 - 10*int(13.9) => 139 - 130 = 9
print(139 % -10)  # 139 - (-10)*(int(139 // -10) => 139 + 10*int(-13.9) => 139 + (-140) = -1
print(-139 % 10)  # -139 - 10*(int(-139 // 10)) => -139 - 10*int(-13.9) => -139 + 140 = 1
print(-139 % -10) # -139 - (-10)*(int(-139 // -10)) => -139 - (-10)*(int(13.9)) => -139 + 130 = -9

运行结果:
[Python] 纯文本查看 复制代码
Python 3.8.1 (tags/v3.8.1:1b293b6, Dec 18 2019, 23:11:46) [MSC v.1916 64 bit (AMD64)] on win32
Type "help", "copyright", "credits" or "license()" for more information.
>>> 
================ RESTART: C:/Users/Administrator/Desktop/fuck.py ===============
9
-1
1
-9
>>> 

依据代码的计算过程,我们可以看出,Python在取余的计算过程中使用int()进行取整数时,计算方法为:向下取整;即正数向0取整,负数向负无穷取整。
C语言代码与Fortran代码计算结果相同
实验程序如下所示:
[C] 纯文本查看 复制代码
#include <stdio.h>
int main(int argc, char *argv[])
{
        printf("%d\n", 139 % 10);
        printf("%d\n", 139 % -10);
        printf("%d\n", -139 % 10);
        printf("%d\n", -139 % -10);
        return 0;
}

运行结果:
[C] 纯文本查看 复制代码
9
9
-9
-9
请按任意键继续. . .

所以C语言在求余数的方法上与Fortran是相同的。
需要注意的地方是:Fortran中使用int( ),C语言中使用(int)( ),Python中使用int( )对正数与负数取整时,均采用的是向0取整;很明显Python在取余过程中,其计算方法在负数方面发生了改变,成为了向负无穷取整。

希望大家多多批评指正!

分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩
菜鸟一枚,请多多指教!

178

帖子

15

主题

0

精华

大宗师

F 币
4973 元
贡献
1152 点
沙发
发表于 2020-8-3 11:00:46 | 只看该作者
python脚本名字是亮点……

3

帖子

1

主题

0

精华

入门

F 币
110 元
贡献
47 点
板凳
 楼主| 发表于 2020-8-3 16:22:40 | 只看该作者
liudy02 发表于 2020-8-3 11:00
python脚本名字是亮点……

可以说阅读相当仔细了!
菜鸟一枚,请多多指教!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-4-24 12:48

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表