JPS58142441A - High radix division system - Google Patents

High radix division system

Info

Publication number
JPS58142441A
JPS58142441A JP2500782A JP2500782A JPS58142441A JP S58142441 A JPS58142441 A JP S58142441A JP 2500782 A JP2500782 A JP 2500782A JP 2500782 A JP2500782 A JP 2500782A JP S58142441 A JPS58142441 A JP S58142441A
Authority
JP
Japan
Prior art keywords
divisor
division
circuit
digit
register
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2500782A
Other languages
Japanese (ja)
Inventor
Sadao Nakamura
中村 定雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Tokyo Shibaura Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Tokyo Shibaura Electric Co Ltd filed Critical Toshiba Corp
Priority to JP2500782A priority Critical patent/JPS58142441A/en
Publication of JPS58142441A publication Critical patent/JPS58142441A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/38Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
    • G06F7/48Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
    • G06F7/52Multiplying; Dividing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)

Abstract

PURPOSE:To realize a division circuit of a high radix with a simple hardware constitution and to attain high speed division, by obtaining a quotinent of a divisor and a dividend obtained from a divisor obtained newly and a partial residual, by one digit each from the upper-order digit sequentially. CONSTITUTION:The upper-order of a divisor (d) on an input line 20 is given to a reciprocal table 22. The (d) is set to a register 25, an output of the table 22, i.e., an approximate value of the reciprocal of the (d) is set to a register 26, the multiplication is executed by taking the content of the registers 25, 26 as dividend and divider respectively and the result is set to pipeline registers 39, 40. The upper-order bits of the registers 39, 40 are summed at a pre-addition circuit 41 to obtain an approximate value of a partial residual pj. A quotinent digit qi+1 is obtained from this approximate value at a quotinent digit generating circuit 42, the product between this value qj+1 and the content of the register 25 is calculated, and a new partial residual pj+1 is obtained from the content of the registers 39, 40 and the partial residual pj. This residual is set to the registers 39, 40 for the operation.

Description

【発明の詳細な説明】 〔発明の技術分野〕 本発明はディジタル除算回路の1つの方式である高基数
除算方式の改良に関する。
DETAILED DESCRIPTION OF THE INVENTION [Technical Field of the Invention] The present invention relates to an improvement in a high radix division method, which is one method of a digital division circuit.

〔発明の技術的背景とその問題点〕[Technical background of the invention and its problems]

従来知られている除算方式には大別して2つの方式があ
る。その1つ祉部分剰余法であシ、他の1つは乗算を用
いた方法(以下乗算法と呼ぶ)である。部分剰余法の最
も簡単なものがよく知られた回復法(r@stor1m
g m@thod )と非回復法(nonrsstor
ing m@thod )である。この回復法によれば
与えられた除数、被除数の商をその上位桁から順に1ピ
ツトづつ求められる。ところで一般の部分剰余法は回復
法や非回復法のアルがリズムを発展させて高基数の除算
を扱えるようにし喪ものである。回復法や非回復法によ
って商を1ピツトづつ求める除算回路方式は最も単純で
あるが、当然その演算実行に多くのステップ数と時間を
要する。高基数法、すなわち一般の部分剰余法によれば
前記回復法や非回復法の場合と同じだけの精度で商を求
めようとするとき、上記回復法や非回復と較べて演算実
行ステップ数が減少し、1ステッグ当りの実行時間は多
少増加するものの全体としての実行時間の減少を図るこ
とが可能である。尚、部分剰余法についての詳しい説明
は後で行うが、一般に部分剰余法による除算回路は次に
述べる乗算法と較べて、その構造及び制御は単純であり
Conventionally known division methods can be broadly classified into two types. One of them is the welfare partial remainder method, and the other is a method using multiplication (hereinafter referred to as the multiplication method). The simplest partial remainder method is the well-known recovery method (r@stor1m
g m@thod ) and non-recovery method (nonrsstor
ing m@thod). According to this recovery method, the quotient of the given divisor and dividend can be found one pit at a time starting from the most significant digit. By the way, the general partial remainder method is different from restorative and non-restorative methods in that it develops a rhythm that allows it to handle divisions with high cardinal numbers. The division circuit method that calculates the quotient one pit at a time by the recovery method or the non-recovery method is the simplest, but it naturally requires a large number of steps and time to execute the calculation. According to the high-radix method, that is, the general partial remainder method, when trying to obtain a quotient with the same precision as the above-mentioned recovery method or non-recovery method, the number of calculation execution steps is smaller than in the above-mentioned recovery method or non-recovery method. Although the execution time per step increases somewhat, it is possible to reduce the overall execution time. A detailed explanation of the partial remainder method will be given later, but in general, a division circuit using the partial remainder method has a simpler structure and control than the multiplication method described below.

壇九基数をあiり大きくとらなければ回路規模も小さく
することが可能であるという利点がある。その反面部分
剰余法による除算回路は実行時間が遅く、九とえ高基数
法によってもハードウェアをあまシ大きくしない範囲で
は高々lステツノ肖り4〜5ビツトまでの除算しか望め
ず、高速化は困難であるという欠点を有する。
There is an advantage that the circuit scale can be made smaller unless the number of bases is too large. On the other hand, the execution time of division circuits based on the partial remainder method is slow, and even if the nine-high radix method is used, as long as the hardware is not made too large, it is only possible to perform divisions of up to 4 to 5 bits, and speed-up is difficult. It has the disadvantage of being difficult.

さて乗算法による除算法の概略は、除数dを分母、被除
数p、を分子としたとき、上記分母、および分子に次々
にある数を乗算することによって上記分母t−1に収束
させ、同時にその結果として上記分子に求めようとして
いる商qを得るものである。この関係は上述した次々に
乗する数七ki (’ =O−1s  2 =・、n)
とすると次式のように表わせる。
Now, the outline of the division method using the multiplication method is, when the divisor d is the denominator and the dividend p is the numerator, the denominator and numerator are multiplied one after another by certain numbers to converge to the denominator t-1, and at the same time, As a result, we obtain the quotient q that we are trying to find for the above numerator. This relationship is expressed as the number 7 ki (' = O-1s 2 = ·, n) which is multiplied one after another as mentioned above.
Then, it can be expressed as the following formula.

この第(1)式においてdkok、・・・kt→lなら
ばp、 1cllk、・・・kn−+Qとして商Qが求
められる。
In this equation (1), if dkok, . . . kt→l, the quotient Q is obtained as p, 1cllk, .

実際には以下に示す手順によって除数d、被(+)  
aの逆数の近似値に、を論理回路又は読み出し専用メモ
リ(以下ROMと書く)によって実現され友逆数テーブ
ルより求める。
Actually, the divisor d, divisor (+) are
An approximate value of the reciprocal of a is obtained from a reciprocal table realized by a logic circuit or a read-only memory (hereinafter referred to as ROM).

このような手順によれば適当な反復回数nの後、ある精
度でもって dn+1:’;1で6 ’) 、Pn+ 
1: Qが求まる。この手順は上記第(1)弐において
kt=(2−aj) (j= 1.2、・・・)と選ん
だことに相当する。ここに示された乗算法による除算方
式は更にいくつかの重要な工夫を加え友上で、実際にハ
ードウェアとして効率的に実現されている。乗算法によ
る除算アルプリズムのハードウェアによる実現例拡、例
えば A山rion、 8.F、山1.          
L“The IBM 8ysleW4/360 Mod
@l 91 : Floatlng−”Po1nt E
xe@ution Unit ”IBM Journa
l、January 1967、PP 34−53゜に
詳しく述べられている。
According to such a procedure, after a suitable number of iterations n, with a certain precision dn+1:'; 1 in 6'), Pn+
1: Find Q. This procedure corresponds to selecting kt=(2-aj) (j=1.2, . . . ) in the second part of (1) above. The multiplication-based division method shown here has been effectively implemented in hardware with several important improvements. Expansion of examples of hardware realization of division al prism using multiplication method, e.g. A-mount rion, 8. F, mountain 1.
L“The IBM 8ysleW4/360 Mod
@l 91: Floatlng-”Po1nt E
xe@ution Unit ”IBM Journal
1, January 1967, PP 34-53°.

しかして上記乗算法による除算方式には以下に示すよう
な利点がある。
However, the division method using the multiplication method described above has the following advantages.

1、)2次収束アルゴリズムであること。即ちlステラ
!ごとに解の精度が2倍づつ向上する。これKよプ少な
いステップ数で高精度の除算が実現できる。
1.) It must be a quadratic convergence algorithm. In other words, Stella! The accuracy of the solution improves by a factor of two for each step. High-precision division can be achieved with a number of steps less than K.

