1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
|
## Copyright (C) 1998-2003 Joao Cardoso.
## Copyright (C) 2004 Rafael Laboissiere
##
## This program is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by the
## Free Software Foundation; either version 2 of the License, or (at your
## option) any later version.
##
## This program is distributed in the hope that it will be useful, but
## WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
## General Public License for more details.
##
## This file is part of plplot_octave.
function __pl_contour(x, y, z, n)
global __pl
strm = plgstrm+1;
unwind_protect
old_empty_list_elements_ok = warning("query","Octave:empty-list-elements");
warning("off","Octave:empty-list-elements");
grid = 0;
if (__pl.grid(strm))
grid = 2;
endif
if (ishold== 0)
xm = min(min(x)); xM = max(max(x));
ym = min(min(y)); yM = max(max(y));
zm = min(min(z)); zM = max(max(z));
if (__pl.axis_st(strm))
## at least x always exist
xm = __pl.axis(strm,1); xM = __pl.axis(strm,2);
z = z(find(x >= xm & x <= xM), :);
if (length(__pl.axis) >= 4)
ym = __pl.axis(strm,3); yM = __pl.axis(strm,4);
z = z(:, find(y >= ym & y <= yM));
else
__pl.axis(strm,3) = ym; __pl.axis(strm,4) = yM;
endif
if (length(__pl.axis) == 6)
zm = __pl.axis(strm,5); zM = __pl.axis(strm,6);
else
__pl.axis(strm,5) = zm; __pl.axis(strm,6) = zM;
endif
else # make axis() return current axis
__pl.axis(strm,1) = xm; __pl.axis(strm,2) = xM;
__pl.axis(strm,3) = ym; __pl.axis(strm,4) = yM;
__pl.axis(strm,5) = zm; __pl.axis(strm,6) = zM;
endif
__pl.type(strm) = 0; ## was -3 ??
__pl.plcol(strm) = 1;
__pl.pllsty(strm) = 1;
__pl.lab_pos(strm) = 1;
__pl.lab_str = "";
plcol0(15); pllsty(1);
__pl_plenv(xm, xM, ym, yM, 0, grid);
else
if (columns(__pl.axis(strm,:)) != 6)
error("You must contour/shade plot something before entering hold mode");
endif
xm = __pl.axis(strm,1); xM = __pl.axis(strm,2);
ym = __pl.axis(strm,3); yM = __pl.axis(strm,4);
zm = __pl.axis(strm,5); zM = __pl.axis(strm,6);
z = z(find(x >= xm & x <= xM), find(y >= ym & y <= yM));
endif
[xlen, ylen] = size(z);
tr = [(xM-xm)/(xlen-1); 0; xm; 0; (yM-ym)/(ylen-1); ym];
if (!isscalar(n))
st = (zM-zm)/length(n)/2;
clevel = n; n = length(n);
else
st = (zM-zm)/n/2;
clevel = (linspace(zm+st, zM-st, n))';
endif
for i=1:n
plcol0(__pl.plcol(strm)); pllsty(__pl.pllsty(strm));
plcont(z, 1, xlen, 1, ylen, clevel(i), tr);
__pl.lab_str = [__pl.lab_str; sprintf("%#+.2G", clevel(i))];
__pl.lab_col(strm,__pl.lab_pos(strm)) = __pl.plcol(strm);
__pl.lab_lsty(strm,__pl.lab_pos(strm)) = __pl.pllsty(strm);
__pl.lab_pos(strm) = __pl.lab_pos(strm) + 1;
__pl.plcol(strm) = rem(__pl.plcol(strm), 15)+1;
if (__pl.line_style(strm))
__pl.pllsty(strm) = rem(__pl.pllsty(strm), 8)+1;
endif
endfor
## to not mess up the legend, print it only after ploting all contour lines
if (__pl.legend(strm))
__pl_draw_legend
endif
plcol0(15);
pllab(tdeblank(__pl.xlabel(strm,:)),
tdeblank(__pl.ylabel(strm,:)),
tdeblank(__pl.tlabel(strm,:)));
plflush;
__pl.items(strm) = 1; # for now!
unwind_protect_cleanup
warning(old_empty_list_elements_ok.state,"Octave:empty-list-elements");
end_unwind_protect
endfunction
|
