From 22356e5e0cf21fdad0029c3afb99a648a3b94f83 Mon Sep 17 00:00:00 2001 From: "Nicolas \"Pixel\" Noble" Date: Mon, 7 Jun 2010 03:57:48 +0200 Subject: Initial commit of the gmp lua plugin. --- src/lgmp.lua | 1639 ++++++++++++++++++++++++++++++++ src/plugin-luagmp.cc | 2532 ++++++++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 4171 insertions(+) create mode 100644 src/lgmp.lua create mode 100644 src/plugin-luagmp.cc (limited to 'src') diff --git a/src/lgmp.lua b/src/lgmp.lua new file mode 100644 index 0000000..2a09639 --- /dev/null +++ b/src/lgmp.lua @@ -0,0 +1,1639 @@ + +-- +-- This file is part of the lgmp package for Lua 5.1. +-- +-- Author: Wim Couwenberg +-- Date : 2007/7/28 +-- +-- The lgmp package is distributed under the MIT license. See the "COPYRIGHT" +-- file that came with the distribution of lgmp for license details. +-- + +-- +-- Adapted for lua-interface by Nicolas "Pixel" Noble +-- + +local type = type +local getmetatable = getmetatable +local assert = assert +local error = error +local match = string.match + +local prv +local aux +local randmeta = {} +local zmeta = {} +local fmeta = {} + +randmeta.__index = randmeta +zmeta.__index = zmeta +fmeta.__index = fmeta + +local function isrand(obj) + return type(obj) == "userdata" and getmetatable(obj) == randmeta +end + +local function isz(obj) + return type(obj) == "userdata" and getmetatable(obj) == zmeta +end + +local function isf(obj) + return type(obj) == "userdata" and getmetatable(obj) == fmeta +end + +local function ntype(obj) + local t = type(obj) + if t == "number" then + if obj%1 == 0 then + if obj >= 0 and obj <= prv.ULONG_MAX then + return "u" + elseif obj >= prv.LONG_MIN and obj < 0 then + return "s" + elseif obj >= -prv.ULONG_MAX and obj < prv.LONG_MIN then + return "n" + end + end + return "d" + elseif t == "userdata" then + local m = getmetatable(obj) + if m == zmeta then + return "z" + elseif m == fmeta then + return "f" + end + end + return "?" +end + +local function checkrand(obj) + assert(isrand(obj), "gmp random state expected") +end + +local function checku(obj) + assert(type(obj) == "number" and obj >= 0 and obj <= prv.ULONG_MAX and obj%1 == 0, "unsigned integer expected") +end + +local function checkn(obj) + assert(type(obj) == "number" and obj >= 0 and obj <= prv.LONG_MAX and obj%1 == 0, "non-negative intger expected") +end + +local function checkbase(obj) + assert(type(obj) == "number" and obj >= 2 and obj <= 62 and obj%1 == 0, "gmp number base must be an integer between 2 and 62") +end + +local function checkz(obj) + assert(isz(obj), "gmp integer expected") +end + +local function checkzopt(obj) + assert(obj == nil or isz(obj), "gmp integer expected") +end + +local function checkf(obj) + assert(isf(obj), "gmp floating point expected") +end + +local function checkfopt(obj) + assert(obj == nil or isf(obj), "gmp floating point expected") +end + +local dtoz = prv.mpz_init_set_d +local dtof = prv.mpf_init_set_d + +function z(value, base) + if value == nil then + return prv.mpz_init() + elseif type(value) == "number" then + return dtoz(value) + elseif type(value) == "string" then + if base == nil then + base = 0 + elseif base ~= 0 then + checkbase(base) + end + local res = prv.mpz_init_set_str(value, base) + if res then return res end + error("not a valid integer constant in base " .. base .. ": " .. value) + elseif isz(value) then + return prv.mpz_init_set(value) + elseif isf(value) then + local res = prv.mpz_init() + prv.mpz_set_f(res, value) + return res + else + error("cannot initialize gmp integer from " .. type(value)) + end +end + +function f(value, base) + if value == nil then + return prv.mpf_init() + elseif type(value) == "number" then + return dtof(value) + elseif type(value) == "string" then + if base == nil then + base = 10 + elseif type(base) == "number" and base < 0 then + checkbase(-base) + else + checkbase(base) + end + local res = prv.mpf_init_set_str(value, base) + if res then return res end + error("not a valid floating point constant in base " .. base .. ": " .. value) + elseif isf(value) then + return prv.mpf_init_set(value) + elseif isz(value) then + local res = prv.mpf_init() + prv.mpf_set_z(res, value) + return res + else + error("cannot initialize gmp floating point from " .. type(value)) + end +end + +function zmeta:__tostring() + checkz(self) + return prv.mpz_get_str(10, self) +end + +function zmeta:__concat(other) + if isz(self) then + return zmeta.get_str(self) .. other + else + return self .. zmeta.get_str(other) + end +end + +function zmeta:__add(other) + if isz(self) then + return zmeta.add(self, other) + else + return zmeta.add(other, self) + end +end + +function zmeta:__sub(other) + if isz(self) then + return zmeta.sub(self, other) + else + return zmeta.rsub(other, self) + end +end + +function zmeta:__mul(other) + if isz(self) then + return zmeta.mul(self, other) + else + return zmeta.mul(other, self) + end +end + +function zmeta:__div(other) + if isz(self) then + return zmeta.fdiv_q(self, other) + elseif type(other) == "number" then + return zmeta.fdiv_q(dtoz(self), other) + else + error("unsupported type") + end +end + +function zmeta:__mod(other) + if isz(self) then + return zmeta.fdiv_r(self, other) + elseif type(self) == "number" then + return zmeta.fdiv_r(dtoz(self), other) + else + error("unsupported type") + end +end + +function zmeta:__pow(exp) + checkz(self) + checku(exp) + return prv.mpz_pow_ui(self, exp) +end + +function zmeta:__unm() + checkz(self) + return prv.mpz_neg(self) +end + +function zmeta:__lt(other) + checkz(self) + checkz(other) + return prv.mpz_cmp(self, other) < 0 +end + +function zmeta:__le(other) + checkz(self) + checkz(other) + return prv.mpz_cmp(self, other) <= 0 +end + +function zmeta:__eq(other) + checkz(self) + checkz(other) + return prv.mpz_cmp(self, other) == 0 +end + +function zmeta:__gc() + checkz(self) + return prv.mpz_clear(self) +end + +function zmeta:abs(res) + checkz(self) + checkzopt(res) + return prv.mpz_abs(self, res) +end + +function zmeta:add(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_add(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_add_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpz_sub_ui(self, -a, res) + elseif t == "z" then + return prv.mpz_add(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:addmul(a1, a2) + checkz(self) + checkz(a1) + local t = ntype(a2) + if t == "d" then + prv.mpz_addmul(self, a1, dtoz(a2)) + elseif t == "u" then + prv.mpz_addmul_ui(self, a1, a2) + elseif t == "s" or t == "n" then + prv.mpz_submul_ui(self, a1, -a2) + elseif t == "z" then + prv.mpz_addmul(self, a1, a2) + else + error("unsupported type") + end +end + +function zmeta:And(a, res) + checkz(self) + checkzopt(res) + if type(a) == "number" then + return prv.mpz_and(self, dtoz(a), res) + elseif isz(a) then + return prv.mpz_and(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:bin(a, res) + checkz(self) + checku(a) + checkzopt(res) + + return prv.mpz_bin_ui(self, a, res) +end + +function bin(a1, a2, res) + checku(a2) + checkzopt(res) + local t = ntype(a1) + if t == "d" then + return prv.mpz_bin_ui(dtoz(a1), a2, res) + elseif t == "u" then + return prv.mpz_bin_uiui(a1, a2, res) + elseif t == "s" or t == "n" then + if a2 <= prv.ULONG_MAX + a1 then + res = prv.mpz_bin_uiui(a2 - a1 - 1, a2, res) + if a2%2 ~= 0 then + prv.mpz_neg(res, res) + end + return res + else + return prv.mpz_bin_ui(dtoz(a1), a2, res) + end + elseif t == "z" then + return prv.mpz_bin_ui(a1, a2, res) + else + error("unsupported type") + end +end + +function zmeta:cdiv_q(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_cdiv_q(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_cdiv_q_ui(self, a, res) + elseif t == "s" or t == "n" then + local r2 + res, r2 = prv.mpz_fdiv_q_ui(self, -a, res) + return prv.mpz_neg(res, res), r2 + elseif t == "z" then + return prv.mpz_cdiv_q(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:cdiv_q_2exp(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_cdiv_q_2exp(self, a, res) +end + +function zmeta:cdiv_qr(a, r1, r2) + checkz(self) + checkzopt(r1) + checkzopt(r2) + local t = ntype(a) + if t == "d" then + return prv.mpz_cdiv_qr(self, dtoz(a), r1, r2) + elseif t == "u" then + return prv.mpz_cdiv_qr_ui(self, a, r1, r2) + elseif t == "s" or t == "n" then + local r3 + r1, r2, r3 = prv.mpz_fdiv_qr_ui(self, -a, r1, r2) + return prv.mpz_neg(r1, r1), r2, r3 + elseif t == "z" then + return prv.mpz_cdiv_qr(self, a, r1, r2) + else + error("unsupported type") + end +end + +function zmeta:cdiv_r(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_cdiv_r(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_cdiv_r_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpz_fdiv_r_ui(self, -a, res) + elseif t == "z" then + return prv.mpz_cdiv_r(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:cdiv_r_2exp(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_cdiv_r_2exp(self, a, res) +end + +function zmeta:cdiv(a) + checkz(self) + local t = ntype(a) + if t == "u" then + return prv.mpz_cdiv_ui(self, a) + elseif t == "s" or t == "n" then + return prv.mpz_cdiv_ui(self, -a) + else + error("unsupported type") + end +end + +function zmeta:clrbit(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_clrbit(self, a, res) +end + +function zmeta:cmp(a) + checkz(self) + local t = ntype(a) + if t == "d" or t == "n" then + return prv.mpz_cmp_d(self, a) + elseif t == "u" then + return prv.mpz_cmp_ui(self, a) + elseif t == "s" then + return prv.mpz_cmp_si(self, a) + elseif t == "z" then + return prv.mpz_cmp(self, a) + else + error("unsupported type") + end +end + +function zmeta:cmpabs(a) + checkz(self) + local t = ntype(a) + if t == "d" then + return prv.mpz_cmpabs_d(self, a) + elseif t == "u" then + return prv.mpz_cmpabs_ui(self, a) + elseif t == "s" or t == "n" then + return prv.mpz_cmpabs_ui(self, -a) + elseif t == "z" then + return prv.mpz_cmpabs(self, a) + else + error("unsupported type") + end +end + +function zmeta:com(res) + checkz(self) + checkzopt(res) + return prv.mpz_com(self, res) +end + +function zmeta:combit(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_combit(self, a, res) +end + +function zmeta:congruent(a1, a2) + checkz(self) + local t1, t2 = ntype(a1), ntype(a2) + if t1 == "z" and t2 == "z" then + return prv.mpz_congruent_p(self, a1, a2) ~= 0 + elseif t1 == "u" and t2 == "u" then + return prv.mpz_congruent_ui_p(self, a1, a2) ~= 0 + else + error("unsupported type") + end +end + +function zmeta:congruent_2exp(a1, a2) + checkz(self) + checkz(a1) + checku(a2) + return prv.mpz_congruent_2exp_p(self, a1, a2) ~= 0 +end + +function zmeta:divexact(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_divexact(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_divexact_ui(self, a, res) + elseif t == "s" or t == "n" then + res = prv.mpz_divexact_ui(self, -a, res) + return prv.mpz_neg(res, res) + elseif t == "z" then + return prv.mpz_divexact(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:divisible(a) + checkz(self) + local t = ntype(a) + if t == "d" then + return prv.mpz_divisible_p(self, dtoz(a)) ~= 0 + elseif t == "u" then + return prv.mpz_divisible_ui_p(self, a) ~= 0 + elseif t == "s" or t == "n" then + return prv.mpz_divisible_ui_p(self, -a) ~= 0 + elseif t == "z" then + return prv.mpz_divisible_p(self, a) ~= 0 + else + error("unsupported type") + end +end + +function zmeta:divisible_2exp(a) + checkz(self) + checku(a) + return prv.mpz_divisible_2exp_p(self, a) ~= 0 +end + +function fac(a, res) + checku(a) + checkzopt(res) + return prv.mpz_fac_ui(a, res) +end + +function zmeta:fdiv_q(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_fdiv_q(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_fdiv_q_ui(self, a, res) + elseif t == "s" or t == "n" then + local r2 + res, r2 = prv.mpz_cdiv_q_ui(self, -a, res) + return prv.mpz_neg(res, res), r2 + elseif t == "z" then + return prv.mpz_fdiv_q(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:fdiv_q_2exp(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_fdiv_q_2exp(self, a, res) +end + +function zmeta:fdiv_qr(a, r1, r2) + checkz(self) + checkzopt(r1) + checkzopt(r2) + local t = ntype(a) + if t == "d" then + return prv.mpz_fdiv_qr(self, dtoz(a), r1, r2) + elseif t == "u" then + return prv.mpz_fdiv_qr_ui(self, a, r1, r2) + elseif t == "s" or t == "n" then + local r3 + r1, r2, r3 = prv.mpz_cdiv_qr_ui(self, -a, r1, r2) + return prv.mpz_neg(r1, r1), r2, r3 + elseif t == "z" then + return prv.mpz_fdiv_qr(self, a, r1, r2) + else + error("unsupported type") + end +end + +function zmeta:fdiv_r(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_fdiv_r(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_fdiv_r_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpz_cdiv_r_ui(self, -a, res) + elseif t == "z" then + return prv.mpz_fdiv_r(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:fdiv_r_2exp(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_fdiv_r_2exp(self, a, res) +end + +function zmeta:fdiv(a) + checkz(self) + local t = ntype(a) + if t == "u" then + return prv.mpz_fdiv_ui(self, a) + elseif t == "s" or t == "n" then + return prv.mpz_fdiv_ui(self, -a) + else + error("unsupported type") + end +end + +function fib(a, res) + checku(a) + checkzopt(res) + return prv.mpz_fib_ui(a, res) +end + +function fib2(a, r1, r2) + checku(a) + checkzopt(r1) + checkzopt(r2) + return prv.mpz_fib2_ui(a, r1, r2) +end + +function zmeta:fits_sint() + checkz(self) + return prv.mpz_fits_sint_p(self) ~= 0 +end + +function zmeta:fits_slong() + checkz(self) + return prv.mpz_fits_slong_p(self) ~= 0 +end + +function zmeta:fits_sshort() + checkz(self) + return prv.mpz_fits_sshort_p(self) ~= 0 +end + +function zmeta:fits_uint() + checkz(self) + return prv.mpz_fits_uint_p(self) ~= 0 +end + +function zmeta:fits_ulong() + checkz(self) + return prv.mpz_fits_ulong_p(self) ~= 0 +end + +function zmeta:fits_ushort() + checkz(self) + return prv.mpz_fits_ushort_p(self) ~= 0 +end + +function zmeta:gcd(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_gcd(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_gcd_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpz_gcd_ui(self, -a, res) + elseif t == "z" then + return prv.mpz_gcd(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:gcdext(a, r1, r2, r3) + checkz(self) + checkz(a) + checkzopt(r1) + checkzopt(r2) + checkzopt(r3) + return prv.mpz_gcdext(self, a, r1, r2, r3) +end + +function zmeta:get_d() + checkz(self) + return prv.mpz_get_d(self) +end + +function zmeta:get_d_2exp() + checkz(self) + return prv.mpz_get_d_2exp(self) +end + +function zmeta:get_str(base) + checkz(self) + if base == nil then base = 10 else checkbase(base) end + return prv.mpz_get_str(base, self) +end + +function zmeta:hamdist(a) + checkz(self) + checkz(a) + return prv.mpz_hamdist(self, a) +end + +function zmeta:invert(a, res) + checkz(self) + checkz(a) + checkzopt(res) + local r2 + res, r2 = prv.mpz_invert(self, a, res) + if r2 ~= 0 then + return res + end +end + +function zmeta:ior(a, res) + checkz(self) + checkzopt(res) + if type(a) == "number" then + return prv.mpz_ior(self, dtoz(a), res) + elseif isz(a) then + return prv.mpz_ior(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:kronecker(a) + checkz(self) + local t = ntype(a) + if t == "d" then + return prv.mpz_kronecker(self, dtoz(a)) + elseif t == "u" then + return prv.mpz_kronecker_ui(self, a) + elseif t == "s" then + return prv.mpz_kronecker_si(self, a) + elseif t == "n" then + return prv.mpz_kronecker(self, dtoz(a)) + elseif t == "z" then + return prv.mpz_kronecker(self, a) + else + error("unsupported type") + end +end + +function zmeta:rkronecker(a) + checkz(self) + local t = ntype(a) + if t == "d" then + return prv.mpz_kronecker(dtoz(a), self) + elseif t == "u" then + return prv.mpz_ui_kronecker(a, self) + elseif t == "s" then + return prv.mpz_si_kronecker(a, self) + elseif t == "n" then + return prv.mpz_kronecker(dtoz(a), self) + elseif t == "z" then + return prv.mpz_kronecker(a, self) + else + error("unsupported type") + end +end + +zmeta.jacobi = zmeta.kronecker +zmeta.rjacobi = zmeta.rkronecker + +zmeta.legendre = zmeta.kronecker +zmeta.rlegendre = zmeta.rkronecker + +function zmeta:lcm(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_lcm(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_lcm_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpz_lcm_ui(self, -a, res) + elseif t == "z" then + return prv.mpz_lcm(self, a, res) + else + error("unsupported type") + end +end + +function lucnum(a, res) + checku(a) + checkzopt(res) + return prv.mpz_lucnum_ui(a, res) +end + +function lucnum2(a, r1, r2) + checku(a) + checkzopt(r1) + checkzopt(r2) + return prv.mpz_lucnum2_ui(a, r1, r2) +end + +function zmeta:mod(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_mod(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_fdiv_r_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpz_cdiv_r_ui(self, -a, res) + elseif t == "z" then + return prv.mpz_mod(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:mul(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_mul(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_mul_ui(self, a, res) + elseif t == "s" then + return prv.mpz_mul_si(self, a, res) + elseif t == "n" then + res = prv.mpz_mul_ui(self, -a, res) + return prv.mpz_neg(res, res) + elseif t == "z" then + return prv.mpz_mul(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:mul_2exp(a, res) + checkz(self) + checku(a) + checkzopt(self) + return prv.mpz_mul_2exp(self, a, res) +end + +function zmeta:neg(res) + checkz(self) + checkzopt(res) + return prv.mpz_neg(self, res) +end + +function zmeta:nextprime(res) + checkz(self) + checkzopt(res) + return prv.mpz_nextprime(self, res) +end + +function zmeta:perfect_power() + checkz(self) + return prv.mpz_perfect_power_p(self) ~= 0 +end + +function zmeta:perfect_square() + checkz(self) + return prv.mpz_perfect_square_p(self) ~= 0 +end + +function zmeta:popcount() + checkz(self) + return prv.mpz_popcount(self) +end + +function zmeta:pow(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_pow_ui(self, a, res) +end + +function zmeta:powm(a1, a2, res) + checkz(self) + checkz(a2) + checkzopt(res) + local t = ntype(a1) + if t == "d" then + return prv.mpz_powm(self, dtoz(a1), a2, res) + elseif t == "u" then + return prv.mpz_powm_ui(self, a1, a2, res) + elseif t == "s" or t == "n" then + return prv.mpz_powm(self, dtoz(a1), a2, res) + elseif t == "z" then + return prv.mpz_powm(self, a1, a2, res) + else + error("unsupported type") + end +end + +function zmeta:probab_prime(a) + checkz(self) + if a == nil then a = 10 else ckechn(a) end + local res = prv.mpz_probab_prime_p(self, a) + return res ~= 0 and res +end + +function zmeta:remove(a, res) + checkz(self) + checkzopt(res) + if type(a) == "number" then + return prv.mpz_remove(self, dtoz(a), res) + elseif isz(a) then + return prv.mpz_rempve(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:root(a, res) + checkz(self) + checku(a) + checkzopt(res) + local r2 + res, r2 = prv.mpz_root(self, a, res) + return res, r2 ~= 0 +end + +function zmeta:rootrem(a, r1, r2) + checkz(self) + checku(a) + checkzopt(r1) + checkzopt(r2) + return prv.mpz_rootrem(self, a, r1, r2) +end + +function zmeta:scan0(a) + checkz(self) + if a == nil then a = 0 else checku(a) end + return prv.mpz_scan0(self, a) +end + +function zmeta:scan1(a) + checkz(self) + if a == nil then a = 0 else checku(a) end + return prv.mpz_scan1(self, a) +end + +function zmeta:set(a, base) + checkz(self) + if type(a) == "number" then + prv.mpz_set_d(self, a) + elseif type(a) == "string" then + if base == nil then + base = 0 + elseif base ~= 0 then + checkbase(base) + end + if prv.mpz_set_str(self, a, base) ~= 0 then + error("not a valid integer constant in base " .. base .. ": " .. a) + end + elseif isz(a) then + prv.mpz_set(self, a) + elseif isf(a) then + prv.mpz_set_f(self, a) + else + error("unsupported type") + end +end + +function zmeta:setbit(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_setbit(self, a, res) +end + +function zmeta:sgn() + checkz(self) + return prv.mpz_sgn(self) +end + +function zmeta:sizeinbase(a) + checkz(self) + if a == nil then a = 10 else checkbase(a) end + return prv.mpz_sizeinbase(self, a) +end + +function zmeta:sqrt(res) + checkz(self) + assert(prv.mpz_sgn(self) >= 0, "taking square of negative number") + checkzopt(res) + return prv.mpz_sqrt(self, res) +end + +function zmeta:sqrtrem(r1, r2) + checkz(self) + assert(prv.mpz_sgn(self) >= 0, "taking square of negative number") + checkzopt(r1) + checkzopt(r2) + return prv.mpz_sqrtrem(self, r1, r2) +end + +function zmeta:sub(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_sub(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_sub_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpz_add_ui(self, -a, res) + elseif t == "z" then + return prv.mpz_sub(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:rsub(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a1) + if t == "d" then + return prv.mpz_sub(dtoz(a), self, res) + elseif t == "u" then + return prv.mpz_ui_sub(a, self, res) + elseif t == "s" or t == "n" then + res = prv.mpz_add_ui(self, -a, res) + return prv.mpz_neg(res, res) + elseif t == "z" then + return prv.mpz_sub(a, self, res) + else + error("unsupported type") + end +end + +function zmeta:submul(a1, a2) + checkz(self) + checkz(a1) + local t = ntype(a2) + if t == "d" then + prv.mpz_submul(self, a1, dtoz(a2)) + elseif t == "u" then + prv.mpz_submul_ui(self, a1, a2) + elseif t == "s" or t == "n" then + prv.mpz_addmul_ui(self, a1, -a2) + elseif t == "z" then + prv.mpz_submul(self, a1, a2) + else + error("unsupported type") + end +end + +function zmeta:tdiv_q(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_tdiv_q(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_tdiv_q_ui(self, a, res) + elseif t == "s" or t == "n" then + local r2 + res, r2 = prv.mpz_tdiv_q_ui(self, -a, res) + return prv.mpz_neg(res, res), r2 + elseif t == "z" then + return prv.mpz_tdiv_q(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:tdiv_q_2exp(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.mpz_tdiv_q_2exp(self, a, res) +end + +function zmeta:tdiv_qr(a, r1, r2) + checkz(self) + checkzopt(r1) + checkzopt(r2) + local t = ntype(a) + if t == "d" then + return prv.mpz_tdiv_qr(self, dtoz(a), r1, r2) + elseif t == "u" then + return prv.mpz_tdiv_qr_ui(self, a, r1, r2) + elseif t == "s" or t == "n" then + local r3 + r1, r2, r3 = prv.mpz_tdiv_qr_ui(self, -a, r1, r2) + return prv.mpz_neg(r1, r1), r2, r3 + elseif t == "z" then + return prv.mpz_tdiv_qr(self, a, r1, r2) + else + error("unsupported type") + end +end + +function zmeta:tdiv_r(a, res) + checkz(self) + checkzopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpz_tdiv_r(self, dtoz(a), res) + elseif t == "u" then + return prv.mpz_tdiv_r_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpz_tdiv_r_ui(self, -a, res) + elseif t == "z" then + return prv.mpz_tdiv_r(self, a, res) + else + error("unsupported type") + end +end + +function zmeta:tdiv_r_2exp(a, res) + checkz(self) + checku(a) + checkzopt(res) + return prv.tdiv_r_2exp(self, a, res) +end + +function zmeta:tdiv(a) + checkz(self) + local t = ntype(a) + if t == "u" then + return prv.mpz_tdiv_ui(self, a) + elseif t == "s" or t == "n" then + return prv.mpz_tdiv_ui(self, -a) + else + error("unsupported type") + end +end + +function zmeta:tstbit(a) + checkz(self) + checku(a) + return prv.mpz_tstbit(self, a) +end + +function pow(a1, a2, res) + checku(a2) + checkzopt(res) + local t = ntype(a1) + if t == "d" then + return prv.mpz_pow_ui(dtoz(a1), a2, res) + elseif t == "u" then + return prv.mpz_ui_pow_ui(a1, a2, res) + elseif t == "s" or t == "n" then + res = prv.mpz_ui_pow_ui(-a1, a2, res) + if a2%2 ~= 0 then + prv.mpz_neg(res, res) + end + return res + elseif t == "z" then + return prv.mpz_pow_ui(a1, a2, res) + else + error("unsupported type") + end + return prv.mpz_ui_pow_ui(a1, a2, res) +end + +function zmeta:xor(a, res) + checkz(self) + checkzopt(res) + if type(a) == "number" then + return prv.mpz_xor(self, dtoz(a), res) + elseif isz(a) then + return prv.mpz_xor(self, a, res) + else + error("unsupported type") + end +end + +function fmeta:abs(res) + checkf(self) + checkfopt(res) + return prv.mpf_abs(self, res) +end + +function fmeta:add(a, res) + checkf(self) + checkfopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpf_add(self, dtof(a), res) + elseif t == "u" then + return prv.mpf_add_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpf_sub_ui(self, -a, res) + elseif t == "f" then + return prv.mpf_add(self, a, res) + else + error("unsupported type") + end +end + +function fmeta:ceil(res) + checkf(self) + checkfopt(res) + return prv.mpf_ceil(self, res) +end + +function fmeta:cmp(a) + checkf(self) + local t = ntype(a) + if t == "d" then + return prv.mpf_cmp_d(self, a) + elseif t == "u" then + return prv.mpf_cmp_ui(self, a) + elseif t == "s" then + return prv.mpf_cmp_si(self, a) + elseif t == "n" then + return prv.mpf_cmp(self, dtof(a)) + elseif t == "f" then + return prv.mpf_cmp(self, a) + else + error("unsupported type") + end +end + +function fmeta:div(a, res) + checkf(self) + checkfopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpf_div(self, dtof(a), res) + elseif t == "u" then + return prv.mpf_div_ui(self, a, res) + elseif t == "s" or t == "n" then + res = prv.mpf_div_ui(self, -a, res) + return prv.mpf_neg(res, res); + elseif t == "f" then + return prv.mpf_div(self, a, res) + else + error("unsupported type") + end +end + +function fmeta:eq(a1, a2) + checkf(self) + checkf(a1) + checku(a2) + return prv.mpf_eq(self, a1, a2) ~= 0 +end + +function fmeta:fits_sint() + checkf(self) + return prv.mpf_fits_sint_p(self) ~= 0 +end + +function fmeta:fits_sshort() + checkf(self) + return prv.mpf_fits_sshort_p(self) ~= 0 +end + +function fmeta:fits_uint() + checkf(self) + return prv.mpf_fits_uint_p(self) ~= 0 +end + +function fmeta:fits_ulong() + checkf(self) + return prv.mpf_fits_ulong_p(self) ~= 0 +end + +function fmeta:fits_ushort() + checkf(self) + return prv.mpf_fits_ushort_p(self) ~= 0 +end + +function fmeta:floor(res) + checkf(self) + checkfopt(res) + return prv.mpf_floor(self, res) +end + +function fmeta:get_d() + checkf(self) + return prv.mpf_get_d(self) +end + +function fmeta:get_d_2exp() + checkf(self) + return prv.mpf_get_d_2exp(self) +end + +get_default_prec = prv.mpf_get_default_prec + +function fmeta:get_prec() + checkf(self) + return prv.mpf_get_prec(self) +end + +function fmeta:get_str(base, size) + checkf(self) + if base == nil then base = 10 else checkbase(base) end + if size == nil then size = 0 else checkn(size) end + return prv.mpf_get_str(base, size, self) +end + +function fmeta:integer() + checkf(self) + return prv.mpf_integer_p(self) ~= 0 +end + +function fmeta:mul(a, res) + checkf(self) + checkfopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpf_mul(self, dtof(a), res) + elseif t == "u" then + return prv.mpf_mul_ui(self, a, res) + elseif t == "s" or t == "n" then + res = prv.mpf_mul_ui(self, -a, res) + return prv.mpf_neg(res, res) + elseif t == "f" then + return prv.mpf_mul(self, a, res) + else + error("unsupported type") + end +end + +function fmeta:mul_2exp(a, res) + checkf(self) + checku(a) + checkfopt(res) + return prv.mpf_mul_2exp(self, a, res) +end + +function fmeta:neg(res) + checkf(self) + checkfopt(res) + return prv.mpf_neg(self, res) +end + +function fmeta:pow(a, res) + checkf(self) + checku(a) + checkfopt(res) + return prv.mpf_pow_ui(self, a, res) +end + +function fmeta:reldiff(a, res) + checkf(self) + checkf(a) + checkfopt(res) + return prv.mpf_reldiff(self, a, res) +end + +function set_default_prec(a) + checku(a) + prv.mpf_set_default_prec(a) +end + +function fmeta:set_prec(a) + checkf(self) + checku(a) + prv.mpf_set_prec(self, a) +end + +function fmeta:set(a, base) + checkf(self) + if type(a) == "number" then + prv.mpf_set_d(self, a) + elseif type(a) == "string" then + if base == nil then + base = 10 + elseif type(base) == "number" and base < 0 then + checkbase(-base) + else + checkbase(base) + end + if prv.mpf_set_str(self, a, base) ~= 0 then + error("not a valid floating point constant in base " .. base .. ": " .. a) + end + elseif isf(a) then + prv.mpf_set(self, a) + elseif isz(a) then + prv.mpf_set_z(self, a) + else + error("unsupported type") + end +end + +function fmeta:sgn() + checkf(self) + return prv.mpf_sgn(self) +end + +function fmeta:sqrt(res) + checkf(self) + checkfopt(res) + return prv.mpf_sqrt(self, res) +end + +function sqrt(a, res) + checku(a) + checkfopt(res) + return prv.mpf_sqrt_ui(a, res) +end + +function fmeta:sub(a, res) + checkf(self) + checkfopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpf_sub(self, dtof(a), res) + elseif t == "u" then + return prv.mpf_sub_ui(self, a, res) + elseif t == "s" or t == "n" then + return prv.mpf_add_ui(self, -a, res) + elseif t == "f" then + return prv.mpf_sub(self, a, res) + else + error("unsupported type") + end +end + +function fmeta:trunc(res) + checkf(self) + checkfopt(res) + return prv.mpf_trunc(self, res) +end + +function fmeta:rdiv(a, res) + checkf(self) + checkfopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpf_div(dtof(a), self, res) + elseif t == "u" then + return prv.mpf_ui_div(a, self, res) + elseif t == "s" or t == "n" then + res = prv.mpf_ui_div(-a, self, res) + return prv.mpf_neg(res, res) + elseif t == "f" then + return prv.mpf_div(a, self, res) + else + error("unsupported type") + end +end + +function fmeta:rsub(a, res) + checkf(self) + checkfopt(res) + local t = ntype(a) + if t == "d" then + return prv.mpf_sub(dtof(a), self, res) + elseif t == "u" then + return prv.mpf_ui_sub(a, self, res) + elseif t == "s" or t == "n" then + res = prv.mpf_add_ui(self, -a, res) + return prv.mpf_neg(res, res) + elseif t == "f" then + return prv.mpf_sub(a, self, res) + else + error("unsupported type") + end +end + +function fmeta:__add(a) + if isf(self) then + return fmeta.add(self, a) + else + return fmeta.add(a, self) + end +end + +function fmeta:__sub(a) + if isf(self) then + return fmeta.sub(self, a) + else + return fmeta.rsub(a, self) + end +end + +function fmeta:__mul(a) + if isf(self) then + return fmeta.mul(self, a) + else + return fmeta.mul(a, self) + end +end + +function fmeta:__div(a) + if isf(self) then + return fmeta.div(self, a) + else + return fmeta.rdiv(a, self) + end +end + +function fmeta:__unm() + checkf(self) + return prv.mpf_neg(self) +end + +function fmeta:__pow(a) + checkf(self) + checku(a) + return prv.mpf_pow_ui(self, a) +end + +function fmeta:__lt(a) + checkf(self) + checkf(a) + return prv.mpf_cmp(self, a) < 0 +end + +function fmeta:__le(a) + checkf(self) + checkf(a) + return prv.mpf_cmp(self, a) <= 0 +end + +function fmeta:__eq(a) + checkf(self) + checkf(a) + return prv.mpf_cmp(self, a) == 0 +end + +function fmeta:__tostring() + return fmeta.format(self, ".g") +end + +function fmeta:__concat(other) + if isf(self) then + return fmeta.__tostring(self) .. other + else + return self .. fmeta.__tostring(other) + end +end + +function fmeta:__gc() + checkf(self) + return prv.mpf_clear(self) +end + +function rand(a) + if a == nil then + return prv.gmp_randinit_default() + else + checkrand(a) + return prv.gmp_randinit_set(a) + end +end + +function randmeta:seed(a) + checkrand(self) + local t = ntype(a) + if a == "u" then + prv.gmp_randseed_ui(self, a) + elseif a == "z" then + prv.gmp_randseed(self, a) + else + error("unsupported type") + end +end + +function randmeta:zbits(a, res) + checkrand(self) + checku(a) + checkzopt(res) + return prv.mpz_urandomb(self, a, res) +end + +function randmeta:z(a, res) + checkrand(self) + checkz(a) + checkzopt(res) + return prv.mpz_urandomm(self, a, res) +end + +function randmeta:fbits(a, res) + checkrand(self) + checku(a) + checkfopt(res) + return prv.mpf_urandomb(self, a, res) +end + +function randmeta:__gc() + checkrand(self) + prv.gmp_randclear(self) +end + +function randmeta:__tostring() + checkrand(self) + return "gmp random state" +end + +function zmeta:format(fmt, p) + checkz(self) + assert(type(fmt) == "string", "gmp integer format string expected") + local fw, prec, conv = + match(fmt, "^%%?([0#+ ]?%d*)(%.?%*?%d*)Z?([dioxX])$") + + if not conv or prec ~= "" and prec ~= ".*" and not match(prec, "^%.%d*$") then + error("invalid format string for gmp integer: " .. fmt) + end + + if prec == ".*" then + checkn(p) + else + assert(p == nil, "precision incorrectly specified") + end + + return prv.mpz_asprintf("%" .. fw .. prec .. "Z" .. conv, self, p) +end + +function fmeta:format(fmt, p) + checkf(self) + assert(type(fmt) == "string", "gmp floating point format string expected") + local fw, prec, conv = + match(fmt, "^%%?([0#+ ]?%d*)(%.?%*?%d*)F?([aAeEfgG])$") + + if not conv or prec ~= "" and prec ~= ".*" and not match(prec, "^%.%d*$") then + error("invalid format string for gmp floating point: " .. fmt) + end + + if prec == ".*" then + checkn(p) + else + assert(p == nil, "precision incorrectly specified") + end + + return prv.mpf_asprintf("%" .. fw .. prec .. "F" .. conv, self, p) +end + +function lgmp_lua_init(l_prv, l_aux) + aux.randmeta = randmeta + aux.zmeta = zmeta + aux.fmeta = fmeta + + prv = l_prv + aux = l_aux + + lgmpversion = prv.version + + lgmp_init = nil +end + diff --git a/src/plugin-luagmp.cc b/src/plugin-luagmp.cc new file mode 100644 index 0000000..f843400 --- /dev/null +++ b/src/plugin-luagmp.cc @@ -0,0 +1,2532 @@ + +/* + * This file is part of the lgmp package for Lua 5.1. + * + * Author: Wim Couwenberg + * Date : 2007/7/28 + * + * The lgmp package is distributed under the MIT license. See the file + * "COPYRIGHT" that came with the lgmp package distribution for license + * details. + */ + +/* + * Adapted for lua-interface by Nicolas "Pixel" Noble + */ + +#include +#include + +#include "BLua.h" +#include "Buffer.h" +#include "gmp.h" + +#ifdef FROM_LUAINTERFACE +#define NO_DLL +#endif + +#ifndef WIN32 +#define WEAK __attribute__ ((weak)) +#else +#define WEAK +#endif + +static const char lgmp_metarand[] = "randmeta"; +static const char lgmp_metaz[] = "zmeta"; +static const char lgmp_metaf[] = "fmeta"; + +static void lgmp_free_str(char *data) +{ + void (*current_free)(void *data, size_t size); + + mp_get_memory_functions(NULL, NULL, ¤t_free); + + if (current_free == NULL) + free(data); + else + (*current_free)(data, strlen(data) + 1); +} + +static gmp_randstate_t *lgmp_rawrand(lua_State *L) +{ + gmp_randstate_t *pr = (gmp_randstate_t *)lua_newuserdata(L, sizeof(gmp_randstate_t)); + lua_getfield(L, LUA_ENVIRONINDEX, lgmp_metarand); + lua_setmetatable(L, -2); + + return pr; +} + +static gmp_randstate_t *lgmp_torand(lua_State *L, int idx) +{ + return (gmp_randstate_t *)lua_touserdata(L, idx); +} + +static mpz_t *lgmp_toz(lua_State *L, int idx) +{ + return (mpz_t *)lua_touserdata(L, idx); +} + +static mpz_t *lgmp_rawz(lua_State *L) +{ + mpz_t *pz = (mpz_t *)lua_newuserdata(L, sizeof(mpz_t)); + lua_getfield(L, LUA_ENVIRONINDEX, lgmp_metaz); + lua_setmetatable(L, -2); + + return pz; +} + +static mpz_t *lgmp_optz(lua_State *L, int idx, int top) +{ + mpz_t *pz; + + if (top > 0) + lua_settop(L, top); + + if (!lua_isnil(L, idx)) + pz = lgmp_toz(L, idx); + else + { + pz = lgmp_rawz(L); + mpz_init(*pz); + lua_replace(L, idx); + } + + return pz; +} + +static mpf_t *lgmp_tof(lua_State *L, int idx) +{ + return (mpf_t *)lua_touserdata(L, idx); +} + +static mpf_t *lgmp_rawf(lua_State *L) +{ + mpf_t *pf = (mpf_t *)lua_newuserdata(L, sizeof(mpf_t)); + lua_getfield(L, LUA_ENVIRONINDEX, lgmp_metaf); + lua_setmetatable(L, -2); + + return pf; +} + +static mpf_t *lgmp_optf(lua_State *L, int idx, int top) +{ + mpf_t *pf; + + if (top > 0) + lua_settop(L, top); + + if (!lua_isnil(L, idx)) + pf = lgmp_tof(L, idx); + else + { + pf = lgmp_rawf(L); + mpf_init(*pf); + lua_replace(L, idx); + } + + return pf; +} + +static int lmpz_asprintf(lua_State *L) +{ + char *str; + const char *a1 = lua_tostring(L, 1); + mpz_t *a2 = lgmp_toz(L, 2); + + if (lua_isnil(L, 3)) + gmp_asprintf(&str, a1, *a2); + else + { + int a3 = lua_tointeger(L, 3); + gmp_asprintf(&str, a1, a3, *a2); + } + + lua_pushstring(L, str); + lgmp_free_str(str); + + return 1; +} + +static int lmpf_asprintf(lua_State *L) +{ + char *str; + const char *a1 = lua_tostring(L, 1); + mpf_t *a2 = lgmp_tof(L, 2); + + if (lua_isnil(L, 3)) + gmp_asprintf(&str, a1, *a2); + else + { + int a3 = lua_tointeger(L, 3); + gmp_asprintf(&str, a1, a3, *a2); + } + + lua_pushstring(L, str); + lgmp_free_str(str); + + return 1; +} + +static int lgmp_randinit_default(lua_State *L) +{ + gmp_randstate_t *res = lgmp_rawrand(L); + + gmp_randinit_default(*res); + + return 1; +} + +static int lgmp_randinit_lc_2exp(lua_State *L) +{ + gmp_randstate_t *res = lgmp_rawrand(L); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + unsigned long arg3 = (unsigned long)lua_tonumber(L, 3); + + gmp_randinit_lc_2exp(*res, *arg1, arg2, arg3); + + return 1; +} + +static int lgmp_randinit_lc_2exp_size(lua_State *L) +{ + gmp_randstate_t *res = lgmp_rawrand(L); + unsigned long arg = (unsigned long)lua_tonumber(L, 1); + + lua_pushinteger(L, gmp_randinit_lc_2exp_size(*res, arg)); + + return 2; +} + +static int lgmp_randinit_mt(lua_State *L) +{ + gmp_randstate_t *res = lgmp_rawrand(L); + + gmp_randinit_mt(*res); + + return 1; +} + +static int lgmp_randinit_set(lua_State *L) +{ + gmp_randstate_t *res = lgmp_rawrand(L); + gmp_randstate_t *arg = lgmp_torand(L, 1); + + gmp_randinit_set(*res, *arg); + + return 1; +} + +static int lgmp_randseed(lua_State *L) +{ + gmp_randstate_t *arg1 = lgmp_torand(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + gmp_randseed(*arg1, *arg2); + + return 0; +} + +static int lgmp_randseed_ui(lua_State *L) +{ + gmp_randstate_t *arg1 = lgmp_torand(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + gmp_randseed_ui(*arg1, arg2); + + return 0; +} + +static int lgmp_randclear(lua_State *L) +{ + gmp_randstate_t *arg = lgmp_torand(L, 1); + + gmp_randclear(*arg); + + return 0; +} + +static int lmpz_abs(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 2, 2); + mpz_t *src = lgmp_toz(L, 1); + + mpz_abs(*res, *src); + + return 1; +} + +static int lmpz_add(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_add(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_add_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_add_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_addmul(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + mpz_t *arg3 = lgmp_toz(L, 3); + + mpz_addmul(*arg1, *arg2, *arg3); + + return 0; +} + +static int lmpz_addmul_ui(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + unsigned long arg3 = (unsigned long)lua_tonumber(L, 3); + + mpz_addmul_ui(*arg1, *arg2, arg3); + + return 0; +} + +static int lmpz_and(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_and(*res, *arg1, *arg2); + + return 1; +} + +#if 0 +#define mpz_array_init __gmpz_array_init +__GMP_DECLSPEC void mpz_array_init __GMP_PROTO ((mpz_ptr, mp_size_t, mp_size_t)); +#endif + +static int lmpz_bin_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_bin_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_bin_uiui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + unsigned long arg1 = (unsigned long)lua_tonumber(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_bin_uiui(*res, arg1, arg2); + + return 1; +} + +static int lmpz_cdiv_q(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_cdiv_q(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_cdiv_q_2exp(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_cdiv_q_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_cdiv_q_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_cdiv_q_ui(*res, *arg1, arg2)); + + return 2; +} + +static int lmpz_cdiv_qr(lua_State *L) +{ + mpz_t *resq = lgmp_optz(L, 3, 4); + mpz_t *resr = lgmp_optz(L, 4, 0); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_cdiv_qr(*resq, *resr, *arg1, *arg2); + + return 2; +} + +static int lmpz_cdiv_qr_ui(lua_State *L) +{ + mpz_t *resq = lgmp_optz(L, 3, 4); + mpz_t *resr = lgmp_optz(L, 4, 0); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_cdiv_qr_ui(*resq, *resr, *arg1, arg2)); + + return 3; +} + +static int lmpz_cdiv_r(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_cdiv_r(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_cdiv_r_2exp(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_cdiv_r_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_cdiv_r_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_cdiv_r_ui(*res, *arg1, arg2)); + + return 2; +} + +static int lmpz_cdiv_ui(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_cdiv_ui(*arg1, arg2)); + + return 1; +} + +static int lmpz_clear(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + mpz_clear(*arg); + + return 0; +} + +static int lmpz_clrbit(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_clrbit(*arg1, arg2); + + return 0; +} + +static int lmpz_cmp(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + lua_pushinteger(L, mpz_cmp(*arg1, *arg2)); + + return 1; +} + +static int lmpz_cmp_d(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + double arg2 = lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_cmp_d(*arg1, arg2)); + + return 1; +} + +static int lmpz_cmp_si(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + signed long arg2 = (signed long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_cmp_si(*arg1, arg2)); + + return 1; +} + +static int lmpz_cmp_ui(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_cmp_ui(*arg1, arg2)); + + return 1; +} + +static int lmpz_cmpabs(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + lua_pushinteger(L, mpz_cmpabs(*arg1, *arg2)); + + return 1; +} + +static int lmpz_cmpabs_d(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + double arg2 = lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_cmpabs_d(*arg1, arg2)); + + return 1; +} + +static int lmpz_cmpabs_ui(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_cmpabs_ui(*arg1, arg2)); + + return 1; +} + +static int lmpz_com(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 2, 2); + mpz_t *arg = lgmp_toz(L, 1); + + mpz_com(*res, *arg); + + return 1; +} + +static int lmpz_combit(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_combit(*arg1, arg2); + + return 0; +} + +static int lmpz_congruent_p(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + mpz_t *arg3 = lgmp_toz(L, 3); + + lua_pushinteger(L, mpz_congruent_p(*arg1, *arg2, *arg3)); + + return 1; +} + +static int lmpz_congruent_2exp_p(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + unsigned long arg3 = (unsigned long)lua_tonumber(L, 3); + + lua_pushinteger(L, mpz_congruent_2exp_p(*arg1, *arg2, arg3)); + + return 1; +} + +static int lmpz_congruent_ui_p(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + unsigned long arg3 = (unsigned long)lua_tonumber(L, 3); + + lua_pushinteger(L, mpz_congruent_ui_p(*arg1, arg2, arg3)); + + return 1; +} + +static int lmpz_divexact(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_divexact(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_divexact_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_divexact_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_divisible_p(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + lua_pushinteger(L, mpz_divisible_p(*arg1, *arg2)); + + return 1; +} + +static int lmpz_divisible_ui_p(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_divisible_ui_p(*arg1, arg2)); + + return 1; +} + +static int lmpz_divisible_2exp_p(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_divisible_2exp_p(*arg1, arg2)); + + return 1; +} + +#if 0 +#define mpz_dump __gmpz_dump +__GMP_DECLSPEC void mpz_dump __GMP_PROTO ((mpz_srcptr)); + +#define mpz_export __gmpz_export +__GMP_DECLSPEC void *mpz_export __GMP_PROTO ((void *, size_t *, int, size_t, int, size_t, mpz_srcptr)); +#endif + +static int lmpz_fac_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 2, 2); + unsigned long arg = (unsigned long)lua_tonumber(L, 1); + + mpz_fac_ui(*res, arg); + + return 1; +} + +static int lmpz_fdiv_q(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_fdiv_q(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_fdiv_q_2exp(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_fdiv_q_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_fdiv_q_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_fdiv_q_ui(*res, *arg1, arg2)); + + return 2; +} + +static int lmpz_fdiv_qr(lua_State *L) +{ + mpz_t *resq = lgmp_optz(L, 3, 4); + mpz_t *resr = lgmp_optz(L, 4, 0); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_fdiv_qr(*resq, *resr, *arg1, *arg2); + + return 2; +} + +static int lmpz_fdiv_qr_ui(lua_State *L) +{ + mpz_t *resq = lgmp_optz(L, 3, 4); + mpz_t *resr = lgmp_optz(L, 4, 0); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_fdiv_qr_ui(*resq, *resr, *arg1, arg2)); + + return 3; +} + +static int lmpz_fdiv_r(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_fdiv_r(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_fdiv_r_2exp(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_fdiv_r_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_fdiv_r_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_fdiv_r_ui(*res, *arg1, arg2)); + + return 2; +} + +static int lmpz_fdiv_ui(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_fdiv_ui(*arg1, arg2)); + + return 1; +} + +static int lmpz_fib_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 2, 2); + unsigned long arg = (unsigned long)lua_tonumber(L, 1); + + mpz_fib_ui(*res, arg); + + return 1; +} + +static int lmpz_fib2_ui(lua_State *L) +{ + mpz_t *res1 = lgmp_optz(L, 2, 3); + mpz_t *res2 = lgmp_optz(L, 3, 0); + unsigned long arg = (unsigned long)lua_tonumber(L, 1); + + mpz_fib2_ui(*res1, *res2, arg); + + return 2; +} + +static int lmpz_fits_sint_p(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_fits_sint_p(*arg)); + + return 1; +} + +static int lmpz_fits_slong_p(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_fits_slong_p(*arg)); + + return 1; +} + +static int lmpz_fits_sshort_p(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_fits_sshort_p(*arg)); + + return 1; +} + +static int lmpz_fits_uint_p(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_fits_uint_p(*arg)); + + return 1; +} + +static int lmpz_fits_ulong_p(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_fits_ulong_p(*arg)); + + return 1; +} + +static int lmpz_fits_ushort_p(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_fits_ushort_p(*arg)); + + return 1; +} + +static int lmpz_gcd(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_gcd(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_gcd_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_gcd_ui(*res, *arg1, arg2)); + + return 2; +} + +static int lmpz_gcdext(lua_State *L) +{ + mpz_t *resg = lgmp_optz(L, 3, 5); + mpz_t *resp = lgmp_optz(L, 4, 0); + mpz_t *resq = lgmp_optz(L, 5, 0); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_gcdext(*resg, *resp, *resq, *arg1, *arg2); + + return 3; +} + +static int lmpz_get_d(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushnumber(L, mpz_get_d(*arg)); + + return 1; +} + +static int lmpz_get_d_2exp(lua_State *L) +{ + signed long resexp; + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushnumber(L, mpz_get_d_2exp(&resexp, *arg)); + lua_pushinteger(L, resexp); + + return 2; +} + +#if 0 +#define mpz_get_si __gmpz_get_si +__GMP_DECLSPEC /* signed */ long int mpz_get_si __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; +#endif + +static int lmpz_get_str(lua_State *L) +{ + int arg1 = lua_tointeger(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + char *res = mpz_get_str(NULL, arg1, *arg2); + lua_pushstring(L, res); + lgmp_free_str(res); + + return 1; +} + +#if 0 +#define mpz_get_ui __gmpz_get_ui +#if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_get_ui) +__GMP_DECLSPEC unsigned long int mpz_get_ui __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; +#endif + +#define mpz_getlimbn __gmpz_getlimbn +#if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_getlimbn) +__GMP_DECLSPEC mp_limb_t mpz_getlimbn __GMP_PROTO ((mpz_srcptr, mp_size_t)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; +#endif +#endif + +static int lmpz_hamdist(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + lua_pushnumber(L, mpz_hamdist(*arg1, *arg2)); + + return 1; +} + +#if 0 +#define mpz_import __gmpz_import +__GMP_DECLSPEC void mpz_import __GMP_PROTO ((mpz_ptr, size_t, int, size_t, int, size_t, __gmp_const void *)); +#endif + +static int lmpz_init(lua_State *L) +{ + mpz_t *res = lgmp_rawz(L); + mpz_init(*res); + + return 1; +} + +#if 0 +#define mpz_init2 __gmpz_init2 +__GMP_DECLSPEC void mpz_init2 __GMP_PROTO ((mpz_ptr, unsigned long)); +#endif + +static int lmpz_init_set(lua_State *L) +{ + mpz_t *res = lgmp_rawz(L); + mpz_t *arg = lgmp_toz(L, 1); + + mpz_init_set(*res, *arg); + + return 1; +} + +static int lmpz_init_set_d(lua_State *L) +{ + mpz_t *res = lgmp_rawz(L); + double arg = lua_tonumber(L, 1); + + mpz_init_set_d(*res, arg); + + return 1; +} + +#if 0 +#define mpz_init_set_si __gmpz_init_set_si +__GMP_DECLSPEC void mpz_init_set_si __GMP_PROTO ((mpz_ptr, signed long int)); +#endif + +static int lmpz_init_set_str(lua_State *L) +{ + mpz_t *res = lgmp_rawz(L); + const char *arg1 = lua_tostring(L, 1); + int arg2 = lua_tointeger(L, 2); + + if (mpz_init_set_str(*res, arg1, arg2) != 0) + lua_pushnil(L); + + return 1; +} + +#if 0 +#define mpz_init_set_ui __gmpz_init_set_ui +__GMP_DECLSPEC void mpz_init_set_ui __GMP_PROTO ((mpz_ptr, unsigned long int)); + +#define mpz_inp_raw __gmpz_inp_raw +#ifdef _GMP_H_HAVE_FILE +__GMP_DECLSPEC size_t mpz_inp_raw __GMP_PROTO ((mpz_ptr, FILE *)); +#endif + +#define mpz_inp_str __gmpz_inp_str +#ifdef _GMP_H_HAVE_FILE +__GMP_DECLSPEC size_t mpz_inp_str __GMP_PROTO ((mpz_ptr, FILE *, int)); +#endif +#endif + +static int lmpz_invert(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + lua_pushinteger(L, mpz_invert(*res, *arg1, *arg2)); + + return 2; +} + +static int lmpz_ior(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_ior(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_kronecker(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + lua_pushinteger(L, mpz_kronecker(*arg1, *arg2)); + + return 1; +} + +static int lmpz_kronecker_si(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + long arg2 = (long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_kronecker_si(*arg1, arg2)); + + return 1; +} + +static int lmpz_kronecker_ui(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_kronecker_ui(*arg1, arg2)); + + return 1; +} + +static int lmpz_si_kronecker(lua_State *L) +{ + long arg1 = (long)lua_tonumber(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + lua_pushinteger(L, mpz_si_kronecker(arg1, *arg2)); + + return 1; +} + +static int lmpz_ui_kronecker(lua_State *L) +{ + unsigned long arg1 = (unsigned long)lua_tonumber(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + lua_pushinteger(L, mpz_ui_kronecker(arg1, *arg2)); + + return 1; +} + +static int lmpz_lcm(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_lcm(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_lcm_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_lcm_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_lucnum_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 2, 2); + unsigned long arg = (unsigned long)lua_tonumber(L, 1); + + mpz_lucnum_ui(*res, arg); + + return 1; +} + +static int lmpz_lucnum2_ui(lua_State *L) +{ + mpz_t *res1 = lgmp_optz(L, 2, 3); + mpz_t *res2 = lgmp_optz(L, 3, 0); + unsigned long arg = (unsigned long)lua_tonumber(L, 1); + + mpz_lucnum2_ui(*res1, *res2, arg); + + return 2; +} + +static int lmpz_millerrabin(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + int arg2 = lua_tointeger(L, 2); + + lua_pushinteger(L, mpz_millerrabin(*arg1, arg2)); + + return 1; +} + +static int lmpz_mod(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_mod(*res, *arg1, *arg2); + + return 1; +} + +#if 0 +#define mpz_mod_ui mpz_fdiv_r_ui /* same as fdiv_r because divisor unsigned */ +#endif + +static int lmpz_mul(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_mul(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_mul_2exp(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_mul_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_mul_si(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + long arg2 = (long)lua_tonumber(L, 2); + + mpz_mul_si(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_mul_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_mul_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_neg(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 2, 2); + mpz_t *arg = lgmp_toz(L, 1); + + mpz_neg(*res, *arg); + + return 1; +} + +static int lmpz_nextprime(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 2, 2); + mpz_t *arg = lgmp_toz(L, 1); + + mpz_nextprime(*res, *arg); + + return 1; +} + +#if 0 +#define mpz_out_raw __gmpz_out_raw +#ifdef _GMP_H_HAVE_FILE +__GMP_DECLSPEC size_t mpz_out_raw __GMP_PROTO ((FILE *, mpz_srcptr)); +#endif + +#define mpz_out_str __gmpz_out_str +#ifdef _GMP_H_HAVE_FILE +__GMP_DECLSPEC size_t mpz_out_str __GMP_PROTO ((FILE *, int, mpz_srcptr)); +#endif +#endif + +static int lmpz_perfect_power_p(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_perfect_power_p(*arg)); + + return 1; +} + +static int lmpz_perfect_square_p(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_perfect_square_p(*arg)); + + return 1; +} + +static int lmpz_popcount(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushnumber(L, mpz_popcount(*arg)); + + return 1; +} + +static int lmpz_pow_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_pow_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_powm(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 4, 4); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + mpz_t *arg3 = lgmp_toz(L, 3); + + mpz_powm(*res, *arg1, *arg2, *arg3); + + return 1; +} + +static int lmpz_powm_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 4, 4); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + mpz_t *arg3 = lgmp_toz(L, 3); + + mpz_powm_ui(*res, *arg1, arg2, *arg3); + + return 1; +} + +static int lmpz_probab_prime_p(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + int arg2 = lua_tointeger(L, 2); + + lua_pushinteger(L, mpz_probab_prime_p(*arg1, arg2)); + + return 1; +} + +#if 0 +#define mpz_random __gmpz_random +__GMP_DECLSPEC void mpz_random __GMP_PROTO ((mpz_ptr, mp_size_t)); + +#define mpz_random2 __gmpz_random2 +__GMP_DECLSPEC void mpz_random2 __GMP_PROTO ((mpz_ptr, mp_size_t)); + +#define mpz_realloc2 __gmpz_realloc2 +__GMP_DECLSPEC void mpz_realloc2 __GMP_PROTO ((mpz_ptr, unsigned long)); +#endif + +static int lmpz_remove(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_remove(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_root(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_root(*res, *arg1, arg2)); + + return 2; +} + +static int lmpz_rootrem(lua_State *L) +{ + mpz_t *resro = lgmp_optz(L, 3, 4); + mpz_t *resre = lgmp_optz(L, 4, 0); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_rootrem(*resro, *resre, *arg1, arg2); + + return 2; +} + +static int lmpz_rrandomb(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + gmp_randstate_t *arg1 = lgmp_torand(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_rrandomb(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_scan0(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_scan0(*arg1, arg2)); + + return 1; +} + +static int lmpz_scan1(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_scan1(*arg1, arg2)); + + return 1; +} + +static int lmpz_set(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_set(*arg1, *arg2); + + return 0; +} + +static int lmpz_set_d(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + double arg2 = lua_tonumber(L, 2); + + mpz_set_d(*arg1, arg2); + + return 0; +} + +static int lmpz_set_f(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpz_set_f(*arg1, *arg2); + + return 0; +} + +#if 0 +#define mpz_set_q __gmpz_set_q +#if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_set_q) +__GMP_DECLSPEC void mpz_set_q __GMP_PROTO ((mpz_ptr, mpq_srcptr)); +#endif + +#define mpz_set_si __gmpz_set_si +__GMP_DECLSPEC void mpz_set_si __GMP_PROTO ((mpz_ptr, signed long int)); +#endif + +static int lmpz_set_str(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + const char *arg2 = lua_tostring(L, 2); + int arg3 = lua_tointeger(L, 3); + + lua_pushinteger(L, mpz_set_str(*arg1, arg2, arg3)); + + return 1; +} + +#if 0 +#define mpz_set_ui __gmpz_set_ui +__GMP_DECLSPEC void mpz_set_ui __GMP_PROTO ((mpz_ptr, unsigned long int)); +#endif + +static int lmpz_setbit(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_setbit(*arg1, arg2); + + return 0; +} + +static int lmpz_sgn(lua_State *L) +{ + mpz_t *arg = lgmp_toz(L, 1); + + lua_pushinteger(L, mpz_sgn(*arg)); + + return 1; +} + +#if 0 +#define mpz_size __gmpz_size +#if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_size) +__GMP_DECLSPEC size_t mpz_size __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; +#endif +#endif + +static int lmpz_sizeinbase(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + int arg2 = lua_tointeger(L, 2); + + lua_pushinteger(L, mpz_sizeinbase(*arg1, arg2)); + + return 1; +} + +static int lmpz_sqrt(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 2, 2); + mpz_t *arg = lgmp_toz(L, 1); + + mpz_sqrt(*res, *arg); + + return 1; +} + +static int lmpz_sqrtrem(lua_State *L) +{ + mpz_t *ress = lgmp_optz(L, 2, 3); + mpz_t *resr = lgmp_optz(L, 3, 0); + mpz_t *arg = lgmp_toz(L, 1); + + mpz_sqrtrem(*ress, *resr, *arg); + + return 2; +} + +static int lmpz_sub(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_sub(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_sub_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_sub_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_ui_sub(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + unsigned long arg1 = (unsigned long)lua_tonumber(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_ui_sub(*res, arg1, *arg2); + + return 1; +} + +static int lmpz_submul(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + mpz_t *arg3 = lgmp_toz(L, 3); + + mpz_submul(*arg1, *arg2, *arg3); + + return 0; +} + +static int lmpz_submul_ui(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + unsigned long arg3 = (unsigned long)lua_tonumber(L, 3); + + mpz_submul_ui(*arg1, *arg2, arg3); + + return 0; +} + +#if 0 +#define mpz_swap __gmpz_swap +__GMP_DECLSPEC void mpz_swap __GMP_PROTO ((mpz_ptr, mpz_ptr)) __GMP_NOTHROW; +#endif + +static int lmpz_tdiv_q(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_tdiv_q(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_tdiv_q_2exp(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_tdiv_q_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_tdiv_q_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_tdiv_q_ui(*res, *arg1, arg2)); + + return 2; +} + +static int lmpz_tdiv_qr(lua_State *L) +{ + mpz_t *resq = lgmp_optz(L, 3, 4); + mpz_t *resr = lgmp_optz(L, 4, 0); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_tdiv_qr(*resq, *resr, *arg1, *arg2); + + return 2; +} + +static int lmpz_tdiv_qr_ui(lua_State *L) +{ + mpz_t *resq = lgmp_optz(L, 3, 4); + mpz_t *resr = lgmp_optz(L, 4, 0); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_tdiv_qr_ui(*resq, *resr, *arg1, arg2)); + + return 3; +} + +static int lmpz_tdiv_r(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_tdiv_r(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_tdiv_r_2exp(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_tdiv_r_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_tdiv_r_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_tdiv_r_ui(*res, *arg1, arg2)); + + return 2; +} + +static int lmpz_tdiv_ui(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushnumber(L, mpz_tdiv_ui(*arg1, arg2)); + + return 1; +} + +static int lmpz_tstbit(lua_State *L) +{ + mpz_t *arg1 = lgmp_toz(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpz_tstbit(*arg1, arg2)); + + return 1; +} + +static int lmpz_ui_pow_ui(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + unsigned long arg1 = (unsigned long)lua_tonumber(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_ui_pow_ui(*res, arg1, arg2); + + return 1; +} + +static int lmpz_urandomb(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + gmp_randstate_t *arg1 = lgmp_torand(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpz_urandomb(*res, *arg1, arg2); + + return 1; +} + +static int lmpz_urandomm(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + gmp_randstate_t *arg1 = lgmp_torand(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_urandomm(*res, *arg1, *arg2); + + return 1; +} + +static int lmpz_xor(lua_State *L) +{ + mpz_t *res = lgmp_optz(L, 3, 3); + mpz_t *arg1 = lgmp_toz(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpz_xor(*res, *arg1, *arg2); + + return 1; +} + +static int lmpf_abs(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 2, 2); + mpf_t *arg = lgmp_tof(L, 1); + + mpf_abs(*res, *arg); + + return 1; +} + +static int lmpf_add(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpf_add(*res, *arg1, *arg2); + + return 1; +} + +static int lmpf_add_ui(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_add_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpf_ceil(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 2, 2); + mpf_t *arg = lgmp_tof(L, 1); + + mpf_ceil(*res, *arg); + + return 1; +} + +static int lmpf_clear(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + mpf_clear(*arg); + + return 0; +} + +static int lmpf_cmp(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + lua_pushinteger(L, mpf_cmp(*arg1, *arg2)); + + return 1; +} + +static int lmpf_cmp_d(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + double arg2 = lua_tonumber(L, 2); + + lua_pushinteger(L, mpf_cmp_d(*arg1, arg2)); + + return 1; +} + +static int lmpf_cmp_si(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + signed long arg2 = (signed long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpf_cmp_si(*arg1, arg2)); + + return 1; +} + +static int lmpf_cmp_ui(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + lua_pushinteger(L, mpf_cmp_ui(*arg1, arg2)); + + return 1; +} + +static int lmpf_div(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpf_div(*res, *arg1, *arg2); + + return 1; +} + +static int lmpf_div_2exp(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_div_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpf_div_ui(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_div_ui(*res, *arg1, arg2); + + return 1; +} + +#if 0 +#define mpf_dump __gmpf_dump +__GMP_DECLSPEC void mpf_dump __GMP_PROTO ((mpf_srcptr)); +#endif + +static int lmpf_eq(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + unsigned long arg3 = (unsigned long)lua_tonumber(L, 3); + + lua_pushinteger(L, mpf_eq(*arg1, *arg2, arg3)); + + return 1; +} + +static int lmpf_fits_sint_p(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushinteger(L, mpf_fits_sint_p(*arg)); + + return 1; +} + +static int lmpf_fits_slong_p(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushinteger(L, mpf_fits_slong_p(*arg)); + + return 1; +} + +static int lmpf_fits_sshort_p(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushinteger(L, mpf_fits_sshort_p(*arg)); + + return 1; +} + +static int lmpf_fits_uint_p(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushinteger(L, mpf_fits_uint_p(*arg)); + + return 1; +} + +static int lmpf_fits_ulong_p(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushinteger(L, mpf_fits_ulong_p(*arg)); + + return 1; +} + +static int lmpf_fits_ushort_p(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushinteger(L, mpf_fits_ushort_p(*arg)); + + return 1; +} + +static int lmpf_floor(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 2, 2); + mpf_t *arg = lgmp_tof(L, 1); + + mpf_floor(*res, *arg); + + return 1; +} + +static int lmpf_get_d(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushnumber(L, mpf_get_d(*arg)); + + return 1; +} + +static int lmpf_get_d_2exp(lua_State *L) +{ + signed long exp; + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushnumber(L, mpf_get_d_2exp(&exp, *arg)); + lua_pushnumber(L, exp); + + return 2; +} + +static int lmpf_get_default_prec(lua_State *L) +{ + lua_pushnumber(L, mpf_get_default_prec()); + + return 1; +} + +static int lmpf_get_prec(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushnumber(L, mpf_get_prec(*arg)); + + return 1; +} + +#if 0 +#define mpf_get_si __gmpf_get_si +__GMP_DECLSPEC long mpf_get_si __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; +#endif + +static int lmpf_get_str(lua_State *L) +{ + mp_exp_t exp; + int arg1 = lua_tointeger(L, 1); + size_t arg2 = lua_tointeger(L, 2); + mpf_t *arg3 = lgmp_tof(L, 3); + + char *res = mpf_get_str(NULL, &exp, arg1, arg2, *arg3); + lua_pushstring(L, res); + lua_pushnumber(L, exp); + lgmp_free_str(res); + + return 2; +} + +#if 0 +#define mpf_get_ui __gmpf_get_ui +__GMP_DECLSPEC unsigned long mpf_get_ui __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; +#endif + +static int lmpf_init(lua_State *L) +{ + mpf_t *res = lgmp_rawf(L); + + mpf_init(*res); + + return 1; +} + +#if 0 +#define mpf_init2 __gmpf_init2 +__GMP_DECLSPEC void mpf_init2 __GMP_PROTO ((mpf_ptr, unsigned long int)); +#endif + +static int lmpf_init_set(lua_State *L) +{ + mpf_t *res = lgmp_rawf(L); + mpf_t *arg = lgmp_tof(L, 1); + + mpf_init_set(*res, *arg); + + return 1; +} + +static int lmpf_init_set_d(lua_State *L) +{ + mpf_t *res = lgmp_rawf(L); + double arg = lua_tonumber(L, 1); + + mpf_init_set_d(*res, arg); + + return 1; +} + +#if 0 +#define mpf_init_set_si __gmpf_init_set_si +__GMP_DECLSPEC void mpf_init_set_si __GMP_PROTO ((mpf_ptr, signed long int)); +#endif + +static int lmpf_init_set_str(lua_State *L) +{ + mpf_t *res = lgmp_rawf(L); + const char *arg1 = lua_tostring(L, 1); + int arg2 = lua_tonumber(L, 2); + + if (mpf_init_set_str(*res, arg1, arg2) != 0) + { + lua_pushnil(L); + lua_replace(L, -2); + } + + return 1; +} + +#if 0 +#define mpf_init_set_ui __gmpf_init_set_ui +__GMP_DECLSPEC void mpf_init_set_ui __GMP_PROTO ((mpf_ptr, unsigned long int)); + +#define mpf_inp_str __gmpf_inp_str +#ifdef _GMP_H_HAVE_FILE +__GMP_DECLSPEC size_t mpf_inp_str __GMP_PROTO ((mpf_ptr, FILE *, int)); +#endif +#endif + +static int lmpf_integer_p(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushinteger(L, mpf_integer_p(*arg)); + + return 1; +} + +static int lmpf_mul(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpf_mul(*res, *arg1, *arg2); + + return 1; +} + +static int lmpf_mul_2exp(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_mul_2exp(*res, *arg1, arg2); + + return 1; +} + +static int lmpf_mul_ui(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_mul_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpf_neg(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 2, 2); + mpf_t *arg = lgmp_tof(L, 1); + + mpf_neg(*res, *arg); + + return 1; +} + +#if 0 +#define mpf_out_str __gmpf_out_str +#ifdef _GMP_H_HAVE_FILE +__GMP_DECLSPEC size_t mpf_out_str __GMP_PROTO ((FILE *, int, size_t, mpf_srcptr)); +#endif +#endif + +static int lmpf_pow_ui(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_pow_ui(*res, *arg1, arg2); + + return 1; +} + +static int lmpf_random2(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mp_size_t arg1 = (mp_size_t)lua_tonumber(L, 1); + mp_exp_t arg2 = (mp_exp_t)lua_tonumber(L, 2); + + mpf_random2(*res, arg1, arg2); + + return 1; +} + +static int lmpf_reldiff(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpf_reldiff(*res, *arg1, *arg2); + + return 1; +} + +static int lmpf_set(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpf_set(*arg1, *arg2); + + return 0; +} + +static int lmpf_set_d(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + double arg2 = lua_tonumber(L, 2); + + mpf_set_d(*arg1, arg2); + + return 0; +} + +static int lmpf_set_default_prec(lua_State *L) +{ + unsigned int arg = (unsigned int)lua_tonumber(L, 1); + + mpf_set_default_prec(arg); + + return 0; +} + +static int lmpf_set_prec(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_set_prec(*arg1, arg2); + + return 0; +} + +#if 0 +#define mpf_set_prec_raw __gmpf_set_prec_raw +__GMP_DECLSPEC void mpf_set_prec_raw __GMP_PROTO ((mpf_ptr, unsigned long int)) __GMP_NOTHROW; + +#define mpf_set_q __gmpf_set_q +__GMP_DECLSPEC void mpf_set_q __GMP_PROTO ((mpf_ptr, mpq_srcptr)); + +#define mpf_set_si __gmpf_set_si +__GMP_DECLSPEC void mpf_set_si __GMP_PROTO ((mpf_ptr, signed long int)); +#endif + +static int lmpf_set_str(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + const char *arg2 = lua_tostring(L, 2); + int arg3 = lua_tointeger(L, 3); + + lua_pushinteger(L, mpf_set_str(*arg1, arg2, arg3)); + + return 1; +} + +#if 0 +#define mpf_set_ui __gmpf_set_ui +__GMP_DECLSPEC void mpf_set_ui __GMP_PROTO ((mpf_ptr, unsigned long int)); +#endif + +static int lmpf_set_z(lua_State *L) +{ + mpf_t *arg1 = lgmp_tof(L, 1); + mpz_t *arg2 = lgmp_toz(L, 2); + + mpf_set_z(*arg1, *arg2); + + return 0; +} + +static int lmpf_sgn(lua_State *L) +{ + mpf_t *arg = lgmp_tof(L, 1); + + lua_pushinteger(L, mpf_sgn(*arg)); + + return 1; +} + +#if 0 +#define mpf_size __gmpf_size +__GMP_DECLSPEC size_t mpf_size __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; +#endif + +static int lmpf_sqrt(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 2, 2); + mpf_t *arg = lgmp_tof(L, 1); + + mpf_sqrt(*res, *arg); + + return 1; +} + +static int lmpf_sqrt_ui(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 2, 2); + unsigned int arg = (unsigned int)lua_tonumber(L, 1); + + mpf_sqrt_ui(*res, arg); + + return 1; +} + +static int lmpf_sub(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpf_sub(*res, *arg1, *arg2); + + return 1; +} + +static int lmpf_sub_ui(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + mpf_t *arg1 = lgmp_tof(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_sub_ui(*res, *arg1, arg2); + + return 1; +} + +#if 0 +#define mpf_swap __gmpf_swap +__GMP_DECLSPEC void mpf_swap __GMP_PROTO ((mpf_ptr, mpf_ptr)) __GMP_NOTHROW; +#endif + +static int lmpf_trunc(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 2, 2); + mpf_t *arg = lgmp_tof(L, 1); + + mpf_trunc(*res, *arg); + + return 1; +} + +static int lmpf_ui_div(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + unsigned long arg1 = (unsigned long)lua_tonumber(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpf_ui_div(*res, arg1, *arg2); + + return 1; +} + +static int lmpf_ui_sub(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + unsigned long arg1 = (unsigned long)lua_tonumber(L, 1); + mpf_t *arg2 = lgmp_tof(L, 2); + + mpf_ui_sub(*res, arg1, *arg2); + + return 1; +} + +static int lmpf_urandomb(lua_State *L) +{ + mpf_t *res = lgmp_optf(L, 3, 3); + gmp_randstate_t *arg1 = lgmp_torand(L, 1); + unsigned long arg2 = (unsigned long)lua_tonumber(L, 2); + + mpf_urandomb(*res, *arg1, arg2); + + return 1; +} + +static const luaL_Reg lgmp_prv[] = +{ + {"gmp_randinit_default", lgmp_randinit_default}, + {"gmp_randinit_lc_2exp", lgmp_randinit_lc_2exp}, + {"gmp_randinit_lc_2exp_size", lgmp_randinit_lc_2exp_size}, + {"gmp_randinit_mt", lgmp_randinit_mt}, + {"gmp_randinit_set", lgmp_randinit_set}, + {"gmp_randseed", lgmp_randseed}, + {"gmp_randseed_ui", lgmp_randseed_ui}, + {"gmp_randclear", lgmp_randclear}, + {"mpz_abs", lmpz_abs}, + {"mpz_add", lmpz_add}, + {"mpz_add_ui", lmpz_add_ui}, + {"mpz_addmul", lmpz_addmul}, + {"mpz_addmul_ui", lmpz_addmul_ui}, + {"mpz_and", lmpz_and}, + {"mpz_asprintf", lmpz_asprintf}, + {"mpz_bin_ui", lmpz_bin_ui}, + {"mpz_bin_uiui", lmpz_bin_uiui}, + {"mpz_cdiv_q", lmpz_cdiv_q}, + {"mpz_cdiv_q_2exp", lmpz_cdiv_q_2exp}, + {"mpz_cdiv_q_ui", lmpz_cdiv_q_ui}, + {"mpz_cdiv_qr", lmpz_cdiv_qr}, + {"mpz_cdiv_qr_ui", lmpz_cdiv_qr_ui}, + {"mpz_cdiv_r", lmpz_cdiv_r}, + {"mpz_cdiv_r_2exp", lmpz_cdiv_r_2exp}, + {"mpz_cdiv_r_ui", lmpz_cdiv_r_ui}, + {"mpz_cdiv_ui", lmpz_cdiv_ui}, + {"mpz_clear", lmpz_clear}, + {"mpz_clrbit", lmpz_clrbit}, + {"mpz_cmp", lmpz_cmp}, + {"mpz_cmp_d", lmpz_cmp_d}, + {"mpz_cmp_si", lmpz_cmp_si}, + {"mpz_cmp_ui", lmpz_cmp_ui}, + {"mpz_cmpabs", lmpz_cmpabs}, + {"mpz_cmpabs_d", lmpz_cmpabs_d}, + {"mpz_cmpabs_ui", lmpz_cmpabs_ui}, + {"mpz_com", lmpz_com}, + {"mpz_combit", lmpz_combit}, + {"mpz_congruent_p", lmpz_congruent_p}, + {"mpz_congruent_2exp_p", lmpz_congruent_2exp_p}, + {"mpz_congruent_ui_p", lmpz_congruent_ui_p}, + {"mpz_divexact", lmpz_divexact}, + {"mpz_divexact_ui", lmpz_divexact_ui}, + {"mpz_divisible_p", lmpz_divisible_p}, + {"mpz_divisible_ui_p", lmpz_divisible_ui_p}, + {"mpz_divisible_2exp_p", lmpz_divisible_2exp_p}, + {"mpz_fac_ui", lmpz_fac_ui}, + {"mpz_fdiv_q", lmpz_fdiv_q}, + {"mpz_fdiv_q_2exp", lmpz_fdiv_q_2exp}, + {"mpz_fdiv_q_ui", lmpz_fdiv_q_ui}, + {"mpz_fdiv_qr", lmpz_fdiv_qr}, + {"mpz_fdiv_qr_ui", lmpz_fdiv_qr_ui}, + {"mpz_fdiv_r", lmpz_fdiv_r}, + {"mpz_fdiv_r_2exp", lmpz_fdiv_r_2exp}, + {"mpz_fdiv_r_ui", lmpz_fdiv_r_ui}, + {"mpz_fdiv_ui", lmpz_fdiv_ui}, + {"mpz_fib_ui", lmpz_fib_ui}, + {"mpz_fib2_ui", lmpz_fib2_ui}, + {"mpz_fits_sint_p", lmpz_fits_sint_p}, + {"mpz_fits_slong_p", lmpz_fits_slong_p}, + {"mpz_fits_sshort_p", lmpz_fits_sshort_p}, + {"mpz_fits_uint_p", lmpz_fits_uint_p}, + {"mpz_fits_ulong_p", lmpz_fits_ulong_p}, + {"mpz_fits_ushort_p", lmpz_fits_ushort_p}, + {"mpz_gcd", lmpz_gcd}, + {"mpz_gcd_ui", lmpz_gcd_ui}, + {"mpz_gcdext", lmpz_gcdext}, + {"mpz_get_d", lmpz_get_d}, + {"mpz_get_d_2exp", lmpz_get_d_2exp}, + {"mpz_get_str", lmpz_get_str}, + {"mpz_hamdist", lmpz_hamdist}, + {"mpz_init", lmpz_init}, + {"mpz_init_set", lmpz_init_set}, + {"mpz_init_set_d", lmpz_init_set_d}, + {"mpz_init_set_str", lmpz_init_set_str}, + {"mpz_invert", lmpz_invert}, + {"mpz_ior", lmpz_ior}, + {"mpz_kronecker", lmpz_kronecker}, + {"mpz_kronecker_si", lmpz_kronecker_si}, + {"mpz_kronecker_ui", lmpz_kronecker_ui}, + {"mpz_si_kronecker", lmpz_si_kronecker}, + {"mpz_ui_kronecker", lmpz_ui_kronecker}, + {"mpz_lcm", lmpz_lcm}, + {"mpz_lcm_ui", lmpz_lcm_ui}, + {"mpz_lucnum_ui", lmpz_lucnum_ui}, + {"mpz_lucnum2_ui", lmpz_lucnum2_ui}, + {"mpz_millerrabin", lmpz_millerrabin}, + {"mpz_mod", lmpz_mod}, + {"mpz_mul", lmpz_mul}, + {"mpz_mul_2exp", lmpz_mul_2exp}, + {"mpz_mul_si", lmpz_mul_si}, + {"mpz_mul_ui", lmpz_mul_ui}, + {"mpz_neg", lmpz_neg}, + {"mpz_nextprime", lmpz_nextprime}, + {"mpz_perfect_power_p", lmpz_perfect_power_p}, + {"mpz_perfect_square_p", lmpz_perfect_square_p}, + {"mpz_popcount", lmpz_popcount}, + {"mpz_pow_ui", lmpz_pow_ui}, + {"mpz_powm", lmpz_powm}, + {"mpz_powm_ui", lmpz_powm_ui}, + {"mpz_probab_prime_p", lmpz_probab_prime_p}, + {"mpz_remove", lmpz_remove}, + {"mpz_root", lmpz_root}, + {"mpz_rootrem", lmpz_rootrem}, + {"mpz_rrandomb", lmpz_rrandomb}, + {"mpz_scan0", lmpz_scan0}, + {"mpz_scan1", lmpz_scan1}, + {"mpz_set", lmpz_set}, + {"mpz_set_d", lmpz_set_d}, + {"mpz_set_f", lmpz_set_f}, + {"mpz_set_str", lmpz_set_str}, + {"mpz_setbit", lmpz_setbit}, + {"mpz_sgn", lmpz_sgn}, + {"mpz_sizeinbase", lmpz_sizeinbase}, + {"mpz_sqrt", lmpz_sqrt}, + {"mpz_sqrtrem", lmpz_sqrtrem}, + {"mpz_sub", lmpz_sub}, + {"mpz_sub_ui", lmpz_sub_ui}, + {"mpz_ui_sub", lmpz_ui_sub}, + {"mpz_submul", lmpz_submul}, + {"mpz_submul_ui", lmpz_submul_ui}, + {"mpz_tdiv_q", lmpz_tdiv_q}, + {"mpz_tdiv_q_2exp", lmpz_tdiv_q_2exp}, + {"mpz_tdiv_q_ui", lmpz_tdiv_q_ui}, + {"mpz_tdiv_qr", lmpz_tdiv_qr}, + {"mpz_tdiv_qr_ui", lmpz_tdiv_qr_ui}, + {"mpz_tdiv_r", lmpz_tdiv_r}, + {"mpz_tdiv_r_2exp", lmpz_tdiv_r_2exp}, + {"mpz_tdiv_r_ui", lmpz_tdiv_r_ui}, + {"mpz_tdiv_ui", lmpz_tdiv_ui}, + {"mpz_tstbit", lmpz_tstbit}, + {"mpz_ui_pow_ui", lmpz_ui_pow_ui}, + {"mpz_urandomb", lmpz_urandomb}, + {"mpz_urandomm", lmpz_urandomm}, + {"mpz_xor", lmpz_xor}, + {"mpf_abs", lmpf_abs}, + {"mpf_add", lmpf_add}, + {"mpf_add_ui", lmpf_add_ui}, + {"mpf_asprintf", lmpf_asprintf}, + {"mpf_ceil", lmpf_ceil}, + {"mpf_clear", lmpf_clear}, + {"mpf_cmp", lmpf_cmp}, + {"mpf_cmp_d", lmpf_cmp_d}, + {"mpf_cmp_si", lmpf_cmp_si}, + {"mpf_cmp_ui", lmpf_cmp_ui}, + {"mpf_div", lmpf_div}, + {"mpf_div_2exp", lmpf_div_2exp}, + {"mpf_div_ui", lmpf_div_ui}, + {"mpf_eq", lmpf_eq}, + {"mpf_fits_sint_p", lmpf_fits_sint_p}, + {"mpf_fits_slong_p", lmpf_fits_slong_p}, + {"mpf_fits_sshort_p", lmpf_fits_sshort_p}, + {"mpf_fits_uint_p", lmpf_fits_uint_p}, + {"mpf_fits_ulong_p", lmpf_fits_ulong_p}, + {"mpf_fits_ushort_p", lmpf_fits_ushort_p}, + {"mpf_floor", lmpf_floor}, + {"mpf_get_d", lmpf_get_d}, + {"mpf_get_d_2exp", lmpf_get_d_2exp}, + {"mpf_get_default_prec", lmpf_get_default_prec}, + {"mpf_get_prec", lmpf_get_prec}, + {"mpf_get_str", lmpf_get_str}, + {"mpf_init", lmpf_init}, + {"mpf_init_set", lmpf_init_set}, + {"mpf_init_set_d", lmpf_init_set_d}, + {"mpf_init_set_str", lmpf_init_set_str}, + {"mpf_integer_p", lmpf_integer_p}, + {"mpf_mul", lmpf_mul}, + {"mpf_mul_2exp", lmpf_mul_2exp}, + {"mpf_mul_ui", lmpf_mul_ui}, + {"mpf_neg", lmpf_neg}, + {"mpf_pow_ui", lmpf_pow_ui}, + {"mpf_random2", lmpf_random2}, + {"mpf_reldiff", lmpf_reldiff}, + {"mpf_set", lmpf_set}, + {"mpf_set_d", lmpf_set_d}, + {"mpf_set_default_prec", lmpf_set_default_prec}, + {"mpf_set_prec", lmpf_set_prec}, + {"mpf_set_str", lmpf_set_str}, + {"mpf_set_z", lmpf_set_z}, + {"mpf_sgn", lmpf_sgn}, + {"mpf_sqrt", lmpf_sqrt}, + {"mpf_sqrt_ui", lmpf_sqrt_ui}, + {"mpf_sub", lmpf_sub}, + {"mpf_sub_ui", lmpf_sub_ui}, + {"mpf_trunc", lmpf_trunc}, + {"mpf_ui_div", lmpf_ui_div}, + {"mpf_ui_sub", lmpf_ui_sub}, + {"mpf_urandomb", lmpf_urandomb}, + {NULL, NULL} +}; + +static int lgmp_initialize(lua_State *L) +{ + lua_pushnumber(L, ULONG_MAX); + lua_setfield(L, 1, "ULONG_MAX"); + + lua_pushnumber(L, LONG_MIN); + lua_setfield(L, 1, "LONG_MIN"); + + lua_pushnumber(L, LONG_MAX); + lua_setfield(L, 1, "LONG_MAX"); + + lua_pushstring(L, gmp_version); + lua_setfield(L, 1, "version"); + + lua_pushvalue(L, 2); + lua_replace(L, LUA_ENVIRONINDEX); + + lua_pushvalue(L, 1); + luaL_register(L, NULL, lgmp_prv); + + return 0; +} + +extern "C" { +extern unsigned int size_lgmplib; +extern unsigned char lgmplib[]; +} + +static int _init_plugin(Lua *L) +{ + L->newtable(); // [1] = prv + L->newtable(); // [2] = aux + L->wrap_open(lgmp_initialize); + + Buffer lgmplib_buff; + lgmplib_buff.write(lgmplib, size_lgmplib); + L->load(&lgmplib_buff); + + L->call("lgmp_lua_init", LUA_GLOBALSINDEX, 2, 0); + + return 0; +} + +extern "C" { + +#ifndef NO_DLL +WEAK void init_plugin(Lua * L) { + _init_plugin(L); +} +#endif + +void luagmp_init(Lua * L) { + _init_plugin(L); +} + +} + -- cgit v1.2.3