From 493d3b1ac0cd61733e809ff8ab439ed91d89d9f6 Mon Sep 17 00:00:00 2001 From: James Stine Date: Tue, 21 Jun 2022 15:07:24 -0500 Subject: [PATCH] Add hex output in bad but okay way --- pipelined/srt/Makefile | 6 +- pipelined/srt/qslc_r4a2b | Bin 0 -> 16064 bytes pipelined/srt/qslc_r4a2b.c | 190 +++++++ pipelined/srt/qslc_r4a2b.tv | 1024 +++++++++++++++++++++++++++++++++++ pipelined/srt/stine/notes | 2 +- 5 files changed, 1220 insertions(+), 2 deletions(-) create mode 100755 pipelined/srt/qslc_r4a2b create mode 100644 pipelined/srt/qslc_r4a2b.c create mode 100644 pipelined/srt/qslc_r4a2b.tv diff --git a/pipelined/srt/Makefile b/pipelined/srt/Makefile index c3c69df8..63146339 100644 --- a/pipelined/srt/Makefile +++ b/pipelined/srt/Makefile @@ -1,4 +1,4 @@ -all: exptestgen testgen qslc_r4a2 +all: exptestgen testgen qslc_r4a2 qslc_r4a2b sqrttestgen: sqrttestgen.c gcc sqrttestgen.c -o sqrttestgen -lm @@ -15,5 +15,9 @@ qslc_r4a2: qslc_r4a2.c gcc qslc_r4a2.c -o qslc_r4a2 -lm ./qslc_r4a2 > qslc_r4a2.sv +qslc_r4a2b: qslc_r4a2b.c + gcc qslc_r4a2b.c -o qslc_r4a2b -lm + ./qslc_r4a2b > qslc_r4a2b.tv + clean: rm -f testgen exptestgen qslc_r4a2 diff --git a/pipelined/srt/qslc_r4a2b b/pipelined/srt/qslc_r4a2b new file mode 100755 index 0000000000000000000000000000000000000000..f719bbf471bfc1094ffe5bdd0d6441d7c2426e9d GIT binary patch literal 16064 zcmeHOeQ;FO6~CLXB*BmkpF;d#p-TBMO9COq4>SZ87AFYApwem|o9rgpHrb848wm~{ zv8k05VXVy9I8*$n(-B*(`GfIeoJx=?GPM@Qahlp0jSkFg@B;`~SB<>>&b{}%ygUa;WcV^8M2bI&=ybMJZQ-uH6vzWa9gt5#;YT!JZE+$czDFI1WeHvE7KLR5$)Vm!_> z#Fb(!)MIc|XtR^|W+?;nT#d^HC%>v#9Zj>AU`)5wCv>T zbzscdZnftnXOLH(EnY(a5Ocd9h4MtPG>a#e4(Ba8FKwE$ABp7Gt^K;SA9KNY%{iZx z6FQdYd^&jvY(&hRevR7C$yX>h!JOJ4y7go-^j+Fpr~Lw6S&eY=8xWG`kIZ#H`ajpNc_jH168`dK)> zMaJryHAW~L4>v>-$#A@8&GM#LG+Yy`Z3+v+XlRZ_jYKjSPa1}h&=G+xvCRP2wFQm( zNHo|KxtDy=vvu0HIT(qGmeyolV=ykDZ>A0|r&66Nt*TmGx!hRjTX?PD;-WIExX`zl z;?6=Zxp4Qn@W*{c1muYy{Yv)c3FSSrl;|rp-cuc8OUTHlYvYIG8xEZAd`XrZ3fPHuKX&{@kwrAm-`tmdssGXE`Av50rTp;V4fL*@Qq{{=>O0(&4;8+i-dvx<17QcPoh@{ ziHCfbL*8JQOFv!RdHlzL&W{2eC(f*^@t5|L9u4eXO7l~5$|4}o*ZVv(elFvIej_x2 z-R^2)*QS$G5c*UqBlJA8&oiT)fQK~)ey+qXo<(f#g+S+-z@g7>3LHvj2V8Fij$cSl zMg=E*s)Agzuii5QakD>NU;9!ri`MxY0v$_Vqhg@*WHLXnYv~&hQhB)UR0yXx+<$`X zx(mL;{QEbACg@%;F;i0b#6P++1fYYDfFk#dA5&#sz8W(Pt75+_3rZbebHO* zIo&hhnb(8x9+Ad5Z`KDq(Vo;x7fchS?L8zOtM@F}(0k9Em#7ow`%iTG&wpHSumnP? z9k>qvdDkxg`Ro0C+fVgB3oX?O7k}!WYjgwt6KJ7qSPQqo6OpY+4S4R55vA!H-s#?mg6BndHCX-CIfM0`?o}ur@eiU_(SJKuDo6ZjS)PhAJm*Ny9M20q z9qAm;1J7cu(b*H&J^6kDVKj{yFT}1fS%I#=?kVYfs9eWx#9(^SeHts*Lw#fG52f>N zxMaKM;1b|dvk`ZH9CZNjwhfRCKp6qt1X!;C&)EQjTRkcZCk6p1QGf?+fO7!QFAS-} z1n>Z$GIfAN`h^EA-f6>+=aBXi-~%k(0Jz=;utJIx;8g`EvH`y0kk%7mp8|Y2f9R^u za!4fv*b0COX(oxy3V+mwFW`^_0jd;WzYSo8^ugBv%vXR1Yyc~yR|zmq0qSjlGaS-B z0`%iSp+ee2;(HL0gmN3+!yz>j;AI7vVgpzqRS}>|0nX1GN`e*Ad;)A%fOl+w&p4!U z1h@?V71GyVL0qNqPucJjIHdOoFk1nBW&^kZkYnj(0=N|*Xal6VcHRbjkB<#7k^Kp?5-Q9d}1O20!1-ql0i9eq1GdR#cWdXA;(G1s?a z_vF`srG(4_j-Hcw$MB@Q0z!vB?dsT;-r?EPhn~XoF_)BmIkJiz6gppn^4$03<2LWu zPIrom{y{hm=6Loz4yTu>?xTlsqKCAcLs+v0GG!=}flLN68OUTHlYvYI{$DdN^u|g> zSa8lXpTtpi-ZYQk_~GZK`6-U}f0^c_FR+ePw8L>d6|T#&ucvpm^e(9Pv}v}0=bkam zF7RpK-QY{WX&y84tZBXpe)|9x&frA!Zj{J%Z?$l>6}m3V&&hoV(>>B(j)Ur)O@I|S z@|D?DDAGJ@6ZEv$6$OPW3#Q%dxny&0ySS<7x_OJPnnh~zzY|9zaMF>KAkrsskiQq4 ziQ0M)$9lwCj`oGWpzs%2%g2w+E_NX#8OyU!H{cq6i>2j|{ZTjBGnC0dCIguaWHOM+ zKqdp33}iBp$v`FpnGAf_4Ddcb-q%M<5+d4LNW?PjB_zsM%-UFp*IJTROjCTimU%zn zWm@KyOOcj&Um>lBh}i!51v5r+H(tsR(GrNL3$Lb#crReDw$rcGM9(NDF40_vBc^DX znRXN1dS4)3dq|RNJd|kNM#OsyFIA#w)D5$KzYdJ!?$NkB&AI(6wahrYyph+3yh;_` zWlH%DT_2XeQs%g@Oq6n^w&PcSDf4yc>i=uZ+6V5au|zL!E9{G02j*g8ZiTY0TE6@` z@0<;_t^m$qNLS0yF`&FtWYxFFSNVbC8i1<=gKO?RurgyzVBGSQ{tys z`*2zRGzufdKiYL2ko8O!>t=EdLwuxqUX%Dsg|pxPBHbog`{%iq_a%Osb^oybJPH$# zhshO6qiKMR zP%u5F@r4fj?=}8$eR1+PBH*9b`iS12Eq_~r{uSscPG>zwp!dS=qz|E=mhZ-I+VGJs4D zgWjEJsx#tc!G*Pm$}m<|FI(d`{I{&2-Tf;z-m+}X>g7=5qX@$v(AL0;YGJIZT3flS z%2>N{<*oi2qh?uUmERC8O+nt&U%|WfX{-LXw*FgI?CHPg_XgomFo`;2n=KD~5FjrY zm6A<;PXJ-~hY3a~5i=Tt(GY#KuzD>}p-9wdO@u=NIl{@PO(e9h{1zb;NwgTX_*Mbw zGYow6z)7I5Ax8YfVfg0_L-En|QksSECAKstgSFtvxZ;gmibU~gMT_u7W67|uA=>I| ziN{*P@#GepsvhK!%%L<2p!=r`p;D~f5NSY{~3y$-?Q=jNt-2>%f!EjgXS#kKiYR&geBs6 zijssVN5Ssz{C!Cw=4}73Fi`vK zAMh$9^Ba_Cot^Ff6BfJwPOpM8f6MYz1s*@FK8XUhj>>zGbiA(Q?=78pmYJW0kH_ZE z^YQ(~RNK&aRnm~m{@=hrW0%)kJWuBDXE`z4jqJyH`~!UH`YnIG?(ALwC1{;u)|vg8 zkA<_{UsPyliQ?Q2SHONO(>$AW+&<3(3$?!!@9e(`x_u4jY6bhxwlu0hu@c!ij~!)E z=Owzr{W`EKpr$D;$HUK`Jm@%f<*Cknv>*dzIV$INm6hh=Bl~|OzpNCZ$l>4+#lHca C!tRp* literal 0 HcmV?d00001 diff --git a/pipelined/srt/qslc_r4a2b.c b/pipelined/srt/qslc_r4a2b.c new file mode 100644 index 00000000..94a3a4cd --- /dev/null +++ b/pipelined/srt/qslc_r4a2b.c @@ -0,0 +1,190 @@ +/* + Program: qslc_r4a2.c + Description: Prints out Quotient Selection Table (assumes CPA is utilized to reduce memory) + User: James E. Stine + +*/ + +#include +#include + +#define DIVISOR_SIZE 3 +#define CARRY_SIZE 7 +#define SUM_SIZE 7 +#define TOT_SIZE 7 + +void disp_binary(double, int, int); + +struct bits { + unsigned int divisor : DIVISOR_SIZE; + int tot : TOT_SIZE; +} pla; + +/* + + Function: disp_binary + Description: This function displays a Double-Precision number into + four 16 bit integers using the global union variable + dp_number + Argument List: double x The value to be converted + int bits_to_left Number of bits left of radix point + int bits_to_right Number of bits right of radix point + Return value: none + +*/ +void disp_binary(double x, int bits_to_left, int bits_to_right) { + int i; + double diff; + + if (fabs(x) < pow(2.0, ((double) -bits_to_right)) ) { + for (i = -bits_to_left + 1; i <= bits_to_right; i++) { + printf("0"); + } + if (i == bits_to_right+1) + ; + + return; + } + + if (x < 0.0) + x = pow(2.0, ((double) bits_to_left)) + x; + + for (i = -bits_to_left + 1; i <= bits_to_right; i++) { + diff = pow(2.0, ((double) -i) ); + if (x < diff) + printf("0"); + else { + printf("1"); + x -= diff; + } + if (i == 0) + ; + + } + +} + +int main() { + int m; + int n; + int o; + pla.divisor = 0; + pla.tot = 0; + for (o=0; o < pow(2.0, DIVISOR_SIZE); o++) { + for (m=0; m < pow(2.0, TOT_SIZE); m++) { + /* + 4 bits for Radix 4 (a=2) + 1000 = +2 + 0100 = +1 + 0000 = 0 + 0010 = -1 + 0001 = -2 + */ + switch (pla.divisor) { + case 0: + if ((pla.tot) >= 12) + printf("8"); + else if ((pla.tot) >= 4) + printf("4"); + else if ((pla.tot) >= -4) + printf("0"); + else if ((pla.tot) >= -13) + printf("2"); + else + printf("1"); + break; + case 1: + if ((pla.tot) >= 14) + printf("8"); + else if ((pla.tot) >= 4) + printf("4"); + else if ((pla.tot) >= -6) + printf("0"); + else if ((pla.tot) >= -15) + printf("2"); + else + printf("1"); + break; + case 2: + if ((pla.tot) >= 15) + printf("8"); + else if ((pla.tot) >= 4) + printf("4"); + else if ((pla.tot) >= -6) + printf("0"); + else if ((pla.tot) >= -16) + printf("2"); + else + printf("1"); + break; + case 3: + if ((pla.tot) >= 16) + printf("8"); + else if ((pla.tot) >= 4) + printf("4"); + else if ((pla.tot) >= -6) + printf("0"); + else if ((pla.tot) >= -18) + printf("2"); + else + printf("1"); + break; + case 4: + if ((pla.tot) >= 18) + printf("8"); + else if ((pla.tot) >= 6) + printf("4"); + else if ((pla.tot) >= -8) + printf("0"); + else if ((pla.tot) >= -20) + printf("2"); + else + printf("1"); + break; + case 5: + if ((pla.tot) >= 20) + printf("8"); + else if ((pla.tot) >= 6) + printf("4"); + else if ((pla.tot) >= -8) + printf("0"); + else if ((pla.tot) >= -20) + printf("2"); + else + printf("1"); + break; + case 6: + if ((pla.tot) >= 20) + printf("8"); + else if ((pla.tot) >= 8) + printf("4"); + else if ((pla.tot) >= -8) + printf("0"); + else if ((pla.tot) >= -22) + printf("2"); + else + printf("1"); + break; + case 7: + if ((pla.tot) >= 24) + printf("8"); + else if ((pla.tot) >= 8) + printf("4"); + else if ((pla.tot) >= -8) + printf("0"); + else if ((pla.tot) >= -24) + printf("2"); + else + printf("1"); + break; + default: printf ("X"); + + } + + printf("\n"); + (pla.tot)++; + } + (pla.divisor)++; + } + +} diff --git a/pipelined/srt/qslc_r4a2b.tv b/pipelined/srt/qslc_r4a2b.tv new file mode 100644 index 00000000..b92d81e8 --- /dev/null +++ b/pipelined/srt/qslc_r4a2b.tv @@ -0,0 +1,1024 @@ +0 +0 +0 +0 +4 +4 +4 +4 +4 +4 +4 +4 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +2 +2 +2 +2 +2 +2 +2 +2 +2 +0 +0 +0 +0 +0 +0 +0 +0 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +2 +2 +2 +2 +2 +2 +2 +2 +2 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +4 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +8 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +2 +0 +0 +0 +0 +0 +0 +0 +0 diff --git a/pipelined/srt/stine/notes b/pipelined/srt/stine/notes index 83be0f20..77311b17 100644 --- a/pipelined/srt/stine/notes +++ b/pipelined/srt/stine/notes @@ -9,7 +9,7 @@ D = 197/256 = 0.76953125 Note: Add lg(r) extra iterations due to shifting of computed q q_{computed} = q / radix -./srt4div 0.68359375 0.76953125 4 10 +./srt4div 0.68359375 0.76953125 8 10 r=2 X = 0.10011111