$result = matrixwise_map(\&f, $mats);

    print dumpf($result), $/;
    # [[f([1,5]),f([2,6])],[f([3,7]),f([4,8])]]