Mercurial > octave-nkf
annotate src/DLD-FUNCTIONS/__voronoi__.cc @ 13871:adf60d6dc1dd
more compatibility fixes for __voronoi__
* __voronoi__.cc (F__voronoi__): Use v option for qhull, not Fv.
Improve argument handling. Call qh_findgood_all to obtain number of
Voronoi vertices directly. Correctly size output values. Use
qh_pointid to place cell indices in the same order as the given points.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Wed, 16 Nov 2011 16:17:13 -0500 |
parents | 6d7e133a4bed |
children | 440d7914cf01 |
rev | line source |
---|---|
6823 | 1 /* |
2 | |
11523 | 3 Copyright (C) 2000-2011 Kai Habel |
6823 | 4 |
5 This file is part of Octave. | |
6 | |
7 Octave is free software; you can redistribute it and/or modify it | |
8 under the terms of the GNU General Public License as published by the | |
7016 | 9 Free Software Foundation; either version 3 of the License, or (at your |
10 option) any later version. | |
6823 | 11 |
12 Octave is distributed in the hope that it will be useful, but WITHOUT | |
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
15 for more details. | |
16 | |
17 You should have received a copy of the GNU General Public License | |
7016 | 18 along with Octave; see the file COPYING. If not, see |
19 <http://www.gnu.org/licenses/>. | |
6823 | 20 |
21 */ | |
22 | |
23 /* | |
24 20. Augiust 2000 - Kai Habel: first release | |
25 */ | |
26 | |
27 /* | |
28 2003-12-14 Rafael Laboissiere <rafael@laboissiere.net> | |
29 Added optional second argument to pass options to the underlying | |
30 qhull command | |
31 */ | |
32 | |
33 #ifdef HAVE_CONFIG_H | |
34 #include <config.h> | |
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 |
9003
0631d397fbe0
replace lo_ieee_isnan by xisnan, add missing includes
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
41 #include "lo-ieee.h" |
9786
2c279308f6ab
fix includes in some src/DLD-FUNCTIONS files
John W. Eaton <jwe@octave.org>
parents:
9003
diff
changeset
|
42 |
6823 | 43 #include "Cell.h" |
9786
2c279308f6ab
fix includes in some src/DLD-FUNCTIONS files
John W. Eaton <jwe@octave.org>
parents:
9003
diff
changeset
|
44 #include "defun-dld.h" |
2c279308f6ab
fix includes in some src/DLD-FUNCTIONS files
John W. Eaton <jwe@octave.org>
parents:
9003
diff
changeset
|
45 #include "error.h" |
2c279308f6ab
fix includes in some src/DLD-FUNCTIONS files
John W. Eaton <jwe@octave.org>
parents:
9003
diff
changeset
|
46 #include "oct-obj.h" |
6823 | 47 |
48 #ifdef HAVE_QHULL | |
49 extern "C" { | |
50 #include <qhull/qhull_a.h> | |
51 } | |
52 | |
53 #ifdef NEED_QHULL_VERSION | |
54 char qh_version[] = "__voronoi__.oct 2007-07-24"; | |
55 #endif | |
56 #endif | |
57 | |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
58 DEFUN_DLD (__voronoi__, args, , |
6823 | 59 "-*- texinfo -*-\n\ |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
60 @deftypefn {Loadable Function} {@var{C}, @var{F} =} __voronoi__ (@var{pts})\n\ |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
61 @deftypefnx {Loadable Function} {@var{C}, @var{F} =} __voronoi__ (@var{pts}, @var{options})\n\ |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
62 @deftypefnx {Loadable Function} {@var{C}, @var{F}, @var{Inf_Pts} =} __voronoi__ (@dots{})\n\ |
6945 | 63 Undocumented internal function.\n\ |
6823 | 64 @end deftypefn") |
65 { | |
66 octave_value_list retval; | |
6945 | 67 |
6823 | 68 #ifdef HAVE_QHULL |
6945 | 69 |
6823 | 70 retval(0) = 0.0; |
71 | |
6945 | 72 int nargin = args.length (); |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11553
diff
changeset
|
73 if (nargin < 1 || nargin > 2) |
6823 | 74 { |
75 print_usage (); | |
76 return retval; | |
77 } | |
78 | |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
79 std::string cmd = "qhull v"; |
6880 | 80 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
81 if (nargin == 2 && ! args(1).is_empty ()) |
6823 | 82 { |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
83 if (args(1).is_string ()) |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
84 cmd += " " + args(1).string_value (); |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
85 else if (args(1).is_cellstr ()) |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
86 { |
13203
b6aba5b4edb1
voronoin: accept options as a cell array of character strings
John W. Eaton <jwe@octave.org>
parents:
12805
diff
changeset
|
87 Array<std::string> tmp = args(1).cellstr_value (); |
b6aba5b4edb1
voronoin: accept options as a cell array of character strings
John W. Eaton <jwe@octave.org>
parents:
12805
diff
changeset
|
88 |
b6aba5b4edb1
voronoin: accept options as a cell array of character strings
John W. Eaton <jwe@octave.org>
parents:
12805
diff
changeset
|
89 for (octave_idx_type i = 0; i < tmp.numel (); i++) |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
90 cmd += " " + tmp(i); |
13203
b6aba5b4edb1
voronoin: accept options as a cell array of character strings
John W. Eaton <jwe@octave.org>
parents:
12805
diff
changeset
|
91 } |
b6aba5b4edb1
voronoin: accept options as a cell array of character strings
John W. Eaton <jwe@octave.org>
parents:
12805
diff
changeset
|
92 else |
b6aba5b4edb1
voronoin: accept options as a cell array of character strings
John W. Eaton <jwe@octave.org>
parents:
12805
diff
changeset
|
93 { |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
94 error ("__voronoi__: OPTIONS argument must be a string, cell array of strings, or empty"); |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
95 return retval; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
96 } |
6823 | 97 } |
98 | |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
99 Matrix points (args(0).matrix_value ()); |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
100 const octave_idx_type dim = points.columns (); |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
101 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
|
102 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
103 points = points.transpose (); |
6945 | 104 |
6823 | 105 boolT ismalloc = false; |
106 | |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
107 // Replace the 0 pointer with stdout for debugging information |
6945 | 108 FILE *outfile = 0; |
6880 | 109 FILE *errfile = stderr; |
110 | |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
111 // Qhull flags and points arguments are not const... |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
112 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
113 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
|
114 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
115 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
|
116 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
117 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
|
118 ismalloc, cmd_str, outfile, errfile); |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
119 if (! exitcode) |
6823 | 120 { |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
121 // Calling findgood_all provides the number of Voronoi vertices |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
122 // (sets qh num_good). |
6823 | 123 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
124 qh_findgood_all (qh facet_list); |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11553
diff
changeset
|
125 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
126 octave_idx_type num_voronoi_regions |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
127 = qh num_vertices - qh_setsize (qh del_vertices); |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
128 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
129 octave_idx_type num_voronoi_vertices = qh num_good; |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
130 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
131 // Find the voronoi centers for all facets. |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
132 |
6945 | 133 qh_setvoronoi_all (); |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
134 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
135 facetT *facet; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
136 vertexT *vertex; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
137 octave_idx_type k; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
138 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
139 // Find the number of Voronoi vertices for each Voronoi cell and |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
140 // store them in NI so we can use them later to set the dimensions |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
141 // of the RowVector objects used to collect them. |
6945 | 142 |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11553
diff
changeset
|
143 FORALLfacets |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
144 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
145 facet->seen = false; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
146 } |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
147 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
148 OCTAVE_LOCAL_BUFFER (octave_idx_type, ni, num_voronoi_regions); |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
149 for (octave_idx_type i = 0; i < num_voronoi_regions; i++) |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
150 ni[i] = 0; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
151 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
152 k = 0; |
6945 | 153 |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11553
diff
changeset
|
154 FORALLvertices |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
155 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
156 if (qh hull_dim == 3) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
157 qh_order_vertexneighbors (vertex); |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
158 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
159 bool infinity_seen = false; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
160 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
161 facetT *neighbor, **neighborp; |
6945 | 162 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
163 FOREACHneighbor_ (vertex) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
164 { |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
165 if (neighbor->upperdelaunay) |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
166 { |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
167 if (! infinity_seen) |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
168 { |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
169 infinity_seen = true; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
170 ni[k]++; |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
171 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
172 } |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
173 else |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
174 { |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
175 neighbor->seen = true; |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
176 ni[k]++; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
177 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
178 } |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
179 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
180 k++; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
181 } |
6823 | 182 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
183 // If Qhull finds fewer regions than points, we will pad the end |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
184 // of the at_inf and C arrays so that they always contain at least |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
185 // as many elements as the given points array. |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
186 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
187 // FIXME -- is it possible (or does it make sense) for |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
188 // num_voronoi_regions to ever be larger than num_points? |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
189 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
190 octave_idx_type nr = (num_points > num_voronoi_regions |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
191 ? num_points : num_voronoi_regions); |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
192 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
193 boolMatrix at_inf (nr, 1, false); |
6823 | 194 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
195 // The list of Voronoi vertices. The first element is always |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
196 // Inf. |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
197 Matrix F (num_voronoi_vertices+1, dim); |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
198 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
199 for (octave_idx_type d = 0; d < dim; d++) |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
200 F(0,d) = octave_Inf; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
201 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
202 // The cell array of vectors of indices into F that represent the |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
203 // vertices of the Voronoi regions (cells). |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
204 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
205 Cell C (nr, 1); |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
206 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
207 // Now loop through the list of vertices again and store the |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
208 // coordinates of the Voronoi vertices and the lists of indices |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
209 // for the cells. |
6945 | 210 |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11553
diff
changeset
|
211 FORALLfacets |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
212 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
213 facet->seen = false; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
214 } |
6945 | 215 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
216 octave_idx_type i = 0; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
217 k = 0; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
218 |
6823 | 219 FORALLvertices |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
220 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
221 if (qh hull_dim == 3) |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
222 qh_order_vertexneighbors (vertex); |
13746
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
223 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
224 bool infinity_seen = false; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
225 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
226 octave_idx_type idx = qh_pointid (vertex->point); |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
227 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
228 octave_idx_type num_vertices = ni[k++]; |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
229 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
230 // Qhull seems to sometimes produces regions with a single |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
231 // vertex. Is that a bug? How can a region have just one |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
232 // vertex? Let's skip it. |
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
233 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
234 if (num_vertices == 1) |
13862
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
235 continue; |
6d7e133a4bed
compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13746
diff
changeset
|
236 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
237 RowVector facet_list (num_vertices); |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
238 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
239 octave_idx_type m = 0; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
240 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
241 facetT *neighbor, **neighborp; |
6945 | 242 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
243 FOREACHneighbor_(vertex) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
244 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
245 if (neighbor->upperdelaunay) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
246 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
247 if (! infinity_seen) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
248 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
249 infinity_seen = true; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
250 facet_list(m++) = 1; |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
251 at_inf(idx) = true; |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
252 } |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11553
diff
changeset
|
253 } |
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11553
diff
changeset
|
254 else |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
255 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
256 if (! neighbor->seen) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
257 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
258 i++; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
259 for (octave_idx_type d = 0; d < dim; d++) |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
260 F(i,d) = neighbor->center[d]; |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
261 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
262 neighbor->seen = true; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
263 neighbor->visitid = i; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
264 } |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
265 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
266 facet_list(m++) = neighbor->visitid + 1; |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
267 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
268 } |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
269 |
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
270 C(idx) = facet_list; |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9786
diff
changeset
|
271 } |
6823 | 272 |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
273 retval(2) = at_inf; |
6823 | 274 retval(1) = C; |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
275 retval(0) = F; |
6823 | 276 } |
6945 | 277 else |
278 error ("__voronoi__: qhull failed"); | |
279 | |
13871
adf60d6dc1dd
more compatibility fixes for __voronoi__
John W. Eaton <jwe@octave.org>
parents:
13862
diff
changeset
|
280 // 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
|
281 qh_freeqhull (! qh_ALL); |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
282 |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
283 int curlong, totlong; |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
284 qh_memfreeshort (&curlong, &totlong); |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
285 |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
286 if (curlong || totlong) |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
287 warning ("__voronoi__: did not free %d bytes of long memory (%d pieces)", |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
288 totlong, curlong); |
7ff0bdc3dc4c
Revamp geometry functions dependent on Qhull (Bug #34604, Bug #33346)
Rik <octave@nomad.inbox5.com>
parents:
13203
diff
changeset
|
289 |
6823 | 290 #else |
291 error ("__voronoi__: not available in this version of Octave"); | |
292 #endif | |
6945 | 293 |
6823 | 294 return retval; |
295 } | |
12805
3641167e5b75
codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents:
11586
diff
changeset
|
296 |
3641167e5b75
codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents:
11586
diff
changeset
|
297 /* |
3641167e5b75
codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents:
11586
diff
changeset
|
298 |
3641167e5b75
codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents:
11586
diff
changeset
|
299 ## 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
|
300 %!assert (1) |
3641167e5b75
codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents:
11586
diff
changeset
|
301 |
3641167e5b75
codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents:
11586
diff
changeset
|
302 */ |