2、)いくつかのハードウェアの工夫及び上の1)の性
質により、非常に高速な除算が可能である。
2.) Very high-speed division is possible due to some hardware innovations and the properties of 1) above.

しかしながらこの除算方式は乗算を基本とする九め、高
速な乗算回路が必要であり、更にいくつかのハードウェ
アの工夫が必要な仁となどによシ、全体の回路規模が大
きくなる欠点がある。この欠点は・・−ドウエアの量を
押さえる必要があるときには重要な問題となる。
However, this division method is based on multiplication, requires a high-speed multiplication circuit, requires some hardware innovation, and has the drawback of increasing the overall circuit scale. . This drawback becomes an important problem when it is necessary to reduce the amount of doware.

〔発明の目的〕[Purpose of the invention]

本発明は上述した従来の除算方式の欠点Kliみてなさ
れたもので、その目的とするとζろは、演算実行のlス
デッグ当シ多くのビット数の商ディジツトを比較的短い
時間で求めることができ、且つ全体として高速な除算回
路を小規模なる。
The present invention was made in view of the above-mentioned drawbacks of the conventional division method, and its purpose is that it is possible to obtain a quotient digit of a large number of bits in a relatively short period of time during calculation execution. , and the overall high-speed division circuit is small-scale.

〔発明の概要〕[Summary of the invention]

本発明は与えられた除数の上位数ビットから選択した適
当な数を上記除数及び与えられた被除数にそれぞれ乗じ
て除数の取シ得る値を狭い範囲に限定した新しい除数と
新しい被除数とを求め、この新しく求められ九被除数を
部分剰余の初期値として、前記新しく求められた除数と
上記部分剰余とを用いて高基数除算の演算を行わしめる
ようにしたものである。
The present invention calculates a new divisor and a new dividend that limit the possible values of the divisor to a narrow range by multiplying the divisor and the given dividend by an appropriate number selected from the upper few bits of the given divisor, respectively. This newly obtained 9 dividend is used as the initial value of the partial remainder, and the newly obtained divisor and the partial remainder are used to perform a high radix division operation.

〔発明の効果〕〔Effect of the invention〕

従って、本発明によれば演算−行に用いられる除数の取
り得る値が狭い範囲に限定されるので、除数および、部
分剰余からの商ディジツトの選択が容易となシ、上記演
算を短時間で実行することか可能となる。しかも、土遥
のように除数の取妙得る値が限定されるので、基数を大
きくしても簡単で小規模なハードウェアによって除算回
路を構成することができる等の効果を奏する。
Therefore, according to the present invention, the possible values of the divisor used in the operation row are limited to a narrow range, so the selection of the divisor and the quotient digit from the partial remainder is easy, and the above operation can be performed in a short time. It is possible to execute it. Moreover, since the possible values of the divisor are limited as in the case of Tsuchiharu, even if the radix is increased, the division circuit can be constructed using simple and small-scale hardware.

〔発明の実施例〕[Embodiments of the invention]

以下、図面を参照して本発明の詳細な説明明する。 Hereinafter, the present invention will be explained in detail with reference to the drawings.

先ず本発明に係る高基数除算方式を明らかにする為に、
部分剰余法のアルゴリズム及び従来の高基数除算回路に
ついて説明する。尚、説明の便宜上、除数dの値の範囲
を次のように仮定する。
First, in order to clarify the high radix division method according to the present invention,
The algorithm of the partial remainder method and the conventional high-radix division circuit will be explained. For convenience of explanation, the value range of the divisor d is assumed as follows.

+≦d<1 部分剰余法においては、除@aと被除数Poの商Qa次
の反復演算によって計算される。
+≦d<1 In the partial remainder method, the quotient Qa of the division @a and the dividend Po is calculated by the following iterative operation.

れる。但し、上式中rは基数である。第(2)式におけ
るセレクト関数S@l@ct(Pj,d)は除数dと第
j部分剰余Pjとから、ある規則によって第(j十+)
商ディジツトqj+1  を求める一のである。この商
ディジツトq1+1  をどのようにして求めるかにつ
いては後ヤ説明する。第(2)式に示される関係から簡
単表導出によって次の第(3)式を求めることができる
It will be done. However, in the above formula, r is a radix. The selection function S@l@ct(Pj, d) in equation (2) selects the (j10+)th by a certain rule from the divisor d and the jth partial remainder Pj.
The first step is to find the quotient digit qj+1. How to obtain this quotient digit q1+1 will be explained later. The following equation (3) can be obtained by deriving a simple table from the relationship shown in equation (2).

1 M,r−i −!! =−!’ r−”     
− − (3)j=l     d   d この第(3)式において、 lPn/dl < 1            = ・
” (4)が成立する壜らば昼q1r−1  はnディ
ジットのj=1 精度で求められる商Q=PO/d に一致することが分
る。そこで次の条件を設定する。
1 M, r-i-! ! =-! 'r-'
− − (3) j=l d d In this equation (3), lPn/dl < 1 = ・
It can be seen that if (4) holds, then the time q1r-1 corresponds to the quotient Q=PO/d, which is obtained with j=1 precision of n digits.Therefore, the following conditions are set.

この第(5)式が成立するならば任意のnについて  
 “Σq1r−1はnディジットの精度で求められる商
j=1 Q = p./d に一致すると言える。第(5)式に
おいてはlよシ小さいkをも許し九が、実際には、ある
下限よ)も大きければ1よシ小さいkであっても除算を
正しく遂行でき、更にkの値は除算の実行において必要
な商デイジツ}0絶対値の最大値と関係することになる
。このことについては後の詳論で明らかになる。
If this formula (5) holds, then for any n
``Σq1r-1 can be said to match the quotient j=1 Q = p./d, which is found with n-digit precision.In equation (5), we allow k that is smaller than l, but in reality, there is If the lower limit) is also large, division can be performed correctly even if k is smaller than 1, and furthermore, the value of k will be related to the maximum value of the absolute value of the quotient digits 0 that is necessary to perform the division. This will become clear in the detailed discussion later.

条件第(5)式によって、除算が実行される為には初期
条件1p./at<k  が成立しなければならない。
According to condition No. (5), in order for division to be executed, initial conditions 1p. /at<k must hold true.

セレクト関数FiIP5/al ( kの仮定のもとで
lPj+1/if<hの成立を保証するような商ディジ
ツトQj+1を生成するものであれによい。
The selection function FiIP5/al may be any function that generates the quotient digit Qj+1 that guarantees that lPj+1/if<h under the assumption of k.

第1図は与えられ九d 、 Pjに対して選択可能な商
デイジツ} Q3+を及びその商ディジツトを選択し九
ときのPj+1  の値を示しえものである。
FIG. 1 shows a selectable quotient digit Q3+ for a given 9d, Pj, and the value of Pj+1 when the quotient digit is selected.

除算の実行は第1図を参照することによって次のように
説明することかで龜る。
The execution of the division can be explained as follows with reference to FIG.

即ち、与えられ九部分剰余rP5+tの値に対応すゐ商
デイジツ} Q3+tと次の部分剰余Pj+tの値を第
1図よシ読みとると、例えばrP5の値が第1図中A点
にあるなら、その時選択可能な商ディジツトはQ3+,
=1又は2として求められる。
That is, the value of Q3+t and the next partial remainder Pj+t corresponding to the value of the given nine partial remainder rP5+t is read in Figure 1.For example, if the value of rP5 is at point A in Figure 1, then The quotient digits that can be selected at that time are Q3+,
= 1 or 2.

qj+1=2を選択したならp3+.の値は図中B点と
々る。仮夛にrP5の値が第1図中C点にあるなら、そ
の時選択可能な商ディジツトはqj+1”−1だけであ
シ、P5+1の値線図中B点となる。このようにしてq
j+1とP3+1が求壕つ九ならj4−j+1として上
の手順を操シ返して実行する。
If qj+1=2 is selected, p3+. The value reaches point B in the figure. Suppose that the value of rP5 is at point C in Figure 1, then the only quotient digit that can be selected is qj+1''-1, which becomes point B in the value diagram of P5+1.In this way, q
If j+1 and P3+1 are 9, repeat the above procedure as j4-j+1.

