Maja Gwozdz

2papers

2 Papers

10.9DSMar 24
No weakly factor-universal cellular automaton

Maja Gwozdz

Hochman asked whether there exists a cellular automaton $F$ such that every cellular automaton is a factor of $F$ in the dynamical sense. In particular, we do not require the factor map to commute with the spatial shifts. We show that no such cellular automaton exists. More generally, if $F$ weakly factors onto the radius-zero $q$-clock automaton $C_q^{(k)}$, then every periodic point of $F$ has period divisible by $q$. For a cellular automaton $F:A^{\mathbb Z^d}\to A^{\mathbb Z^d}$, define $φ_F:A\to A$ by $F(\underline a)=\underline{φ_F(a)}$, and let $g_F$ be the greatest common divisor of the cycle lengths of $φ_F$. We prove that if $C_q^{(k)}$ is a weak factor of $F$, then $q\mid g_F$ holds. It follows that the action of $F$ on constant configurations yields an explicit divisibility obstruction to clock weak factors.

CLMay 9, 2021
High-performance symbolic-numerics via multiple dispatch

Shashi Gowda, Yingbo Ma, Alessandro Cheli et al.

As mathematical computing becomes more democratized in high-level languages, high-performance symbolic-numeric systems are necessary for domain scientists and engineers to get the best performance out of their machine without deep knowledge of code optimization. Naturally, users need different term types either to have different algebraic properties for them, or to use efficient data structures. To this end, we developed Symbolics.jl, an extendable symbolic system which uses dynamic multiple dispatch to change behavior depending on the domain needs. In this work we detail an underlying abstract term interface which allows for speed without sacrificing generality. We show that by formalizing a generic API on actions independent of implementation, we can retroactively add optimized data structures to our system without changing the pre-existing term rewriters. We showcase how this can be used to optimize term construction and give a 113x acceleration on general symbolic transformations. Further, we show that such a generic API allows for complementary term-rewriting implementations. We demonstrate the ability to swap between classical term-rewriting simplifiers and e-graph-based term-rewriting simplifiers. We showcase an e-graph ruleset which minimizes the number of CPU cycles during expression evaluation, and demonstrate how it simplifies a real-world reaction-network simulation to halve the runtime. Additionally, we show a reaction-diffusion partial differential equation solver which is able to be automatically converted into symbolic expressions via multiple dispatch tracing, which is subsequently accelerated and parallelized to give a 157x simulation speedup. Together, this presents Symbolics.jl as a next-generation symbolic-numeric computing environment geared towards modeling and simulation.