octave: libinterp/dldfcn/__voronoi_

annotate libinterp/dldfcn/voronoi.cc @ 21200:fcac5dbbf9ed

maint: Indent #ifdef blocks in libinterp. * builtins.h, Cell.cc, __contourc__.cc, __dispatch__.cc, __dsearchn__.cc, __ichol__.cc, __ilu__.cc, __lin_interpn__.cc, __pchip_deriv__.cc, __qp__.cc, balance.cc, besselj.cc, betainc.cc, bitfcns.cc, bsxfun.cc, c-file-ptr-stream.cc, c-file-ptr-stream.h, cellfun.cc, colloc.cc, comment-list.cc, conv2.cc, daspk.cc, dasrt.cc, dassl.cc, data.cc, debug.cc, defaults.cc, defaults.in.h, defun-dld.h, defun.cc, defun.h, det.cc, dirfns.cc, display.cc, dlmread.cc, dot.cc, dynamic-ld.cc, eig.cc, ellipj.cc, error.cc, errwarn.cc, event-queue.cc, fft.cc, fft2.cc, fftn.cc, file-io.cc, filter.cc, find.cc, gammainc.cc, gcd.cc, getgrent.cc, getpwent.cc, getrusage.cc, givens.cc, gl-render.cc, gl2ps-print.cc, graphics.cc, graphics.in.h, gripes.cc, hash.cc, help.cc, hess.cc, hex2num.cc, input.cc, inv.cc, jit-ir.cc, jit-typeinfo.cc, jit-util.cc, jit-util.h, kron.cc, load-path.cc, load-save.cc, lookup.cc, ls-ascii-helper.cc, ls-hdf5.cc, ls-mat-ascii.cc, ls-mat4.cc, ls-mat5.cc, ls-oct-binary.cc, ls-oct-text.cc, ls-oct-text.h, ls-utils.cc, ls-utils.h, lsode.cc, lu.cc, luinc.cc, mappers.cc, matrix_type.cc, max.cc, mex.h, mexproto.h, mgorth.cc, nproc.cc, oct-errno.in.cc, oct-fstrm.cc, oct-hdf5-types.cc, oct-hdf5.h, oct-hist.cc, oct-iostrm.cc, oct-lvalue.cc, oct-map.cc, oct-prcstrm.cc, oct-procbuf.cc, oct-stream.cc, oct-strstrm.cc, octave-link.cc, ordschur.cc, pager.cc, pinv.cc, pr-output.cc, procstream.cc, profiler.cc, psi.cc, pt-jit.cc, quad.cc, quadcc.cc, qz.cc, rand.cc, rcond.cc, regexp.cc, schur.cc, sighandlers.cc, sparse-xdiv.cc, sparse-xpow.cc, sparse.cc, spparms.cc, sqrtm.cc, str2double.cc, strfind.cc, strfns.cc, sub2ind.cc, svd.cc, sylvester.cc, symtab.cc, syscalls.cc, sysdep.cc, sysdep.h, time.cc, toplev.cc, tril.cc, tsearch.cc, txt-eng-ft.cc, txt-eng.cc, typecast.cc, urlwrite.cc, utils.cc, variables.cc, xdiv.cc, xnorm.cc, xpow.cc, zfstream.cc, __delaunayn__.cc, __eigs__.cc, __fltk_uigetfile__.cc, __glpk__.cc, __init_fltk__.cc, __init_gnuplot__.cc, __magick_read__.cc, __osmesa_print__.cc, __voronoi__.cc, amd.cc, audiodevinfo.cc, audioread.cc, ccolamd.cc, chol.cc, colamd.cc, convhulln.cc, dmperm.cc, fftw.cc, oct-qhull.h, qr.cc, symbfact.cc, symrcm.cc, oct-conf.in.cc, ov-base-diag.cc, ov-base-int.cc, ov-base-mat.cc, ov-base-scalar.cc, ov-base-sparse.cc, ov-base.cc, ov-bool-mat.cc, ov-bool-sparse.cc, ov-bool.cc, ov-builtin.cc, ov-cell.cc, ov-ch-mat.cc, ov-class.cc, ov-classdef.cc, ov-colon.cc, ov-complex.cc, ov-cs-list.cc, ov-cx-diag.cc, ov-cx-mat.cc, ov-cx-sparse.cc, ov-dld-fcn.cc, ov-fcn-handle.cc, ov-fcn-inline.cc, ov-fcn.cc, ov-float.cc, ov-flt-complex.cc, ov-flt-cx-diag.cc, ov-flt-cx-mat.cc, ov-flt-re-diag.cc, ov-flt-re-mat.cc, ov-int16.cc, ov-int32.cc, ov-int64.cc, ov-int8.cc, ov-java.cc, ov-lazy-idx.cc, ov-mex-fcn.cc, ov-null-mat.cc, ov-oncleanup.cc, ov-perm.cc, ov-range.cc, ov-re-diag.cc, ov-re-mat.cc, ov-re-sparse.cc, ov-scalar.cc, ov-str-mat.cc, ov-struct.cc, ov-typeinfo.cc, ov-uint16.cc, ov-uint32.cc, ov-uint64.cc, ov-uint8.cc, ov-usr-fcn.cc, ov.cc, ovl.cc, octave.cc, op-b-b.cc, op-b-bm.cc, op-b-sbm.cc, op-bm-b.cc, op-bm-bm.cc, op-bm-sbm.cc, op-cdm-cdm.cc, op-cell.cc, op-chm.cc, op-class.cc, op-cm-cm.cc, op-cm-cs.cc, op-cm-m.cc, op-cm-s.cc, op-cm-scm.cc, op-cm-sm.cc, op-cs-cm.cc, op-cs-cs.cc, op-cs-m.cc, op-cs-s.cc, op-cs-scm.cc, op-cs-sm.cc, op-dm-dm.cc, op-dm-scm.cc, op-dm-sm.cc, op-dm-template.cc, op-dms-template.cc, op-double-conv.cc, op-fcdm-fcdm.cc, op-fcdm-fdm.cc, op-fcm-fcm.cc, op-fcm-fcs.cc, op-fcm-fm.cc, op-fcm-fs.cc, op-fcn.cc, op-fcs-fcm.cc, op-fcs-fcs.cc, op-fcs-fm.cc, op-fcs-fs.cc, op-fdm-fdm.cc, op-float-conv.cc, op-fm-fcm.cc, op-fm-fcs.cc, op-fm-fm.cc, op-fm-fs.cc, op-fs-fcm.cc, op-fs-fcs.cc, op-fs-fm.cc, op-fs-fs.cc, op-i16-i16.cc, op-i32-i32.cc, op-i64-i64.cc, op-i8-i8.cc, op-int-concat.cc, op-int-conv.cc, op-m-cm.cc, op-m-cs.cc, op-m-m.cc, op-m-s.cc, op-m-scm.cc, op-m-sm.cc, op-pm-pm.cc, op-pm-scm.cc, op-pm-sm.cc, op-pm-template.cc, op-range.cc, op-s-cm.cc, op-s-cs.cc, op-s-m.cc, op-s-s.cc, op-s-scm.cc, op-s-sm.cc, op-sbm-b.cc, op-sbm-bm.cc, op-sbm-sbm.cc, op-scm-cm.cc, op-scm-cs.cc, op-scm-m.cc, op-scm-s.cc, op-scm-scm.cc, op-scm-sm.cc, op-sm-cm.cc, op-sm-cs.cc, op-sm-m.cc, op-sm-s.cc, op-sm-scm.cc, op-sm-sm.cc, op-str-m.cc, op-str-s.cc, op-str-str.cc, op-struct.cc, op-ui16-ui16.cc, op-ui32-ui32.cc, op-ui64-ui64.cc, op-ui8-ui8.cc, pt-arg-list.cc, pt-array-list.cc, pt-assign.cc, pt-binop.cc, pt-bp.cc, pt-cbinop.cc, pt-cell.cc, pt-check.cc, pt-classdef.cc, pt-cmd.cc, pt-colon.cc, pt-colon.h, pt-const.cc, pt-decl.cc, pt-eval.cc, pt-except.cc, pt-exp.cc, pt-fcn-handle.cc, pt-funcall.cc, pt-id.cc, pt-idx.cc, pt-jump.cc, pt-loop.cc, pt-mat.cc, pt-misc.cc, pt-pr-code.cc, pt-select.cc, pt-stmt.cc, pt-unop.cc, pt.cc, token.cc, Array-jit.cc, Array-os.cc, Array-sym.cc, Array-tc.cc, version.cc: Indent #ifdef blocks in libinterp.

