Mercurial > octave
diff scripts/general/nextpow2.m @ 3426:f8dde1807dee
[project @ 2000-01-13 08:40:00 by jwe]
author | jwe |
---|---|
date | Thu, 13 Jan 2000 08:40:53 +0000 |
parents | 5e0a0b1cba43 |
children | 38c61cbf086c |
line wrap: on
line diff
--- a/scripts/general/nextpow2.m Thu Jan 13 08:32:16 2000 +0000 +++ b/scripts/general/nextpow2.m Thu Jan 13 08:40:53 2000 +0000 @@ -1,15 +1,15 @@ ## Copyright (C) 1995, 1996 Kurt Hornik -## +## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by ## the Free Software Foundation; either version 2, or (at your option) ## any later version. -## +## ## This program is distributed in the hope that it will be useful, but ## WITHOUT ANY WARRANTY; without even the implied warranty of ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU -## General Public License for more details. -## +## General Public License for more details. +## ## You should have received a copy of the GNU General Public License ## along with this file. If not, write to the Free Software Foundation, ## 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. @@ -25,7 +25,7 @@ ## @ifinfo ## 2^n >= abs (x). ## @end ifinfo -## +## ## If @var{x} is a vector, return @code{nextpow2 (length (@var{x}))}. ## @end deftypefn ## @seealso{pow2} @@ -35,7 +35,7 @@ ## Adapted-By: jwe function n = nextpow2 (x) - + if (nargin != 1) usage ("nextpow2 (x)"); endif @@ -48,10 +48,10 @@ if (t > 1) x = t; endif - + [f, n] = log2 (abs (x)); if (f == 0.5) n = n - 1; endif - + endfunction