MIPS32的DIV和DIVU实现（除法指令）

《自己动手写CPU》一书的7.11节到7.14节实现了DIV和DIVU指令。

书中通过“试商法”实现除法，并在原有的流水线结构之外另加了状态机进行计算。

照抄书上的实现方法需要另外添加个.v，我实在有点懒，不想在运算指令实现上再另外加个文件。

而且按照书上的实现，DIV和DIVU指令跟其他的运算指令差别很大，不符合我对结构的审美。

那么就必须想想办法自己写一个别的实现方法。

 

于是找到了一个结构化比较好的实现方法：

http://www.cnblogs.com/AndyJee/p/4575152.html

 

参照这个方法，使用verilog实现需要解决2个问题：

 

（1）文章是无符号数的除法，对应DIVU，而DIV是有符号数除法。

这里又要分两种情况讨论。

a、DIV指令的被除数和除数是有符号数

有符号数运算时，输入的被除数和除数都是补码，做除法时需要先将补码转成原码。

正数的补码就是原码，那么只需要在负数时转成原码即可。

b、DIV指令的商和余数正负符号判断

按照初等数学的定义，商的正负取决于被除数和除数的正负，而余数取决于除数的符号。

而计算机语言中，商的正负没有变化，余数符号取决于被除数。

该问题的讨论：

http://blog.sina.com.cn/s/blog_956c84bb0101j31f.html

DIV指令下，被除数/除数为负数，则转换成原码：

　　assign div_dividend = ((aluop_i == `EXE_DIV_OP) && (reg1_i[31] == 1'b1)) ? (~reg1_i + 1) : reg1_i;
　　assign div_divisor = ((aluop_i == `EXE_DIV_OP) && (reg2_i[31] == 1'b1)) ? (~reg2_i + 1) : reg2_i;

DIV指令下，将商和余数转回为补码：

　　assign div_quotient = ((aluop_i == `EXE_DIV_OP) && (reg1_i[31] ^ reg2_i[31] == 1'b1)) ? {1'b1, (~div_quo_o[30:0] + 1)} : div_quo_o;
　　assign div_remainder = ((aluop_i == `EXE_DIV_OP) && (reg1_i[31] == 1'b1)) ? {1'b1, (~div_rem_o[30:0] + 1)} : div_rem_o;

 

 

 

（2）参考的方法通过高级语言描述，需要转化成适合verilog实现的方式。

该方法的核心是通过减法来实现除法，使用中文通俗描述的步骤如下：

基本方法

a、比较被除数和除数的大小；

b、如果被除数比较大，将被除数减去除数，商加1；

c、循环步骤a和b，直到被除数比除数小，这时候被除数剩下的就是余数，商就是结果（初值为0）。

基本方法比较“笨”，效率是最低的。

加速方法

a、比较被除数和n*除数的大小；

b、如果被除数比较大，将被除数减去n*除数，商加n；

c、循环步骤a和b，直到被除数比除数小，这时候被除数剩下的就是余数，商就是结果（初值为0）。

 

加速方法适合使用移位来实现（verilog中，左移1位数值乘于2）。

但是上面文章中描述的方法不适合直接用verilog实现（2层while循环，verilog实现很麻烦）。

可以通过以下几个小技巧避免在verilog中出现类似2层while循环的结构：

a、考虑DIV和DIVU指令的除数位数有限，那么可以先判断最大需要移位几次进行加速方法步骤a（数一数除数第一个1之前有几个0）；

b、在最大移位范围内，按照移位次数从大到小的顺序进行操作（使n*除数中的n尽量大），以尽快收敛结果。

最后实现时，在流水线的寄存器之间进行的运算如下（核心代码）：

                //core method
                if (div_rem_i < (div_divisor << div_shift_cnt_i))
                begin
                    div_quo_o = div_quo_i;
                    div_rem_o = div_rem_i;
                end
                else
                begin
                    div_quo_o = div_quo_i + (1 << div_shift_cnt_i);
                    div_rem_o = div_rem_i + (~(div_divisor << div_shift_cnt_i)) + 1;
                end
            end

这样DIV和DIVU运算所需的周期受除数影响，除数越大则运算周期越短，除数越小则运算周期越长（移位操作次数多）。

例如，当除数是1时，需要31次移位操作（31个时钟）；除数是10'b10_0000_0000时，需要22次移位（22个时钟）。

这种实现方法最大运算周期为32*时钟周期，而《自己动手写CPU》中给的例子对于32位的除法，至少需要32个时钟周期。

另外，这种实现方法与累乘加和累乘减（MADD/MADDU/MSUB/MSUBU）的结构类似，不需要另外再添加.v文件，符合懒人“初始结构以外，最好不做任何结构变动”的习惯。

 

最后，这种方法可以简单实现加速，即通过比较被除数和除数从bit31开始0的数量，可以在被除数和除数绝对值都较小的情况下节省很多运算周期。