# HG changeset patch # User Rik # Date 1282414544 25200 # Node ID 9f0f1a89c7048b6731f5cd258c78f82fedabb552 # Parent 9a64e02e2aad2a08d39794b8539ed7b987b11367 Add new %!tests for negative integer arguments to gamma, lgamma. diff --git a/src/ChangeLog b/src/ChangeLog --- a/src/ChangeLog +++ b/src/ChangeLog @@ -1,3 +1,8 @@ +2010-08-17 Rik + + * mappers.cc (gamma,lgamma): Add %!tests for negative integer + input arguments. + 2010-08-17 Rik * bitfcns.cc (bitmax): Add additional calling argument, class, and diff --git a/src/mappers.cc b/src/mappers.cc --- a/src/mappers.cc +++ b/src/mappers.cc @@ -916,6 +916,12 @@ %! v = single([sqrt(pi), 1, .5*sqrt(pi), 1, 2, 6, 24]); %! assert(gamma(x), v, sqrt(eps('single'))) +%!test +%! x = [-1, 0, 1, Inf]; +%! v = [NaN, Inf, 1, Inf]; +%! assert (gamma(x), v); +%! assert (gamma(single (x)), single (v)); + %!error gamma(); %!error gamma(1,2); @@ -1294,6 +1300,12 @@ %! v = single([sqrt(pi), 1, .5*sqrt(pi), 1, 2, 6, 24]); %! assert(lgamma(x), log(v), sqrt(eps ('single'))) +%!test +%! x = [-1, 0, 1, Inf]; +%! v = [Inf, Inf, 0, Inf]; +%! assert (lgamma(x), v); +%! assert (lgamma(single (x)), single(v)); + %!error lgamma(); %!error lgamma(1,2);