[go: up one dir, main page]

JPH02118725A - 対数関数演算装置 - Google Patents

対数関数演算装置

Info

Publication number
JPH02118725A
JPH02118725A JP63271013A JP27101388A JPH02118725A JP H02118725 A JPH02118725 A JP H02118725A JP 63271013 A JP63271013 A JP 63271013A JP 27101388 A JP27101388 A JP 27101388A JP H02118725 A JPH02118725 A JP H02118725A
Authority
JP
Japan
Prior art keywords
register
pseudo
adder
stl
multiplication
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.)
Granted
Application number
JP63271013A
Other languages
English (en)
Other versions
JP2822399B2 (ja
Inventor
Misayo Nakayama
中山 美小夜
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.)
NEC Corp
Original Assignee
NEC Corp
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 NEC Corp filed Critical NEC Corp
Priority to JP63271013A priority Critical patent/JP2822399B2/ja
Priority to US07/428,301 priority patent/US5041999A/en
Priority to DE68926563T priority patent/DE68926563T2/de
Priority to EP89120128A priority patent/EP0366155B1/en
Publication of JPH02118725A publication Critical patent/JPH02118725A/ja
Application granted granted Critical
Publication of JP2822399B2 publication Critical patent/JP2822399B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F1/00Details not covered by groups G06F3/00 - G06F13/00 and G06F21/00
    • G06F1/02Digital function generators
    • G06F1/03Digital function generators working, at least partly, by table look-up
    • G06F1/0307Logarithmic or exponential functions
    • 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/544Methods 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 for evaluating functions by calculation
    • G06F7/556Logarithmic or exponential functions
    • 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/483Computations with numbers represented by a non-linear combination of denominational numbers, e.g. rational numbers, logarithmic number system or floating-point numbers
    • G06F7/4833Logarithmic number system

Landscapes

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

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は計算機における対数関数演算装置に関するもの
である。
〔従来の技術〕
対数関数は、科学技術計算を取り扱う計算機にとって、
是非とも備えなければならない機能の一つである。
この関数を演算する方式としてはティラー級数展開 lo g (1+x )=x−x 2/2+x 3/3
−− +x ”/n −(1)や、連分数展開や、チエ
ビシエフ級数展開などによる有理関数近似によって求め
る方式がある。しかし、いずれの方式も、多くの乗算あ
るいは除算を必要とするため、演算時間が長くなシ、シ
かも演算精度も満足すべきものが得られないという欠点
があった。
また、マイクロプログラム制御の計算機に適した対数関
数演算方式としてSTL (5equentialTa
ble Look up )がある(−松信著:初等関
数の数値計算、教育出版)。STLは、加算、減算及び
右シフトによって演算できるため、高速な乗算器を持た
ない計算機では有効である。
ここでSTLの演算の原理及び゛算法について一般的な
説明をする。
CSTLの演算原理〕 はじめに2進法で1桁の精度で対数関数tog。
(1+x)を求めるためのSTLの演算原理を以下に示
す。入力数Xは定数r、と数列(a、)を用いると。
以下のように表わせる。
臥下弦日 X士(1+2−’ ) 0X (1+2−2)al・・
・×(1+2−(n−1))a(n−1)×ε    
    (2)ak= (+1 、 o )     
       (3)これよシ。
