Add a rand function to mpn_mod

Author: dlesnoff

doc/source/mpn_mod.rst

@@ -109,6 +109,7 @@ Basic operations and arithmetic
               int mpn_mod_set_mpn(nn_ptr res, nn_srcptr x, slong xn, gr_ctx_t ctx)
               int mpn_mod_set_fmpz(nn_ptr res, const fmpz_t x, gr_ctx_t ctx)
               int mpn_mod_set_other(nn_ptr res, gr_ptr v, gr_ctx_t v_ctx, gr_ctx_t ctx)
+              int mpn_mod_rand(nn_ptr res, flint_rand_t state, gr_ctx_t ctx)
               int mpn_mod_randtest(nn_ptr res, flint_rand_t state, gr_ctx_t ctx)
               int mpn_mod_write(gr_stream_t out, nn_srcptr x, gr_ctx_t ctx)
               int mpn_mod_get_fmpz(fmpz_t res, nn_srcptr x, gr_ctx_t ctx)

src/mpn_mod/ring.c

@@ -202,6 +202,18 @@ mpn_mod_set_other(nn_ptr res, gr_ptr v, gr_ctx_t v_ctx, gr_ctx_t ctx)
     return gr_generic_set_other(res, v, v_ctx, ctx);
 }
 
+int
+mpn_mod_rand(nn_ptr res, flint_rand_t state, gr_ctx_t ctx)
+{
+    fmpz_t t;
+    fmpz_init(t);
+    fmpz_set_ui_array(t, MPN_MOD_CTX_MODULUS(ctx), MPN_MOD_CTX_NLIMBS(ctx));
+    fmpz_randm(t, state, t);
+    GR_IGNORE(mpn_mod_set_fmpz(res, t, ctx));
+    fmpz_clear(t);
+    return GR_SUCCESS;
+}
+
 int
 mpn_mod_randtest(nn_ptr res, flint_rand_t state, gr_ctx_t ctx)
 {