octave-nkf: src/DLD-FUNCTIONS/convhulln.cc annotate

annotate src/DLD-FUNCTIONS/convhulln.cc @ 8920:eb63fbe60fab

update copyright notices

author	John W. Eaton <jwe@octave.org>
date	Sat, 07 Mar 2009 10:41:27 -0500
parents	d0bc587fce55
children	2c279308f6ab

rev	line source
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	1 /*
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	2
8920 eb63fbe60fab update copyright notices John W. Eaton <jwe@octave.org> parents: 8788 diff changeset	3 Copyright (C) 2000, 2007, 2008, 2009 Kai Habel
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	4
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	5 This file is part of Octave.
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	6
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	7 Octave is free software; you can redistribute it and/or modify it
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	8 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6959 diff changeset	9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6959 diff changeset	10 option) any later version.
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	11
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	12 Octave is distributed in the hope that it will be useful, but WITHOUT
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	15 for more details.
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	16
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	17 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6959 diff changeset	18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6959 diff changeset	19 <http://www.gnu.org/licenses/>.
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	20
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	21 */
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	22
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	23 /*
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	24 29. July 2000 - Kai Habel: first release
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	25 2002-04-22 Paul Kienzle
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	26 * Use warning(...) function rather than writing to cerr
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	27 2006-05-01 Tom Holroyd
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	28 * add support for consistent winding in all dimensions; output is
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	29 * guaranteed to be simplicial.
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	30 */
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	31
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	32 #include <sstream>
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	33
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	34 #ifdef HAVE_CONFIG_H
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	35 #include <config.h>
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	36 #endif
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	37 #include "oct.h"
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	38 #include "Cell.h"
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	39
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	40 #ifdef HAVE_QHULL
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	41 #if defined(HAVE__SNPRINTF) && !defined(HAVE_SNPRINTF)
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	42 #define snprintf _snprintf
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	43 #endif
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	44
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	45 extern "C" {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	46 #include <qhull/qhull_a.h>
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	47 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	48
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	49 #ifdef NEED_QHULL_VERSION
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	50 char qh_version[] = "convhulln.oct 2007-07-24";
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	51 #endif
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	52 #endif
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	53
8786 dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	54 DEFUN_DLD (convhulln, args, nargout,
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	55 "-- texinfo --\n\
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	56 @deftypefn {Loadable Function} {@var{h} =} convhulln (@var{p})\n\
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	57 @deftypefnx {Loadable Function} {@var{h} =} convhulln (@var{p}, @var{opt})\n\
8788 d0bc587fce55 doc fix John W. Eaton <jwe@octave.org> parents: 8786 diff changeset	58 @deftypefnx {Loadable Function} {[@var{h}, @var{v}] =} convhulln (@dots{})\n\
8786 dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	59 Return an index vector to the points of the enclosing convex hull.\n\
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	60 The input matrix of size [n, dim] contains n points of dimension dim.\n\n\
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	61 If a second optional argument is given, it must be a string or cell array\n\
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	62 of strings containing options for the underlying qhull command. (See\n\
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	63 the Qhull documentation for the available options.) The default options\n\
8786 dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	64 are \"s Qci Tcv\".\n\
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	65 If the second output @var{V} is requested the volume of the convex hull is\n\
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	66 calculated.\n\n\
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	67 @seealso{convhull, delaunayn}\n\
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	68 @end deftypefn")
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	69 {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	70 octave_value_list retval;
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	71
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	72 #ifdef HAVE_QHULL
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	73 std::string options;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	74
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	75 int nargin = args.length ();
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	76 if (nargin < 1 \|\| nargin > 2)
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	77 {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	78 print_usage ();
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	79 return retval;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	80 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	81
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	82 if (nargin == 2)
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	83 {
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	84 if (args (1).is_string ())
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	85 options = args(1).string_value ();
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	86 else if (args(1).is_cell ())
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	87 {
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	88 Cell c = args(1).cell_value ();
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	89 options = "";
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	90 for (octave_idx_type i = 0; i < c.numel (); i++)
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	91 {
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	92 if (! c.elem(i).is_string ())
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	93 {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	94 error ("convhulln: second argument must be a string or cell array of strings");
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	95 return retval;
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	96 }
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	97
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	98 options = options + c.elem(i).string_value() + " ";
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	99 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	100 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	101 else
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	102 {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	103 error ("convhulln: second argument must be a string or cell array of strings");
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	104 return retval;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	105 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	106 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	107 else
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	108 // turn on some consistency checks
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	109 options = "s Qci Tcv";
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	110
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	111 Matrix p (args(0).matrix_value ());
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	112
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	113 const octave_idx_type dim = p.columns ();
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	114 const octave_idx_type n = p.rows ();
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	115 p = p.transpose ();
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	116
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	117 double *pt_array = p.fortran_vec ();
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	118
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	119 boolT ismalloc = False;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	120
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	121 OCTAVE_LOCAL_BUFFER (char, flags, 250);
6880 e00a8f661f06 [project @ 2007-09-07 21:39:55 by dbateman] dbateman parents: 6823 diff changeset	122
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	123 // hmm, lots of options for qhull here
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	124 // QJ guarantees that the output will be triangles
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	125 snprintf (flags, 250, "qhull QJ %s", options.c_str ());
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	126
7520 b166043585a8 use 0 instead of NULL John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	127 if (! qh_new_qhull (dim, n, pt_array, ismalloc, flags, 0, stderr))
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	128 {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	129 // If you want some debugging information replace the NULL
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	130 // pointer with stdout
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	131
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	132 vertexT vertex, *vertexp;
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	133 facetT *facet;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	134 setT *vertices;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	135 unsigned int nf = qh num_facets;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	136
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	137 Matrix idx (nf, dim);
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	138
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	139 octave_idx_type j, i = 0;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	140 FORALLfacets
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	141 {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	142 j = 0;
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	143 if (! facet->simplicial)
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	144 // should never happen with QJ
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	145 error ("convhulln: non-simplicial facet");
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	146
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	147 if (dim == 3)
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	148 {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	149 vertices = qh_facet3vertex (facet);
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	150 FOREACHvertex_ (vertices)
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	151 idx(i, j++) = 1 + qh_pointid(vertex->point);
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	152 qh_settempfree (&vertices);
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	153 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	154 else
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	155 {
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	156 if (facet->toporient ^ qh_ORIENTclock)
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	157 {
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	158 FOREACHvertex_ (facet->vertices)
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	159 idx(i, j++) = 1 + qh_pointid(vertex->point);
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	160 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	161 else
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	162 {
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	163 FOREACHvertexreverse12_ (facet->vertices)
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	164 idx(i, j++) = 1 + qh_pointid(vertex->point);
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	165 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	166 }
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	167 if (j < dim)
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	168 // likewise but less fatal
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	169 warning ("facet %d only has %d vertices", i, j);
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	170 i++;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	171 }
8786 dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	172
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	173 if (nargout == 2)
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	174 // calculate volume of convex hull
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	175 // taken from qhull src/geom2.c
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	176 {
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	177 realT area;
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	178 realT dist;
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	179
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	180 FORALLfacets
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	181 {
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	182 if (! facet->normal)
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	183 continue;
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	184
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	185 if (facet->upperdelaunay && qh ATinfinity)
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	186 continue;
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	187
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	188 facet->f.area = area = qh_facetarea (facet);
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	189 facet->isarea = True;
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	190
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	191 if (qh DELAUNAY)
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	192 {
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	193 if (facet->upperdelaunay == qh UPPERdelaunay)
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	194 qh totarea += area;
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	195 }
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	196 else
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	197 {
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	198 qh totarea += area;
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	199 qh_distplane (qh interior_point, facet, &dist);
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	200 qh totvol += -dist * area/ qh hull_dim;
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	201 }
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	202 }
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	203
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	204 retval(1) = octave_value (qh totvol);
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	205 }
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	206
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	207 retval(0) = octave_value (idx);
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	208 }
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	209
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	210 // free long memory
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	211 qh_freeqhull (! qh_ALL);
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	212
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	213 // free short memory and memory allocator
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	214 int curlong, totlong;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	215 qh_memfreeshort (&curlong, &totlong);
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	216
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	217 if (curlong \|\| totlong)
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	218 warning ("convhulln: did not free %d bytes of long memory (%d pieces)",
47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	219 totlong, curlong);
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	220 #else
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	221 error ("convhulln: not available in this version of Octave");
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	222 #endif
6959 47f4f4e88166 [project @ 2007-10-04 20:43:32 by jwe] jwe parents: 6882 diff changeset	223
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	224 return retval;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	225 }
8786 dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	226
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	227 /*
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	228 %!test
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	229 %! cube = [0 0 0;1 0 0;1 1 0;0 1 0;0 0 1;1 0 1;1 1 1;0 1 1];
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	230 %! [h, v] = convhulln(cube,'Pp');
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	231 %! assert (v, 1.0, 1e6*eps);
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	232 %!test
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	233 %! tetrahedron = [1 1 1;-1 -1 1;-1 1 -1;1 -1 -1];
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	234 %! [h, v] = convhulln(tetrahedron,'Pp');
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	235 %! assert (v, 8/3, 1e6*eps);
dbd428efbf56 Add calculation convex hull volume Kai Habel <kai.habel@gmx.de> parents: 7520 diff changeset	236 */

Mercurial > octave-nkf

annotate src/DLD-FUNCTIONS/convhulln.cc @ 8920:eb63fbe60fab