curv/arithmetic/
samplable.rs1use rand::{rngs::OsRng, RngCore};
2
3use super::traits::{BitManipulation, Converter, Samplable, Zero};
4use super::BigInt;
5
6impl Samplable for BigInt {
7 fn sample_below(upper: &Self) -> Self {
8 assert!(*upper > Self::zero());
9
10 let bits = upper.bit_length();
11 loop {
12 let n = Self::sample(bits);
13 if n < *upper {
14 return n;
15 }
16 }
17 }
18
19 fn sample_range(lower: &Self, upper: &Self) -> Self {
20 assert!(upper > lower);
21 lower + Self::sample_below(&(upper - lower))
22 }
23
24 fn strict_sample_range(lower: &Self, upper: &Self) -> Self {
25 assert!(upper > lower);
26 loop {
27 let n = lower + Self::sample_below(&(upper - lower));
28 if n > *lower && n < *upper {
29 return n;
30 }
31 }
32 }
33
34 fn sample(bit_size: usize) -> Self {
35 if bit_size == 0 {
36 return BigInt::zero();
37 }
38 let mut rng = OsRng;
39 let bytes = (bit_size - 1) / 8 + 1;
40 let mut buf: Vec<u8> = vec![0; bytes];
41 rng.fill_bytes(&mut buf);
42 BigInt::from_bytes(&*buf) >> (bytes * 8 - bit_size)
43 }
44
45 fn strict_sample(bit_size: usize) -> Self {
46 if bit_size == 0 {
47 return BigInt::zero();
48 }
49 loop {
50 let n = Self::sample(bit_size);
51 if n.bit_length() == bit_size {
52 return n;
53 }
54 }
55 }
56}