C 这个MIPS程序需要这么长时间吗？（近似pi的程序）

其目的是编写一个使用该算法的MIPS程序。 是的。它可以工作，只需大约30分钟即可完成delta=1E-2。 C程序（使用gcc编译）使用该增量大约需要1分半钟。 我在C程序上尝试了delta=1E-3，但在超过2小时后不得不中止

我只是想知道：这应该发生吗？结果对我来说足够精确（3.13909200，δ=1E-2）。我做错什么了吗

我知道这个算法可能不是最有效的，MIPS或MARS（我用于MIPS）也不是

MIPS代码：

    .data

l_cubo:     .double     1.0
delta:      .double     1E-2
zero:       .double     0.0
dois:       .double     2.0
six:        .double     6.0


    .text
    .globl main

main:
    la  $a0,l_cubo
    l.d $f20,0($a0) #l_cubo
    la  $a0,dois

    l.d $f4,0($a0)
    div.d   $f22,$f20,$f4   #r_esfera
    la  $a0,delta
    l.d $f24,0($a0) #delta
    la  $a0,zero
    l.d     $f26,0($a0)     #v_cubo ou v_total
    l.d $f28,0($a0)     #v_esfera

    l.d $f4,0($a0)      #x
    l.d $f6,0($a0)      #y
    l.d $f8,0($a0)      #z

loop_x:

    c.lt.d  $f4,$f20
    bc1f    end_loop_x
    l.d $f6,0($a0)
    loop_y:

        c.lt.d  $f6,$f20
        bc1f    end_loop_y
        l.d $f8,0($a0)
        loop_z:

            c.lt.d  $f8,$f20
            bc1f    end_loop_z
            add.d   $f26,$f26,$f24
            mov.d   $f12,$f4
            mov.d   $f14,$f6
            mov.d   $f30,$f8

            jal in_esfera


            l.d $f10,0($a0)

            beqz    $v0,continue

                add.d   $f28,$f28,$f24
            continue:
                add.d   $f8,$f8,$f24
                j loop_z
            end_loop_z:
        add.d   $f6,$f6,$f24
        j loop_y
        end_loop_y:
        add.d   $f4,$f4,$f24
        j loop_x
end_loop_x:

mul.d   $f24,$f24,$f24
mul.d   $f28,$f28,$f24
mul.d   $f26,$f26,$f24

div.d   $f28,$f28,$f26

la  $a0,six
l.d $f10,0($a0)
mul.d   $f28,$f28,$f10

li  $v0,3       #
mov.d   $f12,$f28   #
syscall         # print pi

li $v0,10       #
syscall         #exit

####################################

    .text
    .globl  in_esfera

in_esfera:

    sub.d   $f12,$f12,$f22
    mul.d   $f12,$f12,$f12
    sub.d   $f14,$f14,$f22
    mul.d   $f14,$f14,$f14
    sub.d   $f30,$f30,$f22
    mul.d   $f30,$f30,$f30
    add.d   $f30,$f12,$f30
    add.d   $f30,$f14,$f30

    mul.d   $f16,$f22,$f22

    li $v0,0
    c.le.d  $f30,$f16
    bc1f    continue2
        li $v0,1
    continue2:
        jr  $ra

---

我只是想知道我的教授将如何纠正一个需要30分钟才能执行的程序。

我假设它使用与

C

版本相同的算法。通过测试立方体中点的三维网格，以查看它们是否在球体内，从而近似Pi值。这是一个

O（N^3）

计算，其中

N

是网格每个维度中的单位数（增量）

所以。。。对预计MIPS代码需要很长时间来计算Pi的精确近似值

• 如果l_cubo是4，delta是1/100，那么您应该执行400 x 400 x 400=64000000次迭代。30分钟似乎太长了

• 如果l_cubo是4，delta是1/1000，那么您应该执行4000 x 4000 x 4000=64000000000次迭代

---

但是，如果您想检查它是否正常，那么当在具有相同参数的相同硬件上运行时，您的MIPs代码应该与C实现的速度一样快（如果不是更快的话）。（注意：如果您在MIPS模拟器上运行MIPS代码，您将无法执行此操作。）

我在模拟器上运行MIPS。30分钟似乎太长了。好吧，我确实停了几次来检查它的运行情况，同时我还有很多其他的东西在运行。那我就把这个寄给我的教授。谢谢