ni...@lysator.liu.se (Niels Möller) writes: > Below is a revised draft interface. Comments appreciated. I think I'll > try to get this in fairly soon.
Another revision, this time as a patch. As you can see, I renamed things a bit more. Compile tested only... Regards, /Niels diff -r d5352c3cce8b gmp-h.in --- a/gmp-h.in Wed Feb 20 22:01:47 2013 +0100 +++ b/gmp-h.in Wed Feb 20 22:42:36 2013 +0100 @@ -1117,6 +1117,22 @@ using std::FILE; #define mpz_eor __gmpz_xor __GMP_DECLSPEC void mpz_xor (mpz_ptr, mpz_srcptr, mpz_srcptr); +#define mpz_limbs_read __gmpz_limbs_read +__GMP_DECLSPEC mp_srcptr mpz_limbs_read (mpz_srcptr); + +#define mpz_limbs_write __gmpz_limbs_write +__GMP_DECLSPEC mp_ptr mpz_limbs_write (mpz_ptr, mp_size_t); + +#define mpz_limbs_modify __gmpz_limbs_modify +__GMP_DECLSPEC mp_ptr mpz_limbs_modify (mpz_ptr, mp_size_t); + +#define mpz_limbs_finish __gmpz_limbs_finish +__GMP_DECLSPEC void mpz_limbs_finish (mp_ptr, mp_size_t); + +#define mpz_roinit_n __gmpz_roinit_n +__GMP_DECLSPEC mpz_srcptr mpz_roinit_n (mpz_ptr, mp_srcptr, mp_size_t); + +#define MPZ_ROINIT_N(xp, xs) {{0, (xs),(xp) }} /**************** Rational (i.e. Q) routines. ****************/ diff -r d5352c3cce8b mpz/Makefile.am --- a/mpz/Makefile.am Wed Feb 20 22:01:47 2013 +0100 +++ b/mpz/Makefile.am Wed Feb 20 22:42:36 2013 +0100 @@ -43,12 +43,13 @@ libmpz_la_SOURCES = aors.h aors_ui.h fit import.c init.c init2.c inits.c inp_raw.c inp_str.c \ invert.c ior.c iset.c iset_d.c iset_si.c iset_str.c iset_ui.c \ jacobi.c kronsz.c kronuz.c kronzs.c kronzu.c \ - lcm.c lcm_ui.c lucnum_ui.c lucnum2_ui.c mfac_uiui.c millerrabin.c \ + lcm.c lcm_ui.c limbs_read.c limbs_write.c limbs_modify.c limbs_finish.c \ + lucnum_ui.c lucnum2_ui.c mfac_uiui.c millerrabin.c \ mod.c mul.c mul_2exp.c mul_si.c mul_ui.c n_pow_ui.c neg.c nextprime.c \ oddfac_1.c \ out_raw.c out_str.c perfpow.c perfsqr.c popcount.c pow_ui.c powm.c \ powm_sec.c powm_ui.c pprime_p.c prodlimbs.c primorial_ui.c random.c random2.c \ - realloc.c realloc2.c remove.c root.c rootrem.c rrandomb.c \ + realloc.c realloc2.c remove.c roinit_n.c root.c rootrem.c rrandomb.c \ scan0.c scan1.c set.c set_d.c set_f.c set_q.c set_si.c set_str.c \ set_ui.c setbit.c size.c sizeinbase.c sqrt.c sqrtrem.c sub.c sub_ui.c \ swap.c tdiv_ui.c tdiv_q.c tdiv_q_2exp.c tdiv_q_ui.c tdiv_qr.c \ diff -r d5352c3cce8b mpz/limbs_finish.c --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mpz/limbs_finish.c Wed Feb 20 22:42:36 2013 +0100 @@ -0,0 +1,29 @@ +/* mpz_finish_limbs -- Update mpz after writing to the limb array. + +Copyright 2013 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library is free software; you can redistribute it and/or modify +it under the terms of the GNU Lesser General Public License as published by +the Free Software Foundation; either version 3 of the License, or (at your +option) any later version. + +The GNU MP Library is distributed in the hope that it will be useful, but +WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY +or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ + +#include "gmp.h" +#include "gmp-impl.h" + +void +mpz_limbs_finish (mpz_ptr x, mp_size_t xs) +{ + mp_size_t xn = ABS(xs); + MPN_NORMALIZE (PTR (x), xn); + SIZ (x) = xs < 0 ? -xn : xn; +} diff -r d5352c3cce8b mpz/limbs_modify.c --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mpz/limbs_modify.c Wed Feb 20 22:42:36 2013 +0100 @@ -0,0 +1,28 @@ +/* mpz_limbs_modify -- Read-and-modify access to the mpn-style limb array. + +Copyright 2013 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library is free software; you can redistribute it and/or modify +it under the terms of the GNU Lesser General Public License as published by +the Free Software Foundation; either version 3 of the License, or (at your +option) any later version. + +The GNU MP Library is distributed in the hope that it will be useful, but +WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY +or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ + +#include "gmp.h" +#include "gmp-impl.h" + +mp_ptr +mpz_limbs_modify (mpz_ptr x, mp_size_t n) +{ + ASSERT (n > 0); + return MPZ_REALLOC (x, n); +} diff -r d5352c3cce8b mpz/limbs_read.c --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mpz/limbs_read.c Wed Feb 20 22:42:36 2013 +0100 @@ -0,0 +1,27 @@ +/* mpz_limbs_read -- Read access to the mpn-style limb array. + +Copyright 2013 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library is free software; you can redistribute it and/or modify +it under the terms of the GNU Lesser General Public License as published by +the Free Software Foundation; either version 3 of the License, or (at your +option) any later version. + +The GNU MP Library is distributed in the hope that it will be useful, but +WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY +or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ + +#include "gmp.h" +#include "gmp-impl.h" + +mp_srcptr +mpz_limbs_read (mpz_srcptr x) +{ + return PTR(x); +} diff -r d5352c3cce8b mpz/limbs_write.c --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mpz/limbs_write.c Wed Feb 20 22:42:36 2013 +0100 @@ -0,0 +1,28 @@ +/* mpz_limbs_write -- Write access to the mpn-style limb array. + +Copyright 2013 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library is free software; you can redistribute it and/or modify +it under the terms of the GNU Lesser General Public License as published by +the Free Software Foundation; either version 3 of the License, or (at your +option) any later version. + +The GNU MP Library is distributed in the hope that it will be useful, but +WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY +or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ + +#include "gmp.h" +#include "gmp-impl.h" + +mp_ptr +mpz_limbs_write (mpz_ptr x, mp_size_t n) +{ + ASSERT (n > 0); + return MPZ_NEWALLOC (x, n); +} diff -r d5352c3cce8b mpz/roinit_n.c --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mpz/roinit_n.c Wed Feb 20 22:42:36 2013 +0100 @@ -0,0 +1,30 @@ +/* mpz_roinit_n -- Initialize mpz with read-only limb array. + +Copyright 2013 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library is free software; you can redistribute it and/or modify +it under the terms of the GNU Lesser General Public License as published by +the Free Software Foundation; either version 3 of the License, or (at your +option) any later version. + +The GNU MP Library is distributed in the hope that it will be useful, but +WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY +or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ + +#include "gmp.h" +#include "gmp-impl.h" + +mpz_srcptr +mpz_roinit_n (mpz_ptr x, mp_srcptr xp, mp_size_t xs) +{ + ALLOC (x) = 0; + SIZ (x) = xs; + PTR (x) = xp; + return x; +} -- Niels Möller. PGP-encrypted email is preferred. Keyid C0B98E26. Internet email is subject to wholesale government surveillance. _______________________________________________ gmp-devel mailing list gmp-devel@gmplib.org http://gmplib.org/mailman/listinfo/gmp-devel