上述した除算手続きが正しく実行される為にはに≧十で
なければならない。何故ならば第1図中、商ディジツト
の重複領域の幅tを求めるとz=(2h−1)aである
から、もしk〈+であると、あるrPjの値に対して選
択可能な商ディジツト’j+1が存在しないという事態
を生じる。
In order for the division procedure described above to be executed correctly, it must be greater than or equal to 10. This is because, in Fig. 1, the width t of the overlapping area of the quotient digits is z = (2h-1)a, so if k〈+, the selectable quotient for a certain value of rPj is A situation arises in which digit 'j+1 does not exist.

従ってkは+≦に≦lなゐ条件を満たさなければならな
い。
Therefore, k must satisfy the condition that +≦and≦l.

さて除算を実行する上で必要な商ディジツトの絶対値の
最大値pFiis冨rk(r−+)lによって求まる。
Now, the maximum value of the absolute value of the quotient digits required to execute the division is determined by pFiis(r-+)l.

ζζにraltfaの天井値、即ちa≦鳳表る整数nの
中で最小なものである。そしてすペての商ディジツトは
ディジット集合 (−戸、−jl+1.・・・I  ’Tot’l・・・
、β−1,ρ)の中から選択される。
ζζ is the ceiling value of raltfa, that is, the minimum value among the integers n expressed as a≦O. And all the quotient digits are the digit set (-to, -jl+1...I 'Tot'l...
, β-1, ρ).

尚、上に述べ友ようにkを先に与えてからρを求めるの
でなく、初めにディジット集合を与えておき、k = 
s/(r−+ )によってkの値を求めて4よい。
Note that instead of giving k first and then finding ρ as mentioned above, we give the digit set first and then calculate k =
Find the value of k using s/(r-+).

以上の説明で明らかなように商ディジツトのは最小冗員
の場合であり、部分剰余rP5から商ディジツトqjl
が一意に定まる。またに=−1(ρ=r−+)は最大冗
長の場合であり、すべての部分剰余rPjに対して選択
可能な商デイジッ) qj+lが2通シ存在する。特に
r=2 、に=1のと奮ディジット集合は(t、O,x
)であるが最大冗長であるので商ディジツト0が不必要
になる。即ち 、Qj+t= (−1、l )によって
除算が遂行できる。
As is clear from the above explanation, the quotient digit is the case of minimum redundancy, and the quotient digit qjl is obtained from the partial remainder rP5.
is uniquely determined. Further, =-1 (ρ=r-+) is the case of maximum redundancy, and there are two selectable quotients qj+l for all partial remainders rPj. In particular, the digit set r=2, ni=1 is (t, O, x
), but since it is maximally redundant, the quotient digit 0 is unnecessary. That is, division can be performed by Qj+t=(-1,l).

このような理論に基づいて実行される除算演算が非回復
法に相当する。
A division operation executed based on such a theory corresponds to a non-recovery method.

さて、第1図に示すrPj−d  平面上において町+
1が一定の領域は2つの直線 rPj=(k+qj++)d t rPj= (−に+
Qj++)dに挾まれた領域である。それぞれのqj+
、の値に対してこの領域を図示し九のが第2図に示すP
D線図である。この第2図は特K r = 4、k==
4、ρ=2の場合のPD線図である。第2図中、斜線領
域線部ディジットの選択に冗長性があることを示し、ま
た斜線のない領域は商ディノットが一意に定まる仁とを
示している。例えば、第2図中領域■はqj+l=0と
して一意に定−tC1領域■はqj+1=0又はqj+
l=1のどちらでもよく、更には領域■はqj十1=1
に一意に定まる。
Now, on the rPj-d plane shown in Figure 1, the town +
The area where 1 is constant has two straight lines rPj = (k + qj++) d t rPj = (- to +
This is the area between Qj++)d. each qj+
This region is illustrated for the value of , and the ninth one is P shown in FIG.
It is a D diagram. This figure 2 is especially K r = 4, k = =
4. It is a PD diagram in the case of ρ=2. In FIG. 2, the shaded area indicates that there is redundancy in the selection of the shaded digits, and the shaded area indicates that the quotient dinot is uniquely determined. For example, the region ■ in FIG. 2 is uniquely defined as qj+l=0, and the region
Either l=1 is fine, and furthermore, the area ■ is qj11=1
is uniquely determined.

このようなPD線図は与えられた部分剰余rPjと除数
dとから選択可能な商ディジツトを示すものであり、換
言すればセレクト関数を図示したものとなる。実際に除
算を実行する上では七 飯しクト関数によって商ディジ
ツトが一意に決ま  :べろようにすることが必要であ
る。その為にFiPD線図上の商ディジツトの冗長領域
、即ち#l2FjA中斜線領域に境界を入れなければな
らない。ところでこの境界は斜線領域内に納まってる限
りどんなものでもよいから、除数d及び部分剰余rP4
を2道表示したとき、区切りのいい点(例えば少ないピ
ット数で表現できる地点)K上記境界を定めることがで
きる。r=4、ρ=2、k=にの場合境界の例を第3図
に示す。この例では除数及び部分剰余の上位4ピツトか
ら正しい商ディジツトを求めることができる。この第3
図のように境界を設けたPD線図のハードウェア化は、
例えば読み出し専用メモリ(ROM’)によって実現で
き、このROMをテーブルとして高速に商ディジツトを
求めることが可能となる。このハードウェア化されたP
D線図を以後PDテーブルと呼ぶ。以上の考え方は従来
、既に知られており、例えば下記の文献 Daniel  E、  Atklns“ H1gh@
r−Radlx  Divison  Using  
Istlmatssof  the  Dlvjser
  mad  Partial  R*majmdsr
s”IEICE Transaetions on C
@mput@r@、 vol、c−17゜A 10,0
ctob@r  1968.PP925−934゜に詳
しく述べられている。
Such a PD diagram shows the quotient digits that can be selected from the given partial remainder rPj and divisor d, in other words, it is a diagram illustrating a selection function. When actually performing division, it is necessary to ensure that the quotient digit is uniquely determined by the digit function. For this purpose, a boundary must be placed in the redundant area of the quotient digits on the FiPD diagram, that is, in the shaded area in #12FjA. By the way, this boundary can be anything as long as it falls within the shaded area, so the divisor d and partial remainder rP4
When two roads are displayed, the boundary can be defined as a point with a good separation (for example, a point that can be expressed with a small number of pits). An example of the boundary in the case of r=4, ρ=2, and k= is shown in FIG. In this example, the correct quotient digit can be found from the top four pits of the divisor and partial remainder. This third
The hardware version of the PD diagram with boundaries as shown in the figure is
For example, it can be realized by a read-only memory (ROM'), and the quotient digit can be obtained at high speed by using this ROM as a table. This hardware P
The D diagram will be referred to as a PD table hereinafter. The above idea has already been known in the past, for example, the following document Daniel E, Atklns “H1gh@
r-Radlx Divison Using
Istlmatssof the Dlvjser
mad Partial R*majmdsr
s”IEICE Translations on C
@mput@r@, vol, c-17°A 10,0
ctob@r 1968. It is described in detail in PP925-934°.

さて次に従来の高基数除算回路の代表的な例を示し、そ
の問題点を明らかにする。
Next, we will show a typical example of a conventional high-radix division circuit and clarify its problems.

第6図は従来の高基数除算回路の一例を示す構成図であ
る。図中、jd除除数管置くレジスタ、2は部分剰余と
して初期設定される被除数P0およびそれ以後の各演算
ステップで部分剰余を置くレジスタである。5は乗算回
路であり、レジスタ1に置かれている商デイゾッ) q
3+1と前記レジスタ3に置かれている除数dを入力し
てその積qj+ldを得るものである。また8は加算回
路であり、前記レジスタ2に置かれている部分剰余rP
jと乗算回路5の出力(lj+、dとから新しい部分剰
余Pj+1を計算するものである。尚、基数rとして祉
2のべき乗r=2を仮定しておシ、従って1”Pjは部
分剰余Pjをnピット左シフトした亀のとなる。即ち加
算回路8は次の計算をする。   ′ PJ+1 = rP6−1’J+l dまた4は前置加
算回路であシ前記レジスタ2に置かれる部分剰余の上位
N、ピットと前記乗算回路5の出力qj+ldの上位N
mピットから次の演算ステツノに用いる部分剰余Pj+
1の近似値P3+1を計算している。この近似値Pj+
 !はP3+1の上位N、/  ピットにほぼ勢しい。
FIG. 6 is a block diagram showing an example of a conventional high radix division circuit. In the figure, the register 2 in which the jd divisor tube is placed is a register in which the dividend P0, which is initially set as a partial remainder, and the partial remainder in each calculation step thereafter are placed. 5 is a multiplication circuit, and the quotient placed in register 1) q
3+1 and the divisor d placed in the register 3 are input to obtain the product qj+ld. Further, 8 is an adder circuit, and the partial remainder rP placed in the register 2 is
A new partial remainder Pj+1 is calculated from j and the output (lj+, d) of the multiplier circuit 5. It is assumed that the base r is a power of 2, r=2, so 1''Pj is the partial remainder It becomes a turtle by shifting Pj to the left by n pits.In other words, the adder circuit 8 performs the following calculation: 'PJ+1 = rP6-1'J+l d or 4 is a pre-adder circuit.The partial remainder placed in the register 2 is The upper N of the pit and the output of the multiplication circuit 5 qj+ld
Partial remainder Pj+ used for the next calculation step from m pit
An approximate value of 1, P3+1, is being calculated. This approximate value Pj+
! is in the top N of P3+1, and is almost in the pit position.

