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