dashu-cmplx 0.5.0

A big arbitrary precision complex number library
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48

//! `rand` 0.9 `Distribution` impls for `CBig` (enable the `rand_v09` feature).

use crate::cbig::CBig;
use crate::third_party::rand::UniformCBig;
use dashu_float::round::Round;
use dashu_float::FBig;
use dashu_int::Word;
use rand_v09::distr::{Distribution, Open01, OpenClosed01, StandardUniform};
use rand_v09::Rng;

fn bridge<R: Rng + ?Sized>(rng: &mut R) -> impl dashu_int::rand::BitRng + '_ {
    dashu_int::rand::bridge_v09(rng)
}

impl<R: Round, const B: Word> Distribution<CBig<R, B>> for UniformCBig<R, B> {
    #[inline]
    fn sample<RNG: Rng + ?Sized>(&self, rng: &mut RNG) -> CBig<R, B> {
        self.sample_cbig(&mut bridge(rng))
    }
}

impl<R: Round, const B: Word> Distribution<CBig<R, B>> for StandardUniform {
    /// Each part uniform in `[0, 1)` → the unit square `[0, 1)²`, at inline precision.
    #[inline]
    fn sample<RNG: Rng + ?Sized>(&self, rng: &mut RNG) -> CBig<R, B> {
        let re: FBig<R, B> = StandardUniform.sample(rng);
        let im: FBig<R, B> = StandardUniform.sample(rng);
        CBig::from_parts(re, im)
    }
}

impl<R: Round, const B: Word> Distribution<CBig<R, B>> for Open01 {
    #[inline]
    fn sample<RNG: Rng + ?Sized>(&self, rng: &mut RNG) -> CBig<R, B> {
        let re: FBig<R, B> = Open01.sample(rng);
        let im: FBig<R, B> = Open01.sample(rng);
        CBig::from_parts(re, im)
    }
}

impl<R: Round, const B: Word> Distribution<CBig<R, B>> for OpenClosed01 {
    #[inline]
    fn sample<RNG: Rng + ?Sized>(&self, rng: &mut RNG) -> CBig<R, B> {
        let re: FBig<R, B> = OpenClosed01.sample(rng);
        let im: FBig<R, B> = OpenClosed01.sample(rng);
        CBig::from_parts(re, im)
    }
}