incorrect test on gamma function ?

[Marco Atzeri]

Hi,
could someone explain me why
we expect:

gamma(0) = Inf 
gamma(-1) = NaN 

>>>>> processing /pub/hg/octave/src/mappers.cc
  ***** test
 x = [-1, 0, 1, Inf];
 v = [NaN, Inf, 1, Inf];
 assert (gamma(x), v);
 assert (gamma(single (x)), single (v));
!!!!! test failed
assert (gamma (x),v) expected
   NaN   Inf     1   Inf
but got
   Inf   Inf     1   Inf
NaNs don't match

considering the graphics
http://en.wikipedia.org/wiki/Gamma_function

I presume the two values should be both Inf (or both NaN)

Marco