tag、 (1+x ) =aoX To+a1Xγ1
+・・・+1n−4xγn−1+ε      (4)
−k(5) ただし、γ□=Log、 (1+2  )となる。
このとき、Xから(輸)を求めることを疑似除算。
(a、)からZo g e (1+ x )を求めるこ
とを疑似乗算という。
W=ψ□−ξk(6) W2Oならば。
ak= + 1                  
 (7)ψに+1=W(8) ξ  =ξ+2−に×ξk(9) k−Hk β=fi +r   (r =Log  (1+2−k
))   00に+1    k    k     
  kw (0ならば。
ak=Oαつ ψl(+1 =ψ、                
  (2)ξに+1=ξk(13 βに+1=βk                  
  a→これらを繰シ返し実行することで、ψえをOに
近づけてゆくとT togo(1+ x ) =β、が
求められる。
このとき1式(6)は(ak)を決定するだめの疑似除
算である。また1式(8) 、 (9) 、 (1(e
 、αり、α3 、 (14)は(楡)からtoge(
1+x )を求めるための疑似乗算である。
次にSTLの算法について説明する。
〔1〕xo=x(1/2≦x〈1)、yo=l、Lo=
0を初期値とする。
[2) k=o + 1.2、、−、 n−tに対して
、〔3〕を反復する。
[3] w=xk−yk              
   HW≧0ならば 輸=+t                 (161
χに+1=W(17) ”k+I = yk+ ykx 2−k(IIICzk
+1−9k+γk               (t
lWくOならば ak=o                     
Hk+1ke1) ”k+1=”k                  
(2)k+1   k               
   りただし、rkは(5)式を満たす定数である。
[4] toge(1+x)=zkが得られる。
従来例 第2図は従来のSTLの算法による対数関数演算装置で
、第3図のフローチャートに示す手順で行われる。
第2図において、207,208,209は3種の2進
数の変数xk、yk、殊を格納するレソスタ。
203は値を任意桁だけ右シフトできるバレル・シフタ
、202は(5)式の定数rkを発生するn語のROM
 、 204は入力Aの内容と入力Bの内容を加算また
は減算する加減算器である。206はSTLのループ回
数を制御するカウンタで、201はバレル・シフタ20
3のシフト桁数の指定と、 ROM202の定数番号に
の指定を行なうカウンタである。
次に、第2図の動作を第3図のSTLの算法に基′づい
て説明する。
〔1〕レジスタ207.208,209にx、y。
z(301,302,303)を、カウンタ201.2
06にO,n(305,304)を初期値として与える
〔2〕 カウンタ206の値がnから1までのnステッ
プの間1次のステップ〔3〕を繰り返す。
〔3〕レソスタ208には(8)、(イ)式の変数へが
、レジスタ207には(9) 、 (13式の変数xk
が、レジスタ209にはα0.α1式のzkが、レジス
タ205には(6)式の変数Wが、カウンタ201には
変数kが、カウンタ206には変数nが格納されている
バレル・シフタ203はバスから出力さレル値を、カウ
ンタ201が示す数だけ右シフトすることにより、2 
倍する。
レジスタ208からyを入力Aに転送し、レジスタ20
7からXを入力Bに転送し、加減算器204で減算し、
演算結果をレジスタ205に転送する(306)。
もし、レジスタ205の符号桁の値が正のときは(30
7)、以下の演算を行う。まず、前記演算結果であるレ
ジスタ205の内容をレジスタ208に転送する(30
9)。次に、レジスタ207からXを入力Aに転送し、
Xをバレル・シフタ203でカウンタ201が示す値だ
け右シフトしたものを入力已に転送し、加減算器204
で加算し、演算結果をし・シフタ207に転送する(3
08)。次に、レジスタ209から2を入力Aに転送し
、カウンタ201が示すROM 202の値を入力Bに
転送し、加減算器204で加算し、演算結果をレジスタ
209に転送する(310)。
もし、レジスタ205の符号桁の値が負のときは(30
7)、何もしない(311、312、313)。
〔4〕上記ステツゾをn回実行した後(316)、レジ
スタ209の内容はzn=Loge(1+y/x )と
なる(317)。
従来は9以上の手順をマイクロプログラムで行っていた
。このとき、上記ステップ〔3〕の1回の処理にαクロ
ックかかるとすると、全体の演算処理にはnXαクロッ
ク要し、膨大な時間がかかつていた。
〔発明が解決しようとする課題〕
上述した従来の対数関数演算装置には以下の問題点があ
る。
問題点1:演算結果の精度が良くない。
(5)式のrkはkが増大するに従って小さくなるため
+ rkの有効数字が減少し、 LSB付近に丸め誤差
が蓄積される。
また1式(8) 、 (9) 、 (10、口2 、 
(13、α→は固定小数点数での算法であるため、 x
 、 log。(1+x )が浮動小数点数の場合、X
を固定小数点数に変換してからlog (1+x)を求
め、浮動小数点数に変換しなければならない。このとき
、Xが小さい数である場合、固定小数点数への変換過程
で有効数字が大幅に減少してしまう。
問題点2:実行時間が長い。
STLの処理(式(6) 、 (8) 、 (9) 、
 (10、(ロ)、@1.α4)をすべてマイクロプロ
グラムで行っていたため。
STL 1回の処理にαクロックかかるとすると、全体
ではn×αクロック要し、膨大な時間がかかっでしまう
本発明は従来のもののこのような課題を解決しようとす
るもので、精度の良い演算を短時間で行うことのできる
対数関数演算装置を提供するものである。
〔課題を解決するための手段〕
本発明による対数関数to ge(1+y/x)を演算
する対数関数演算装置は、 2 X log。(1+2
律)またはtoge (1+2  ) (ただしに=0
.1.・−・、m−1)をに=m−1からに=oまで逐
次発生できる定数発生器(104)と、第1のレジスタ
(110)と。
第2のレジスタ(115)と、前記第1のレジスタの内
容をに桁(ただしに==1.2.−、m−2,m −1
、m:)’r:けシフトができる1個のバレル・シフタ
(105)と、前記第1のレジスタの内容に、前記バレ
ル・シフタの内容か、0か、または前記定数発生器の内
容を加算または減算し、前記第1のしシフタに出力する
第1の加減算器(108)と、前記第2のレジスタの内
容に、0か、または前記第1のレジスタの内容を加算ま
たは減算し、前記第2のレジスタに出力する第2の加減
算器(113)と、前記第2のレジスタの符号を入力で
き、かつ第1の加減算器と第2の加減算器が行う演算が
加算か転送か減算かを制御できるm段のFirst−I
nLast−Out (以下FILOと略す)方式のス
タック(116)と、除算器(117)とを含んで構成
される。
すなわち従来の対数関数演算装置は、1個のバレル・シ
フタと1個の加減算器を用いて、 STLのループの処
理をマイクロプログラムで求めるのに。
対し2本発明は1個のバレル・シフタと2個の加減算器
とFILO方式のスタックと除算器を用いて。
疑似除算と疑似乗算を区別して行うことにより。
従来例より高速で、かつ従来例より高い精度で。
Log (1+y/x)を求めるものである。
〔本発明の演算の原理と算法〕
はじめに数式を用いて2本発明の対数関数の演算原理を
説明する。
本発明の演算原理はSTLと同様であυ22本発明算法
はSTL算法の変形と言える。
STLでは疑似除算と疑似乗算を同時に実行しているが
1本発明ではまず疑似除算を行い、その後。
疑似乗算を行う。
STLでは(2〕式のεを無視し、Yo=1を初期値と
して疑似乗算を行っていたため、2進n桁の精度を得る
にはnステップの演算が必要である。ここで、(2)式
でε<2nであるため、2進数で2n桁の精度で求める
ならば(1)式よ#)togo(1+ε)嬌εと近似で
きることを利用し、これを初期値として疑似乗算を行う
そのため、疑似除算、疑似乗算合わせてのステップ数は
STL同様で約nステップ必要である。
本発明では、精度良く使用するために、疑似除算で1ス
テツプごとにXkを1桁下位の桁にずらしながら演算す
る。また、疑似乗算は疑似除算とは逆の順序で、Y、を
1桁上位の桁にずらしながら行う。
とし、に=mからに=i+14で反復される。
Yk−1=Y’に+ ak             
  (財)このように、Xと1が浮動小数点数でありy
/x<1のときはe log (1+y/x ) # 
y/xであるため、疑似除算のステップを途中から始め
、疑似乗算のステップを途中で終わらせる。その結果。
桁合わせKよる有効数字の大幅減少を防ぎ、また性能も
向上できる。
算法 次に演算の算法について説明する演算結果の精度を2進
数でn(=2m)桁とする。
CD Log (1+ y/ x )を求めるに当り、
Xとy(0≦y<x<+■)を入力する。
(2) x =2−’ XX (1≦X(2、lは整)
、Y=2−JXY(1≦Y(2,jは整数)としてi 
  j、X、Yを決定し、x、=X+ yj=Y+ 1
=j−tを初期値とする。
(3) k=i 、 i+1 r i+2、−+ m−
1に対して、〔4〕全反復する(疑似除算を行う)。
[4) W=2Yk−X、             
    (ハ)W2Oならば ak=+1      
       (ト)Yk+1=Xk+2 XXk  
     @Yk+1:WtA w<oならば ak=0              
翰Xk+、=Xk             (7)Y
k+1=2XYkOめ 〔5〕Ym=Ym/Xmを疑似乗算の初期値とする。 
  0■[6] k=m、m−1,m−2,−−−、j
+1に対して、〔7〕を反復する(疑似乗算を行う)。
〔7〕殊=+1ならば、Yk+、=(Yk+Fk)/2
(33ak−0ならば、Yk+、=Yk/2(ロ)(8
) Lage(1+y/x )=Yiが得られる。
〔実施例〕
次に1本発明を図面を参照して説明する。
第1図は本発明に従って対数関数togo(1+y/x
)を演算する演算装置の第1の実施例の構成図である。
第1図において、110,115は2種の変数X、yを
格納するレジスタ、108は式に)、@。
(7)、0■で加算または転送を行う加減算器、105
はレジスタ110の内容を任意の桁数だけ右シフトでき
るバレル・シフタ、109は加減算器108の出力をL
/2倍してレジスタ110に出力するシフタ、114は
加減算器113の出力を2倍してレジスタ115に出力
するシフタである。104は(至)式の定数Fkを格納
するm語の定数発生器。
103はバレル・シフタ105のシフト数およびROM
 104のアドレスを制御し、対数部の演算を行う対数
部演算器、116は数列(a、)を蓄積でき、かつ加減
算器108で行う演算が加算か転送かを制御するか、ま
たは、加減算器113で行う演算が加算か転送か減算か
を制御するためのFILO方式のスタックである。10
1は入出力データおよび演算の中間結果を転送するため
のパス、102はレジスタ110の値を転送するための
パスである。
次に、第1図の動作を説明する。
〔1〕2進浮動小数点数で表現されたx、y(0≦y〈
x<+Oo)の値をパス101から入力し、f=2−1
XX(1≦X(2、iは整数) r y=2’″jxy
(t≦Y<2゜jは整数)となるような指数部の補数1
==j−4を対数部演算器103で求める。
〔2〕xの仮数部Xをレジスタ110に、yの仮数部Y
をレジスタ115に入力する。
〔3〕対数部演算器103からの出力値kをi、i+1
、i+2.・・・2m−1とインクリメントしながら。
〔4〕を反復する。
〔4〕まず、レジスタ110のXkの値をバス102を
介してレジスタ111の入力Cに転送し、レジスタ11
5のY、の値をバス101を介してレジスタ112の入
力りに転送する。加減算器113は入力りから入力Cを
減算し、レジスタ115に書き込む。
次に、レジスタ115の符号桁の反転値をスタック11
6に入力(7°ツシユ)する。
次に、レジスタ110の内容をパス102を介してレジ
スタ106に転送する。また、レジスタ110の内容を
バス101を介してバレル−シフタ105に転送し、対
数部演算器103の示す値だけシフトしたものをレジス
タ107に転送する。
もし、レジスタ115の符号桁が正ならば。
レジスタ106の内容とレジスタ107の内容を加算し
、レジスタ110に転送する。
もし、レジスタ115の符号桁が負ならば。
レジスタ112の内容と0を加算したものをシフタ11
4で2倍にしてレジスタ115に転送する。
〔5〕  レジスタ115のYmをレジスタ110のX
mで、除算器117を使って除算した値(Ym=Yrl
y/′Xm)を疑似乗算の初期値Ymとして、レジスタ
110に転送する。そして、対数部演算器103の出力
値kをに=m 、 m −1、m−2、−1+1とデク
リメントしながら、〔6〕を反復する。
〔6〕 レジスタ110のYkの値をバス101を介し
てレジスタ106の入力Aに転送し、同時に。
対数部演算器103が示す定数発生器104の内容をレ
ジスタ107の入力Bに転送する。
スタック116からポツプされた値が正ならば、加減算
器108は入力Aと入力Bを加算し。
シック109で1/2倍してからレジスタ110に書き
込む。
スタック116からポツプされた値が負ならば、加減算
器108は入力Aと0を加算し、シフタ109で1/2
倍してからレジスタ110に書き込む。
〔7〕レジスタ110にtage(1+y/x )の仮
数部が得られる。
このように、ステップ〔4〕と〔6〕の演算は1クロツ
クで処理でき、全体の処理時間は2(m−1)≦ロクロ
ック要する。
以上述べたように本実施例では、バレル・シフタ1個と
加減算器2個とスタックと除算器を用いて、高速かつ高
精度にtoge(1+y/ x )を演算できる。
上述した第1の実施例では浮動小数点を扱うが。
同一のハードウェアで固定小数点数も扱える。変数1を
1−〇に固定すればよい。
算法は1次のように変形できる。
(1) X、=x 、 Yi=y (0≦y < x 
<+” )を初期値とすす る。(x、 、 y、を固定小数点数にする)〔2〕変
数iをi = Oとする。
〔3〜6〕 (浮動小数点の算法と同様)[7] to
ge (1+y/x )が得られる。
〔発明の効果〕
以上説明したように1本発明による対数関数演算装置は
1個の/’Pラレル・シフタと、2個の加減算器と、 
FILO方式のスタックと、除算器を用いて。
疑似除算と疑似乗算を区別して行うことによシ。
次の2つの効果を有する。
〔効果1〕 演算結果の精度が良い。すなわち。
本発明の算法は、常に有効数字が最大になるように桁合
わせを行いながら演算するために、演算精度が良く、ま
た+ x + y + toge(1+y/x )が浮
動小数点数の場合、y/Xが小さい数であっても、有効
数字が減少しない。
〔効果2〕 演算時間を短縮できる。すなわち。
本発明の演算装置はSTLの処理をノ・−ドウエアで実
現している。このため、マイクロプログラムで制御する
のに比べ、1/(STLの1回のループにかかった時間
)の時間で演算ができる。
【図面の簡単な説明】
第1図は本発明の一実施例の構成図、第2図は従来のも
のの一例の構成図、第3図は第2図に示すもののフロー
チャートである。 107.110,111,112,115・・・レジス
タ、109,114・・・シフタ、116・・・スタッ
ク、117・・・除算器。 第3Il!! 手続補 正  書 (自発) 平成2年7 月:22日 特許庁長官 吉 1)文 毅 殿 1、事件の表示 昭和63年特許願第271013号 2、発明の名称 函 対数(2)数演算装置 3、補正をする者 事件との関係    特許出願人 名称  (423)日本電気株式会社 4、代理人〒105 住所 東京都港区西新橋1丁目4番10号5、補正の対
象 6、補正の内容 明細書の記載全欠のとおり訂正する。 1)第2頁末行を次のように変更する lr togo(1+y )=y−y2/2+y3/3
−・・・十y n/n・・・ (1)J2)第3頁15
行目のr (STLの演算原理〕」を削除する。 3)同頁17行目及び18行目のrxJをryJに変え
る。(2箇所) 4)第4頁1行目、5行目、9行目及び10行目のrx
JをryJに変える。(4箇所)5)同頁12行目〜第
5頁8行目:「W=φ2−ξ5・・・疑似乗算である。 」を削除する。 6)第5頁9行目〜第6頁6行目:「次にSTLの乗算
・・・が得られる。」を次のように変更する。 「次にSTL法によりLog (1+y/x) k求め
る場合について説明する。 [1] yO= y(1/2≦yく1)、xo=1.z
o二〇を初期値とする。 [2] k=o 、 1 、2.−・−、n−1に対し
て、〔3〕のル−2′を反復する。 [3] W””、ykXk(6) W2Oならば ak=: +1                  
(7)Xk十、=y                
     (8)x   ==x  +x  X2−k k+1   k  k               
   (9)”k+1=zk+γk         
   α1w (Oならば a k=O(11) )’に+1=’lk             (2)
k++    k(13) k+1   k                  
  o4)ただし、γ2は(5)、式を満たす定数であ
る。 [4] tog、#、(]+y)=殊が得られる。」7
)第7頁3−4行目:rx+y+zJをr XOgy、
3 r zOJに変える。 8)第15頁1行目:[幀よ「αυ」に変える。 9)第15頁1行目から第12行目までを次のように変
更する。 rc4〕W= 2Yk−Xk(16) W2Oならばak:= +I           Q
7)Yk+1”Xk+2xX、、    08)2Yk
+1=2W(19) w<oならば輸=0            (20)
Xk+、=XkQ1) 2Yk+、=4 XY、、       (2の1:5
) Y、=Y、/Xmを疑似乗算の初期値とする。■)
[6:] k二m g m−1* m  2 g ・・
・* F + 1に対して。 〔7〕を反復する(疑似乗算を行う)。 〔7〕ak=+1ならば、Yk+、=(Yk十Fk)/
2(24)輸=0ならばt Yy+ 1 = Yi2 
     (25)[8) Z o g (1+y/x
 ) =Y tが得られる。」10)同頁末行から第1
5頁1行の「(ハ)、@、(ト)。 0■」を「0Q、α呻、Qρ、に)」に変える。 手続補正書(方式) %式% 1、事件の表示 昭和63年特許願第271013号 2、発明の名称 対数関数演算装置装置 3、補正をする者 事件との関係    特許出願人 氏   名  (423)日本電気株式会社4、代理人
 〒105 住 所  東京都港区西新橋1丁目4番10号手続補正
