Mercurial > octave
view scripts/general/polyarea.m @ 31198:863730dd0f83 stable
nextpow2: Fix for input between 0.5 and 1 (bug #62947).
* scripts/general/nextpow2.m: Switch to a naïve implementation using log2 with
a single output argument and ceil.
author | Markus Mützel <markus.muetzel@gmx.de> |
---|---|
date | Wed, 24 Aug 2022 17:15:34 +0200 |
parents | 796f54d4ddbf |
children | 397d29f7135c |
line wrap: on
line source
######################################################################## ## ## Copyright (C) 1999-2022 The Octave Project Developers ## ## See the file COPYRIGHT.md in the top-level directory of this ## distribution or <https://octave.org/copyright/>. ## ## This file is part of Octave. ## ## Octave 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 3 of the License, or ## (at your option) any later version. ## ## Octave 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. ## ## You should have received a copy of the GNU General Public License ## along with Octave; see the file COPYING. If not, see ## <https://www.gnu.org/licenses/>. ## ######################################################################## ## -*- texinfo -*- ## @deftypefn {} {} polyarea (@var{x}, @var{y}) ## @deftypefnx {} {} polyarea (@var{x}, @var{y}, @var{dim}) ## ## Determine area of a polygon by triangle method. ## ## The variables @var{x} and @var{y} define the vertex pairs, and must ## therefore have the same shape. They can be either vectors or arrays. If ## they are arrays then the columns of @var{x} and @var{y} are treated ## separately and an area returned for each. ## ## If the optional @var{dim} argument is given, then @code{polyarea} works ## along this dimension of the arrays @var{x} and @var{y}. ## ## @end deftypefn ## FIXME: Add moments for centroid, etc. ## ## Bugs and limitations: ## Probably ought to be an optional check to make sure that ## traversing the vertices doesn't make any sides cross ## (Is simple closed curve the technical definition of this?). function a = polyarea (x, y, dim) if (nargin < 2) print_usage (); elseif (! size_equal (x, y)) error ("polyarea: X and Y must have the same shape"); endif if (nargin == 2) a = abs (sum (x .* (shift (y, -1) - shift (y, 1)))) / 2; else a = abs (sum (x .* (shift (y, -1, dim) - shift (y, 1, dim)), dim)) / 2; endif endfunction %!shared x, y %! x = [1;1;3;3;1]; %! y = [1;3;3;1;1]; %!assert (polyarea (x,y), 4, eps) %!assert (polyarea ([x,x],[y,y]), [4,4], eps) %!assert (polyarea ([x,x],[y,y],1), [4,4], eps) %!assert (polyarea ([x,x]',[y,y]',2), [4;4], eps)