/usr/lib/perl5/PDL/Demos/General.pm is in pdl 1:2.007-2build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | # Copyright (C) 1998 Tuomas J. Lukka.
# All rights reserved, except redistribution
# with PDL under the PDL License permitted.
package PDL::Demos::General;
use PDL;
PDL::Demos::Routines->import();
sub comment($);
sub act($);
sub output;
sub run {
comment q|
Welcome to a short tour of PDL's capabilities.
This tour shows some of the main selling points
of PDL. However, because we want this script to
run everywhere, some modules which require external
modules for use are explicitly excluded, namely
- PDL::Graphics::TriD (3D Graphics) [*]
- PDL::Graphics::PGPLOT (PGPLOT graphics)
- PDL::IO::FlexRaw (flexible raw input/output)
[*]: this module has its separate demos in a subdirectory.
Note that your own scripts must start with
use PDL;
to work properly, so that you can simply say
perl script.pl
or you can just try some of the commands illustrated
in the demos by just retyping them at the perldl or pdl
'pdl>' command prompt.
|;
act q|
$a = zeroes 5,5; # 5x5 matrix
output $a;
|;
act q|
# Now, don't think that the number of dimensions is limited
# to two:
$m = zeroes(3,2,2); # 3x2x2 cube
output $m;
|;
act q|
$a ++; # Operators like increment work..
output $a;
|;
act q|
# xvals and yvals (yes, there is also zvals...)
# give you piddles which give the coordinate value.
$b = xvals $a;
output $b;
|;
act q|
# So you can do things like
$b = $a + 0.1 * xvals($a) + 0.01 * yvals($a);
output $b;
|;
act q|
# Arithmetic operations work:
$x = xvals(10) / 5;
output $x,"\n";
output ((sin $x),"\n");
|;
act q|
# You can also take slices:
output $b;
output $b->slice(":,2:3"); # rows 2 and 3
|;
act q|
output $b->slice("2:3,:"); # or columns 2 and 3
|;
act q|
output $b;
output $b->diagonal(0,1),"\n"; # 0 and 1 are the dimensions
|;
act q|
# One of the really nifty features is that the
# slices are actually references back to the original
# piddle:
$diag = $b->diagonal(0,1);
output $b;
output $diag,"\n";
$diag+=100;
output "AFTER:\n";
output $diag,"\n";
output "Now, guess what \$b looks like?\n";
|;
act q|
# Yes, it has changed:
output $b;
|;
act q|
# Another example (we only modify elements 0,2 and 4 of
# each row):
$t = $b->slice("0:4:2"); $t += 50;
output $b;
|;
act q|
# There are lots of useful functions in e.g. PDL::Primitive
# and PDL::Slices - we can't show you all but here are some
# examples:
output $b;
output $b->sum, "\n";
output $b->sumover,"\n"; # Only over first dim.
|;
act q|
output $b->xchg(0,1);
output $b->minimum,"\n"; # over first dim.
output $b->min,"\n";
|;
act q|
output $b->random;
|;
act q|
# Here are some more advanced tricks for selecting
# parts of 1-D vectors:
$a = (xvals 12)/3;
$i = which(sin($a) > 0.5); # Indices of those sines > 0.5
output $a,"\n";
output $i,"\n";
output $a->index($i),"\n";
# and we can have the effect of the last command in one
# go using 'where' instead of 'which' and 'index' as in
output $a->where(sin($a) > 0.5),"\n";
# and finally take the sin of these elements
# (to show that these are indeed the correct ones)
output sin($a->index($i)),"\n";
|;
comment q|
We hope you enjoyed these demos illustrating some
of the basic capabilities of PDL.
We encourage you to play with these commands in
the perldl or pdl2 shell and use its online help support
to find out more about these and other commands and
features of PDL.
Just type 'help' to get started.
|;
}
1;
|