arnaucube
arnaucube/poulpyPublic Self-hosted repository
poulpy/base2k/examples/vmp.rs
base2k/examples/vmp.rsRaw
 1// use base2k::{
 2//     Encoding, FFT64, MatZnxDft, MatZnxDftOps, Module, VecZnx, VecZnxBig, VecZnxBigOps, VecZnxDft, VecZnxDftOps, VecZnxOps,
 3//     ZnxInfos, ZnxLayout, alloc_aligned,
 4// };
 5
 6fn main() {
 7    // let log_n: i32 = 5;
 8    // let n: usize = 1 << log_n;
 9
10    // let module: Module<FFT64> = Module::<FFT64>::new(n);
11    // let log_base2k: usize = 15;
12
13    // let a_cols: usize = 2;
14    // let a_size: usize = 5;
15
16    // let log_k: usize = log_base2k * a_size - 5;
17
18    // let mat_rows: usize = a_size;
19    // let mat_cols_in: usize = a_cols;
20    // let mat_cols_out: usize = 2;
21    // let mat_size: usize = a_size + 1;
22
23    // let mut tmp_bytes_vmp: Vec<u8> = alloc_aligned(
24    //     module.vmp_prepare_row_tmp_bytes(mat_cols_out, mat_size)
25    //         | module.vmp_apply_dft_tmp_bytes(
26    //             a_size,
27    //             a_size,
28    //             mat_rows,
29    //             mat_cols_in,
30    //             mat_cols_out,
31    //             mat_size,
32    //         ),
33    // );
34
35    // let mut tmp_bytes_dft: Vec<u8> = alloc_aligned(module.bytes_of_vec_znx_dft(mat_cols_out, mat_size));
36
37    // let mut a: VecZnx = module.new_vec_znx(a_cols, a_size);
38
39    // (0..a_cols).for_each(|i| {
40    //     let mut values: Vec<i64> = vec![i64::default(); n];
41    //     values[1 + i] = (1 << log_base2k) + 1;
42    //     a.encode_vec_i64(i, log_base2k, log_k, &values, 32);
43    //     a.normalize(log_base2k, i, &mut tmp_bytes_vmp);
44    //     a.print(n, i);
45    //     println!();
46    // });
47
48    // let mut mat_znx_dft: MatZnxDft<FFT64> = module.new_mat_znx_dft(mat_rows, mat_cols_in, mat_cols_out, mat_size);
49
50    // (0..a.size()).for_each(|row_i| {
51    //     let mut tmp: VecZnx = module.new_vec_znx(mat_cols_out, mat_size);
52    //     (0..mat_cols_out).for_each(|j| {
53    //         tmp.at_mut(j, row_i)[1 + j] = 1 as i64;
54    //     });
55    //     (0..mat_cols_in).for_each(|j| {
56    //         module.vmp_prepare_row(&mut mat_znx_dft, row_i, j, &tmp, &mut tmp_bytes_vmp);
57    //     })
58    // });
59
60    // let mut c_dft: VecZnxDft<FFT64> = module.new_vec_znx_dft_from_bytes_borrow(mat_cols_out, mat_size, &mut tmp_bytes_dft);
61    // module.vmp_apply_dft(&mut c_dft, &a, &mat_znx_dft, &mut tmp_bytes_vmp);
62
63    // let mut res: VecZnx = module.new_vec_znx(mat_cols_out, a_size);
64    // let mut c_big: VecZnxBig<FFT64> = c_dft.alias_as_vec_znx_big();
65    // (0..mat_cols_out).for_each(|i| {
66    //     module.vec_znx_idft_tmp_a(&mut c_big, i, &mut c_dft, i);
67    //     module.vec_znx_big_normalize(log_base2k, &mut res, i, &c_big, i, &mut tmp_bytes_vmp);
68
69    //     let mut values_res: Vec<i64> = vec![i64::default(); n];
70    //     res.decode_vec_i64(i, log_base2k, log_k, &mut values_res);
71    //     res.print(n, i);
72    //     println!();
73    //     println!("{:?}", values_res);
74    //     println!();
75    // });
76
77    // module.free();
78}