Re: Błędny epsilon - this is not a bug, this is ? - Grupy dyskusyjne w eGospodarka.pl

eGospodarka.pl › Grupy › pl.comp.programming › Błędny epsilon - this is not a bug, this is ? › Re: Błędny epsilon - this is not a bug, this is ?

Path: news-archive.icm.edu.pl!agh.edu.pl!news.agh.edu.pl!news.cyf-kr.edu.pl!news.nask
.pl!news.nask.org.pl!news.unit0.net!feeder.erje.net!eu.feeder.erje.net!proxad.n
et!feeder1-2.proxad.net!209.85.216.88.MISMATCH!c7no10651041qap.0!news-out.googl
e.com!gf5ni18188417qab.0!nntp.google.com!c7no10651038qap.0!postnews.google.com!
glegroupsg2000goo.googlegroups.com!not-for-mail
Newsgroups: pl.comp.programming
Date: Wed, 7 Nov 2012 13:03:03 -0800 (PST)
In-Reply-To: <5...@g...com>
Complaints-To: g...@g...com
Injection-Info: glegroupsg2000goo.googlegroups.com; posting-host=31.60.190.247;
posting-account=Sb6m8goAAABbWsBL7gouk3bfLsuxwMgN
NNTP-Posting-Host: 31.60.190.247
References: <50924bb1$0$1308$65785112@news.neostrada.pl>
<k72sqt$gq$1@news.dialog.net.pl>
<509505f2$0$1317$65785112@news.neostrada.pl>
<k734vm$43g$1@news.dialog.net.pl>
<509533d1$0$26687$65785112@news.neostrada.pl>
<k797v2$t98$1@news.dialog.net.pl>
<509848c6$0$1208$65785112@news.neostrada.pl>
<k79hhh$26k$1@news.dialog.net.pl> <k7am9f$n0p$1@zeus.man.szczecin.pl>
<k7apte$eiu$1@news.task.gda.pl>
<0...@g...com>
<k7b4tp$9ke$1@news.task.gda.pl> <k7b8b5$nqq$1@node1.news.atman.pl>
<k7bbgo$pth$1@news.task.gda.pl>
<50996f4d$0$1297$65785112@news.neostrada.pl>
<k7c0hi$jp4$1@news.task.gda.pl>
<50999d69$0$1215$65785112@news.neostrada.pl>
<k7c774$3tq$1@news.task.gda.pl> <k7dic5$ssp$1@node2.news.atman.pl>
<5...@g...com>
User-Agent: G2/1.0
MIME-Version: 1.0
Message-ID: <e...@g...com>
Subject: Re: Błędny epsilon - this is not a bug, this is ?
From: kenobi <p...@g...com>
Injection-Date: Wed, 07 Nov 2012 21:03:04 +0000
Content-Type: text/plain; charset=ISO-8859-2
Content-Transfer-Encoding: quoted-printable
Lines: 52
Xref: news-archive.icm.edu.pl pl.comp.programming:200626
[ ukryj nagłówki ]
W dniu środa, 7 listopada 2012 19:50:41 UTC+1 użytkownik e...@g...com
napisał:
> W dniu środa, 7 listopada 2012 06:57:57 UTC-5 użytkownik bartekltg napisał:
>
> > Rozmawiasz z facetem, który przez tydzień, mimo wysiłków
>
> >
>
> > grupy, nie dał sobie wytłumaczyć i uparcie twierdził, że
>
> >
>
> > "jakość" całki numerycznej nie zależy od rzędu kwadratury
>
> >
>
> > (Newtona-Cotesa, czyli tych zwykłych, kawałkami wielomianów).
>
> >
>
> > Wszytko rzędu wyższego niż "trapezy" miało być przereklamowanymi
>
> >
>
> > akademickimi bredniami nie dającymi poprawy wyniku względem
>
> >
>
> > metody trapezów:)
>
>
>
> A moze tak przy okazji ktos (moze OP?) wytlumaczy, jak sie
>
> liczy blad maszynowy wyrazenia sin(b/a) dla a=100 i powiedzmy
>
> b=549.755813887? Najlepiej na postawie tego DBL_EPS.
>
>
ciekawe pytanie, jak sie sumuje bledy w takich
wyrazeniach ja mialem na miernictwie w technikum i na pracowni fizycznej na studiach
- pewnie stosuje sie te wzorki

niby mozna zaryzykowac twierdzenie ze to moze byc wlasnie DBL_EPS, ale trzeba
sprawdzic czy blad z przeciwdziedziny sinusa nie jest wiekszy, choc chyba
5.49755813887 miesci sie w dublu z niedokladnoscia DBL_EPSILON , blad
przeciwdziedziny sinusa zdaje sie nie jest wiekszy niz to wiec chyba ten blad wynosi
wlasnie DBL_EPSILON