この前置加算回路4社上記部分剰余PJ+ 1の近似値
Pj+1をより速く求める九めに置かれた4のである。
These four pre-addition circuits are placed in the ninth position to more quickly obtain the approximate value Pj+1 of the partial remainder PJ+1.

6は商デイジツト発生回路、即ち前述したPDテーブル
であり、前記レジスタ3に置かれている除数dの近似7
(to上位Ndピット)と前記前置加算回路4の出力”
jj’3 + 1とから商デイジツ)J+2を生成する
。このようにして除11kdの近似値dと部分剰余Pj
の近似値p5かり正しい商が得られるのは、#IK述べ
たようにd t P、5+1の上位数ピットでPDテー
ブルCを引くことができるようKし九九めである。但し
、上の説明中のNP、 Nm、 N、’、 Ndの値社
商ディジットの最大値ρと基数r=2とを与える仁とに
よって定まる。詳しくは後述するが、例えばr=2’ 
、I=2’−1の場合には、Np=11ビット、Nm 
= 11ビツト、NP’= 6ビツト、Nd=5ピット
あれば十分である。従ってPDテーブル6紘入力NP’
 + N、1 = 11ビツト、出力4ピツトOR6M
を用いることによって実現できる。えだしPDテーブル
6の出力qJ +Iの符号は別に設けられ九回路によっ
て発生される。
6 is a quotient digit generation circuit, that is, the PD table mentioned above, which calculates the approximation 7 of the divisor d placed in the register 3.
(to upper Nd pit) and the output of the pre-adder circuit 4”
From jj'3 + 1, the quotient digit) J+2 is generated. In this way, the approximate value d of the division 11kd and the partial remainder Pj
The correct quotient can be obtained from the approximate value p5 of d t P, as mentioned above, in order to be able to draw the PD table C using the upper number pits of 5+1. However, the values of NP, Nm, N,', and Nd in the above explanation are determined by the maximum value ρ of the company digit and the value that gives the base r=2. The details will be described later, but for example, r=2'
, I=2'-1, Np=11 bits, Nm
= 11 bits, NP' = 6 bits, and Nd = 5 pits are sufficient. Therefore, PD table 6 Hiro input NP'
+N, 1 = 11 bits, output 4 pits OR6M
This can be achieved by using The sign of the output qJ +I of the starting PD table 6 is separately provided and generated by nine circuits.

第7図はこのように構成された従来の除算回路の動作手
順を示すタイムシーケンス図である。
FIG. 7 is a time sequence diagram showing the operating procedure of the conventional division circuit configured as described above.

このタイムシーケンス図を参照してその除算実行手順を
説明すると、次のようになる。
The division execution procedure will be explained as follows with reference to this time sequence diagram.

(1)初期設定、除数、被除数をそれぞれレジスタS、
レジスタ2にセットするとともにレジスタ1をクリアす
る。(j4−0 )第7図中タイミングtoに相当する
(1) Initial settings, divisor, and dividend are set in register S, respectively.
Set in register 2 and clear register 1. (j4-0) Corresponds to timing to in FIG.

(2)  レジスタ2の上位ピットとレジスタ3の上位
ピットからPDチー!ル6によって第1の商ディジツト
q1を生成してレジスタ1にセットする。このとき乗算
回路5の出力  jは0であるから、レジスタ2の上位
ピット   □゛は前置加算回路4をそのまま通過する
(2) PD chi! from the upper pit of register 2 and the upper pit of register 3! A first quotient digit q1 is generated by the register 6 and set in the register 1. At this time, since the output j of the multiplication circuit 5 is 0, the upper pit □゛ of the register 2 passes through the pre-addition circuit 4 as is.

(j+11第7図中タイミングt、に相当する。(j+11 corresponds to timing t in FIG. 7).

(3)  レジスタSの内容とレジスタ1の内容の積が
乗算回路5によって計算され、この積の結果とレジスタ
2の内容を入力して加算回路1は次のステツノに対する
部分剰余PJ+1を計算して出力する。tfc上記動作
と並行して乗算回路5の出力の上位ビットとレジスタ2
の上位ピットから前置加算回路4Fi部分剰余の近似値
Pj+ttf¥i、9. PDテーブル6は前記Pj+
1とレジスタ3の上位ビットから商ディジツトqj+l
を生成する。
(3) The product of the contents of register S and the contents of register 1 is calculated by multiplier circuit 5, and by inputting the result of this product and the contents of register 2, adder circuit 1 calculates the partial remainder PJ+1 for the next step. Output. tfc In parallel with the above operation, the upper bit of the output of the multiplication circuit 5 and the register 2
Approximate value Pj+ttf\i of the partial remainder of the pre-addition circuit 4Fi from the upper pit of 9. PD table 6 is the Pj+
1 and the upper bit of register 3 to the quotient digit qj+l
generate.

(4)加算回路8の出力及びPDテーブル6の出力をそ
れぞれレジスタ2、レジスタ1にセットしてステツノ(
3)へ戻る。(14−j++)第7図中タイtングLt
 、t、*・・・に相当する。
(4) Set the output of the adder circuit 8 and the output of the PD table 6 to register 2 and register 1, respectively.
Return to 3). (14-j++) Timing Lt in Figure 7
, t, *...

このようにしてステラf (3) (4)の1回の反復
ごとに商ディジツトが1デイジツトづつ求t)、m回の
反復の結果商社Σqtr−” なる形として1=1 求まる。丸だしこの例における除算回路の生成する商デ
ィジツトは符号つきなので、これによって求められる商
は符号付ディジット数値表現(81g5*d dlgl
t numb@r 5yst@m )となる。この符号
付ディジット!!!現とは例えば次のような轄値表現を
意味する。
In this way, the quotient digit is calculated by one digit for each iteration of Stella f (3) (4) (t), and as a result of m iterations, 1=1 is found in the form ``Shosha Σqtr-''. Since the quotient digit generated by the division circuit in is signed, the quotient obtained by this is expressed as a signed digit numerical representation (81g5*d dlgl
t number@r 5yst@m). This signed digit! ! ! For example, ``present'' means the following expression of value.

Nm2百102・・・ (基数4) この例で示さ゛れるように、ディジット上の横棒記号が
、そのディジットが負であることを示している。このよ
うな符号付ディジット表現を普通の表現、例えば2の補
数表現か或いは符号絶対値表現に変換する為に蝶、実際
に商ディジツトの加算仝q1r−4を行う必要がある。
Nm2102... (Radix 4) As shown in this example, the horizontal bar symbol on the digit indicates that the digit is negative. In order to convert such a signed digit representation into an ordinary representation, such as a two's complement representation or a signed absolute value representation, it is necessary to actually perform the addition of quotient digits (q1r-4).

上の餌1=1 では例えば次のようにして行われる。Upper bait 1 = 1 For example, this is done as follows.

Nm2 X 4−1−3X4”−” +lX4−s+0
X4−’−2X4−’ +−・・特に符号付ディジット
表現に0を含まないとき上記変換において?ローが2桁
以上伝播しないのでその変換を容易に行うことができる
。首九、最大冗長(k=x)の場合には商ディジッ  
−トの冗長性によって商ディジツトに0を用いないで除
算を実行することが可能である。
Nm2 X 4-1-3X4"-" +lX4-s+0
X4-'-2X4-' +-... Especially in the above conversion when the signed digit representation does not include 0? Since rows do not propagate more than two orders of magnitude, the conversion is easy. In the case of maximum redundancy (k=x), the quotient digit
- It is possible to perform the division without using zero in the quotient digit due to the redundancy of the bits.

