kai/pounder_test
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

atan2: test min and max angle inputs

Browse Source
This commit is contained in:
Matt Huszagh 2020-12-17 11:34:39 -08:00
parent 6ffc42021e
commit 9c5e68ceea
1 changed files with 25 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
dsp/src/trig.rs
View File

@ -179,10 +179,10 @@ mod tests {
test_vals[i] = val_bounds.0 + i as f64 * val_delta; test_vals[i] = val_bounds.0 + i as f64 * val_delta;
} }
let atol: f64 = 4e-5;
let rtol: f64 = 0.127;
for &x in test_vals.iter() { for &x in test_vals.iter() {
for &y in test_vals.iter() { for &y in test_vals.iter() {
let atol: f64 = 4e-5;
let rtol: f64 = 0.127;
let actual = (y.atan2(x) as f64 * i16::MAX as f64).round() let actual = (y.atan2(x) as f64 * i16::MAX as f64).round()
/ i16::MAX as f64; / i16::MAX as f64;
let tol = atol + rtol * angle_to_axis(actual).abs(); let tol = atol + rtol * angle_to_axis(actual).abs();
@ -202,6 +202,29 @@ mod tests {
} }
} }
} }
// test min and max explicitly
for (x, y) in [
((i16::MIN as i32 + 1) << 16, -(1 << 16) as i32),
((i16::MIN as i32 + 1) << 16, (1 << 16) as i32),
]
.iter()
{
let yf = *y as f64 / ((i16::MAX as i32) << 16) as f64;
let xf = *x as f64 / ((i16::MAX as i32) << 16) as f64;
let actual =
(yf.atan2(xf) * i16::MAX as f64).round() / i16::MAX as f64;
let computed = (atan2(*y, *x) >> 16) as f64 / i16::MAX as f64 * PI;
let tol = atol + rtol * angle_to_axis(actual).abs();
if !isclose(computed, actual, 0., tol) {
println!("(x, y) : {}, {}", *x, *y);
println!("actual : {}", actual);
println!("computed : {}", computed);
println!("tolerance: {}\n", tol);
assert!(false);
}
}
} }
#[test] #[test]