--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/toolbox/chop.m	Thu May 07 18:36:24 2015 +0200
@@ -0,0 +1,18 @@
+function c = chop(x, t)
+%CHOP    Round matrix elements.
+%        CHOP(X, t) is the matrix obtained by rounding the elements of X
+%        to t significant binary places.
+%        Default is t = 24, corresponding to IEEE single precision.
+
+if nargin < 2, t = 24; end
+[m, n] = size(x);
+
+%  Use the representation:
+%  x(i,j) = 2^e(i,j) * .d(1)d(2)...d(s) * sign(x(i,j))
+
+%  On the next line `+(x==0)' avoids passing a zero argument to LOG, which
+%  would cause a warning message to be generated.
+
+y = abs(x) + (x==0);
+e = floor(log2(y) + 1);
+c = pow2(round( pow2(x, t-e) ), e-t);