以上示したように従来の高基数除算回路は、その基数が
小さい場合にはPDテーブル6が小さくなるので簡潔な
ハードウェア構造によって実現できる。特に基数がよシ
小さい場合には、第5図における前置加算路4や乗算回
路5やPDテーブル6が非常に単純表ものとなる。その
極限が基数2に対応する非回復法による除算回路に和尚
する。ところが、従来の高基数除算回路の基数rは24
程度が限界である。何故ならもし基数をr=24よシ大
きくしようとするとPDテーブルCの大きさは非常に膨
大なハードウェア構成となシ、ハードウェアの他の部分
とのバランスから見て、現実的でなくなる。この様子に
ついて、次Kl!際に必要なRAMテーブル(PDテー
ブルC)の大きさにつき、解析してみる。
As shown above, the conventional high radix division circuit can be realized with a simple hardware structure because the PD table 6 becomes small when the radix is small. Particularly when the radix is very small, the pre-adder circuit 4, multiplier circuit 5 and PD table 6 in FIG. 5 become very simple tables. The limit corresponds to a division circuit using a non-recovery method corresponding to radix 2. However, the radix r of the conventional high radix division circuit is 24.
The extent is the limit. This is because if the radix were to be made larger than r=24, the size of the PD table C would become an extremely large hardware configuration, which would be unrealistic from the standpoint of balance with other parts of the hardware. About this situation, next Kl! Let's analyze the size of the RAM table (PD table C) required in this case.

第3図に示したPD線図につき再び着目すると、図中書
き込まれた階段の1ステツプの最小の幅はQj+、 =
=βとq!+1−i −+の重複領域の最小除数d=十
の近辺に現れる。従って、この近辺OPDPD線図べる
と必要なPDテーブル6の大きさを知ることができる。
Focusing again on the PD diagram shown in Figure 3, the minimum width of one step of the stairs drawn in the diagram is Qj+, =
= β and q! It appears near the minimum divisor d=10 in the overlapping region of +1−i −+. Therefore, if the OPDPD diagram in this vicinity can be obtained, the necessary size of the PD table 6 can be known.

第4図はこの領域を拡大したものである。鮪4図中のb
 、 dy 、 Pyの値を求めると次のようになる。
FIG. 4 shows an enlarged view of this area. b in tuna 4 diagram
, dy , and Py are calculated as follows.

