use criterion::{criterion_group, criterion_main, Criterion};

extern crate blas_src;
extern crate lapack_src;

use rlst::prelude::*;

use green_kernels::helmholtz_3d::Helmholtz3dKernel;
use green_kernels::traits::Kernel;
use green_kernels::types::GreenKernelEvalType;

use rand::SeedableRng;

const NPOINTS: usize = 1000;

pub fn helmholtz_c64_test_standard(c: &mut Criterion) {
    let mut rng = rand_chacha::ChaCha8Rng::seed_from_u64(0);

    let mut sources = rlst_dynamic_array2!(f64, [3, NPOINTS]);
    let mut targets = rlst_dynamic_array2!(f64, [3, NPOINTS]);

    let mut charges = rlst_dynamic_array1!(c64, [NPOINTS]);

    let mut result = rlst_dynamic_array1!(c64, [NPOINTS]);

    sources.fill_from_equally_distributed(&mut rng);
    targets.fill_from(sources.view());

    charges.fill_from_standard_normal(&mut rng);

    c.bench_function("Helmholtz evaluate c64", |b| {
        b.iter(|| {
            Helmholtz3dKernel::<c64>::new(1.0).evaluate_st(
                GreenKernelEvalType::Value,
                sources.data(),
                targets.data(),
                charges.data(),
                result.data_mut(),
            );
        })
    });
}

criterion_group!(benches, helmholtz_c64_test_standard,);
criterion_main!(benches);