Mercurial > hg > octave-max
diff src/mappers.cc @ 7667:fb3a6c53c2b2
Allow negative zero imaginary part to be treated as zero for erf, erfc, gamma and lgamma mapper function
| author | David Bateman <dbateman@free.fr> |
|---|---|
| date | Fri, 28 Mar 2008 18:26:29 +0100 |
| parents | 2df457529cfa |
| children | 87865ed7405f |
line wrap: on
line diff
--- a/src/mappers.cc +++ b/src/mappers.cc @@ -323,6 +323,14 @@ return retval; } +/* + +%!test +%! a = -1i*sqrt(-1/(6.4187*6.4187)); +%! assert (erf(a), erf(real(a))); + +*/ + DEFUN (erfc, args, , "-*- texinfo -*-\n\ @deftypefn {Mapping Function} {} erfc (@var{z})\n\ @@ -347,6 +355,14 @@ return retval; } +/* + +%!test +%! a = -1i*sqrt(-1/(6.4187*6.4187)); +%! assert (erfc(a), erfc(real(a))); + +*/ + DEFUN (exp, args, , "-*- texinfo -*-\n\ @deftypefn {Mapping Function} {} exp (@var{x})\n\ @@ -467,6 +483,14 @@ return retval; } +/* + +%!test +%! a = -1i*sqrt(-1/(6.4187*6.4187)); +%! assert (gamma(a), gamma(real(a))); + +*/ + DEFUN (imag, args, , "-*- texinfo -*-\n\ @deftypefn {Mapping Function} {} imag (@var{z})\n\ @@ -753,6 +777,14 @@ return retval; } +/* + +%!test +%! a = -1i*sqrt(-1/(6.4187*6.4187)); +%! assert (lgamma(a), lgamma(real(a))); + +*/ + DEFUN (log, args, , "-*- texinfo -*-\n\ @deftypefn {Mapping Function} {} log (@var{x})\n\