#include <math.h>
#include <stddef.h>
#include "ts.h"

bool ts_test_acosh( void )
{
	fp_t const x0[] = { 1.000000, 1.500000, 2.000000, 3.000000, 5.000000 };
	fp_t const e0[] = { 0.000000, 0.962424, 1.316958, 1.762747, 2.292432 };

	fp_t const * const x[1] = { x0 };
	fp_t const * const e[1] = { e0 };

	bool result = ts_test(ts_acosh, ts_acosh_offset, x, e, sizeof( x0 ) / sizeof( *x0 ), (ts_constant_t *)NULL, &ts_acosh_p);

	fp_t const xb[] = { 0.000000, -1.000000, 0.500000 };

	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_acosh(yi, 1, xi, 1, 1, NULL);
		result &= isnan( tmp );
	}

	return result;
}