Mercurial > hg > octave-lyh

annotate libinterp/corefcn/__contourc__.cc @ 17535:c12c688a35ed default tip lyh

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix warnings
author LYH <lyh.kernel@gmail.com>
date Fri, 27 Sep 2013 17:43:27 +0800
parents 704e15f8fecd
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
1 /* Contour lines for function evaluated on a grid.
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
2
14138
72c96de7a403 maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents: 12805
diff changeset
3 Copyright (C) 2007-2012 Kai Habel
7017
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe]
jwe
parents: 7016
diff changeset
4 Copyright (C) 2004, 2007 Shai Ayal
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
5
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
6 Adapted to an oct file from the stand alone contourl by Victro Munoz
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
7 Copyright (C) 2004 Victor Munoz
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
8
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
9 Based on contour plot routine (plcont.c) in PLPlot package
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
10 http://plplot.org/
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
11
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
12 Copyright (C) 1995, 2000, 2001 Maurice LeBrun
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
13 Copyright (C) 2000, 2002 Joao Cardoso
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
14 Copyright (C) 2000, 2001, 2002, 2004 Alan W. Irwin
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
15 Copyright (C) 2004 Andrew Ross
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
16
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
17 This file is part of Octave.
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
18
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
19 Octave is free software; you can redistribute it and/or modify it
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
20 under the terms of the GNU General Public License as published by the
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
21 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
22 option) any later version.
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
23
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
24 Octave is distributed in the hope that it will be useful, but WITHOUT
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
25 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
26 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
27 for more details.
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
28
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
29 You should have received a copy of the GNU General Public License
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
30 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
31 <http://www.gnu.org/licenses/>.
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
32
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
33 */
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
34
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
35 #ifdef HAVE_CONFIG_H
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
36 #include <config.h>
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
37 #endif
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
38
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
39 #include <cfloat>
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
40
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
41 #include "quit.h"
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
42
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14501
diff changeset
43 #include "defun.h"
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
44 #include "error.h"
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
45 #include "oct-obj.h"
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
46
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
47 static Matrix this_contour;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
48 static Matrix contourc;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
49 static int elem;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
50
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
51 // This is the quanta in which we increase this_contour.
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
52 #define CONTOUR_QUANT 50
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
53
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
54 // Add a coordinate point (x,y) to this_contour.
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
55
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
56 static void
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
57 add_point (double x, double y)
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
58 {
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
59 if (elem % CONTOUR_QUANT == 0)
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
60 this_contour = this_contour.append (Matrix (2, CONTOUR_QUANT, 0));
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
61
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
62 this_contour (0, elem) = x;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
63 this_contour (1, elem) = y;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
64 elem++;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
65 }
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
66
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
67 // Add contents of current contour to contourc.
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
68 // this_contour.cols () - 1;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
69
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
70 static void
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
71 end_contour (void)
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
72 {
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
73 if (elem > 2)
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
74 {
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
75 this_contour (1, 0) = elem - 1;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
76 contourc = contourc.append (this_contour.extract_n (0, 0, 2, elem));
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
77 }
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
78
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
79 this_contour = Matrix ();
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
80 elem = 0;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
81 }
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
82
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
83 // Start a new contour, and add contents of current one to contourc.
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
84
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
85 static void
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
86 start_contour (double lvl, double x, double y)
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
87 {
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
88 end_contour ();
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
89 this_contour.resize (2, 0);
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
90 add_point (lvl, 0);
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
91 add_point (x, y);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
92 }
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
93
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
94 static void
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
95 drawcn (const RowVector& X, const RowVector& Y, const Matrix& Z,
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
96 double lvl, int r, int c, double ct_x, double ct_y,
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
97 unsigned int start_edge, bool first, charMatrix& mark)
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
98 {
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
99 double px[4], py[4], pz[4], tmp;
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
100 unsigned int stop_edge, pt[2];
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
101
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
102 // Continue while next facet is not done yet.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
103 while (r >= 0 && c >= 0 && r < mark.rows () && c < mark.cols ()
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
104 && mark(r, c) > 0)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
105 {
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
106
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
107 //get x, y, and z - lvl for current facet
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
108 px[0] = px[3] = X(c);
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
109 px[1] = px[2] = X(c+1);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
110
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
111 py[0] = py[1] = Y(r);
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
112 py[2] = py[3] = Y(r+1);
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
113
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
114 pz[3] = Z(r+1, c) - lvl;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
115 pz[2] = Z(r+1, c + 1) - lvl;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
116 pz[1] = Z(r, c+1) - lvl;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
117 pz[0] = Z(r, c) - lvl;
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
118
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
119 // Facet edge and point naming assignment.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
120 //
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
121 // 0-----1 .-0-.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
122 // | | | |
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
123 // | | 3 1
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
124 // | | | |
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
125 // 3-----2 .-2-.
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
126
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
127 // Get mark value of current facet.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
128 char id = static_cast<char> (mark(r, c));
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
129
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
130 // Check startedge s.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
131 if (start_edge == 255)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
132 {
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
133 // Find start edge.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
134 for (unsigned int k = 0; k < 4; k++)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
135 if (static_cast<char> (1 << k) & id)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
136 start_edge = k;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
137 }
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
138
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
139 if (start_edge == 255)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
140 break;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
141
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
142 // Decrease mark value of current facet for start edge.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
143 mark(r, c) -= static_cast<char> (1 << start_edge);
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
144
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
145 // Next point (clockwise).
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
146 pt[0] = start_edge;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
147 pt[1] = (pt[0] + 1) % 4;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
148
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
149 // Calculate contour segment start if first of contour.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
150 if (first)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
151 {
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
152 tmp = fabs (pz[pt[1]]) / fabs (pz[pt[0]]);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
153
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
154 if (xisnan (tmp))
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
155 ct_x = ct_y = 0.5;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
156 else
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
157 {
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
158 ct_x = px[pt[0]] + (px[pt[1]] - px[pt[0]])/(1 + tmp);
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
159 ct_y = py[pt[0]] + (py[pt[1]] - py[pt[0]])/(1 + tmp);
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
160 }
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
161
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
162 start_contour (lvl, ct_x, ct_y);
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
163 first = false;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
164 }
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
165
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
166 // Find stop edge.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
167 // FIXME -- perhaps this should use a while loop?
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
168 for (unsigned int k = 1; k <= 4; k++)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
169 {
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
170 if (start_edge == 0 || start_edge == 2)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
171 stop_edge = (start_edge + k) % 4;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
172 else
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
173 stop_edge = (start_edge - k) % 4;
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
174
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
175 if (static_cast<char> (1 << stop_edge) & id)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
176 break;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
177 }
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
178
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
179 pt[0] = stop_edge;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
180 pt[1] = (pt[0] + 1) % 4;
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
181 tmp = fabs (pz[pt[1]]) / fabs (pz[pt[0]]);
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
182
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
183 if (xisnan (tmp))
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
184 ct_x = ct_y = 0.5;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
185 else
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
186 {
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
187 ct_x = px[pt[0]] + (px[pt[1]] - px[pt[0]])/(1 + tmp);
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
188 ct_y = py[pt[0]] + (py[pt[1]] - py[pt[0]])/(1 + tmp);
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
189 }
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
190
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
191 // Add point to contour.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
192 add_point (ct_x, ct_y);
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
193
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
194 // Decrease id value of current facet for start edge.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
195 mark(r, c) -= static_cast<char> (1 << stop_edge);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
196
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
197 // Find next facet.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
198 if (stop_edge == 0)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
199 r--;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
200 else if (stop_edge == 1)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
201 c++;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
202 else if (stop_edge == 2)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
203 r++;
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
204 else if (stop_edge == 3)
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
205 c--;
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
206
15867
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
207 // Go to next facet.
704e15f8fecd Modify contourc recursion to a loop to avoid stack overflow (bug #37891)
Mike Miller <mtmiller@ieee.org>
parents: 15220
diff changeset
208 start_edge = (stop_edge + 2) % 4;
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
209
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
210 }
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
211 }
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
212
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
213 static void
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
214 mark_facets (const Matrix& Z, charMatrix& mark, double lvl)
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
215 {
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
216 unsigned int nr = mark.rows ();
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
217 unsigned int nc = mark.cols ();
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
218
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
219 double f[4];
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
220
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
221 for (unsigned int c = 0; c < nc; c++)
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
222 for (unsigned int r = 0; r < nr; r++)
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
223 {
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
224 f[0] = Z(r, c) - lvl;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
225 f[1] = Z(r, c+1) - lvl;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
226 f[3] = Z(r+1, c) - lvl;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
227 f[2] = Z(r+1, c+1) - lvl;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
228
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
229 for (unsigned int i = 0; i < 4; i++)
15220
61822c866ba1 use std::numeric_limits<T>::epsilon in C++ code
John W. Eaton <jwe@octave.org>
parents: 15195
diff changeset
230 if (fabs(f[i]) < std::numeric_limits<double>::epsilon ())
61822c866ba1 use std::numeric_limits<T>::epsilon in C++ code
John W. Eaton <jwe@octave.org>
parents: 15195
diff changeset
231 f[i] = std::numeric_limits<double>::epsilon ();
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
232
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
233 if (f[1] * f[2] < 0)
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
234 mark(r, c) += 2;
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
235
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
236 if (f[0] * f[3] < 0)
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
237 mark(r, c) += 8;
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
238 }
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
239
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
240 for (unsigned int r = 0; r < nr; r++)
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
241 for (unsigned int c = 0; c < nc; c++)
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
242 {
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
243 f[0] = Z(r, c) - lvl;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
244 f[1] = Z(r, c+1) - lvl;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
245 f[3] = Z(r+1, c) - lvl;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
246 f[2] = Z(r+1, c+1) - lvl;
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
247
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
248 for (unsigned int i = 0; i < 4; i++)
15220
61822c866ba1 use std::numeric_limits<T>::epsilon in C++ code
John W. Eaton <jwe@octave.org>
parents: 15195
diff changeset
249 if (fabs(f[i]) < std::numeric_limits<double>::epsilon ())
61822c866ba1 use std::numeric_limits<T>::epsilon in C++ code
John W. Eaton <jwe@octave.org>
parents: 15195
diff changeset
250 f[i] = std::numeric_limits<double>::epsilon ();
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
251
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
252 if (f[0] * f[1] < 0)
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
253 mark(r, c) += 1;
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
254
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
255 if (f[2] * f[3] < 0)
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
256 mark(r, c) += 4;
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
257 }
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
258 }
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
259
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
260 static void
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
261 cntr (const RowVector& X, const RowVector& Y, const Matrix& Z, double lvl)
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
262 {
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
263 unsigned int nr = Z.rows ();
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
264 unsigned int nc = Z.cols ();
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
265
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
266 charMatrix mark (nr - 1, nc - 1, 0);
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
267
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
268 mark_facets (Z, mark, lvl);
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
269
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
270 // Find contours that start at a domain edge.
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
271
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
272 for (unsigned int c = 0; c < nc - 1; c++)
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
273 {
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
274 // Top.
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
275 if (mark(0, c) & 1)
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
276 drawcn (X, Y, Z, lvl, 0, c, 0.0, 0.0, 0, true, mark);
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
277
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
278 // Bottom.
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
279 if (mark(nr - 2, c) & 4)
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
280 drawcn (X, Y, Z, lvl, nr - 2, c, 0.0, 0.0, 2, true, mark);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
281 }
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
282
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
283 for (unsigned int r = 0; r < nr - 1; r++)
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
284 {
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
285 // Left.
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
286 if (mark(r, 0) & 8)
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
287 drawcn (X, Y, Z, lvl, r, 0, 0.0, 0.0, 3, true, mark);
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
288
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
289 // Right.
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
290 if (mark(r, nc - 2) & 2)
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
291 drawcn (X, Y, Z, lvl, r, nc - 2, 0.0, 0.0, 1, true, mark);
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
292 }
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
293
7070
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
294 for (unsigned int r = 0; r < nr - 1; r++)
7593f8e83a2e [project @ 2007-10-25 20:41:16 by jwe]
jwe
parents: 7042
diff changeset
295 for (unsigned int c = 0; c < nc - 1; c++)
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
296 if (mark (r, c) > 0)
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
297 drawcn (X, Y, Z, lvl, r, c, 0.0, 0.0, 255, true, mark);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
298 }
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
299
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14501
diff changeset
300 DEFUN (__contourc__, args, ,
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
301 "-*- texinfo -*-\n\
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14501
diff changeset
302 @deftypefn {Built-in Function} {} __contourc__ (@var{x}, @var{y}, @var{z}, @var{levels})\n\
6945
6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe]
jwe
parents: 6257
diff changeset
303 Undocumented internal function.\n\
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
304 @end deftypefn")
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
305 {
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
306 octave_value retval;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
307
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
308 if (args.length () == 4)
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
309 {
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
310 RowVector X = args (0).row_vector_value ();
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
311 RowVector Y = args (1).row_vector_value ();
7042
e54cc03d53f6 [project @ 2007-10-19 20:43:32 by jwe]
jwe
parents: 7017
diff changeset
312 Matrix Z = args (2).matrix_value ();
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
313 RowVector L = args (3).row_vector_value ();
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
314
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
315 if (! error_state)
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
316 {
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
317 contourc.resize (2, 0);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
318
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
319 for (int i = 0; i < L.length (); i++)
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
320 cntr (X, Y, Z, L (i));
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
321
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
322 end_contour ();
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
323
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
324 retval = contourc;
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
325 }
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
326 else
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
327 error ("__contourc__: invalid argument values");
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
328 }
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
329 else
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
330 print_usage ();
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
331
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
332 return retval;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents:
diff changeset
333 }
12805
3641167e5b75 codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents: 11586
diff changeset
334
3641167e5b75 codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents: 11586
diff changeset
335 /*
3641167e5b75 codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents: 11586
diff changeset
336 ## 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
337 %!assert (1)
3641167e5b75 codesprint: *.cc helper functions do not need tests
Rik <octave@nomad.inbox5.com>
parents: 11586
diff changeset
338 */