書(自発) 平成2年1月22日 特許庁長官 吉 1)文 毅 殿 1、事件の表示 昭和63年特許願第271013号 2、発明の名称 対数関数演算装置 3、補正をする者 事件との関係    特許出願人 氏   名  (423)日本電気株式会社4、代理人
 〒105 住 所  東京都港区西新橋1丁目4番lo号5、補正
命令の日付 自発補正 6、補正の対象 ■平成2年1月22日付手続補正占の発明の名称の欄 7 補正の内容 5、補正の対象 明細書の発明の詳細な説明の欄

Claims (1)

  1. 【特許請求の範囲】 1、計算機における対数関数log_e(1+y/x)
    を演算する装置であって、 2^k×log_e(1+2^−^k)またはlog_
    e(1+2^−^k)(ただしk=0、1、…、m−1
    )をk=m−1からk=0まで逐次発生できる定数発生
    器と、 第1のレジスタ及び第2のレジスタと、 前記第1のレジスタの内容をk桁(ただしk=1、2、
    …、m−2、m−1、m)だけシフトができる1個のバ
    レル・シフタと、 前記第1のレジスタの内容に、前記バレル・シフタの内
    容か、0か、または前記定数発生器の内容を加算または
    減算し、前記第1のレジスタに出力する第1の加減算器
    と、 前記第2のレジスタの内容に、0か、または前記第1の
    レジスタの内容を加算または減算し、前記第2のレジス
    タに出力する第2の加減算器と、前記第2のレジスタの
    符号を入力でき、かつ第1の加減算器と第2の加減算器
    が行う演算が加算か転送か減算かを制御できるm段のF
    ILO(First−In Last−Out)方式の
    スタックと、除最器とを含むことを特徴とする対数関数
    演算装置。
