ブースの乗算アルゴリズム

ブースの乗算アルゴリズムは、2の補数表現のふたつの符号付整数の乗算の手法である。

このアルゴリズムは1950年ごろアンドリュー・ドナルド・ブースがロンドン大学バークベック・カレッジで結晶学を研究しているときに発明したものである。ブースは卓上計算機を使って研究していて、計算速度を向上させるために乗算を高速化する方法を探していてこれを発明した。彼の時代のマシンではシフトは加算よりも高速であり、ある種の数値では彼のアルゴリズムは高速であった。しかも負数についてもこのアルゴリズムは機能した。ブースのアルゴリズムはコンピュータ・アーキテクチャの研究において興味深い。

1. アルゴリズム
---------

ブースのアルゴリズムは、符号付きの2の補数表現の//N//ビットの乗数 //Y// において、最下位ビットよりさらに下に //y//-1 = 0 というビットを暗黙のうちに補って隣接する2つのビットを調べる。//Y//の各ビット //y////i//について、//y////i// と //y////i//-1 を調べる。それら2ビットが同じ場合、積 //P// は変化しない。//y////i// = 0 かつ //y////i//-1 = 1 の場合、被乗数に 2//i// をかけたものを //P// に加算する。//y//i = 1 かつ //y//i-1 = 0 の場合、被乗数に 2//i// をかけたものを //P// から減算する。このようにして得られた //P// の最終的な値が符号付きの積となっている。

被乗数と積の表現は特に指定されていない。一般にそれらも乗数と同様に2の補数表現とするが、加減算が可能な任意の数値表現が可能である。また、乗数をどちらから調べるかも規定されていない。通常はLSBからMSBへと進め、//i// = 0 からスタートする。この場合 2//i// をかける演算は //P// をステップごとに1ビットずつシフトする操作に置き換え可能となる。シフトではみ出たビットは捨てられ、その後の加減算は//P//の上位//N//ビットについてのみ行う。実際には様々な派生や最適化が存在する。

2. 典型的実装
--------

ブースのアルゴリズムは、2つの事前に決定された値 //A// および //S// のうちの1つを積 //P// に繰り返し加算し、//P// を右方向に算術シフトするという形で実装できる。ここで、m を被乗数、r を乗数とし、m のビット数を //x//、rのビット数を //y// とする。

//A// と //S// の値を決定し、//P// を初期状態にする。いずれも長さは  ビットとする。
A: 最上位側のビット列に m の値を格納する。残る  ビットは全て0にする。
S: 最上位側のビット列に  の2の補数表現を格納する。残る  ビットは全て0にする。
P: 最上位側の //x// ビットを全て0にする。その右のビット列に r の値を格納する。最下位ビットは0とする。
//P// の最下位2ビットを調べる。
その中身が 01 の場合、//P//   +   //A// を計算し、オーバーフローは無視する。
その中身が 10 の場合、//P//   +   //S// を計算し、オーバーフローは無視する。
その中身が 00 の場合、何もしない。現在の //P// をそのまま次ステップで使用する。
その中身が 11 の場合、何もしない。現在の //P// をそのまま次ステップで使用する。
上のステップで得られた値を右に1ビット算術シフトする。その値を //P// の新たな値とする。
ステップ2と3を //y// 回反復する。
//P//の最下位ビットを除去する。これが m と r の積である。

3. 実施例
------

3 × -4 を実行する。m = 3、r = −4、//x// = 4、//y// = 4 である。

m = 0011, -m = 1101, r = 1100
A = 0011 0000 0
S = 1101 0000 0
P = 0000 1100 0
ループを 4回実行する。
P = 0000 1100 0. 最後の 2ビットは 00.
P = 0000 0110 0. 右シフト.
P = 0000 0110 0. 最後の 2ビットは 00.
P = 0000 0011 0. 右シフト.
P = 0000 0011 0. 最後の 2ビットは 10.
P = 1101 0011 0. P = P + S.
P = 1110 1001 1. 右シフト.
P = 1110 1001 1. 最後の 2ビットは 11.
P = 1111 0100 1. 右シフト.
積は 1111 0100 となり、 -12 である。2の補数表現では、表現可能な最小の負数の2の補数を計算すると本来のビット数では表現できない。例えば、被乗数が4ビットなら -8  という値の2の補数は 8  となる。つまり A と S で必要とするビット数が異なることになり、上述の実装ではうまく計算できない。1つの解決策として、A、S、P のビット数を1ビット増やせばよい。以下は被乗数も乗数も4ビットのとき −8 × 2 を計算する例である。

A = 1 1000 0000 0
S = 0 1000 0000 0
P = 0 0000 0010 0
ループを 4回実行する。
P = 0 0000 0010 0. 最後の 2ビットは 00.
P = 0 0000 0001 0. 右シフト.
P = 0 0000 0001 0. 最後の 2ビットは 10.
P = 0 1000 0001 0. P = P + S.
P = 0 0100 0000 1. 右シフト.
P = 0 0100 0000 1. 最後の 2ビットは 01.
P = 1 1100 0000 1. P = P + A.
P = 1 1110 0000 0. 右シフト.
P = 1 1110 0000 0. 最後の 2ビットは 00.
P = 1 1111 0000 0. 右シフト.
積は 11110000 となり、−16 である。

4. 動作原理
-------

1の前後を 0で囲んだ正の乗数を考える。例えば、00111110 である。積は以下のようになる。

ここで M は被乗数である。演算回数は以下のように書き換えることで 2回にまで減らすことができる。

積は被乗数の 1回の加算と 1回の減算で生成される。この手法を 1 のかたまり毎に実施することでいかなる乗数にも適用できる。したがって以下のようになる。

ブースのアルゴリズムはこの手法に従って、1のかたまりの最初に当たったら加算を実行し、1のかたまりの最後に当たったら減算を実行するものである。この手法は負数の乗算でもうまく機能する。乗数の中の 1 が連続していれば、ブースのアルゴリズムでは一般的な乗算アルゴリズムよりも加算や減算が少なくて済む。逆に 1 が2つ以上連続していない場合、一般的なアルゴリズムよりも加算や減算の回数が多くなってしまう。

5. 脚注
-----

6. 参考文献
-------

Andrew D. Booth. [1]. The Quarterly Journal of Mechanics and Applied Mathematics, Volume IV, Pt. 2
Collin, Andrew. [2]. //Resurrection//, Issue 5, Spring 1993. London: Computer Conservation Society.
Patterson, David and John Hennessy. //Computer Organization and Design: The Hardware/Software Interface, Second Edition//. ISBN 1-55860-428-6. San Francisco, California: Morgan Kaufmann Publishers. 1998.
Stallings, William. //Computer Organization and Architecture: Designing for performance, Fifth Edition//. ISBN 0-13-081294-3. New Jersey: Prentice-Hall, Inc.. 2000.

7. 関連項目
-------

乗算器

8. 外部リンク
--------

[3]
[4] in [5]
[6]
[7]
[8]
[9]
[10]