Article View: pl.comp.os.linux.programowanie
Article #72449Re: DBL_EPSILON i Calculated Machine epsilon
From: Maciej Witkowiak
Date: Sat, 28 Apr 2012 12:39
Date: Sat, 28 Apr 2012 12:39
41 lines
2478 bytes
2478 bytes
Adam Majewski napisał: > Witam. > > Napisałem program w c ( kod poniżej). Oblicza ile potrzeba kroków ( > iteracji tu LastIteration) aby punkt z1 uciekł z koła o promieniu ER, > czyli sprawdza abs(z1)<ER (...) > Jak to się ma do : > > DBL_EPSILON = 2.220446e-16 = > 0.00000000000000022204460492503130808473 ( z float.h ) > ? Przy wywołaniu: gcc -lm ... oraz gcc -lm -mfpmath=sse ... dostałem identyczny wynik, jak Ty. Ale po kompilacji przez: gcc -lm -fpmath=387 ... pojawiły się nowe wiersze wyniku: n= 27 distance(z1,zp) = 7.450581e-09 = 0.00000000745058059692382812500000000000 LastIteration = 134217761 n= 28 distance(z1,zp) = 3.725290e-09 = 0.00000000372529029846191406250000000000 LastIteration = 268435874 n= 29 distance(z1,zp) = 1.862645e-09 = 0.00000000186264514923095703125000000000 LastIteration = 536883630 n= 30 distance(z1,zp) = 9.313226e-10 = 0.00000000093132257461547851562500000000 LastIteration = 1074147225 n= 31 distance(z1,zp) = 4.656613e-10 = 0.00000000046566128730773925781250000000 LastIteration = 2160053374 n= 32 distance(z1,zp) = 2.328306e-10 = 0.00000000023283064365386962890625000000 LastIteration = 4245834220 Tak na pierwszy rzut oka - może lepiej przepisać ten program tak, żeby używał GMP? I. -- Najlepsza sygnatura to brak sygnatury. http://bossstation.dnsalias.org/
Technical Details
Message-ID:
<slrnjpnp87.ldn.ytm@salyut.bossstation.dnsalias.org>Path:
mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.icm.edu.pl!.POSTED!not-for-mailReferences:
<jmehj3$mkj$1@news.icm.edu.pl>Show Full Headers
Complete Article Headers
Path: mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.icm.edu.pl!.POSTED!not-for-mail From: Maciej Witkowiak <ytm@elysium.pl.andremowe.me> Newsgroups: pl.comp.os.linux.programowanie Subject: Re: DBL_EPSILON i Calculated Machine epsilon Date: Sat, 28 Apr 2012 12:39:52 +0000 (UTC) Organization: Zakon Braci Paradygmatow Lines: 41 Message-ID: <slrnjpnp87.ldn.ytm@salyut.bossstation.dnsalias.org> References: <jmehj3$mkj$1@news.icm.edu.pl> Reply-To: ytm@elysiumremowe.pl NNTP-Posting-Host: 77-45-90-10.sta.asta-net.com.pl Mime-Version: 1.0 X-Trace: news.icm.edu.pl 1335616792 26635 77.45.90.10 (28 Apr 2012 12:39:52 GMT) X-Complaints-To: usenet@news.icm.edu.pl NNTP-Posting-Date: Sat, 28 Apr 2012 12:39:52 +0000 (UTC) User-Agent: slrn/0.9.9p1 (Linux 2.6.37-5/x86_64) X-GeekCode-V3.1: GCS/M d+@ s-:- a- C++(+++) ULO++ P++ L++++ E--- W++(--) N++(+++) o? K w++ O? M? V? PS+ PE Y++> PGP++ t++ 5-- X- R+> tv- b++ DI D++ G++ e+++ h r++ y? X-University: UAM X-Registered-Linux-User: 150987 X-Paszport-Polsatu: 81062721 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit