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
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
|
#!/usr/bin/perl
use lib 'lib';
use PDL;
use PDL::NiceSlice;
use PDL::Graphics::Gnuplot qw(plot plot3d gpwin);
@windows = ();
use feature qw(say);
# data I use for 2D testing
my $x = sequence(21) - 10;
# data I use for 3D testing
my($th,$ph,$x_3d,$y_3d,$z_3d);
$th = zeros(30)->xlinvals( 0, 3.14159*2)->transpose;
$ph = zeros(30)->xlinvals( -3.14159/2, 3.14159/2);
$x_3d = (cos($ph)*cos($th))->flat->glue(0, (cos($ph)*cos($th))->transpose->flat);
$y_3d = (cos($ph)*sin($th))->flat->glue(0, (cos($ph)*sin($th))->transpose->flat);
$z_3d = (sin($ph)*$th->ones)->flat->glue(0,(sin($ph)*$th->ones)->transpose->flat);
$rv = rvals(21,21);
$zv = cos($rv) / (3 + $rv);
#################################
# Now the tests!
#################################
# first, some very basic stuff. Testing implicit domains, multiple curves in
# arguments, packed in ndarrays, etc
sub prompt {
print $_[0]. " (Press <RETURN> to continue)";
$a=<>;
}
$w=gpwin(x11);
$w->plot($x**2);
prompt("A simple parabola");
$w->plot(-$x, $x**3);
prompt("A cubic");
$w->plot(-$x, $x**3,{},
$x, $x**2);
prompt("Parabola and cubic");
$w->plot(PDL::cat($x**2, $x**3));
prompt("Another way to plot a parabola and cubic");
$w->plot(-$x,
PDL::cat($x**2, $x**3));
prompt("Yet another way");
push(@windows,$w);
# some more varied plotting, using the object-oriented interface
$w->options(globalwith => 'linespoints', xmin => -10);
$w->plot( title => 'Error bars and other things',
y2tics=>10,
{with => 'lines', 'lw' => 4},
legend => ['Parabola A','Parabola B','Parabola C'],
axes => 'x1y2',
PDL::cat($x, $x*2, $x*3), $x**2 - 300,
with => 'xyerrorbars',
axes=>'x1y1',
$x**2 * 10, $x**2/40, $x**2/2, # implicit domain
{with => 'line', legend => 'cubic', tuplesize => 2},
{legend => ['shifted cubic A','shifted cubic B']},
$x, PDL::cat($x**3, $x**3 - 100) );
prompt("Error bars and other things");
# a way to control the point size
$w->plot(cbmin => -600, cbmax => 600, title=>"Variable pointsize",
{with => qw(points pointtype 7 pointsize variable)},
$x, $x/2, (10-abs($x))/2);
prompt("Variable pointsize");
## filled curve
$w->plot(
{ with => 'filledcurves', fc => 'red', below => 'y=0' },
-$x, (-$x)**3,
{ with => 'filledcurves', fillcolor => 'green', above => 'y=0' },
$x, ($x)**3
);
prompt("Area curve with filled colors");
################################
# some 3d stuff
################################
$w->reset;
# plot a gridded surface
$w->plot3d(
with=>'linespoints',
xvals($zv), yvals($zv),
$zv
);
prompt("A gridded surface (feed in a 2-D column)\n (Incidentally, if you're using the X11 or wxt device you can\n change the view by dragging your mouse around the window)\n");
# Plot a collection of lines
$w->plot3d(
with=>'linespoints',
cdim=>1,
xvals($zv), yvals($zv),
$zv
);
prompt("A collection of lines (same data, treated as a collection of 1-D columns)");
# plot a sphere
$w->plot3d( j=>1, zrange=>[-1,1],yrange=>[-1,1],xrange=>[-1,1],
{with=>'linespoints'},
$x_3d, $y_3d, $z_3d,
);
prompt("A sphere");
# sphere, ellipse together
$w->plot3d(
j=>1,
{legend => ['sphere', 'ellipsoid'],cd=>1},
$x_3d->cat($x_3d*2),
$y_3d->cat($y_3d*2), $z_3d );
prompt("A sphere and an ellipsoid");
# some paraboloids plotted on an implicit 2D domain
{
my($xy,$z,$xy_half,$z_half);
$xy = zeros(21,21)->ndcoords - pdl(10,10);
$z = inner($xy, $xy);
$xy_half = zeros(11,11)->ndcoords;
$z_half = inner($xy_half, $xy_half);
$w->plot3d(
globalwith => 'points', title => 'gridded paraboloids',tuplesize=>1,
{legend => ['zplus','zminus']}, $z->cat(-$z),
{legend => 'zplus2'}, $z*2);
prompt("Some paraboloids");
$w->reset;
}
# 3d, variable color, variable pointsize
{
my($pi,$theta,$z);
$pi = 3.14159;
$theta = zeros(200)->xlinvals(0, 6*$pi);
$z = zeros(200)->xlinvals(0, 5);
$w->plot3d(
title => 'double helix',
{ with => 'linespoints', pointsize => 'variable', pointtype => 7, palette => 1, tuplesize => 5,
legend => ['spiral 1','spiral 2'],cdim=>1},
# 2 sets of x, y, z:
cos($theta)->cat(-cos($theta)),
sin($theta)->cat(-sin($theta)),
$z,
# pointsize, color
0.5 + abs(cos($theta)), sin(2*$theta) );
prompt("A double helix");
}
# implicit domain heat map
{
my $xy;
$xy = zeros(21,21)->ndcoords - pdl(10,10);
$w->plot3d(
title => 'Paraboloid heat map',
extracmds => 'set view 0,0',
zrange => [-1,1],
with => 'image', inner($xy, $xy));
prompt("An image");
}
say STDERR 'I should complain about an invalid "with":';
say STDERR "=================================";
eval( <<'EOM' );
plot(with => 'bogusstyle', $x);
EOM
print STDERR $@ if $@;
say STDERR "=================================\n\n";
say STDERR 'PDL::Graphics::Gnuplot can detect I/O hangs. Here I ask for a delay, so I should detect this and quit after a few seconds:';
say STDERR "=================================";
eval( <<'EOM' );
my $xy = zeros(21,21)->ndcoords - pdl(10,10);
plot( extracmds => 'pause 10',
sequence(5));
EOM
print STDERR $@ if $@;
say STDERR "=================================\n\n";
say STDERR "This concludes the gnuplot demo! Enjoy using PDL::Graphics::Gnuplot!\n";
|
