1 files changed, 35 insertions, 0 deletions
diff --git a/test/atanh.c b/test/atanh.c
new file mode 100644
index 0000000..91a9169
--- /dev/null
+++ b/test/atanh.c
@@ -0,0 +1,35 @@
+#include <math.h>
+#include <stddef.h>
+#include "ts.h"
+
+bool ts_test_atanh( void )
+{
+	fp_t const x0[] = { -0.900000, -0.500000, 0.000000, 0.500000, 0.900000 };
+	fp_t const e0[] = { -1.472219, -0.549306, 0.000000, 0.549306, 1.472219 };
+
+	fp_t const * const x[1] = { x0 };
+	fp_t const * const e[1] = { e0 };
+
+	bool result = ts_test(ts_atanh, ts_atanh_offset, x, e, sizeof( x0 ) / sizeof( *x0 ), (ts_constant_t *)NULL, &ts_atanh_p);
+
+	fp_t const xb[] = { -1.000000, 1.000000, -1.500000, 1.500000 };
+	fp_t const eb[] = { -INFINITY, INFINITY, NAN, NAN };
+
+	for( uint_t i = 0; i < sizeof( xb ) / sizeof( *xb ); ++i )
+	{
+		fp_t tmp = 0;
+		fp_t const * const xi[1] = { &xb[i] };
+		fp_t * const yi[1] = { &tmp };
+		ts_atanh(yi, 1, xi, 1, 1, NULL);
+		if( isinf( eb[i] ) )
+		{
+			result &= isinf( tmp ) && signbit( tmp ) == signbit( eb[i] );
+		}
+		else if( isnan( eb[i] ) )
+		{
+			result &= isnan( tmp );
+		}
+	}
+
+	return result;
+}