SSE2: Функция журнала с двойной точностью

Мне нужна реализация функции журнала с открытым исходным кодом (без ограничений по лицензии), что-то с подписью

__m128d _mm_log_pd(__m128d);

Это доступно в Intel Short Vector Math Библиотека (часть ICC), но ICC не является ни бесплатным, ни открытым исходным кодом. Я ищу реализацию, использующую только встроенные функции.

Она должна использовать специальные приближения рациональных функций. Мне нужно что-то почти такое же точное, как журнал cmath, скажем, 9-10 десятичных цифр, но быстрее.