octave: libinterp/dldfcn/__voronoi_

annotate libinterp/dldfcn/voronoi.cc @ 31605:e88a07dec498 stable

maint: Use macros to begin/end C++ namespaces. * oct-conf-post-public.in.h: Define two macros (OCTAVE_BEGIN_NAMESPACE, OCTAVE_END_NAMESPACE) that can be used to start/end a namespace. * mk-opts.pl, build-env.h, build-env.in.cc, __betainc__.cc, __contourc__.cc, __dsearchn__.cc, __eigs__.cc, __expint__.cc, __ftp__.cc, __gammainc__.cc, __ichol__.cc, __ilu__.cc, __isprimelarge__.cc, __lin_interpn__.cc, __magick_read__.cc, __pchip_deriv__.cc, __qp__.cc, amd.cc, auto-shlib.cc, auto-shlib.h, balance.cc, base-text-renderer.cc, base-text-renderer.h, besselj.cc, bitfcns.cc, bsxfun.cc, c-file-ptr-stream.cc, c-file-ptr-stream.h, call-stack.cc, call-stack.h, ccolamd.cc, cellfun.cc, chol.cc, colamd.cc, colloc.cc, conv2.cc, daspk.cc, dasrt.cc, dassl.cc, data.cc, data.h, debug.cc, defaults.cc, defaults.h, defun-int.h, defun.cc, det.cc, dirfns.cc, display.cc, display.h, dlmread.cc, dmperm.cc, dot.cc, dynamic-ld.cc, dynamic-ld.h, eig.cc, ellipj.cc, environment.cc, environment.h, error.cc, error.h, errwarn.h, event-manager.cc, event-manager.h, event-queue.cc, event-queue.h, fcn-info.cc, fcn-info.h, fft.cc, fft2.cc, fftn.cc, file-io.cc, filter.cc, find.cc, ft-text-renderer.cc, ft-text-renderer.h, gcd.cc, getgrent.cc, getpwent.cc, getrusage.cc, givens.cc, gl-render.cc, gl-render.h, gl2ps-print.cc, gl2ps-print.h, graphics-toolkit.cc, graphics-toolkit.h, graphics.cc, graphics.in.h, gsvd.cc, gtk-manager.cc, gtk-manager.h, hash.cc, help.cc, help.h, hess.cc, hex2num.cc, hook-fcn.cc, hook-fcn.h, input.cc, input.h, interpreter-private.cc, interpreter-private.h, interpreter.cc, interpreter.h, inv.cc, jsondecode.cc, jsonencode.cc, kron.cc, latex-text-renderer.cc, latex-text-renderer.h, load-path.cc, load-path.h, load-save.cc, load-save.h, lookup.cc, ls-ascii-helper.cc, ls-ascii-helper.h, ls-oct-text.cc, ls-utils.cc, ls-utils.h, lsode.cc, lu.cc, mappers.cc, matrix_type.cc, max.cc, mex-private.h, mex.cc, mgorth.cc, nproc.cc, oct-fstrm.cc, oct-fstrm.h, oct-hdf5-types.cc, oct-hdf5-types.h, oct-hist.cc, oct-hist.h, oct-iostrm.cc, oct-iostrm.h, oct-opengl.h, oct-prcstrm.cc, oct-prcstrm.h, oct-procbuf.cc, oct-procbuf.h, oct-process.cc, oct-process.h, oct-stdstrm.h, oct-stream.cc, oct-stream.h, oct-strstrm.cc, oct-strstrm.h, oct-tex-lexer.in.ll, oct-tex-parser.yy, ordqz.cc, ordschur.cc, pager.cc, pager.h, pinv.cc, pow2.cc, pr-flt-fmt.cc, pr-output.cc, procstream.cc, procstream.h, psi.cc, qr.cc, quad.cc, quadcc.cc, qz.cc, rand.cc, rcond.cc, regexp.cc, schur.cc, settings.cc, settings.h, sighandlers.cc, sighandlers.h, sparse-xdiv.cc, sparse-xdiv.h, sparse-xpow.cc, sparse-xpow.h, sparse.cc, spparms.cc, sqrtm.cc, stack-frame.cc, stack-frame.h, stream-euler.cc, strfind.cc, strfns.cc, sub2ind.cc, svd.cc, sylvester.cc, symbfact.cc, syminfo.cc, syminfo.h, symrcm.cc, symrec.cc, symrec.h, symscope.cc, symscope.h, symtab.cc, symtab.h, syscalls.cc, sysdep.cc, sysdep.h, text-engine.cc, text-engine.h, text-renderer.cc, text-renderer.h, time.cc, toplev.cc, tril.cc, tsearch.cc, typecast.cc, url-handle-manager.cc, url-handle-manager.h, urlwrite.cc, utils.cc, utils.h, variables.cc, variables.h, xdiv.cc, xdiv.h, xnorm.cc, xnorm.h, xpow.cc, xpow.h, __delaunayn__.cc, __fltk_uigetfile__.cc, __glpk__.cc, __init_fltk__.cc, __init_gnuplot__.cc, __ode15__.cc, __voronoi__.cc, audiodevinfo.cc, audioread.cc, convhulln.cc, fftw.cc, gzip.cc, mk-build-env-features.sh, mk-builtins.pl, cdef-class.cc, cdef-class.h, cdef-fwd.h, cdef-manager.cc, cdef-manager.h, cdef-method.cc, cdef-method.h, cdef-object.cc, cdef-object.h, cdef-package.cc, cdef-package.h, cdef-property.cc, cdef-property.h, cdef-utils.cc, cdef-utils.h, ov-base.cc, ov-base.h, ov-bool-mat.cc, ov-builtin.h, ov-cell.cc, ov-class.cc, ov-class.h, ov-classdef.cc, ov-classdef.h, ov-complex.cc, ov-fcn-handle.cc, ov-fcn-handle.h, ov-fcn.h, ov-java.cc, ov-java.h, ov-mex-fcn.h, ov-null-mat.cc, ov-oncleanup.cc, ov-struct.cc, ov-typeinfo.cc, ov-typeinfo.h, ov-usr-fcn.cc, ov-usr-fcn.h, ov.cc, ov.h, octave.cc, octave.h, mk-ops.sh, 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-fcdm-fcdm.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-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-m-cm.cc, op-m-cs.cc, op-m-m.cc, op-m-s.cc, op-m-scm.cc, op-m-sm.cc, op-mi.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, ops.h, anon-fcn-validator.cc, anon-fcn-validator.h, bp-table.cc, bp-table.h, comment-list.cc, comment-list.h, filepos.h, lex.h, lex.ll, oct-lvalue.cc, oct-lvalue.h, oct-parse.yy, parse.h, profiler.cc, profiler.h, pt-anon-scopes.cc, pt-anon-scopes.h, pt-arg-list.cc, pt-arg-list.h, pt-args-block.cc, pt-args-block.h, pt-array-list.cc, pt-array-list.h, pt-assign.cc, pt-assign.h, pt-binop.cc, pt-binop.h, pt-bp.cc, pt-bp.h, pt-cbinop.cc, pt-cbinop.h, pt-cell.cc, pt-cell.h, pt-check.cc, pt-check.h, pt-classdef.cc, pt-classdef.h, pt-cmd.h, pt-colon.cc, pt-colon.h, pt-const.cc, pt-const.h, pt-decl.cc, pt-decl.h, pt-eval.cc, pt-eval.h, pt-except.cc, pt-except.h, pt-exp.cc, pt-exp.h, pt-fcn-handle.cc, pt-fcn-handle.h, pt-id.cc, pt-id.h, pt-idx.cc, pt-idx.h, pt-jump.h, pt-loop.cc, pt-loop.h, pt-mat.cc, pt-mat.h, pt-misc.cc, pt-misc.h, pt-pr-code.cc, pt-pr-code.h, pt-select.cc, pt-select.h, pt-spmd.cc, pt-spmd.h, pt-stmt.cc, pt-stmt.h, pt-tm-const.cc, pt-tm-const.h, pt-unop.cc, pt-unop.h, pt-vm-eval.cc, pt-walk.cc, pt-walk.h, pt.cc, pt.h, token.cc, token.h, Range.cc, Range.h, idx-vector.cc, idx-vector.h, range-fwd.h, CollocWt.cc, CollocWt.h, aepbalance.cc, aepbalance.h, chol.cc, chol.h, gepbalance.cc, gepbalance.h, gsvd.cc, gsvd.h, hess.cc, hess.h, lo-mappers.cc, lo-mappers.h, lo-specfun.cc, lo-specfun.h, lu.cc, lu.h, oct-convn.cc, oct-convn.h, oct-fftw.cc, oct-fftw.h, oct-norm.cc, oct-norm.h, oct-rand.cc, oct-rand.h, oct-spparms.cc, oct-spparms.h, qr.cc, qr.h, qrp.cc, qrp.h, randgamma.cc, randgamma.h, randmtzig.cc, randmtzig.h, randpoisson.cc, randpoisson.h, schur.cc, schur.h, sparse-chol.cc, sparse-chol.h, sparse-lu.cc, sparse-lu.h, sparse-qr.cc, sparse-qr.h, svd.cc, svd.h, child-list.cc, child-list.h, dir-ops.cc, dir-ops.h, file-ops.cc, file-ops.h, file-stat.cc, file-stat.h, lo-sysdep.cc, lo-sysdep.h, lo-sysinfo.cc, lo-sysinfo.h, mach-info.cc, mach-info.h, oct-env.cc, oct-env.h, oct-group.cc, oct-group.h, oct-password.cc, oct-password.h, oct-syscalls.cc, oct-syscalls.h, oct-time.cc, oct-time.h, oct-uname.cc, oct-uname.h, action-container.cc, action-container.h, base-list.h, cmd-edit.cc, cmd-edit.h, cmd-hist.cc, cmd-hist.h, f77-fcn.h, file-info.cc, file-info.h, lo-array-errwarn.cc, lo-array-errwarn.h, lo-hash.cc, lo-hash.h, lo-ieee.h, lo-regexp.cc, lo-regexp.h, lo-utils.cc, lo-utils.h, oct-base64.cc, oct-base64.h, oct-glob.cc, oct-glob.h, oct-inttypes.h, oct-mutex.cc, oct-mutex.h, oct-refcount.h, oct-shlib.cc, oct-shlib.h, oct-sparse.cc, oct-sparse.h, oct-string.h, octave-preserve-stream-state.h, pathsearch.cc, pathsearch.h, quit.cc, quit.h, unwind-prot.cc, unwind-prot.h, url-transfer.cc, url-transfer.h : Use new macros to begin/end C++ namespaces.

