-
Data: 2012-11-03 12:54:25
Temat: Re: Błędny epsilon - this is not a bug, this is ?
Od: "slawek" <h...@s...pl> szukaj wiadomości tego autora
[ pokaż wszystkie nagłówki ]Użytkownik "Tomasz Sowa" napisał w wiadomości grup
dyskusyjnych:k72sqt$gq$...@n...dialog.net.pl...
>Na wikipedii jest 2.220446e-16 (zjedź na dół do przykładu)
A ty przeczytaj tabelę u góry (na tej samej stronie Wikipedii).
Jak na to nie patrzeć - trzeba poprawić - albo u góry, albo u dołu, albo w
obu miejscach.
>A po co taki naiwny algorytm? definicja maszynowego epsilon chyba jest
>jasna?
A po co lepszy? Skoro i tak liczy się w mniej niż sekundę? (Czyli o parę
rzędów wielkości krócej, niż trwało np. napisanie przez ciebie tekstu?)
>I jest to prawidłowa wartość.
Możliwe. Pod warunkiem, że przyjmiemy inną definicję niż podana w samym
pliku float.h.
Wniosek - albo definicja w pliku float.h (z pakietu MSVS) jest błędna, albo
wartość w pliku float.h (ibidem) jest błędna, albo obie są błędne.
I to właśnie (podobnie jak hasło w Wikipedii) wymagałoby wyprostowania.
>Pokaż ten program.
Nie będę ciebie obrażał podejrzeniem, że nie potrafisz. Do sprawdzenia, że
1.0+1.5E-16 > 1.0 to wystarczy ci jedna linijka zaczynająca się od cout.
>Nie wiem co tu jest do dzielenia, aby obliczyć maszynowe epsilon nic nie
>trzeba dzielić, przykład:
Owszem, ale ty nie obliczasz epsilona jako inf { x in R : op(1,x) > 1 },
czyli jako najmniejszej liczby, która w wyniku "operacji dodawania
maszynowego" do liczby 1 daje wynik większy niż 1. Zajrzyj sobie do
Teukolsky'ego - oczywiście masz prawo twierdzić, że są tam bzdury. Ale w
takim razie warto napisać do Teukolsky'ego - i pouczyć go jak ma zmienić
fragment rozdziału.
Prawdą jest (w opisanych warunkach), że (double)1.0 + (double)2.22E-16 >
(double)1.0 , ale prawdą jest też, że np. (double)1.0 + (double)1.5E-16 >
(double)1.0E-16.
Z mojego punktu widzenia, jeżeli mam macierz 1000x1000 a chcę przekształcić
ja na macierz rzadką przeglądając jej elementy i uznając za zerowe te, które
są mniejsze niż eps, to lepiej jeżeli wartości równe DLB_EPSILON zostaną jak
są, bo jeżeli odrzucę pół miliona takich - to może polecieć mi w porywach
pięć miejsc znaczących. Przy epsilon = 1.11E-16 czy uznam za zero, czy
zostawię jak są - wynik sumowania z 1.0 będzie jednakowy i nic nie stracę
jeżeli liczby były 64-bitowe z 53-bitową mantysą.
Następne wpisy z tego wątku
- 03.11.12 14:07 Tomasz Sowa
- 03.11.12 16:10 slawek
- 03.11.12 17:59 Michoo
- 03.11.12 22:22 slawek
- 04.11.12 15:42 kenobi
- 04.11.12 22:50 AK
- 04.11.12 23:00 AK
- 05.11.12 08:26 Adam Wysocki
- 05.11.12 11:07 Michoo
- 05.11.12 11:52 slawek
- 05.11.12 12:14 slawek
- 05.11.12 12:48 slawek
- 05.11.12 13:04 slawek
- 05.11.12 13:10 Michoo
- 05.11.12 16:17 AK
Najnowsze wątki z tej grupy
- Can you activate BMW 48V 10Ah Li-Ion battery, connecting to CAN-USB laptop interface ?
- We Wrocławiu ruszyła Odra 5, pierwszy w Polsce komputer kwantowy z nadprzewodzącymi kubitami
- Ada-Europe - AEiC 2025 early registration deadline imminent
- John Carmack twierdzi, że gdyby gry były optymalizowane, to wystarczyły by stare kompy
- Ada-Europe Int.Conf. Reliable Software Technologies, AEiC 2025
- Linuks od wer. 6.15 przestanie wspierać procesory 486 i będzie wymagać min. Pentium
- ,,Polski przemysł jest w stanie agonalnym" - podkreślił dobitnie, wskazując na brak zamówień.
- Rewolucja w debugowaniu!!! SI analizuje zrzuty pamięci systemu M$ Windows!!!
- Brednie w wiki - hasło Dehomag
- Perfidne ataki krakerów z KRLD na skrypciarzy JS i Pajton
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- U nas propagują modę na SI, a w Chinach naukowcy SI po kolei umierają w wieku 40-50lat
- C++. Podróż Po Języku - komentarz
Najnowsze wątki
- 2025-07-12 Warszawa => PC Hardware Expert / Specjalista PC <=
- 2025-07-12 Warszawa => Account Manager - Usługi rekrutacyjne <=
- 2025-07-12 Warszawa => Administrator IT <=
- 2025-07-12 Warszawa => IT Administrator <=
- 2025-07-12 Warszawa => Asystent/tka ds. Administracji <=
- 2025-07-12 Warszawa => Specjalista/stka ds. Organizacji <=
- 2025-07-12 Warszawa => MENA New Business Manager <=
- 2025-07-12 Gdynia => Controlling systems Consultant <=
- 2025-07-12 Warszawa => Developer Microsoft Dynamics 365 Finance & Operations (D36
- 2025-07-12 Warszawa => Programista Microsoft Dynamics 365 Finance & Operations (D
- 2025-07-12 Warszawa => Dyrektor IT <=
- 2025-07-12 Warszawa => IT Director <=
- 2025-07-12 Czy wypowiedź Kaczyńskiego o Braunie jest skarżalna? ["działa z OBCEJ inspiracji"]
- 2025-07-11 Rejestrator temperatur - termopara, siec
- 2025-07-11 DPD, przeniesienie numerów z a2mobile i z Orange