| ... | @@ -95,7 +95,7 @@ fn exp2_32(x: f32) f32 { | ... | @@ -95,7 +95,7 @@ fn exp2_32(x: f32) f32 { |
| 95 | uf -= redux; | 95 | uf -= redux; |
| 96 | | 96 | |
| 97 | const z: f64 = x - uf; | 97 | const z: f64 = x - uf; |
| 98 | var r: f64 = exp2ft[i_0]; | 98 | var r: f64 = exp2ft[@intCast(usize, i_0)]; |
| 99 | const t: f64 = r * z; | 99 | const t: f64 = r * z; |
| 100 | r = r + t * (P1 + z * P2) + t * (z * z) * (P3 + z * P4); | 100 | r = r + t * (P1 + z * P2) + t * (z * z) * (P3 + z * P4); |
| 101 | return @floatCast(f32, r * uk); | 101 | return @floatCast(f32, r * uk); |
| ... | @@ -418,8 +418,8 @@ fn exp2_64(x: f64) f64 { | ... | @@ -418,8 +418,8 @@ fn exp2_64(x: f64) f64 { |
| 418 | | 418 | |
| 419 | // r = exp2(y) = exp2t[i_0] * p(z - eps[i]) | 419 | // r = exp2(y) = exp2t[i_0] * p(z - eps[i]) |
| 420 | var z = x - uf; | 420 | var z = x - uf; |
| 421 | const t = exp2dt[2 * i_0]; | 421 | const t = exp2dt[@intCast(usize, 2 * i_0)]; |
| 422 | z -= exp2dt[2 * i_0 + 1]; | 422 | z -= exp2dt[@intCast(usize, 2 * i_0 + 1)]; |
| 423 | const r = t + t * z * (P1 + z * (P2 + z * (P3 + z * (P4 + z * P5)))); | 423 | const r = t + t * z * (P1 + z * (P2 + z * (P3 + z * (P4 + z * P5)))); |
| 424 | | 424 | |
| 425 | return math.scalbn(r, ik); | 425 | return math.scalbn(r, ik); |