sub transpose_gmch {
    my ($matrix) = shift ;

    my @result = () ;
    my $lc = $#{$matrix->[0]} ;

    for my $col (0..$lc) {
      push @result, [map $_->[$col], @$matrix] ;
    } ;

    return( \@result );
}