只是想知道 MIPS 中的一个函数，它可以在我的 MIPS 程序中以任何数字为根。
【问题讨论】：
pow(x, 1.0/3)
你可以使用简化版的牛顿法求整数的根
x=N
iterate 20 times:
x'=(x+N/x) /2
x=x'
Mips 实现
.data
.text
main:
li $t0,25 #N
move $t1,$t0 #x
li $t4,0 #loop variable
sqrLoop:
#Newton Formula
div $t3, $t0, $t1 # N/x
add $t1, $t3, $t1 # x + N/x
div $t1, $t1, 2 # (x + N/x)/2
#loop
add $t4, $t4, 1
blt $t4, 20, sqrLoop
end:
li $v0,10
syscall
【讨论】：
srl 除以 1 来除以 2，而不是另一个慢速 div。
add $t4, $t4, -1 / bne $zero, $t4 sqrLoop 以将 $t4 向下 计数为零。 （两个寄存器之间的 BNE 是真正的硬件指令；针对立即数的 BNE 只是伪指令，必须先将 20 放入寄存器中）。
试试这个。
#DATA
.data
square: .asciiz "Enter the number you wish to find the square root for: "
answer: .asciiz "The answer is: "
newline: .asciiz "\n"
#Text
.text
.globl main
main:
li $v0, 4 #Prompt user for input
la $a0, square
syscall
li $v0, 5 #Receive said input
syscall
move $a0, $v0
move $t4, $zero #Move variables to t registers
move $t1, $a0
addi $t0, $zero, 1 #Set $t0 to 1
sll $t0, $t0, 30 #Bit Shift $t0 left by 30
#For loop
loop1:
slt $t2, $t1, $t0
beq $t2, $zero, loop2
nop
srl $t0, $t0, 2 #Shift $t0 right by 2
j loop1
loop2:
beq $t0, $zero, return
nop
add $t3, $t4, $t0 #if $t0 != zero add t0 and t4 into t3
slt $t2, $t1, $t3
beq $t2, $zero, else1
nop
srl $t4, $t4, 1 #shift $t4 right by 1
j loopEnd
else1:
sub $t1, $t1, $t3 #Decrement $t1 by $t3
srl $t4, $t4, 1 #Shift $t4 right by 1
add $t4, $t4, $t0 #then add $t0 to that
loopEnd:
srl $t0, $t0, 2 #shift $t0 to the right
j loop2
return:
li $v0, 4 #print out the answer then exit
la $a0, answer
syscall
li $v0, 1
move $a0, $t4
syscall
li $v0, 10
syscall
【讨论】：
b 不是j 指令之后？）也许只是翻译相同的伪代码/C 的巧合。
.数据
.文本
.globl 主要
.ent main
平方：
移动 $v1, $a1 # $v0 = x = N li $t0, 0 # 计数器
sqr循环：
div $t8, $a1, $v1 # N/x 添加 $v1, $t8, $v1 # x + N/x div $v1, $v1, 2 # (x + N/x)/2 添加 $t0, $t0, 1 blt $t0, 20, sqrLoop
jr $ra
.end Sqrt
这是一个函数，调用时会取出数字的平方根。只需在必要或需要的地方编写“jal Sqrt”即可调用该函数。
【讨论】：
你可以试试这个算法，它给出的整数小于或等于你的数字的平方根。
假设你想要 n 的平方根。然后不断重复以下计算：
x = (x + n/x) / 2
选择 x = n 开始并不断重复直到 x 停止变化。
这是您可以在程序中添加的以下 MIPS 程序库。
#SquareRoot.s
#DATA
.data
square: .asciiz "Enter the number you wish to find the square root for: "
answer: .asciiz "The answer is: "
newline: .asciiz "\n"
#Text
.text
.globl main
main:
li $v0, 4 #Prompt user for input
la $a0, square
syscall
li $v0, 5 #Receive said input
syscall
move $a0, $v0
move $t4, $zero #Move variables to t registers
move $t1, $a0
addi $t0, $zero, 1 #Set $t0 to 1
sll $t0, $t0, 30 #Bit Shift $t0 left by 30
#For loop
loop1:
slt $t2, $t1, $t0
beq $t2, $zero, loop2
nop
srl $t0, $t0, 2 #Shift $t0 right by 2
j loop1
loop2:
beq $t0, $zero, return
nop
add $t3, $t4, $t0 #if $t0 != zero add t0 and t4 into t3
slt $t2, $t1, $t3
beq $t2, $zero, else1
nop
srl $t4, $t4, 1 #shift $t4 right by 1
j loopEnd
else1:
sub $t1, $t1, $t3 #Decrement $t1 by $t3
srl $t4, $t4, 1 #Shift $t4 right by 1
add $t4, $t4, $t0 #then add $t0 to that
loopEnd:
srl $t0, $t0, 2 #shift $t0 to the right
j loop2
return:
li $v0, 4 #print out the answer then exit
la $a0, answer
syscall
li $v0, 1
move $a0, $t4
syscall
li $v0, 10
syscall
【讨论】：
x = (x + n/x) / 2 Newton 迭代不相同的二进制逐位算法。 （算法看起来像Fastest Integer Square Root in the least amount of instructions 中的算法之一）。所以它可能是抄袭的，并且错误地描述了它实现的算法。