On 22.07.2014 16:16, firr wrote:
>>>
>>>
>>> Zerknąłem do tego, co robi oryginalny log. Też głownie
>>> mnoży i dodaje, ma jednak nieco więcej skoków.

> powody sa lub ich nie ma.. (co do tego to moze zalezy od platformy

Przestań "dresować".

> a jaki to jest 'oryginalny log'? jest to jakies zrodło w c w asmie
> czy cos takiego? ztcw to w x87 sa dwie funkcje

Napisałem o tym w ostatniej linijce posta. Kompilator nie użył
koprocesora, tylko liczy jakiś szereg używając sse,
po drodze używając paru porównań.

Wymuszenie użycia koprocesora przez -mfpmath=387 nic nie daje,
bo treść log2 się linkuje w wersji sse.

Jak bardzo chcesz źródło, to masz.
http://pastebin.com/BZpVhHGb
najpierw std, potem to co wypluł kompilator z f.wymiernej + frexp.
[bardzo ładnie sam przeplata liczenie licznika i mianownika]


> FYL2X - liczy y*lg_2(x) (jesli y=lg_b(2) => liczy lg_b(x) )
> FYL2XP1 - y*lg_2(x+1) "more precise than lg_2(x) if x is close to
> zero" (acz tego troche nie rozumiem - to jak sie tego uzywa?)

Jeśli masz liczbę postaci 1+dx to logartym (naturalny dla
uproszczenia) tegobędzie z grubsza wynosił dx. Ale precyzja
1+dx to 16 cyfr, jeśli dx jesst na poziomie 10^-10 to
dx ma tylko 6 cyfr znaczących. I tyle ma też wynik.

A logartym w tym punkcie jest przydatny. Zwłaszcza, ze
dx może być równe 10^-40 ;-)
log1p (x) = log(1+x) tyle, że gdy x jest małe, znacznie dokładniej.

pzdr
bartekltg