use strict;
use warinings;
use PDL;
use PDL::NiceSlice;

sub submatrix_sum {
    my $m = pdl @_;

    $m->range(ndcoords($m(1:,1:)), 2)->reorder(2, 3, 0, 1)
        ->clump(2)->sumover;
}