ntt.rs 884 B

1234567891011121314151617181920212223242526
  1. use criterion::{black_box, criterion_group, criterion_main, Criterion};
  2. use rand::Rng;
  3. use spiral_rs::ntt::*;
  4. use spiral_rs::util::*;
  5. fn criterion_benchmark(c: &mut Criterion) {
  6. let params = get_test_params();
  7. let mut v1 = vec![0; params.crt_count * params.poly_len];
  8. let mut rng = rand::thread_rng();
  9. for i in 0..params.crt_count {
  10. for j in 0..params.poly_len {
  11. let idx = calc_index(&[i, j], &[params.crt_count, params.poly_len]);
  12. let val: u64 = rng.gen();
  13. v1[idx] = val % params.moduli[i];
  14. }
  15. }
  16. c.bench_function("nttf 2048", |b| {
  17. b.iter(|| ntt_forward(black_box(&params), black_box(&mut v1)))
  18. });
  19. c.bench_function("ntti 2048", |b| {
  20. b.iter(|| ntt_inverse(black_box(&params), black_box(&mut v1)))
  21. });
  22. }
  23. criterion_group!(benches, criterion_benchmark);
  24. criterion_main!(benches);