JP63271013A 1988-10-28 1988-10-28 対数関数演算装置 Expired - Fee Related JP2822399B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP63271013A JP2822399B2 (ja) 1988-10-28 1988-10-28 対数関数演算装置
US07/428,301 US5041999A (en) 1988-10-28 1989-10-30 Logarithmic function arithmetic unit including means for separately processing pseudo division and multiplication
DE68926563T DE68926563T2 (de) 1988-10-28 1989-10-30 Arithmetische Einheit für logarithmische Funktionen, versehen mit Mitteln zur gesonderten Verarbeitung von Pseudodivision und -multiplikation
EP89120128A EP0366155B1 (en) 1988-10-28 1989-10-30 Logarithmic function arithmetic unit including means for separately processing pseudo division and multiplication

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP63271013A JP2822399B2 (ja) 1988-10-28 1988-10-28 対数関数演算装置

Publications (2)

Publication Number Publication Date
JPH02118725A true JPH02118725A (ja) 1990-05-07
JP2822399B2 JP2822399B2 (ja) 1998-11-11

Family

ID=17494194

Family Applications (1)

Application Number Title Priority Date Filing Date
JP63271013A Expired - Fee Related JP2822399B2 (ja) 1988-10-28 1988-10-28 対数関数演算装置

Country Status (4)

