use strict;
use warnings;
use PDL;

my $x = pdl([1.0, 0.0, 3.0, 8.0],
            [2.0, 0.0, 1.0, 1.0],
            [1.5, 2.5, 1.0, 0.0]);
my $n = 4;

my $mask = ( $x == $x-> rint )-> andover
            &
           ( $x-> sumover == $n );

print $x-> transpose
        -> whereND( $mask )
        -> transpose;

##</code><code>##

print $z-> info( '%M' );

##</code><code>##

my $a = ones(3);
my $b = 2 * ones(3);

$b += $a;
$a /= -2;
$b *= $a;

print $b;

##</code><code>##

my $z = random( 10 );
$z( $z-> maximum_ind ) .= 0;
print $z;

##</code><code>##

my $z = zeroes( 10 );
my $i = pdl( 1, 3, 5, 3, 1 );
indadd( 1, $i, $z );
print $z;

##</code><code>##

my $z = 1 + sequence 14;
my $len = 4;

print $z-> lags( 0, 1, 1 + $z-> nelem - $len )
        -> slice( '','-1:0' );

##</code><code>##

my $x = sequence 16, 16;

print $x-> lags( 1, 4, 4 )
        -> slice( '', '', '-1:0' )
        -> xchg( 0, 1 )
        -> sumover
        -> lags( 0, 4, 4 )
        -> slice( '', '-1:0' )
        -> sumover;

##</code><code>##

my $x = sequence 16, 16;

print $x-> reshape( 4, 4, 4, 4 )
        -> reorder( 0, 2, 1, 3 )
        -> clump( 2 )
        -> sumover