1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
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;
}
|