Relay synchronization in complex networks is characterized by the synchronization of remote parts of the network due to their interaction via a relay. In multilayer networks, distant layers that are not connected directly can synchronize due to signal propagation via relay layers. In this work, we investigate relay synchronization of partial synchronization patterns like chimera states in three-layer networks of interacting FitzHugh-Nagumo oscillators. We demonstrate that the phenomenon of relay synchronization is robust to topological random inhomogeneities of small-world type in the layer networks. We show that including randomness in the connectivity structure either of the remote network layers or of the relay layer increases the range of interlayer coupling strength where relay synchronization can be observed.

NetworkDynamics.jl is an easy-to-use and computationally efficient package for simulating heterogeneous dynamical systems on complex networks, written in Julia, a high-level, high-performance, dynamic programming language. By combining state-of-the-art solver algorithms from DifferentialEquations.jl with efficient data structures, NetworkDynamics.jl achieves top performance while supporting advanced features such as events, algebraic constraints, time delays, noise terms, and automatic differentiation.