Function gcd [src]

rma may alias x or y. x and y may alias each other. Asserts that rma has enough limbs to store the result. Upper bound is @min(x.limbs.len, y.limbs.len). limbs_buffer is used for temporary storage during the operation. When this function returns, it will have the same length as it had when the function was called.

Prototype

pub fn gcd(rma: *Mutable, x: Const, y: Const, limbs_buffer: *std.ArrayList(Limb)) !void

Parameters

rma: *Mutablex: Consty: Constlimbs_buffer: *std.ArrayList(Limb)

Source

  pub fn gcd(rma: *Mutable, x: Const, y: Const, limbs_buffer: *std.ArrayList(Limb)) !void {
    const prev_len = limbs_buffer.items.len;
    defer limbs_buffer.shrinkRetainingCapacity(prev_len);
    const x_copy = if (rma.limbs.ptr == x.limbs.ptr) blk: {
        const start = limbs_buffer.items.len;
        try limbs_buffer.appendSlice(x.limbs);
        break :blk x.toMutable(limbs_buffer.items[start..]).toConst();
    } else x;
    const y_copy = if (rma.limbs.ptr == y.limbs.ptr) blk: {
        const start = limbs_buffer.items.len;
        try limbs_buffer.appendSlice(y.limbs);
        break :blk y.toMutable(limbs_buffer.items[start..]).toConst();
    } else y;

    return gcdLehmer(rma, x_copy, y_copy, limbs_buffer);
}  