# 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)