lgamma(-n) for n = 0, 1, 2, .. does not raise FE_DIVBYZERO.

[kargl]

Tests with olibm: lgammaf( 1) = 0. Got no exception. Expected 0 without exception. lgammaf( 2) = 0. Got no exception. Expected 0 without exception. lgammaf( 0) = inf. Got no exception. Expected inf with FE_DIVBYZERO. lgammaf(-0) = inf. Got no exception. Expected inf with FE_DIVBYZERO. lgammaf(-1) = inf. Got no exception. Expected inf with FE_DIVBYZERO.

Tests with FreeBSD libm. lgammaf( 1) = 0. Got no exception. Expected 0 without exception. lgammaf( 2) = 0. Got no exception. Expected 0 without exception. lgammaf( 0) = inf. Got FE_DIVBYZERO. Expected inf with FE_DIVBYZERO. lgammaf(-0) = inf. Got FE_DIVBYZERO. Expected inf with FE_DIVBYZERO. lgammaf(-1) = inf. Got FE_DIVBYZERO. Expected inf with FE_DIVBYZERO.

[kargl]

Note, lgammal(-n) correctly raises the FE_DIVBYZERO exception. lgammaf and lgamma have issues.

[oscardssmith]

This hasn't yet been fixed