author	Rik <rik@octave.org>
date	Thu, 01 Dec 2022 14:23:45 -0800
parents	796f54d4ddbf
children	aac27ad79be6

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 //
30564 796f54d4ddbf update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 30390 diff changeset	3 // Copyright (C) 2000-2022 The Octave Project Developers
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	4 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	5 // See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	6 // distribution or <https://octave.org/copyright/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	7 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	8 // This file is part of Octave.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	9 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	10 // Octave is free software: you can redistribute it and/or modify it
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	11 // under the terms of the GNU General Public License as published by
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	12 // the Free Software Foundation, either version 3 of the License, or
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	13 // (at your option) any later version.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	14 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	15 // Octave is distributed in the hope that it will be useful, but
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	16 // WITHOUT ANY WARRANTY; without even the implied warranty of
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	18 // GNU General Public License for more details.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	19 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	20 // You should have received a copy of the GNU General Public License
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	21 // along with Octave; see the file COPYING. If not, see
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	22 // <https://www.gnu.org/licenses/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	23 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	24 ////////////////////////////////////////////////////////////////////////
6823 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 20. Augiust 2000 - Kai Habel: first release
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	28 */
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	29
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 2003-12-14 Rafael Laboissiere <rafael@laboissiere.net>
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	32 Added optional second argument to pass options to the underlying
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	33 qhull command
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	34 */
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	35
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21723 diff changeset	36 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21200 diff changeset	37 # include "config.h"
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	38 #endif
9786 2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	39
2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	40 #include <cstdio>
2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	41
23024 a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	42 #include <limits>
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	43 #include <string>
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	44
23024 a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	45 #include "Array.h"
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	46 #include "boolMatrix.h"
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	47 #include "dMatrix.h"
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	48 #include "dRowVector.h"
19269 65554f5847ac don't include oct-locbuf.h in header files unnecessarily John W. Eaton <jwe@octave.org> parents: 18077 diff changeset	49 #include "oct-locbuf.h"
23024 a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	50 #include "unwind-prot.h"
9786 2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	51
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	52 #include "Cell.h"
9786 2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	53 #include "defun-dld.h"
2c279308f6ab fix includes in some src/DLD-FUNCTIONS files John W. Eaton <jwe@octave.org> parents: 9003 diff changeset	54 #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	55 #include "errwarn.h"
23024 a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	56 #include "ov.h"
20940 48b2ad5ee801 maint: Rename oct-obj.[cc\|h] to ovl.[cc\|h] for clarity. Rik <rik@octave.org> parents: 20939 diff changeset	57 #include "ovl.h"
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	58
14043 f913363318e0 handle new names and locations of qhull include files (bug #33712) John W. Eaton <jwe@octave.org> parents: 13903 diff changeset	59 #if defined (HAVE_QHULL)
21200 fcac5dbbf9ed maint: Indent #ifdef blocks in libinterp. Rik <rik@octave.org> parents: 21126 diff changeset	60 # include "oct-qhull.h"
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	61 # if defined (NEED_QHULL_R_VERSION)
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	62 char qh_version[] = "__voronoi__.oct 2007-07-24";
21200 fcac5dbbf9ed maint: Indent #ifdef blocks in libinterp. Rik <rik@octave.org> parents: 21126 diff changeset	63 # endif
29969 7a4f16bff8fd use fewer OCTAVE_NAMESPACE_BEGIN/END blocks where possible John W. Eaton <jwe@octave.org> parents: 29964 diff changeset	64 #endif
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	65
31605 e88a07dec498 maint: Use macros to begin/end C++ namespaces. Rik <rik@octave.org> parents: 30564 diff changeset	66 OCTAVE_BEGIN_NAMESPACE(octave)
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29629 diff changeset	67
29969 7a4f16bff8fd use fewer OCTAVE_NAMESPACE_BEGIN/END blocks where possible John W. Eaton <jwe@octave.org> parents: 29964 diff changeset	68 #if defined (HAVE_QHULL)
7a4f16bff8fd use fewer OCTAVE_NAMESPACE_BEGIN/END blocks where possible John W. Eaton <jwe@octave.org> parents: 29964 diff changeset	69
14339 3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	70 static void
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	71 free_qhull_memory (qhT *qh)
24260 0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	72 {
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	73 qh_freeqhull (qh, ! qh_ALL);
24260 0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	74
0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	75 int curlong, totlong;
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	76 qh_memfreeshort (qh, &curlong, &totlong);
24260 0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	77
0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	78 if (curlong \|\| totlong)
0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	79 warning ("__voronoi__: did not free %d bytes of long memory (%d pieces)",
0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	80 totlong, curlong);
0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	81 }
0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	82
18077 ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	83 static bool
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	84 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	85 {
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	86 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	87 {
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	88 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	89
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	90 if (dim > maxval \|\| n > maxval)
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	91 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	92 }
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	93
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	94 return true;
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	95 }
ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	96
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	97 #endif
263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	98
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	99 DEFUN_DLD (__voronoi__, args, ,
21966 112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21942 diff changeset	100 doc: /* -- texinfo --
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21942 diff changeset	101 @deftypefn {} {@var{C}, @var{F} =} __voronoi__ (@var{caller}, @var{pts})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21942 diff changeset	102 @deftypefnx {} {@var{C}, @var{F} =} __voronoi__ (@var{caller}, @var{pts}, @var{options})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21942 diff changeset	103 @deftypefnx {} {@var{C}, @var{F}, @var{Inf_Pts} =} __voronoi__ (@dots{})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21942 diff changeset	104 Undocumented internal function.
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21942 diff changeset	105 @end deftypefn */)
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	106 {
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	107 #if defined (HAVE_QHULL)
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	108
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	109 int nargin = args.length ();
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	110
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	111 if (nargin < 2 \|\| nargin > 3)
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	112 print_usage ();
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	113
20700 68e3a747ca02 rename octave_value value extractors that accept error message args John W. Eaton <jwe@octave.org> parents: 20579 diff changeset	114 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	115
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	116 octave_value_list retval;
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	117
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	118 Matrix points = args(1).matrix_value ();
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	119 const octave_idx_type dim = points.columns ();
adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	120 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	121
18077 ac74b0c4c564 avoid overflow when passing problem dimensions to qhull with --enable-64 John W. Eaton <jwe@octave.org> parents: 17787 diff changeset	122 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	123 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	124
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	125 points = points.transpose ();
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	126
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	127 std::string options;
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	128
16343 dda043ccad7c __voronoi__.cc: Change QHull arguments for 4D inputs to match Matlab. Rik <rik@octave.org> parents: 15195 diff changeset	129 if (dim <= 3)
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	130 options = " Qbb";
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	131 else
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	132 options = " Qbb Qx";
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	133
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	134 if (nargin == 3)
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	135 {
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	136 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	137
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	138 if (opt_arg.is_string ())
23807 336f89b6208b Use character literals 'c' rather than string literals "c" when possible. Rik <rik@octave.org> parents: 23577 diff changeset	139 options = ' ' + opt_arg.string_value ();
23577 80c42f4cca13 maint: Deprecate is_empty and replace with isempty. Rik <rik@octave.org> parents: 23575 diff changeset	140 else if (opt_arg.isempty ())
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	141 ; // Use default options.
23575 e95738a119da maint: Deprecate is_cellstr and replace with iscellstr. Rik <rik@octave.org> parents: 23220 diff changeset	142 else if (opt_arg.iscellstr ())
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	143 {
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	144 options = "";
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	145
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	146 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	147
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	148 for (octave_idx_type i = 0; i < tmp.numel (); i++)
23807 336f89b6208b Use character literals 'c' rather than string literals "c" when possible. Rik <rik@octave.org> parents: 23577 diff changeset	149 options += ' ' + tmp(i);
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	150 }
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	151 else
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	152 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	153 caller.c_str ());
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	154 }
440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	155
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	156 boolT ismalloc = false;
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	157
29517 78ccd8bf439c Suppress extraneous error messages from Qhull (bug #57727). Rik <rik@octave.org> parents: 29359 diff changeset	158 // Set the outfile pointer to stdout for status information.
78ccd8bf439c Suppress extraneous error messages from Qhull (bug #57727). Rik <rik@octave.org> parents: 29359 diff changeset	159 FILE *outfile = nullptr;
6880 e00a8f661f06 [project @ 2007-09-07 21:39:55 by dbateman] dbateman parents: 6823 diff changeset	160 FILE *errfile = stderr;
14339 3e4350f09a55 close temporary files opened for Qhull John W. Eaton <jwe@octave.org> parents: 14309 diff changeset	161
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	162 qhT context = { };
93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	163 qhT *qh = &context;
93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	164
13879 440d7914cf01 fix regression in __voronoi__ and convhulln option processing John W. Eaton <jwe@octave.org> parents: 13871 diff changeset	165 std::string cmd = "qhull v" + options;
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	166
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	167 int exitcode = qh_new_qhull (qh, dim, num_points, points.fortran_vec (),
93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	168 ismalloc, &cmd[0], outfile, errfile);
13203 b6aba5b4edb1 voronoin: accept options as a cell array of character strings John W. Eaton <jwe@octave.org> parents: 12805 diff changeset	169
29961 7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	170 unwind_action free_memory ([qh] () { free_qhull_memory (qh); });
24260 0f2dc8d6c34d Eliminate possible segfaults related to not cleaning up Qhull workspace. Rik <rik@octave.org> parents: 23219 diff changeset	171
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	172 if (exitcode)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	173 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	174
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	175 // Calling findgood_all provides the number of Voronoi vertices
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	176 // (sets qh->num_good).
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	177
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	178 qh_findgood_all (qh, qh->facet_list);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	179
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	180 octave_idx_type num_voronoi_regions
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	181 = qh->num_vertices - qh_setsize (qh, qh->del_vertices);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	182
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	183 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	184
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	185 // Find the voronoi centers for all facets.
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	186
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	187 qh_setvoronoi_all (qh);
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	188
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	189 facetT *facet;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	190 vertexT *vertex;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	191 octave_idx_type k;
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	192
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	193 // 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	194 // 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	195 // 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	196
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	197 FORALLfacets
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 facet->seen = false;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	200 }
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	201
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	202 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	203 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	204 ni[i] = 0;
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 18077 diff changeset	205
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	206 k = 0;
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	207
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	208 FORALLvertices
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	209 {
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	210 if (qh->hull_dim == 3)
93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	211 qh_order_vertexneighbors (qh, vertex);
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	212
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	213 bool infinity_seen = false;
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	214
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	215 facetT neighbor, *neighborp;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	216
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	217 FOREACHneighbor_ (vertex)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	218 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	219 if (neighbor->upperdelaunay)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	220 {
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	221 if (! infinity_seen)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	222 {
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	223 infinity_seen = true;
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	224 ni[k]++;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	225 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	226 }
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	227 else
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	228 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	229 neighbor->seen = true;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	230 ni[k]++;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	231 }
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	232 }
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	233
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	234 k++;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	235 }
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	236
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	237 // 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	238 // 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	239 // 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	240
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	241 // 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	242 // 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	243
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	244 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	245 ? num_points : num_voronoi_regions);
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	246
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	247 boolMatrix at_inf (nr, 1, false);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	248
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	249 // 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	250 // Inf.
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	251 Matrix F (num_voronoi_vertices+1, dim);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	252
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	253 for (octave_idx_type d = 0; d < dim; d++)
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 29969 diff changeset	254 F(0, d) = numeric_limits<double>::Inf ();
20866 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 // 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	257 // vertices of the Voronoi regions (cells).
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	258
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	259 Cell C (nr, 1);
13862 6d7e133a4bed compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13746 diff changeset	260
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	261 // 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	262 // 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	263 // for the cells.
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	264
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	265 FORALLfacets
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 facet->seen = false;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	268 }
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	269
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	270 octave_idx_type i = 0;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	271 k = 0;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	272
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	273 FORALLvertices
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	274 {
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	275 if (qh->hull_dim == 3)
93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	276 qh_order_vertexneighbors (qh, vertex);
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	277
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	278 bool infinity_seen = false;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	279
29629 93c8df989ea0 qhull: Use reentrant libqhull_r (bug #60016). Stefan Brüns <stefan.bruens@rwth-aachen.de> parents: 29517 diff changeset	280 octave_idx_type idx = qh_pointid (qh, vertex->point);
20866 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 num_vertices = ni[k++];
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	283
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	284 // 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	285 // 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	286 // vertex? Let's skip it.
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 (num_vertices == 1)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	289 continue;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	290
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	291 RowVector facet_list (num_vertices);
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	292
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	293 octave_idx_type m = 0;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	294
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	295 facetT neighbor, *neighborp;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	296
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	297 FOREACHneighbor_(vertex)
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->upperdelaunay)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	300 {
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	301 if (! infinity_seen)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	302 {
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	303 infinity_seen = true;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	304 facet_list(m++) = 1;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	305 at_inf(idx) = true;
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	306 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	307 }
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	308 else
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	309 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	310 if (! neighbor->seen)
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	311 {
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	312 i++;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	313 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	314 F(i,d) = neighbor->center[d];
13871 adf60d6dc1dd more compatibility fixes for __voronoi__ John W. Eaton <jwe@octave.org> parents: 13862 diff changeset	315
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	316 neighbor->seen = true;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	317 neighbor->visitid = i;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	318 }
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	319
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	320 facet_list(m++) = neighbor->visitid + 1;
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	321 }
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9786 diff changeset	322 }
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	323
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	324 C(idx) = facet_list;
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	325 }
20866 2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	326
2bd3b13e2c8e 2015 Code Sprint: __voronoi__.cc: use ovl(). Rik <rik@octave.org> parents: 20853 diff changeset	327 retval = ovl (F, C, at_inf);
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6880 diff changeset	328
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	329 return retval;
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20866 diff changeset	330
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	331 #else
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	332
263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	333 octave_unused_parameter (args);
21937 55f7de37b618 maint: Strip trailing whitespace from source files. John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	334
21126 ba0a4b4f021d include errwarn.h in files that use err_disabled_feature conditionally John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	335 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	336 = (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	337 ? args(0).xstring_value ("__voronoi__: CALLER must be a string")
23829 01899bdd2a3a Eliminate unnecessary std::string ("...") constructor calls when "..." suffices. Rik <rik@octave.org> parents: 23807 diff changeset	338 : "__voronoi__");
21126 ba0a4b4f021d include errwarn.h in files that use err_disabled_feature conditionally John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	339
21109 bd1752782e56 Use err_disabled_feature, warn_disabled_feature throughout code base. Rik <rik@octave.org> parents: 20940 diff changeset	340 err_disabled_feature (caller, "Qhull");
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	341
6823 9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	342 #endif
9fddcc586065 [project @ 2007-08-24 08:27:27 by dbateman] dbateman parents: diff changeset	343 }
12805 3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	344
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	345 /*
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	346 ## 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	347 %!assert (1)
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	348 */
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29629 diff changeset	349
31605 e88a07dec498 maint: Use macros to begin/end C++ namespaces. Rik <rik@octave.org> parents: 30564 diff changeset	350 OCTAVE_END_NAMESPACE(octave)

Mercurial > octave

annotate libinterp/dldfcn/__voronoi__.cc @ 31605:e88a07dec498 stable

annotate libinterp/dldfcn/voronoi.cc @ 31605:e88a07dec498 stable