coupled_map_lattice()

MAYAFLUX_API double MayaFlux::Kinesis::Discrete::coupled_map_lattice ( std::span< const double >  window, std::span< double >  coefs  )

noexcept

Ring of logistic maps with nearest neighbour coupling.

Each call maps every cell, f(x) = r x (1 - x), then mixes each cell with its two neighbours (periodic): x' = (1 - eps) f(x) + eps / 2 (f(left) + f(right)). Returns the mean cell minus one half. Uniform cells stay uniform; different cells give the pattern formation the coupling allows.

Definition at line 561 of file Kernels.cpp.

{
    if (coefs.size() < 3)
        return 0.0;
 
    const auto n = static_cast<size_t>(std::max(coefs.front(), 0.0));
    if (n == 0 || coefs.size() < 3 + 2 * n)
        return 0.0;
 
    const double r = coefs[1];
    const double epsilon = coefs[2];
    const auto cells = coefs.subspan(3, n);
    const auto mapped = coefs.subspan(3 + n, n);
 
    logistic_map(cells.data(), mapped.data(), r, n);
    blend_neighbours(cells.data(), mapped.data(), epsilon, n);
    return sum_of(cells) / static_cast<double>(n) - 0.5;
}