add a test for fast_acos()

1 files changed, 27 insertions, 0 deletions

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)