qudit-tensor 0.3.3

Accelerated and Extensible Quantum 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
49
50
51
52
53

use qudit_core::{ComplexScalar, memory::MemoryBuffer};
use qudit_expr::{DifferentiationLevel, WriteFunc};

use super::super::buffer::SizedTensorBuffer;

pub struct WriteStruct<C: ComplexScalar, const D: DifferentiationLevel> {
    write_fns: [WriteFunc<C::R>; D],
    param_map: Vec<u64>,
    output_map: Vec<u64>,
    const_map: Vec<bool>,
    buffer: SizedTensorBuffer<C>,
}

impl<C: ComplexScalar, const D: DifferentiationLevel> WriteStruct<C, D> {
    pub fn new(
        write_fns: [WriteFunc<C::R>; D],
        param_map: Vec<u64>,
        output_map: Vec<u64>,
        const_map: Vec<bool>,
        buffer: SizedTensorBuffer<C>,
    ) -> Self {
        Self {
            write_fns,
            param_map,
            output_map,
            const_map,
            buffer,
        }
    }

    #[inline(always)]
    pub fn evaluate<const E: DifferentiationLevel>(
        &self,
        params: &[C::R],
        memory: &mut MemoryBuffer<C>,
    ) {
        if E > D {
            panic!("Unsafe TNVM evaluation.");
        }

        unsafe {
            let ptr = self.buffer.as_ptr_mut(memory) as *mut C::R;
            (self.write_fns[E - 1])(
                params.as_ptr(),
                ptr,
                self.param_map.as_ptr(),
                self.output_map.as_ptr(),
                (self.buffer.unit_memory_size() * 2) as u64,
                self.const_map.as_ptr(),
            );
        }
    }
}