#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;
}