Lôgarit với SSE, hoặc chuyển sang FPU?

Question

Tôi đang thực hiện một số tính toán thống kê. Tôi cần chúng nhanh, vì vậy tôi viết lại hầu hết nó để sử dụng SSE. Tôi khá ngạc nhiên với nó, vì vậy tôi đã tự hỏi cách tiếp cận đúng ở đây là:

Theo hiểu biết của tôi, không có hàm log2 hoặc ln trong SSE, ít nhất là không lên đến 4.1, là phiên bản mới nhất được hỗ trợ bởi phần cứng tôi sử dụng.

Là nó tốt hơn để:

1. chiết xuất 4 phao, và làm các phép tính FPU vào chúng để xác định enthropy - Tôi sẽ không cần phải tải bất kỳ của những giá trị trở lại vào thanh ghi SSE, chỉ tổng hợp họ lên khác nổi
2. tìm thấy một chức năng cho SSE mà không log2

Answer 1

Dường như có một vài triển khai SSE log2 xung quanh, ví dụ: this one.

Ngoài ra còn có Intel Approximate Maths Library có chức năng log2 trong số những người khác - cũ (2000) nhưng là SSE2 và vẫn hoạt động tốt.

---

Xem thêm:

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

Answer 2

không có hướng dẫn SSE mà thực hiện một hàm logarit. Tuy nhiên, cũng không có lệnh x86 đơn nào thực hiện lôgarit chung. Nếu bạn đang nghĩ đến việc sử dụng hàm logarit như log hoặc log10 từ thư viện chuẩn C, bạn nên xem xét triển khai được sử dụng trong thư viện nguồn mở như libc. Bạn có thể dễ dàng cuộn xấp xỉ logarit của riêng bạn hoạt động trên tất cả các phần tử trong thanh ghi SSE.

Hàm như vậy thường được thực hiện bằng cách sử dụng phép tính xấp xỉ đa thức hợp lệ trong một số đặc điểm chính xác về một vùng nhất định của các đối số đầu vào, chẳng hạn như chuỗi Taylor. Sau đó, bạn có thể tận dụng các thuộc tính lôgarit để bọc một đối số đầu vào chung vào phạm vi đầu vào có thể chấp nhận được cho quy trình logarit của bạn. Bên cạnh đó, bạn có thể parameterize cơ sở của logarit bằng cách tận dụng tài sản:

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

đâu a là cơ sở của các thói quen logarit mà bạn đã tạo.