Logarytm z SSE lub przejście na FPU?

Question

Robię kilka obliczeń statystycznych. Potrzebuję ich, aby były szybkie, więc przepisałem większość z nich na SSE. Jestem całkiem nowy, więc zastanawiałem się, jakie jest właściwe podejście:

Według mojej wiedzy, nie ma funkcji log2 ani ln w SSE, przynajmniej nie do wersji 4.1, która jest najnowszą wersją obsługiwane przez sprzęt, którego używam.

Czy lepiej:

1. ekstraktu 4 pływaków i zrobić na nich obliczeń FPU do określenia entropii - I nie trzeba będzie ładować żadnej z tych wartości z powrotem do rejestrów SSE, tylko zsumować je do inny pływak
2. znaleźć funkcję dla wiatru, który nie log2

Answer 1

Wydaje się, że istnieje kilka implementacji SSE log2 wokół, np. this one.

Istnieje również Intel Approximate Maths Library, który ma między innymi funkcję log2 - jest stara (2000), ale jest to SSE2 i powinna nadal działać dość dobrze.

---

Zobacz także:

• sse_mathfun - SSE vector math library
• avx_mathfun - AVX vector math library
• libmvec - vector math library added in glibc 2.22

Answer 2

brak instrukcji SSE, który implementuje funkcję logarytmu. Jednak nie ma też jednej instrukcji x86, która wykonuje również ogólny logarytm. Jeśli myślisz o użyciu funkcji logarytmu, takiej jak log lub log10 ze standardowej biblioteki C, warto rzucić okiem na implementację, która jest używana w bibliotece o otwartym kodzie źródłowym, takiej jak libc. Możesz łatwo przetworzyć własne przybliżenie logarytmu, które działa we wszystkich elementach rejestru SSE.

Taka funkcja jest często realizowana za pomocą wielomianowej aproksymacji, która jest ważna w pewnych specyfikacjach dokładności w pewnym regionie argumentów wejściowych, takich jak seria Taylora. Następnie można skorzystać z właściwości logarytmu, aby zawrzeć ogólny argument wejściowy w dopuszczalnym zakresie wejściowym dla procedury logarytmicznej. Ponadto, można sparametryzować podstawę logarytmu wykorzystując właściwości:

log_y(x) = log_a(x)/log_a(y) 

Gdzie a jest podstawa logarytmu rutyny, który został utworzony.