Intel汇编-部分余数

# An example of using the FPREM1 instruction.
#
# 部分余数
#
# 部分余数是浮点数除法中难于处理的一个部分，它的概念与如何执行浮点除法有关。
# 除法操作的余数由被除数对除数的一系列减法决定。在每次减法迭代时，中间余数称为部分余数(partial remainder)
# 当部分余数小于除数时(不能再执行减法操作，否则就会产生负数)，迭代停止。
# 在除法操作结束时，最终的答案是代表减法迭代次数的整数值(称为商(quotient))，
# 以及表示最终的部分余数的浮点值(称为余数(remainder))
#
# 根据执行除法需要多少次迭代，可能有很多部分余数。迭代的次数取决于被除数和除数的指数值之间的差值。
# 每次减法不能使被除数的指数值的减少量超过63。
#
# FPREM 和 FPREM1 指令都计算浮点除法的余数值。
#
# 确定除法余数的基本方法是确定被除数和除数的除法的浮点商，然后把这个值舍入到最近的整数。
# 那么，余数就是除数和商相乘的结果与被除数之间的差值。
#
# 例:
# 20.65 / 3.97 = 5.201511335, 舍入到5(这是商)
# 5 * 3.97 = 19.85
# 20.65 - 19.85 = 0.8 (这是余数)
#
# 困难的部分在于舍入过程。
#

    .section .data
value1:
    .float 20.65
value2:
    .float 3.97

    .section .bss
.lcomm result, 4

    .section .text
    .globl main
main:
    nop
    finit
    flds value2
    flds value1
loop:
    fprem1
    fstsw %ax
    testb $4, %ah
    jnz loop

    fsts result
    movl $1, %eax
    movl $0, %ebx
    int $0x80

# gcc -g -o 07 07-partial-remainder.s -m32
#
# st0            0.799999475479125976562 (raw 0x3ffeccccc40000000000)
# st1            3.97000002861022949219 (raw 0x4000fe147b0000000000)
# st2            0                   (raw 0x00000000000000000000)
# st3            0                   (raw 0x00000000000000000000)
# st4            0                   (raw 0x00000000000000000000)
# st5            0                   (raw 0x00000000000000000000)
# st6            0                   (raw 0x00000000000000000000)
# st7            0                   (raw 0x00000000000000000000)