SummationByParts/benches/euler.rs

57 lines
1.7 KiB
Rust
Raw Normal View History

2020-01-27 20:00:44 +00:00
use criterion::{black_box, criterion_group, criterion_main, Criterion};
2020-01-28 18:31:14 +00:00
use sbp::euler::System;
2020-01-27 20:20:34 +00:00
use sbp::operators::{SbpOperator, Upwind4, UpwindOperator, SBP4};
2020-01-27 20:00:44 +00:00
2020-01-28 18:31:14 +00:00
fn advance_system<SBP: SbpOperator>(universe: &mut System<SBP>, n: usize) {
2020-01-27 20:00:44 +00:00
for _ in 0..n {
universe.advance(1.0 / 40.0 * 0.2);
}
}
2020-01-28 18:31:14 +00:00
fn advance_system_upwind<UO: UpwindOperator>(universe: &mut System<UO>, n: usize) {
2020-01-27 20:00:44 +00:00
for _ in 0..n {
universe.advance_upwind(1.0 / 40.0 * 0.2);
}
}
fn performance_benchmark(c: &mut Criterion) {
let mut group = c.benchmark_group("EulerSystem");
group.sample_size(25);
let w = 40;
let h = 26;
let x = ndarray::Array1::linspace(-10.0, 10.0, w);
let x = x.broadcast((h, w)).unwrap();
let y = ndarray::Array1::linspace(-10.0, 10.0, h);
let y = y.broadcast((w, h)).unwrap().reversed_axes();
2020-01-28 18:31:14 +00:00
let mut universe = System::<Upwind4>::new(x.into_owned(), y.into_owned());
2020-01-27 20:00:44 +00:00
group.bench_function("advance", |b| {
b.iter(|| {
2020-01-28 18:31:14 +00:00
universe.init_with_vortex(0.0, 0.0);
2020-01-27 20:00:44 +00:00
advance_system(&mut universe, black_box(20))
})
});
2020-01-28 18:31:14 +00:00
let mut universe = System::<Upwind4>::new(x.into_owned(), y.into_owned());
2020-01-27 20:00:44 +00:00
group.bench_function("advance_upwind", |b| {
b.iter(|| {
2020-01-28 18:31:14 +00:00
universe.init_with_vortex(0.0, 0.0);
2020-01-27 20:00:44 +00:00
advance_system_upwind(&mut universe, black_box(20))
})
});
2020-01-28 18:31:14 +00:00
let mut universe = System::<SBP4>::new(x.into_owned(), y.into_owned());
2020-01-27 20:00:44 +00:00
group.bench_function("advance_trad4", |b| {
b.iter(|| {
2020-01-28 18:31:14 +00:00
universe.init_with_vortex(0.0, 0.0);
2020-01-27 20:00:44 +00:00
advance_system(&mut universe, black_box(20))
})
});
group.finish();
}
criterion_group!(benches, performance_benchmark);
criterion_main!(benches);