Date: Fri, 14 Sep 2012 17:34:08 -0700 From: Steve Kargl <sgk@troutmask.apl.washington.edu> To: Stephen Montgomery-Smith <stephen@missouri.edu> Cc: freebsd-numerics@freebsd.org, Bruce Evans <brde@optusnet.com.au> Subject: Re: cexp error Message-ID: <20120915003408.GA70269@troutmask.apl.washington.edu> In-Reply-To: <5053CB4F.3030709@missouri.edu> References: <5048D00B.8010401@missouri.edu> <504D3CCD.2050006@missouri.edu> <504FF726.9060001@missouri.edu> <20120912191556.F1078@besplex.bde.org> <20120912225847.J1771@besplex.bde.org> <50511B40.3070009@missouri.edu> <20120913204808.T1964@besplex.bde.org> <5052A923.9030806@missouri.edu> <20120914143553.X870@besplex.bde.org> <5053CB4F.3030709@missouri.edu>
index | next in thread | previous in thread | raw e-mail
[-- Attachment #1 --]
On Fri, Sep 14, 2012 at 07:26:55PM -0500, Stephen Montgomery-Smith wrote:
> On 09/14/2012 01:05 AM, Bruce Evans wrote:
>
> >My tests also determined the exact minimum for all multiples up to the
> >thresholds for "large" multiples in e_rem_pio2.c, except for ld128,
> >to verify that there are enough bits in the special approximations for
> >Pi/2 there, except for ld128. The maximum multiples handled there are
> >2**28*Pi/2, except for ld128 they are 2**45*Pi/2. 2**45 is too many
> >to check exhaustively.
>
> But presumably one could prove a result to this effect using continued
> fractions, right?
>
See attached. It shows a method for determining the number
of needed bits.
--
Steve
[-- Attachment #2 --]
%PDF-1.4
%쏢
17 0 obj
<</Length 18 0 R/Filter /FlateDecode>>
stream
x]Y%Gqqo/Ӫ}!$$@`9BhfҌάʪ:Gse }k=*}{T'^~<~~2N~ۻ0kWz~Cx_sA~~'I
7z'ӯ9mC><+5{EQ+M/{w zQ*e|4۠;Yhxo,?
ZH==8xaA 1M(w18|558M;z'eTXwoݳ{WD4sBMG6 >ݝTKsѝ>\a"ry%JoiL%5P"gπfK?ICNOT5N0=7=npPx;[gg dQz
oʅ3!ڇ`K2*<'/p?j"a.:k!
+<.*5g7heڳt65=?_!
O%lŀ<P}3Dqn@Dz -][>~\H ~F:}_A:;G(
2I"zl
FJv0gLNNg
B:q\_SX>QN) <zE7poqs2h!Q 15 &"q/V0APN=\e(HMzR=59 T#kTQ^"EE|#T0`C <(`f&YE*k0D@8F6ӛ]2uF2XdQBpM4*0EwUUM!@&)A!;AP#DQ#-`T>N,J=`!<kWOd筜~ )a?8'C!}`]Keqnɨ|rIs@$:t1 |Ѡ&V+``pcpf:)a*i"a`!ZD0E V{.Υ|%{E>/XOÉ
:H{5 3*f&d"j&5H -oJ"Ք.ߠy@~տ: 5MhQa4@7Av&ZH$ұPkY
TܚTTvN]lW7jmtd:g`h.
;#Jȳ4ޛd<q6J?!`.4ogwLݜ*<)xU{6`D֥~YD:+7'p)TCdK3Ėou+A2}#
QU4U
YYO^U7rlFb Ecapi;WdmPd+#"ܐln>~XpIuSy$yu+8]n}|*u|ֈ|$ huvƘHGd9Erȝ0HDپ1
AQ#(6Hkkpg`|ug5F`z
tF,a!L3ZR{nLS/c1SrÃ!k| #\dIH2D_?I6Y%JOޕ͐ K5-FTWͣNq5!WZ?W5K
mՐАSjX&BV' [n=`C*|I UK8ꃁ\u.z4M&,t#Cښ*5y%&I_c24Yaj^1(lb.X:&JLÄ(QqzBan@
LL Qy?|Y
K0k2DΑ )Zz,ڿ$Ūk#qDfX9T:4/A$5`CYI6%fIYAX٦Ș5S9x j?,
]H$a g&TՊ;%̾jܤn!>Uw^qyXڐ$[/D>XY16ƤYQz
eEPTc<Fq'%H@3|yxݷ'lҵ0
lg<iD"zpN/ݛݛzw>8RsrC(?^lWB?
8(kjulhGHu
h92T+lNsn]D<2 QşMLeVK6QdNM,]\8V0XbɫM2/G>G=^j; x J
3T&sm&m8nȠ44d Uc h2)Z9J.=;DrQ'DF@dNU^'~hʘ7Xw)Ebw!
lPhO]B?A'Y(_03Cye<E嚈mǻW&@N%3*HEN]C*kU935pA흏k]XH+f+'\'Da26h QOw H>S "cXs@#38{|I6F8>1Q*YN6
<"!Xn
..Y;1Rʀ{)AZȎ(ɚN ̴Nן_yKjv
^E@4{Zƣ 7Cѐd
klQ~9_4<Fxj3Heui.#%Lc qX
bh6Kj~#taA=3^U7HڰHdvYT7oj!$ԃbQIB!~T?U1GJ4;
Y)웒r.m?v5Rc-Ƹ
$BbԲ_/לx[
6&CFgyXu'Q=GXW)5U@ ϿhtilRZ<Đgk@#Cqwo]²2 !]h"/bR{Z=YWJKa}1'T~mzxvkeܞ*ʛʛ %^wxx kict!OJxХ}UÔ}sӲ9{vg?(XQX0y68>/xǂm|Kݒ~6}MRsOKF9X #`gh@k
jBaBz<.nв]46{5axd`x۸.3A/KF84!r?A)I/Sp!Q"U({J\i@zO%LǧITc2*N'Q2ձ'p ~2?L誧"kMٲ1e
bv:#*tgďAf$w_ߢL2@k 6a^$oZ`
VBX7ZŊ(6E0X%I"l.E"HNa)u-=f%wzk%'Y"Wb79M M*u]T{Pǒz2HCvQi6N)oOA2SqnAZ=X&HB ?`Q&{rOFI@ceѤ5jR9ohO_]6<j pv`s/[l/&cUw*V_m'NUKuQLlAȄ,V ,^fV`Ti{Ƞ yK9
mOe͚tqXfZG:0Y i۸)3r [4pTu/]N;ƈ!Jd7$&r
%kzHʳ(jW~B=HKS@X,Iadi1!ӢxEqy"HOv~#4:u.-Siчd7u
tx+To.vHOZDuHޮ
kDMvrs=7<Ѻmlhx툺aRD]Y"F"Ik
"% #vWa6]:G"qr4'Ŝ ;P)
#s aȫ%rg̜AjoɍȥRW#.ŶTB'Gmy
-"ŹY[G5^:m)eܡ&]~#tl!Dc{o-`C
>ﺯz?z96j=u"=.-vkXD|CpP.uT]c1uRm9T ^x lunob!((g; eŔvj1v +mVfF&ݎm4Nqq{f͂0Hsw< )GWUl_
(nN\kڭ,%3D&>b=IO6e#h:@UO|:o?mo=(r)ZGIH[^lɘ/E6o6Ni8+$bB6DFhRU 9qxʀģ9Q$WECwfu0دc@gx@Y0YU( u%=(7p&Ee"NDV,ZPvA&*dY}oONPL
;t@smqԡ߶D(3ĆeXD|kP[SxpK皼7)b$i|h\ÇU46
Q398hr?_a~8upv$j:x`ƒ~ElW
FEiѭr.\F>)_ n{{ޓ,|U&S4{=wķaFr tHҩ9#CSۮRHgqv-`)&]M
. 6NۏX&J@vBD@V>T9#mzݯ_d0#֡b9[*w2zYn`;6=s;*!㊍,~c)DW d|?`p%9{OIac:\6=b-9VCyr#ƑNS0v9\RZvsM/Gr Z{kBKnAZMoHixؽLZպmyJ y"4
l;+
MG-j
{YdG4:{(9B*?{1F:w)MY40Sѹηs/1l,I
JMo 'Z~c9?3ڜrg19̃X SeC}Xill
zX2-`{=^LH_6g(<.4dׯ?bzց%XIndl٦蛣ژTz=vDAØ.Oo"kU]ZݔLˤgT((O^=,E8jC#k~\3T$H??4ֺendstream
endobj
18 0 obj
5827
endobj
39 0 obj
<</Length 40 0 R/Filter /FlateDecode>>
stream
x\[%c4)#KD qg&X!a}lw.>3"iv|U>Y/^<P7\?o.~
˧ׯ~%RB;}BɃsT̬WO/>=IQ*tiGrZ.Px.u|r3x\O&W_K=;~9YJӧ,<m_#_.>}Hr8\b1r(]<x_]< /o=z@D$e,
KD+T0ǟKd_@D
$)IK2IVJM$E/hw:"/{c/RJ`u|wU)G| G̲w.2vˏuߦ`q#}W,}ɂK;`6J/QqucF:q9T!(:^ e0kvT\A OULbWN \ ZU>|ώrIED
Ѽek>z,lkUIU`Cz]wt궮\UPǧGCӵx6
@~"綟hG*ʿyx^8~ySt=ҺT.:Yʌj
D h]a"hsX', S_5AK't <FJv7٪A>_Ƀm[!
cc0:-r+K
c1_ZBF܈1"hԷdNsvfa *V{ 10T|\#a{QDR
lCtBBbx\!vǒ~
UDW1$leUe95z:kf&L(7E@d"X -Ig^G,ۈs0eTd`XU /.(O] ˂l[-Pit-DBc MTi-fuD<K=m@iag8E;ġ!;HcP4<7!Xm1ˁQft_w7.a #*% y ,
S\#+S`'9Xm<pWMѩS+/{iT^y8{l^pFƷ9Ŀ^`~Xܗ.SZ#f)Y-"U<ukh
_C)<