I need open source (no restriction on license) implementation of log function, something with signature

```
__m128d _mm_log_pd(__m128d);
```

It is available in Intel Short Vector Math Library (part of ICC), but ICC is neither free nor open source. I am looking for implementation using intrinsics only.

It should use special rational function approximations. I need something almost as accurate as cmath log, say 9-10 decimal digits, but faster.