Function simulate

tools/streamline.rs:261–293

(
    geometries: G,
    field: &'v VectorField,
    timestep: f64,
    num_timesteps: usize,
    random_timesteps: bool,
    rng: &'v mut StdRng,
    streamline_kind: StreamlineKind,
    record_strea

#[allow(clippy::too_many_arguments)]
fn simulate<'v, G>(
    geometries: G,
    field: &'v VectorField,
    timestep: f64,
    num_timesteps: usize,
    random_timesteps: bool,
    rng: &'v mut StdRng,
    streamline_kind: StreamlineKind,
    record_streamlines: bool,
) -> impl Iterator<Item = (Geometry, Vec<LineString>)> + 'v
where
    G: IntoIterator<Item = Geometry>,
    G: 'v,
{
    geometries.into_iter().map(move |g| {
        let num_timesteps = if random_timesteps {
            let n = num_timesteps * 2; // changes mean
            let p = 0.5; // changes skew
            let dist = Binomial::new(n as u64, p).unwrap();
            dist.sample(rng) as usize
        } else {
            num_timesteps
        };
        simulate_geometry(
            g,
            field,
            timestep,
            num_timesteps,
            streamline_kind,
            record_streamlines,
        )
    })
}

fn simulate_geometry(
    geometry: Geometry,