Country Link
US (1) US5041999A (ja)
EP (1) EP0366155B1 (ja)
JP (1) JP2822399B2 (ja)
DE (1) DE68926563T2 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0554035A2 (en) * 1992-01-27 1993-08-04 Mitsubishi Denki Kabushiki Kaisha Solid state color video camera
US6788342B1 (en) 1992-01-27 2004-09-07 Mitsubishi Denki Kabushiki Kaisha Color video camera for generating a luminance signal with unattenuated harmonics
JP2007023418A (ja) * 2005-07-15 2007-02-01 Shoei Co Ltd フルフェイス型ヘルメット
WO2020090025A1 (ja) * 2018-10-31 2020-05-07 富士通株式会社 演算処理装置および演算処理装置の制御方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6055553A (en) * 1997-02-25 2000-04-25 Kantabutra; Vitit Apparatus for computing exponential and trigonometric functions
US6567832B1 (en) * 1999-03-15 2003-05-20 Matsushita Electric Industrial Co., Ltd. Device, method, and storage medium for exponentiation and elliptic curve exponentiation
US6976043B2 (en) * 2001-07-30 2005-12-13 Ati Technologies Inc. Technique for approximating functions based on lagrange polynomials
US7509363B2 (en) * 2001-07-30 2009-03-24 Ati Technologies Ulc Method and system for approximating sine and cosine functions
US7606850B2 (en) * 2005-03-30 2009-10-20 Lockheed Martin Corporation Method and apparatus for providing a base-2 logarithm approximation to a binary number
US8510360B2 (en) * 2010-06-04 2013-08-13 International Business Machines Corporation Calculating large precision common logarithms
US10838718B2 (en) * 2017-09-28 2020-11-17 Fujitsu Limited Processing device, arithmetic unit, and control method of processing device

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4089060A (en) * 1976-10-15 1978-05-09 Mitchell Donald K Digital logarithmic apparatus
JPS6051733B2 (ja) * 1977-11-19 1985-11-15 日本電気株式会社 指数関数演算装置

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0554035A2 (en) * 1992-01-27 1993-08-04 Mitsubishi Denki Kabushiki Kaisha Solid state color video camera
EP0554035A3 (ja) * 1992-01-27 1995-02-15 Mitsubishi Electric Corp
US5450124A (en) * 1992-01-27 1995-09-12 Mitsubishi Denki Kabushiki Kaisha Color video camera for generating a luminance signal with unattenuated harmonics
US5652620A (en) * 1992-01-27 1997-07-29 Mitsubishi Denki Kabushiki Kaisha Color video camera
US5966170A (en) * 1992-01-27 1999-10-12 Mitsubishi Denki Kabushiki Color video camera for generating a Luminance signal with unattenuated harmonics
US6108038A (en) * 1992-01-27 2000-08-22 Mitsubishi Denki Kabushiki Kaisha Color video camera for generating a luminance signal with unattenuated harmonics
US6788342B1 (en) 1992-01-27 2004-09-07 Mitsubishi Denki Kabushiki Kaisha Color video camera for generating a luminance signal with unattenuated harmonics
JP2007023418A (ja) * 2005-07-15 2007-02-01 Shoei Co Ltd フルフェイス型ヘルメット
WO2020090025A1 (ja) * 2018-10-31 2020-05-07 富士通株式会社 演算処理装置および演算処理装置の制御方法
JPWO2020090025A1 (ja) * 2018-10-31 2021-09-24 富士通株式会社 演算処理装置および演算処理装置の制御方法

