clfft 0.3.2

Bindings for clFFT, a FFT library for OpenCL.
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
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
111
112
113
114
115
116
117
118
119
120
121
122
123
124

extern crate clfft;
extern crate ocl;

use clfft::*;
use ocl::{MemFlags, ProQue, Buffer, EventList};


// Our kernel source code:
static KERNEL_SRC: &'static str = r#"
    #define mulc64(a,b) ((double2)((double)mad(-(a).y, (b).y, (a).x * (b).x), (double)mad((a).y, (b).x, (a).x * (b).y)))

    __kernel void multiply_vector(
                __global double2 const* const coeff,
                __global double2* const srcres)
    {
        uint const idx = get_global_id(0);
        srcres[idx] = mulc64(srcres[idx], coeff[idx]);
    }
"#;

fn main() {
    // Prepare some data
    let mut source = vec![0.0; 100];
    for i in 0..source.len()/2 {
        let x = std::f64::consts::PI * 4.0 * (i as f64 / 2.0 / source.len() as f64);
        source[2 * i] = x.sin();
        source[2 * i + 1] = x.cos();
    }
    
    let orig = source.clone();
    
    let mut triang = vec![0.0; source.len()];
    for i in 0..source.len() / 2 {
        let half = source.len() / 4;
        triang[2 * i] = if i < half {
            (half - i) as f64 / half as f64
        } else {
            (i - half) as f64 / half as f64
        };
    }
    
    println!("Frequency response: {:?}", triang);
    
    // Build ocl ProQue
    let ocl_pq = ProQue::builder()
        .src(KERNEL_SRC)
        .dims([source.len()])
        .build()
        .expect("Building ProQue");
        
    // Create buffers
    let mut in_buffer = unsafe { Buffer::builder()
                .queue(ocl_pq.queue().clone())
                .flags(MemFlags::new().read_write().copy_host_ptr())
                .len(ocl_pq.dims().clone())
                .use_host_slice(&source)
                .build().expect("Failed to create GPU input buffer") };
                
    let coef_buffer = unsafe { Buffer::builder()
                .queue(ocl_pq.queue().clone())
                .flags(MemFlags::new().read_only().copy_host_ptr())
                .len(ocl_pq.dims().clone())
                .use_host_slice(&triang)
                .build().expect("Failed to create GPU input buffer") };
    
    // Use events to schedule our kernels.
    // When `fft_finish_event` is signaled 
    // then `start_mul_event` gets triggered.
    // Also when `mul_finish_event` is signaled
    // then `start_ifft_event` gets triggered.
    // That leads to a schedule where first the FFT 
    // is executed, then the multiplication and afterwards
    // the IFFT.
    let mut fft_finish_event = EventList::new();
    let start_mul_event = fft_finish_event.clone();
    let mut mul_finish_event = EventList::new();
    let start_ifft_event = mul_finish_event.clone();
    // Make a plan
    let plan = 
        builder::<f64>()
        .precision(Precision::Precise)
        .dims([source.len() / 2])
        .input_layout(Layout::ComplexInterleaved)
        .output_layout(Layout::ComplexInterleaved)
        .bake_inplace_plan(&ocl_pq).unwrap();
        
    // Execute plan
    let mut plan = plan.enew(&mut fft_finish_event);
    plan
        .enq(Direction::Forward, &mut in_buffer)
        .expect("Enq FFT");
        
    let mul = ocl_pq.kernel_builder("multiply_vector")
        .arg_named("coef", Some(&coef_buffer))
        .arg_named("srcres", Some(&in_buffer))
		.build()
		.unwrap();
	unsafe {
		mul.cmd()
			.ewait(&start_mul_event)
			.enew(&mut mul_finish_event)
			.global_work_size([source.len() / 2])
			.enq()
			.expect("Enq Mul")
	};
    
    let mut no_events = EventList::new();
    plan
        .enew(&mut no_events) // Remove the `enew` event by passing an empty list
        .ewait(&start_ifft_event)
        .enq(Direction::Backward, &mut in_buffer)
        .expect("Enq IFFT");
    
    // Wait for calculation to finish and read results
    in_buffer.cmd()
        .read(&mut source)
        .enq()
        .expect("Transferring result vector from the GPU back to memory failed");
        
        
    ocl_pq.queue().finish().unwrap();
    println!("Input data (time): {:?}", orig);
    println!("Output data (time): {:?}", source);
}