From 02174afc5f27a6b255762f52b7e30af41941be0d Mon Sep 17 00:00:00 2001 From: tlatorre Date: Sun, 2 Jun 2019 13:00:29 -0400 Subject: add a test for fast_acos() --- src/test.c | 27 +++++++++++++++++++++++++++ 1 file changed, 27 insertions(+) diff --git a/src/test.c b/src/test.c index 5ab4fba..b1e0426 100644 --- a/src/test.c +++ b/src/test.c @@ -2074,6 +2074,32 @@ err: return 1; } +int test_fast_acos(char *err) +{ + /* Tests that the fast_acos() function returns values within 0.1% of acos(). */ + size_t i; + double x, result, expected; + + init_genrand(0); + + for (i = 0; i < 100; i++) { + x = genrand_real2()*2 - 1; + + result = fast_acos(x); + expected = acos(x); + + if (!isclose(result, expected, 0, 1e-3)) { + sprintf(err, "fast_acos() returned %.5g, but expected %.5g", result, expected); + goto err; + } + } + + return 0; + +err: + return 1; +} + struct tests { testFunction *test; char *name; @@ -2125,6 +2151,7 @@ struct tests { {test_argmax, "test_argmax"}, {test_argmin, "test_argmin"}, {test_electron_get_angular_pdf_norm, "test_electron_get_angular_pdf_norm"}, + {test_fast_acos, "test_fast_acos"}, }; int main(int argc, char **argv) -- cgit