1
1
# From openlibm/test/libm-test-ulps.h, openlibm/test/libm-test.c
2
2
3
- using SpecialFunctions: _lgamma_r, _loggamma_r, _lgammaf_r, _loggammaf_r
3
+ using SpecialFunctions: _lgamma_r, _lgammaf_r
4
4
5
5
# lgamma_test block
6
- for (T, lgamma, labsgamma) in ((Float64, _loggamma_r, _lgamma_r),
7
- (Float32, _loggammaf_r, _lgammaf_r))
6
+ for (T, lgamma) in ((Float64, _lgamma_r), (Float32, _lgammaf_r))
8
7
@testset " lgamma_test, $T " begin
9
- @test lgamma (T (Inf )) === T (Inf )
10
- @test lgamma (T (0 )) === T (Inf )
11
- @test lgamma (T (NaN )) === T (NaN )
8
+ @test lgamma (T (Inf ))[ 1 ] === T (Inf )
9
+ @test lgamma (T (0 ))[ 1 ] === T (Inf )
10
+ @test lgamma (T (NaN ))[ 1 ] === T (NaN )
12
11
13
- @test lgamma (T (- 3 )) === T (Inf )
14
- @test lgamma (T (- Inf )) === T (Inf )
12
+ @test lgamma (T (- 3 ))[ 1 ] === T (Inf )
13
+ @test lgamma (T (- Inf ))[ 1 ] === T (Inf )
15
14
16
15
# lgamma(1) == 0, lgamma (1) sets signgam to 1
17
16
y, signgam = labsgamma (T (1 ))
@@ -32,7 +31,6 @@ for (T, lgamma, labsgamma) in ((Float64, _loggamma_r, _lgamma_r),
32
31
y, signgam = labsgamma (T (- 0.5 ))
33
32
@test y === T (0.5 log (4 π))
34
33
@test signgam == - 1
35
- @test_throws DomainError lgamma (T (- 0.5 ))
36
34
37
35
# In the two "broken" tests, an exact match not possible, even
38
36
# in Float64, thus, we check for as close a tolerance as
0 commit comments