raku_rmultinom

sub rmultinom(int32, CArray[num64], int32, CArray[int32])
    is native( RMATH ) { * };

sub raku_rmultinom($n , $size, @prob is copy) is export {
    @prob.append(0 xx $size - @prob);
    @prob = @prob.map: {$_.Num};
    my $prob = CArray[num64].new(@prob);
    my $ints = CArray[int32].allocate($size);
    rmultinom($n, $prob, $size, $ints);
    return $ints.list
}
my @prob = [0.1,0.2,1.3,0.2,0.3];
set_seed(123,456);
say raku_rmultinom(100,12, @prob);


# Outputs
# (0 0 -1 0 0 0 0 0 0 0 0 0)