oat_rust 0.2.0

User-friendly tools for applied topology
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

//! Functions to generate point clouds


use std::{f64::consts::PI, ops::Range};

use rand::Rng;

/// Return `m` points evenly spaced on the unit circle
pub fn unit_circle( m: usize, noise_range: Option<Range<f64>> ) -> Vec< Vec< f64 > > {
    let mut rng = rand::thread_rng();

    let circpoint = |k: usize | { 
        let theta = k as f64 * 2.0 * PI / m as f64; 
        let vec = vec![ theta.cos(), theta.sin() ]; 
        vec          
    };
    let mut vec: Vec<_> = (0..m).map(circpoint ).collect();
    if let Some( r ) = noise_range {
        for row in vec.iter_mut() {
            for t in row.iter_mut() { *t += rng.gen_range( r.clone()) }
        }
    }           
    vec
}