Mercurial > hg > octave-lyh

comparison scripts/polynomial/polygcd.m @ 7795:df9519e9990c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Handle single precision eps values
author David Bateman <dbateman@free.fr>
date Mon, 12 May 2008 22:57:11 +0200
parents 120f3135952f
children 6f2d95255911
comparison
equal deleted inserted replaced
7794:2b458dfe31ae 7795:df9519e9990c
43 43
44 function x = polygcd (b, a, tol) 44 function x = polygcd (b, a, tol)
45 45
46 if (nargin == 2 || nargin == 3) 46 if (nargin == 2 || nargin == 3)
47 if (nargin == 2) 47 if (nargin == 2)
48 tol = sqrt (eps); 48 if (isa (a, "single") || isa (b, "single"))
49 tol = sqrt (eps ("single"));
50 else
51 tol = sqrt (eps);
52 endif
49 endif 53 endif
50 if (length (a) == 1 || length (b) == 1) 54 if (length (a) == 1 || length (b) == 1)
51 if (a == 0) 55 if (a == 0)
52 x = b; 56 x = b;
53 elseif (b == 0) 57 elseif (b == 0)