cumath 0.2.7

Cuda-based matrix/vector computations
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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110


use std::ptr;

use super::{cuda::*, curand_ffi::*};
use cuvector::CuVectorDeref;
use meta::result::*;

pub use super::curand_ffi::CurandRngType;



pub struct CurandGenerator {
    handle: *mut StructCurandGenerator,
}

impl Drop for CurandGenerator {
    fn drop(&mut self) {
        curand_destroy_generator(self.handle);
    }
}

impl CurandGenerator {
    pub fn new(rng_type: CurandRngType) -> CumathResult<CurandGenerator> {
        let mut handle = ptr::null_mut();
        match curand_create_generator(&mut handle, rng_type) {
            Some(err) => Err(CumathError::new(err)),
            None => Ok(CurandGenerator { handle }),
        }
    }

    pub fn set_stream(&mut self, stream: &CudaStream) {
        curand_set_stream(self.handle, stream.stream);
    }

    pub fn generate_uniform(&mut self, output: &mut CuVectorDeref<f32>) {
        curand_generate_uniform(self.handle, output.as_mut_ptr(), output.len());
    }

    pub fn generate_uniform_range(&mut self, output: &mut CuVectorDeref<f32>, min: f32, max: f32, stream: &CudaStream) {
        #[cfg(not(feature = "disable_checks"))] {
            assert!(min <= max, "min > max");
        }
            curand_generate_uniform(self.handle, output.as_mut_ptr(), output.len());
        unsafe {
            ::kernel::VectorPacked_aypb_f32(max-min, output.as_mut_ptr(), min, output.len() as i32, stream.stream);
        }
    }

    pub fn generate_normal(&mut self, output: &mut CuVectorDeref<f32>, mean: f32, stddev: f32) {
        #[cfg(not(feature = "disable_checks"))] {
            assert!(stddev >= 0.0, "stddev < 0.0");
        }
        curand_generate_normal(self.handle, output.as_mut_ptr(), output.len(), mean, stddev);
    }

    pub fn generate_lognormal(&mut self, output: &mut CuVectorDeref<f32>, mean: f32, stddev: f32) {
        #[cfg(not(feature = "disable_checks"))] {
            assert!(stddev >= 0.0, "stddev < 0.0");
        }
        curand_generate_lognormal(self.handle, output.as_mut_ptr(), output.len(), mean, stddev);
    }

    pub fn generate_poisson(&mut self, output: &mut CuVectorDeref<f32>, lambda: f32) {
        #[cfg(not(feature = "disable_checks"))] {
            assert!(lambda >= 0.0, "lambda < 0.0");
        }
        curand_generate_poisson(self.handle, output.as_mut_ptr(), output.len(), lambda);
    }
}


#[cfg(test)]
mod tests {
    use super::*;
    use cuvector::*;

    #[test]
    fn curand_generate_uniform() {
        let mut generator = CurandGenerator::new(CurandRngType::PseudoDefault).unwrap();

        let mut vector = CuVector::<f32>::new(0.0, 10);
        generator.generate_uniform(&mut vector);

        let mut buffer = vec![0.0; 10];
        vector.clone_to_host(&mut buffer);

        buffer.iter().for_each(|x| {
            assert!(x >= &0.0 && x <= &1.0);
        });
    }

    #[test]
    fn curand_generate_uniform_range() {
        let min = -5.0;
        let max = 15.0;

        let mut generator = CurandGenerator::new(CurandRngType::PseudoDefault).unwrap();

        let mut vector = CuVector::<f32>::new(0.0, 10);
        generator.generate_uniform_range(&mut vector, min, max, &super::DEFAULT_STREAM);

        let mut buffer = vec![0.0; 10];
        vector.clone_to_host(&mut buffer);

        buffer.iter().for_each(|x| {
            assert!(*x >= min && *x <= max);
        });
    }

}