1 --- boost/math/special_functions/detail/lanczos_sse2.hpp
2 +++ boost/math/special_functions/detail/lanczos_sse2.hpp
4 static_cast<double>(23531376880.41075968857200767445163675473L),
5 static_cast<double>(0u)
7 - register __m128d vx = _mm_load1_pd(&x);
8 - register __m128d sum_even = _mm_load_pd(coeff);
9 - register __m128d sum_odd = _mm_load_pd(coeff+2);
10 - register __m128d nc_odd, nc_even;
11 - register __m128d vx2 = _mm_mul_pd(vx, vx);
12 + __m128d vx = _mm_load1_pd(&x);
13 + __m128d sum_even = _mm_load_pd(coeff);
14 + __m128d sum_odd = _mm_load_pd(coeff+2);
15 + __m128d nc_odd, nc_even;
16 + __m128d vx2 = _mm_mul_pd(vx, vx);
18 sum_even = _mm_mul_pd(sum_even, vx2);
19 nc_even = _mm_load_pd(coeff + 4);
21 static_cast<double>(56906521.91347156388090791033559122686859L),
22 static_cast<double>(0u)
24 - register __m128d vx = _mm_load1_pd(&x);
25 - register __m128d sum_even = _mm_load_pd(coeff);
26 - register __m128d sum_odd = _mm_load_pd(coeff+2);
27 - register __m128d nc_odd, nc_even;
28 - register __m128d vx2 = _mm_mul_pd(vx, vx);
29 + __m128d vx = _mm_load1_pd(&x);
30 + __m128d sum_even = _mm_load_pd(coeff);
31 + __m128d sum_odd = _mm_load_pd(coeff+2);
32 + __m128d nc_odd, nc_even;
33 + __m128d vx2 = _mm_mul_pd(vx, vx);
35 sum_even = _mm_mul_pd(sum_even, vx2);
36 nc_even = _mm_load_pd(coeff + 4);