最小の除数の間隔及び部分剰余の間隔はそれぞれdw/
!2 t PW/’Iと考えれば十分にその解析を行い
得る。除数d、部分剰余rPjの範囲はそれぞれ+≦d
(1、O≦trPj1(rであるから、除数、部分剰余
に必要なビット数Nd 、 N、は次のように求められ
る。
The minimum divisor interval and partial remainder interval are respectively dw/
! 2 t PW/'I can be sufficiently analyzed. The range of divisor d and partial remainder rPj is +≦d, respectively.
(1, O≦trPj1(r), so the number of bits Nd and N required for the divisor and partial remainder can be found as follows.

この場合特に基数r=21、最大冗長に=1のときFi
                         
賑Nd=1+tピット、  Np=2+tピットとなJ
)、 k=1 、 r=2’とすればPDテーブル  
′6の大きさは入力11ビツト、出力4ビツトとして現
実的な太きさとして求められる。しかしにツ1.r=2
’のときはPDテーブル6の大きさは入力15ピツト、
出力6ピツトと非常に大きくなる。但し、ここではPD
テーブル6の出力の符号は別の回路によって求めるもの
として上記PDテーブル6の大きさには含めてない。
In this case, especially when the base r=21 and the maximum redundancy=1, Fi

Lively Nd=1+t pit, Np=2+t pit and J
), k=1, r=2', PD table
The size of '6 is determined as a realistic thickness assuming 11 input bits and 4 output bits. However, 1. r=2
', the size of PD table 6 is input 15 pits,
It has a very large output of 6 pits. However, here, PD
The sign of the output of the table 6 is not included in the size of the PD table 6 as it is determined by a separate circuit.

このように従来の高基数除算方式では、その高速化のた
めに基数を大きくしようとすると、PDテーブル6の大
きさが大きくなりすぎるという問題があつ九。この点本
発明においてはPDテーブル6の大きさに関する間−を
次のようにして取シ除いている。本発明に係る基本原理
は除数、被除数にそれぞれ同一のある数を掛けることに
よ如除数をlに近づけることである。これにより除数の
値を取p得る範囲が限定され、この結果、PDテーブル
6の大きさを減小させることができる。
As described above, the conventional high-radix division method has a problem in that when the radix is increased to increase the speed, the size of the PD table 6 becomes too large. In this regard, in the present invention, the gap regarding the size of the PD table 6 is removed as follows. The basic principle of the present invention is to bring the divisor closer to l by multiplying the divisor and dividend by the same number. This limits the range from which the divisor value can be obtained, and as a result, the size of the PD table 6 can be reduced.

このことについて、更に詳しく説明すると次の通勤であ
る。             −今、何んらかの方法
によって除数−の値の取シ得る範囲か1≦dく1と表っ
たとする。、この時のPDテーブル6の大きさは前述し
た従来例の場合と同様な手法によって次のように求める
ことができる。
I will explain this in more detail in my next commute. - Now, suppose that the possible range of the value of the divisor - is expressed as 1≦d×1 by some method. , the size of the PD table 6 at this time can be determined as follows using the same method as in the conventional example described above.

ここで上式におけるaが1に近づけばNP * Ndは
減小する。この効果は特に除数のピット数Ndにおいて
著しい。従ってNdがある程度小さくなれば(1〜3ピ
ツト)、被除数の上位Npビットが表わす値に除数の上
位Ndピットに応じたある修正を施すという考え方に基
づくことによって、PDテーブル6を比較的簡単外論理
回路によって実現することが可能となる。更に鳳を1に
近づければ商ディジツトの選択において部分剰余の上位
ピットのみを用い、除数についての情報が不必要になる
。即ちこの場合、商ディジツトを部分剰余の上位ビット
から極めて簡単に作ることができる。
Here, as a in the above equation approaches 1, NP*Nd decreases. This effect is particularly remarkable in the number of divisor pits Nd. Therefore, if Nd becomes small to a certain extent (1 to 3 pits), the value represented by the upper Np bits of the dividend is modified in accordance with the upper Nd pits of the divisor, and PD table 6 can be modified relatively easily. This can be realized using logic circuits. Furthermore, if Otori is brought closer to 1, only the upper pits of the partial remainder are used in selecting the quotient digit, and information about the divisor becomes unnecessary. That is, in this case, the quotient digit can be created extremely easily from the upper bits of the partial remainder.

ちなみに商ディジツトを求めるとき、その除数について
の情報を不必要にする為の1の値の下限を導出するには
次のようにすればよい。即ち、第5図に示すように1を
選び商ディジツトρとll−1の重複領域に境界ムを入
れると、商デーイジツトの選択は部分剰余のみに依存し
、除数が関与しなくなる。仁のときの1の値を求めると
、 となる。従って係数ζを適当に選んで なる条件が満九されるようKすればよい。
Incidentally, when calculating the quotient digit, the lower limit of the value of 1 can be derived as follows to make information about the divisor unnecessary. That is, as shown in FIG. 5, if 1 is selected and a boundary is placed in the overlapping region of the quotient digit ρ and ll-1, the selection of the quotient digit depends only on the partial remainder and the divisor is not involved. When we find the value of 1 when it is jin, we get the following. Therefore, it is only necessary to select the coefficient ζ appropriately and set K so that the condition formed by satisfying 9 is satisfied.

本発明はこのような理論に立脚してなされ九もので、除
数、被除数に掛ける定数ζを除数dO逆数の近似値とし
て論理回路又はROMテーブルによって求め、これによ
って除数の値を1に近づけるようにしたものである。仁
のハードウェア化され九逆数テーブルの存在によってP
Dテーブルの大きさが減小せしめ、あるいはPDテーブ
ルを不要とすることができる。このPDテーブルの減小
の効果は逆数テーブルの付加によるハードウェアの増加
よシずつと大きい。これを証明する丸めに以下において
本発明方式を実施するに際して必要とされる逆数テーブ
ルの大きさについて説明する。
The present invention has been made based on this theory, and the constant ζ multiplied by the divisor and dividend is determined as an approximate value of the reciprocal of the divisor dO using a logic circuit or a ROM table, thereby bringing the value of the divisor closer to 1. This is what I did. Due to the existence of Jin's hardware nine reciprocal table, P
The size of the D table can be reduced or the PD table can be made unnecessary. The effect of reducing the PD table is greater than the increase in hardware due to the addition of the reciprocal table. To prove this, the size of the reciprocal table required when implementing the method of the present invention will be explained below.

この設定された問題を明確にすると、a (1であって
1に近い定@aが与えられたとき、条件 1≦ζdく1        ・・・・・・alを満足
するような数ζを除数dの近似値から発生しようとする
とき、必要なdの近似値の精度と数ζ自身のピット長を
求めることである。
To clarify this set problem, a (1 and a constant @a that is close to 1 is given, then the condition 1≦ζd×1 ...... al is the number ζ that satisfies the divisor. When attempting to generate from an approximate value of d, the required accuracy of the approximate value of d and the pit length of the number ζ itself are determined.

今、a = 1−2  とし、且つ−の上位tピットを
d(+とし、d(1)の正確な逆数をζ0.ζeの上位
Uピットをζ(1ゝとすると、次の関係が成立する。 
   (d=d”’+2” d”   、+≦a<l 
、  a(1)(x 、a■〈□   tlζ。=ζ(
+)+2−u+1 (■、  1<(’。≦2 、 (
”)< 2 、 ζ’<1これらの関係から積ζ(1)
dを計算するとζ”’d=  1+2−tζd■−2−
H+l  (■dとな夛、従って次の関係式が得られる
Now, if a = 1-2, and the upper t pit of - is d(+, and the exact reciprocal of d(1) is ζ0.The upper U pit of ζe is ζ(1゜), the following relationship holds. do.
(d=d"'+2"d",+≦a<l
, a(1)(x , a■〈□ tlζ.=ζ(
+)+2-u+1 (■, 1<('.≦2, (
”) < 2, ζ'< 1 From these relationships, the product ζ (1)
Calculating d is ζ"'d= 1+2-tζd■-2-
H+l (■d), so the following relational expression is obtained.

1−2−”’<ζ(1)a (1+ 2− t ” 1
ここでζ=ζ(+1−2−t+2  とお砂ばζ、−ζ
(1)d2−t+24≦ζ(j)4 2−t+1(1即
ち 1−21+1−2−”<ζd〈1 とl)簡単のために一=tであるとすれば、1−2−1
+2<ζd<1 となる。この関係からu = t = s + 2が求
められる仁とになる。
1-2-”'<ζ(1)a (1+ 2- t ” 1
Here ζ=ζ(+1-2-t+2 and Osunaba ζ, -ζ
(1) d2-t+24≦ζ(j)4 2-t+1 (1, that is, 1-21+1-2-"<ζd<1 and l) For simplicity, if 1=t, then 1-2-1
+2<ζd<1. From this relationship, u = t = s + 2 becomes the required value.

以上まとめると除数−の上位 +2ビツトを取〉出して
その正確な逆数の上位 +2ビツトから21を引い走数
をことすると、このζは前述し九条件第輔式を満足する
。従って逆数テーブルは入力8−2ビツト、出力8+2
ビツトあれば十分に構成することができる。更に”A<
a<1を考慮すると逆数テーブルの入力および出力とも
に1ピツトづつ減小できる。即ち逆数テーブルの大きさ
は 入力=8+1ビット、出カニS+1ピット ・・・ (
ロ)であればよいことが示される。
To summarize the above, if we take the upper +2 bits of the divisor - and subtract 21 from the upper +2 bits of its exact reciprocal, then this ζ satisfies the ninth condition (the fifth equation) mentioned above. Therefore, the reciprocal table has an input of 8-2 bits and an output of 8+2 bits.
Bits are sufficient for configuration. Furthermore, “A<
Considering a<1, both the input and output of the reciprocal table can be reduced by one pit. In other words, the size of the reciprocal table is input = 8 + 1 bits, output crab S + 1 pit... (
b) indicates that it is acceptable.

このような逆数テーブルの設置によって達成される本発
明の効果は例えば次の通りである。
The effects of the present invention achieved by installing such a reciprocal table are as follows, for example.

即ち、第1の例として与えられた除数d(+≦d<1)
に対して1−2″″4≦ζdくlとする為には第(2)
式によって入力5ピツト出力5ピツトの逆数テーブルを
設置すればよいことが判る。このとき、基数r=2一完
全冗長に≦1の除算を考えると、必要なPDテーブルの
大きさは第(8)式および第(9)式によシ入力N、 
十Nd = 12ピツト、出力6ピツトである。但し、
第(8)、(9)式は必要なピット数を多めに見積っで
あるので実際に必要なPDテーブルの入力ビット数は1
ピツトあるいは2ピツト減小するととがあるかもしれな
い。ともか(PDテーブルの大きさは現実的なハードウ
ェア量である。これに反し従来の方式では同じ除算に対
してPDテーブル6の大きさは入力15ビツト、出力6
ピツトが必要でToり、ハードウエア化に際して、非現
実的である。
That is, the divisor d (+≦d<1) given as the first example
In order to satisfy 1-2″″4≦ζd×l, use (2)
It can be seen from the equation that a reciprocal number table with 5 input pits and 5 output pits should be installed. At this time, considering the division of radix r = 2 - completely redundant ≦1, the size of the required PD table is determined by the input N, according to equations (8) and (9).
10 Nd = 12 pits, output 6 pits. however,
Equations (8) and (9) estimate the number of pits required, so the actual number of input bits for the PD table required is 1.
There may be a decrease in pits or two pits. (The size of the PD table is a realistic amount of hardware. On the other hand, in the conventional method, the size of the PD table 6 for the same division is 15 bits for input and 6 bits for output.
This requires a large pit and is unrealistic in terms of hardware implementation.

まえ別の例を考えると、商ディジツトを生成するのに除
数d4Cついての情報が不必要になる為の条件は第(至
)式が成立することである。基数r工2−斃[k= 1
の除算を行うものとすると、この条件は次のように々る
Considering another example, the condition for the information about the divisor d4C to become unnecessary for generating the quotient digit is that the following equation holds true. Radix r work 2-[k= 1
If we perform the division of , then this condition is as follows.

1−2″″aくζd<1 この場合、与えられた除数dに対して上式を満足させる
ような−の逆数の近似値ζを発生するのに必要な逆数テ
ーブルの大きさは第(2)式によって求められ、入力9
ビツト、出力9ビツトと々る。これに反し従来の高基数
除算方式では同じ除算に対して必要なPDテーブル6の
大きさ拡入力19ビット、出力8ピツFであ如、膨大な
ハードウェア量となる。
1-2""a x ζd<1 In this case, the size of the reciprocal table necessary to generate an approximate value ζ of the reciprocal of - that satisfies the above equation for a given divisor d is the size of the reciprocal table ( 2) Calculated by the formula, input 9
The output reaches 9 bits. On the other hand, in the conventional high-radix division method, the size of the PD table 6 required for the same division is 19 bits for expanded input and 8 bits for output, resulting in an enormous amount of hardware.

次に、上述した演算理論に立脚し九本発明方式の一実施
例につき説明する。
Next, one embodiment of the present invention method will be described based on the above-mentioned arithmetic theory.

第8図は本発明の実施例方式を適用した除算回路O構成
図であ夛、高速表乗算・除算回路を実現するものである
。図中、23.24はセレクタ、25.26は入力レジ
スタ、28,29゜30.31,811は前置乗算回路
、33,34゜ss、sr;、sr、saはキャリー保
留加算回路、19.40はノfイブラインレジスタ、4
6はキャリー伝播加算回路でToシ、これらによって乗
算回路が構成されている。更に、22は逆数テーブル、
41は前置加算回路、42社商デイジツト発生回路、4
6はノ4イブラインレジスタであシ、これらと前記乗算
回路とKよって除算回路が構成されている。
FIG. 8 is a block diagram of a division circuit O to which an embodiment of the present invention is applied, which realizes a high-speed table multiplication/division circuit. In the figure, 23.24 is a selector, 25.26 is an input register, 28, 29° 30.31, 811 is a premultiplier circuit, 33, 34° ss, sr;, sr, sa are carry-pending addition circuits, 19 .40 is nof bline register, 4
6 is a carry propagation adder circuit; these constitute a multiplication circuit. Furthermore, 22 is a reciprocal table,
41 is a pre-addition circuit; 42 is a commercial digit generation circuit; 4
Reference numeral 6 denotes a four-line register, and these, the multiplication circuit, and K constitute a division circuit.

ここでは、上記乗算回路は32ピツ)XS)ピットの演
算能力を持っている。この乗算の原理は次式によって表
わせる。
Here, the multiplier circuit has an arithmetic capacity of 32 pits)XS) pits. The principle of this multiplication can be expressed by the following equation.

X:被乗数、Y:乗数 = (−2yo+ys+yt)2−” + (−2yt
+ys +74 )2−”+(2F4+F!J+76)
2−”+(−276+yy+ys)2−丁+(−2ya
 +ys)2−・          j積X Y =
2−”(−276+yl +7m )X+2−”(2)
’*+3’t+74)x   ’+2””(−2y++
ys+ys)X+2イ(2FM+3’? +FI )X
+2′(−2ya +y* )X 図中、zg、ip、・・・、J2はそれぞれ上式におけ
る各項(Zyo +Fr +F冨)Xs (2yt+F
s+y4)X+・・・* (2yi ” F@ )X 
 の演算をそれぞれ実現するもので、Yのピットノ4タ
ーンに応じてXの0倍、+1倍、+2倍を生成する。4
1KKの負数倍は例えばまずXを正数倍し、次にこの結
果の1の補数をと多、更に補正入力27の適尚な位置K
lを加えて2の補数を得るととくよってなされる。
X: Multiplicand, Y: Multiplier = (-2yo+ys+yt)2-” + (-2yt
+ys +74)2-”+(2F4+F!J+76)
2-”+(-276+yy+ys)2-cho+(-2ya
+ys)2−・j product X Y =
2-”(-276+yl +7m)X+2-”(2)
'*+3't+74)x '+2""(-2y++
ys+ys)X+2i(2FM+3'? +FI)X
+2'(-2ya +y*)X In the figure, zg, ip,..., J2 are each term (Zyo +Fr +F)
s+y4)X+...* (2yi ”F@)X
It realizes the calculations respectively, and generates 0 times, +1 times, and +2 times of X according to the four turns of the pit nozzle of Y. 4
To multiply 1KK by a negative number, for example, first multiply X by a positive number, then multiply the 1's complement of this result, and then select the appropriate position K of the correction input 27.
This is done by adding l to obtain the two's complement number.

この実施例では32ピツ)X9ビツトの積の結果はサム
出力およびキャリ出力の2っに分離してそれぞれノ譬イ
ブラインレジスタ41.14にセットされる。完全外積
の結果を得ゐ為には2つの/ヤイ!ラインレジスタ43
.44の出力をキャリー伝播加算Fm*45によって加
算することが必要である。
In this embodiment, the result of the product of 32 bits (32 bits) x 9 bits is separated into two outputs, a sum output and a carry output, and is set in the parallel register 41.14, respectively. To get a perfect cross product result, you need two /Yay! line register 43
.. It is necessary to add the outputs of 44 by carry propagation addition Fm*45.

本実施例Kj?いて32ビツト×32ピツトの乗算は、
32ピツト×8ピツトの乗算を4回行い最後Ktサム出
力キャリー出力とをキャリ伝播加算N路46によって加
算することによって実行される。第9図はこのような乗
算処理実行過穆を示すタイムシーケンス図である。
This example Kj? The multiplication of 32 bits x 32 pits is
This is executed by performing 32-pit x 8-pit multiplication four times and adding the final Kt sum output and carry output by the carry propagation addition N path 46. FIG. 9 is a time sequence diagram showing an overexecution of such multiplication processing.

次に本実施例における除算の実行について第1θ図に示
す除算実行タイムシーケンス図を用いて説明する。仁の
除算は以下の手順で行われる。即ち、 (1)  除数dを入力線ja上に乗せる。この人力線
20の上に乗せられえ除数dの上位9ビツトは逆数テー
ブル22に与えられる。第10図中タイミングtに相当
する。
Next, the execution of division in this embodiment will be explained using the division execution time sequence diagram shown in FIG. 1θ. The division of jin is performed in the following steps. That is, (1) Place the divisor d on the input line ja. The upper 9 bits of the divisor d that can be placed on the human force line 20 are given to the reciprocal table 22. This corresponds to timing t in FIG.

(2)  入力線20上のデータ、即ち除数dをレジス
タ25にセットし、同時に逆数テーブル22の出力、即
ちdの逆数の近似値をレジスタ26にセットする。第1
0図中タイミングt1に相当する。
(2) Set the data on the input line 20, ie, the divisor d, in the register 25, and at the same time set the output of the reciprocal table 22, ie, the approximate value of the reciprocal of d, in the register 26. 1st
0 corresponds to timing t1 in the figure.

(3)  レジスタzs、zgの内容をそれぞれ被除数
、除数として乗算を実行し、その結果を・譬イ!ライン
レジスタ39及び4oにセットする。
(3) Execute multiplication using the contents of registers zs and zg as the dividend and divisor, respectively, and calculate the result. Set in line registers 39 and 4o.

この動作と並行して入力線:lO,,):に被除数P、
が乗せられて前記レジスタ25にセットされる。嬉10
図中タイ建ングt、に相当する。
In parallel with this operation, the input line: lO,, ): receives the dividend P,
is loaded and set in the register 25. Happy 10
This corresponds to the tie construction t in the figure.

(4)  レジスタ25とレジスタ26の内容が相互に
乗算され、その結果がノ々イグラインレジスタ39.4
0にそれぞれセットされる。この動作と並行してノダイ
グラインレジスタ39゜40の内容が加算回路45によ
って加算され、その結果がレジスタ25にセットされる
。第10図中タイミング1.に相当する。
(4) The contents of register 25 and register 26 are multiplied together, and the result is sent to non-ignition register 39.4.
Each is set to 0. In parallel with this operation, the contents of the no-dig line registers 39 and 40 are added by the adder circuit 45, and the result is set in the register 25. Timing 1 in Figure 10. corresponds to

以上によシ除算の初期設定が(j4−0)完了し、続い
て以下のように除算の実行が行われる。
The initial setting of the division is completed (j4-0) as described above, and then the division is executed as follows.

(5)  レジスタ39.40のそれぞれ上位数ヒツト
を前置加算回路41によって加算して部分剰余pro近
似をP、1を求める。このPjから商デイジツト発生回
路42によって商ディジツト”j+1を得る。前述し九
乗算回路によって”j+1とレジスタ25の内容との積
が計算され、レジスタ19及び400部分剰余から上記
積の値を引き算して新しい部分剰余Pj+1を求める。
(5) Add the high-order hits of the registers 39 and 40 by the pre-adder circuit 41 to obtain the partial remainder pro approximation P,1. The quotient digit "j+1" is obtained from this Pj by the quotient digit generating circuit 42.The product of "j+1" and the contents of the register 25 is calculated by the aforementioned nine multiplication circuit, and the value of the above product is subtracted from the partial remainders of the registers 19 and 400. to find a new partial remainder Pj+1.

(6)演算結果即ち新しい部分剰余P、1+1をレジス
タ3m、40にセットし、上記ステラf(5)へ戻る。
(6) Set the operation result, that is, the new partial remainder P, 1+1, in the registers 3m, 40, and return to the above Stella f(5).

[j+−j+13第10図中部分剰余のタイミングはj
4t  jl+  461  t、に相当する。
[j+-j+13 In Figure 10, the timing of the partial remainder is j
This corresponds to 4t jl+461 t.

このようにして除算はlステツブ当や8ピツトづつ行わ
れ、32ピツトと32ピツトの除算は上記ステラf (
5) I (6)を反復して5回縁秒返す。
In this way, division is performed per l step or 8 pits at a time, and division between 32 pits and 32 pits is performed using the above Stellar f (
5) Repeat I (6) to return edge seconds 5 times.

このとき各ステップごとに符号付商ディジツトが生成さ
れるが、前述したように正しい商を得る為にはこの商デ
ィジツトの加算を行なうことが必要となる。しかしてこ
の商ディジツトの加算は上記ステップ(5) I (6
)の反復と並行して加算回路46によって行なわれる。
At this time, signed quotient digits are generated for each step, but as described above, it is necessary to add these quotient digits in order to obtain a correct quotient. However, addition of this quotient digit is done in step (5) I (6
) is performed by the adder circuit 46 in parallel with the repetition of .

この商ディジツトの加算は第10図中キャリ伝播加算の
タイミングtswLslt丁、t、に相当する。
This addition of the quotient digits corresponds to the timing tswLsltd,t of the carry propagation addition in FIG.

このようにして本実施例では高基数除算はl求めたよう
に入力9ビツト出力9ピツトであれ   ゛ばよい。従
ってこの逆数テーブル22を附加することによって従来
不可能であった基数21の高基数除算回路が小規模なハ
ードウェアによって実現できることになる。
In this way, in this embodiment, high radix division can be performed as long as the input is 9 bits and the output is 9 pits, as determined by l. Therefore, by adding this reciprocal table 22, a high radix division circuit of radix 21, which was previously impossible, can be realized with small-scale hardware.

以上述べ喪ように本発明によれば従来の高基数除算方式
では達成できないような高基数の除算回路を簡単なハー
ドウェア構造によって1!現でき、高速な除算が可能に
なる等の絶大なる効果を奏し、その実用的利点は絶大で
ある。
As stated above, according to the present invention, a high-radix division circuit that cannot be achieved with conventional high-radix division methods can be achieved using a simple hardware structure. It has great effects, such as being able to express and perform high-speed division, and its practical advantages are enormous.

尚、本発明は上記実施例に限定される亀のではない。例
えば除数dおよび被除数p、のビット数は、演算仕様に
応じて定めればよい亀のであり、要するに本発明はその
要旨を逸脱しない籠囲で種々変形して実施することがで
きる。
It should be noted that the present invention is not limited to the above embodiments. For example, the number of bits of the divisor d and the dividend p may be determined according to the calculation specifications, and in short, the present invention can be implemented with various modifications without departing from the gist thereof.

【図面の簡単な説明】[Brief explanation of drawings]

第1図は除数及び現在の部分剰、余と、選択可能な商デ
ィジツト及びそれを選択し九ときの新しい部分剰余との
関係を示す図、第2開拡除数と部分剰余とから選択可能
な商ディジツトを与える為のr=4.#x20場合のF
D−図、第3mlは第2図に示すPD−図における商テ
イジツ)0重横領域の境界を示す図、第4開拡第3図に
示すPD線図の一部を拡大して示した図、第5図はPD
テーブルの減少を説明する為のPD線図、第6開拡従来
の高基数除算回路の一例を示す構成図、第7開拡第5図
に示す回路の動作タイミング図、第8図は本発明の一実
施例を適用した高基数除算回路の構成図、第9図は第8
図に示す回路の乗算動作タイミング図、第1θ図は第8
図に示す回路の除算動作タイきング図である。 1・・・セレクタ、ill、3.7・・・レジスタ、4
・・・前置加算回路、5・・・乗算回路、6・・・PD
テーブル、8・・・加算回路、23,24.4’l、4
3゜44−・セレクタ、25.16.3B、40.46
・・・レジスタ、33,34.:16,36,37゜3
8・・・キャリー保留加算回路、211.29,30゜
31.32・・・前置乗算回路、42・・・前置加算回
路、45・・・キャリー伝播加算回路、42・・・商デ
イジツト発生回路、22・・・逆数テーブル。 出願人代理人 弁理士 鈴 江 武 彦第5図 第6図 qj争1 第7図 t□     tl     t2     t3第9
図 I4傳 キャリーイλ− タイS〉2′5占3゛ 第10図 靭νU惺 MP*2金 1ζ−)−434に−iπi哨1 タイS〉りf九3
Figure 1 shows the relationship between the divisor, the current partial remainder, the remainder, the selectable quotient digits, and the new partial remainder when they are selected. r=4 to give the quotient digit. F for #x20
Figure D, 3ml is a diagram showing the boundary of the 0-fold horizontal area in the PD diagram shown in Figure 2, and a part of the PD diagram shown in Figure 3 is enlarged. Figure 5 is PD
A PD diagram for explaining the reduction in the number of tables, a block diagram showing an example of a conventional high-radix division circuit in the sixth expansion, an operation timing diagram of the circuit shown in FIG. 5 in the seventh expansion, and FIG. FIG. 9 is a block diagram of a high radix division circuit to which an embodiment of the invention is applied.
The multiplication operation timing diagram of the circuit shown in the figure, 1θ is the 8th
FIG. 3 is a division operation timing diagram of the circuit shown in the figure. 1...Selector, ill, 3.7...Register, 4
...Pre-addition circuit, 5...Multiplication circuit, 6...PD
Table, 8...addition circuit, 23, 24.4'l, 4
3゜44-・Selector, 25.16.3B, 40.46
...Register, 33, 34. :16,36,37゜3
8... Carry hold addition circuit, 211.29, 30°31.32... Prefix multiplication circuit, 42... Prefix addition circuit, 45... Carry propagation addition circuit, 42... Quotient digit Generation circuit, 22... reciprocal table. Applicant's representative Patent attorney Takehiko Suzue Figure 5 Figure 6 qj dispute 1 Figure 7 t□ tl t2 t3 No. 9
Figure I4 Carry λ- Tie S〉2'5 Fortune 3゛Figure 10 teu νU惺MP*2金1ζ-)-434-iπi 1 Tie S>ri f93

Claims (1)

【特許請求の範囲】[Claims] 与えられた除数の上位数ピットよシ選択した値を上記除
数および与えられ九被除数にそれぞれ乗じて前記除数が
取シ得る値の範囲を狭く限定してなる新しい除数と被#
数とを求め、この新しく求められた被除数を部分剰余の
初期値として、前記新しく求められた除数と上記部分剰
余とから前記与えられた除数および被除数の商をその上
位桁から順次1桁づつ求めてなる仁とを特徴とする高基
数除算方式。
A new divisor and dividend are obtained by multiplying the above divisor and the given dividend by the value selected by the upper number pit of the given divisor to narrowly limit the range of values that the divisor can take.
Using this newly found dividend as the initial value of the partial remainder, calculate the quotient of the given divisor and dividend from the newly found divisor and the partial remainder one by one, starting from the most significant digit. A high-radix division method that is characterized by
JP2500782A 1982-02-18 1982-02-18 High radix division system Pending JPS58142441A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2500782A JPS58142441A (en) 1982-02-18 1982-02-18 High radix division system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2500782A JPS58142441A (en) 1982-02-18 1982-02-18 High radix division system

Publications (1)

Publication Number Publication Date
JPS58142441A true JPS58142441A (en) 1983-08-24

Family

ID=12153874

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2500782A Pending JPS58142441A (en) 1982-02-18 1982-02-18 High radix division system

Country Status (1)

Country Link
JP (1) JPS58142441A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0158530A2 (en) * 1984-04-09 1985-10-16 Fujitsu Limited Nonrestoring divider
JPH01233520A (en) * 1988-03-14 1989-09-19 Fujitsu Ltd High cardinal number irreparable type dividing device
JPH0273422A (en) * 1988-07-18 1990-03-13 Lsi Logic Corp Digital dividing method and divisor conversion means
JPH04500574A (en) * 1989-02-17 1992-01-30 ディジタル イクイプメント コーポレイション Optimized division circuit

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0158530A2 (en) * 1984-04-09 1985-10-16 Fujitsu Limited Nonrestoring divider
JPH01233520A (en) * 1988-03-14 1989-09-19 Fujitsu Ltd High cardinal number irreparable type dividing device
JPH0273422A (en) * 1988-07-18 1990-03-13 Lsi Logic Corp Digital dividing method and divisor conversion means
JPH04500574A (en) * 1989-02-17 1992-01-30 ディジタル イクイプメント コーポレイション Optimized division circuit

Similar Documents

Publication Publication Date Title
JP3541066B2 (en) Method and apparatus for performing division and square root calculations in a computer
US7921149B2 (en) Division and square root arithmetic unit
JP2000259394A (en) Floating point multiplier
US9785407B2 (en) Data processing apparatus having combined divide-square root circuitry
Kaur et al. Vhdl implementation of non restoring division algorithm using high speed adder/subtractor
JP2972498B2 (en) Automatic logic circuit design method, system and device thereof, and multiplier
JPH0833817B2 (en) Radix 16 divider
JPH10124298A (en) Constant multiplier, method and device for automatically generating constant multiplier, and storage medium stored with automatic constant multiplier generating program
JPH05204611A (en) Method and device for executing pre-scale type division
JPS58142441A (en) High radix division system
JPH0273422A (en) Digital dividing method and divisor conversion means
Antelo et al. Computation of/spl radic/(x/d) in a very high radix combined division/square-root unit with scaling and selection by rounding
JP2857505B2 (en) Division device
US7702715B2 (en) Division arithmatic unit of variable radix
JP3660075B2 (en) Dividing device
JP2518532B2 (en) Subtractor shift type divider
JP3517162B2 (en) Division and square root arithmetic unit
JP4293665B2 (en) Remainder multiplier
JPH0368415B2 (en)
JPS6345633A (en) Multiplication circuit
JP2995721B2 (en) Division device and division method
JPS63254525A (en) Dividing device
JPH08249161A (en) Cubic root arithmetic unit
JPS60247736A (en) Dividing circuit
King Prime appearance in elliptic divisibility sequences