Also Published As

Publication number Publication date
DE68926563D1 (de) 1996-07-04
DE68926563T2 (de) 1997-01-23
JP2822399B2 (ja) 1998-11-11
EP0366155A2 (en) 1990-05-02
US5041999A (en) 1991-08-20
EP0366155B1 (en) 1996-05-29
EP0366155A3 (en) 1991-11-13

Similar Documents

Publication Publication Date Title
USRE39385E1 (en) Method and apparatus for performing mathematical functions using polynomial approximation and a rectangular aspect ratio multiplier
US7395304B2 (en) Method and apparatus for performing single-cycle addition or subtraction and comparison in redundant form arithmetic
US5457804A (en) Accumulating multiplication circuit executing a double-precision multiplication at a high speed
US5241493A (en) Floating point arithmetic unit with size efficient pipelined multiply-add architecture
GB2220090A (en) Cordic complex multiplier
JPH03204720A (ja) 初等関数演算装置
US5184318A (en) Rectangular array signed digit multiplier
JPH02118725A (ja) 対数関数演算装置
EP0297588B1 (en) Trigonometric function arithmetic processor using pseudo-division
JP2504102B2 (ja) 逆三角関数演算装置
Ananda Mohan Specialized residue number systems
JP2508784B2 (ja) 指数関数演算装置
GB2549153B (en) Apparatus and method for supporting a conversion instruction
US6349317B1 (en) Efficient radix-4 CORDIC vector rotators and computers of sine and cosine functions
Kostin et al. Enveloping of Riemann’s zeta function values and curious approximation
US6055553A (en) Apparatus for computing exponential and trigonometric functions
JP2972326B2 (ja) 平方根計算装置
Kuhlmann et al. A novel CORDIC rotation method for generalized coordinate systems
Sathyavathi et al. Single-Precision Floating-Point Addition Under HUB
JPH06105421B2 (ja) 逆三角関数演算装置
JPH0414173A (ja) 固定小数点演算器
JPH08263269A (ja) 浮動小数点演算回路
Jaiswal et al. Dual-mode double precision division architecture
JPH02266427A (ja) 乗算器
JP2546916B2 (ja) 三角関数演算装置

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees