Function rbits

src/lib/rng/utils.ts:25–37

(bits: number, unif_rand: IRNG)

}

function rbits(bits: number, unif_rand: IRNG) {
    let v = 0n;

    for (let n = 0; n <= bits; n += 16) {
        // FIXME:  should this not be "*65535" ?, its ok, 65534 would be an error!!! (as some part of the range would not be reachable)
        //    but 65536 not matter so much
        const ran = unif_rand.random();
        const v1 = new Int32Array([trunc(ran * 65536)]);
        v = 65536n * v + BigInt(v1[0]); //
    }
    // mask out the bits in the result that are not needed
    return v & ((1n << BigInt(bits)) - 1n);
}

function R_unif_index_0(dn: number, unif_rand: IRNG): number {
    const cut = unif_rand.cut;