author	Rik <rik@octave.org>
date	Fri, 05 Feb 2016 16:29:08 -0800
parents	ba0a4b4f021d
children	40de9f8f23a6

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
19697 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 19596 diff changeset	3 Copyright (C) 2000-2015 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: 6945 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: 6945 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: 6945 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: 6945 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 20. Augiust 2000 - Kai Habel: first release
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	25 */
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	26
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	27 /*
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	28 2003-12-14 Rafael Laboissiere <rafael@laboissiere.net>
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	29 Added optional second argument to pass options to the underlying
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	30 qhull command
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
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	33 #ifdef HAVE_CONFIG_H
21200 fcac5dbbf9ed maint: Indent #ifdef blocks in libinterp. Rik <rik@octave.org> parents: 21126 diff changeset	34 # include <config.h>
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	35 #endif
9786 2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	36
2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	37 #include <cstdio>
2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	38
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	39 #include <list>
6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	40
19269 65554f5847ac don't include oct-locbuf.h in header files unnecessarily John W. Eaton <jwe@octave.org> parents: 18077 diff changeset	41 #include "oct-locbuf.h"
9003 0631d397fbe0 replace lo_ieee_isnan by xisnan, add missing includes Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	42 #include "lo-ieee.h"
9786 2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	43
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	44 #include "Cell.h"
9786 2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	45 #include "defun-dld.h"
2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	46 #include "error.h"
21126 ba0a4b4f021d include errwarn.h in files that use err_disabled_feature conditionally John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	47 #include "errwarn.h"
20940 48b2ad5ee801 maint: Rename oct-obj.[cc\|h] to ovl.[cc\|h] for clarity. Rik <rik@octave.org> parents: 20939 diff changeset	48 #include "ovl.h"
14339 3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	49 #include "unwind-prot.h"
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	50
14043 f913363318e0 handle new names and locations of qhull include files (bug #33712) John W. Eaton <jwe@octave.org> parents: 13903 diff changeset	51 #if defined (HAVE_QHULL)
21200 fcac5dbbf9ed maint: Indent #ifdef blocks in libinterp. Rik <rik@octave.org> parents: 21126 diff changeset	52 # include "oct-qhull.h"
fcac5dbbf9ed maint: Indent #ifdef blocks in libinterp. Rik <rik@octave.org> parents: 21126 diff changeset	53 # if defined (NEED_QHULL_VERSION)
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	54 char qh_version[] = "__voronoi__.oct 2007-07-24";
21200 fcac5dbbf9ed maint: Indent #ifdef blocks in libinterp. Rik <rik@octave.org> parents: 21126 diff changeset	55 # endif
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	56 #endif
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	57
14339 3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	58 static void
3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	59 close_fcn (FILE *f)
3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	60 {
3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	61 gnulib::fclose (f);
3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	62 }
3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	63
18077 ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	64 static bool
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	65 octave_qhull_dims_ok (octave_idx_type dim, octave_idx_type n, const char *who)
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	66 {
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	67 if (sizeof (octave_idx_type) > sizeof (int))
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	68 {
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	69 int maxval = std::numeric_limits<int>::max ();
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	70
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	71 if (dim > maxval \|\| n > maxval)
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	72 error ("%s: dimension too large for Qhull", who);
18077 ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	73 }
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	74
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	75 return true;
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	76 }
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	77
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	78 DEFUN_DLD (__voronoi__, args, ,
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	79 "-- texinfo --\n\
20853 1142cf6abc0d 2015 Code Sprint: remove class of function from docstring for all C++ files. Rik <rik@octave.org> parents: 20825 diff changeset	80 @deftypefn {} {@var{C}, @var{F} =} __voronoi__ (@var{caller}, @var{pts})\n\
1142cf6abc0d 2015 Code Sprint: remove class of function from docstring for all C++ files. Rik <rik@octave.org> parents: 20825 diff changeset	81 @deftypefnx {} {@var{C}, @var{F} =} __voronoi__ (@var{caller}, @var{pts}, @var{options})\n\
1142cf6abc0d 2015 Code Sprint: remove class of function from docstring for all C++ files. Rik <rik@octave.org> parents: 20825 diff changeset	82 @deftypefnx {} {@var{C}, @var{F}, @var{Inf_Pts} =} __voronoi__ (@dots{})\n\
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	83 Undocumented internal function.\n\
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	84 @end deftypefn")
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	85 {
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	86 #if defined (HAVE_QHULL)
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	87 int nargin = args.length ();
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	88
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	89 if (nargin < 2 \|\| nargin > 3)
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	90 print_usage ();
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	91
20700 68e3a747ca02 rename octave_value value extractors that accept error message args John W. Eaton <jwe@octave.org> parents: 20579 diff changeset	92 std::string caller = args(0).xstring_value ("__voronoi__: CALLER must be a string");
13903 7b5da9754554 __voronoi__: avoid undeclared variable warning John W. Eaton <jwe@octave.org> parents: 13880 diff changeset	93
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	94 octave_value_list retval;
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	95
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	96 Matrix points = args(1).matrix_value ();
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	97 const octave_idx_type dim = points.columns ();
adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	98 const octave_idx_type num_points = points.rows ();
13746 7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	99
18077 ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	100 if (! octave_qhull_dims_ok (dim, num_points, "__voronoi__"))
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	101 return ovl (0.0);
18077 ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	102
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	103 points = points.transpose ();
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	104
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	105 std::string options;
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	106
16343 dda043ccad7c __voronoi__.cc: Change QHull arguments for 4D inputs to match Matlab. Rik <rik@octave.org> parents: 15195 diff changeset	107 if (dim <= 3)
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	108 options = " Qbb";
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	109 else
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	110 options = " Qbb Qx";
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	111
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	112 if (nargin == 3)
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	113 {
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	114 octave_value opt_arg = args(2);
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	115
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	116 if (opt_arg.is_string ())
13880 1dc120c06ce0 finish previous change to __voronoi__ John W. Eaton <jwe@octave.org> parents: 13879 diff changeset	117 options = " " + opt_arg.string_value ();
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	118 else if (opt_arg.is_empty ())
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	119 ; // Use default options.
13880 1dc120c06ce0 finish previous change to __voronoi__ John W. Eaton <jwe@octave.org> parents: 13879 diff changeset	120 else if (opt_arg.is_cellstr ())
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	121 {
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	122 options = "";
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	123
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	124 Array<std::string> tmp = opt_arg.cellstr_value ();
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	125
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	126 for (octave_idx_type i = 0; i < tmp.numel (); i++)
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	127 options += " " + tmp(i);
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	128 }
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	129 else
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	130 error ("%s: OPTIONS must be a string, cell array of strings, or empty",
66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	131 caller.c_str ());
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	132 }
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	133
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	134 boolT ismalloc = false;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	135
14339 3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	136 unwind_protect frame;
3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	137
14309 824e5d362aba Fix Qhull calling convention by passing true file pointer to qh_new_qhull Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	138 // Replace the outfile pointer with stdout for debugging information.
824e5d362aba Fix Qhull calling convention by passing true file pointer to qh_new_qhull Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	139 #if defined (OCTAVE_HAVE_WINDOWS_FILESYSTEM) && ! defined (OCTAVE_HAVE_POSIX_FILESYSTEM)
824e5d362aba Fix Qhull calling convention by passing true file pointer to qh_new_qhull Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	140 FILE *outfile = gnulib::fopen ("NUL", "w");
824e5d362aba Fix Qhull calling convention by passing true file pointer to qh_new_qhull Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	141 #else
824e5d362aba Fix Qhull calling convention by passing true file pointer to qh_new_qhull Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	142 FILE *outfile = gnulib::fopen ("/dev/null", "w");
824e5d362aba Fix Qhull calling convention by passing true file pointer to qh_new_qhull Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	143 #endif
6880 e00a8f661f06 [project @ 2007-09-07 21:39:55 by dbateman] dbateman parents: 6823 diff changeset	144 FILE *errfile = stderr;
14339 3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	145
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	146 if (! outfile)
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	147 error ("__voronoi__: unable to create temporary file for output");
6880 e00a8f661f06 [project @ 2007-09-07 21:39:55 by dbateman] dbateman parents: 6823 diff changeset	148
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	149 frame.add_fcn (close_fcn, outfile);
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	150
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	151 // qh_new_qhull command and points arguments are not const...
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	152
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	153 std::string cmd = "qhull v" + options;
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	154
adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	155 OCTAVE_LOCAL_BUFFER (char, cmd_str, cmd.length () + 1);
13203 b6aba5b4edb1 voronoin: accept options as a cell array of character strings John W. Eaton <jwe@octave.org> parents: 12805 diff changeset	156
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	157 strcpy (cmd_str, cmd.c_str ());
13203 b6aba5b4edb1 voronoin: accept options as a cell array of character strings John W. Eaton <jwe@octave.org> parents: 12805 diff changeset	158
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	159 int exitcode = qh_new_qhull (dim, num_points, points.fortran_vec (),
adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	160 ismalloc, cmd_str, outfile, errfile);
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	161 if (exitcode)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	162 error ("%s: qhull failed", caller.c_str ());
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11553 diff changeset	163
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	164 // Calling findgood_all provides the number of Voronoi vertices
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	165 // (sets qh num_good).
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	166
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	167 qh_findgood_all (qh facet_list);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	168
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	169 octave_idx_type num_voronoi_regions
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	170 = qh num_vertices - qh_setsize (qh del_vertices);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	171
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	172 octave_idx_type num_voronoi_vertices = qh num_good;
13746 7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	173
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	174 // Find the voronoi centers for all facets.
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	175
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	176 qh_setvoronoi_all ();
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	177
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	178 facetT *facet;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	179 vertexT *vertex;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	180 octave_idx_type k;
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	181
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	182 // Find the number of Voronoi vertices for each Voronoi cell and
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	183 // store them in NI so we can use them later to set the dimensions
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	184 // of the RowVector objects used to collect them.
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 18077 diff changeset	185
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	186 FORALLfacets
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	187 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	188 facet->seen = false;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	189 }
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	190
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	191 OCTAVE_LOCAL_BUFFER (octave_idx_type, ni, num_voronoi_regions);
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	192 for (octave_idx_type i = 0; i < num_voronoi_regions; i++)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	193 ni[i] = 0;
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 18077 diff changeset	194
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	195 k = 0;
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	196
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	197 FORALLvertices
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	198 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	199 if (qh hull_dim == 3)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	200 qh_order_vertexneighbors (vertex);
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	201
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	202 bool infinity_seen = false;
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	203
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	204 facetT neighbor, *neighborp;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	205
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	206 FOREACHneighbor_ (vertex)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	207 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	208 if (neighbor->upperdelaunay)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	209 {
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	210 if (! infinity_seen)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	211 {
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	212 infinity_seen = true;
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	213 ni[k]++;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	214 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	215 }
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	216 else
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	217 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	218 neighbor->seen = true;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	219 ni[k]++;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	220 }
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	221 }
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	222
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	223 k++;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	224 }
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	225
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	226 // If Qhull finds fewer regions than points, we will pad the end
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	227 // of the at_inf and C arrays so that they always contain at least
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	228 // as many elements as the given points array.
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	229
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	230 // FIXME: is it possible (or does it make sense) for
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	231 // num_voronoi_regions to ever be larger than num_points?
13746 7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	232
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	233 octave_idx_type nr = (num_points > num_voronoi_regions
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	234 ? num_points : num_voronoi_regions);
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	235
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	236 boolMatrix at_inf (nr, 1, false);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	237
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	238 // The list of Voronoi vertices. The first element is always
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	239 // Inf.
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	240 Matrix F (num_voronoi_vertices+1, dim);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	241
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	242 for (octave_idx_type d = 0; d < dim; d++)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	243 F(0,d) = octave_Inf;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	244
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	245 // The cell array of vectors of indices into F that represent the
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	246 // vertices of the Voronoi regions (cells).
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	247
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	248 Cell C (nr, 1);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	249
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	250 // Now loop through the list of vertices again and store the
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	251 // coordinates of the Voronoi vertices and the lists of indices
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	252 // for the cells.
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	253
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	254 FORALLfacets
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	255 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	256 facet->seen = false;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	257 }
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	258
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	259 octave_idx_type i = 0;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	260 k = 0;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	261
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	262 FORALLvertices
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	263 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	264 if (qh hull_dim == 3)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	265 qh_order_vertexneighbors (vertex);
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	266
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	267 bool infinity_seen = false;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	268
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	269 octave_idx_type idx = qh_pointid (vertex->point);
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	270
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	271 octave_idx_type num_vertices = ni[k++];
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	272
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	273 // Qhull seems to sometimes produces regions with a single
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	274 // vertex. Is that a bug? How can a region have just one
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	275 // vertex? Let's skip it.
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	276
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	277 if (num_vertices == 1)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	278 continue;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	279
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	280 RowVector facet_list (num_vertices);
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	281
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	282 octave_idx_type m = 0;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	283
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	284 facetT neighbor, *neighborp;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	285
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	286 FOREACHneighbor_(vertex)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	287 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	288 if (neighbor->upperdelaunay)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	289 {
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	290 if (! infinity_seen)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	291 {
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	292 infinity_seen = true;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	293 facet_list(m++) = 1;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	294 at_inf(idx) = true;
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	295 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	296 }
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	297 else
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	298 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	299 if (! neighbor->seen)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	300 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	301 i++;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	302 for (octave_idx_type d = 0; d < dim; d++)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	303 F(i,d) = neighbor->center[d];
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	304
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	305 neighbor->seen = true;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	306 neighbor->visitid = i;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	307 }
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	308
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	309 facet_list(m++) = neighbor->visitid + 1;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	310 }
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	311 }
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	312
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	313 C(idx) = facet_list;
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	314 }
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	315
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	316 retval = ovl (F, C, at_inf);
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	317
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	318 // Free memory from Qhull
13746 7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	319 qh_freeqhull (! qh_ALL);
7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	320
7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	321 int curlong, totlong;
7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	322 qh_memfreeshort (&curlong, &totlong);
7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	323
7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	324 if (curlong \|\| totlong)
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	325 warning ("%s: qhull did not free %d bytes of long memory (%d pieces)",
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	326 caller.c_str (), totlong, curlong);
13746 7ff0bdc3dc4c Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346) Rik <octave@nomad.inbox5.com> parents: 13203 diff changeset	327
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	328 return retval;
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	329
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	330 #else
21126 ba0a4b4f021d include errwarn.h in files that use err_disabled_feature conditionally John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	331 std::string caller
ba0a4b4f021d include errwarn.h in files that use err_disabled_feature conditionally John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	332 = (args.length () > 0
ba0a4b4f021d include errwarn.h in files that use err_disabled_feature conditionally John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	333 ? args(0).xstring_value ("__voronoi__: CALLER must be a string")
ba0a4b4f021d include errwarn.h in files that use err_disabled_feature conditionally John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	334 : std::string ("__voronoi__"));
ba0a4b4f021d include errwarn.h in files that use err_disabled_feature conditionally John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	335
21109 bd1752782e56 Use err_disabled_feature, warn_disabled_feature throughout code base. Rik <rik@octave.org> parents: 20940 diff changeset	336 err_disabled_feature (caller, "Qhull");
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	337 #endif
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	338 }
12805 3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	339
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	340 /*
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	341 ## No test needed for internal helper function.
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	342 %!assert (1)
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	343 */

Mercurial > octave

annotate libinterp/dldfcn/__voronoi__.cc @ 21200:fcac5dbbf9ed

annotate libinterp/dldfcn/voronoi.cc @ 21200:fcac5dbbf9ed