GECCO 2026

Asynchronous Neuro-Evolutionary Symbolic Regression with Maturity-Based Replacement

1 CIIRC, CTU in Prague, Czech Republic   2 Dept. of Cognitive Robotics, TU Delft, The Netherlands

Neural network representation used in ANSR. The figure shows a feedforward network with two hidden layers, learnable units (white), copy units (green), active learnable links (blue), inactive learnable links (blue dotted), active skip connections (red), and inactive skip connections (red dotted).

Neural network representing the expression f(x₀, x₁) = c₀ · (c₁ · (c₂ · x₀) · (c₃ · x₁)) / (c₄ · x₀ + c₅ · x₁).
Blue edges are trainable links; red edges are skip connections.
Bold edges remain enabled after training and form the analytic formula.


Abstract

We consider a neuro-evolutionary symbolic regression, an approach in which mathematical formulas are internally represented by feedforward neural networks. Gradient-free global exploration of the network topology space is combined with gradient-based parameter learning.

From the perspective of evolutionary algorithms, this approach has a potential weakness: newly generated networks produced by genetic operators usually require intensive parameter tuning. If this optimization has not yet converged, such networks may be prematurely discarded due to poor performance, thereby losing the opportunity to contribute their true potential to the evolutionary process. Such a waste of innovative, potentially promising topologies can reduce the algorithm's efficiency.

To address this limitation, we propose a neuro-evolutionary algorithm that relies on asynchronous refinement of network parameters and a maturity-based replacement strategy. Thus, newly created networks are allowed to coexist and develop alongside their parental networks, while being compared only with peers of similar maturity. The proposed method was evaluated on six benchmarks and compared to the original evolutionary approach.

Key Result: ANSR consistently outperforms the standard neuro-evolutionary baseline (SEASR) in solution quality and convergence efficiency across all six benchmark problems — including three synthetic and three physics-based test cases.

Model evolution for the parallel resistors problem

Model evolution for the Nguyen-9c problem

Model evolution for f(x) = sin(2.7 x) + 0.37 x^2 + 0.4

Model evolution for f(x) = sin(2.7 x) + 0.37 x^2 + 0.4

Solution obtained using the approximation of sin(2.7 x)  =  2 sin(1.35 x) cos(1.35 x)