Source author record

Henryk Fukś

Henryk Fukś appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2026arXiv

Approximating dynamics of a number-conserving cellular automaton by a finite-dimensional dynamical system

The local structure theory for cellular automata (CA) can be viewed as an finite-dimensional approximation of infinitely-dimensional system. While it is well known that this approximation works surprisingly well for some cellular automata, it is still not clear why it is the case, and which CA rules have this property. In order to shed some light on this problem, we present an example of a four input CA for which probabilities of occurrence of short blocks of symbols can be computed exactly. This rule is number conserving and possesses a blocking word. Its local structure approximation correctly predicts steady-state probabilities of small length blocks, and we present a rigorous proof of this fact, without resorting to numerical simulations. We conjecture that the number-conserving property together with the existence of the blocking word are responsible for the observed perfect agreement between the finite-dimensional approximation and the actual infinite-dimensional dynamical system.

preprint2026arXiv

Mathematics in the liturgical books of the Catholic Church: phases of the ecclesiastical moon

We use contemporary mathematical notation to describe the method for determining the age of the ecclesiastical moon as mandated by pope Gregory XIII and elaborated in the book of Christopher Clavius \emph{Romani calendarii explicatio}. The algorithm is first introduced by using the tabular method employed by liturgical books such as the Roman Missal, Breviary and Martyrology. Then we construct the recurrence equation for the epacts, derive its solution, and give a simple expression for the age of the moon on a given day of the year. We also consider the problems which can occur at the transition from December 31 to January 1 of the next year, when there could be a ``jump'' in moon's age ("saltus lunae") in years when epact corrections are applied. We propose a simple solution which fixes these problems. A summary of the formulae and listing of the implementation of relevant functions in Python is provided in the last section.

preprint2026arXiv

Ternary cellular automata induced by semigroups of order 3 are solvable

The minimal number of inputs in the local function of a non-trivial cellular automaton is two. Such a function can be viewed as as a kind of binary operation. If this operation is associative, it forms, together with the set of states, a semigroup. There are 18 semigroups of order 3 up to equivalence, and they define 18 cellular automata rules with three states. We investigate these rules with respect to solvability and show that all of them are solvable, meaning that the state of a given cell after $n$ iterations can be expressed by an explicit formula. We derive the relevant formulae for all 18 rules using some additional properties possessed by particular semigroups of order 3, such as commutativity and idempotence.

preprint2025arXiv

Solving the initial value problem for cellular automata by pattern decomposition

For many cellular automata, it is possible to express the state of a given cell after $n$ iterations as an explicit function of the initial configuration. We say that for such rules the solution of the initial value problem can be obtained. In some cases, one can construct the solution formula for the initial value problem by analyzing the spatiotemporal pattern generated by the rule and decomposing it into simpler segments which one can then describe algebraically. We show an example of a rule when such approach is successful, namely elementary rule 156. Solution of the initial value problem for this rule is constructed and then used to compute the density of ones after $n$ iterations, starting from a random initial condition. We also show how to obtain probabilities of occurrence of longer blocks of symbols.

preprint2020arXiv

Dynamics of large scale networks following a merger

We study the dynamic network of relationships among avatars in the massively multiplayer online game Planetside 2. In the spring of 2014, two separate servers of this game were merged, and as a result, two previously distinct networks were combined into one. We observed the evolution of this network in the seven month period following the merger and report our observations. We found that some structures of original networks persist in the combined network for a long time after the merger. As the original avatars are gradually removed, these structures slowly dissolve, but they remain observable for a surprisingly long time. We present a number of visualizations illustrating the post-merger dynamics and discuss time evolution of selected quantities characterizing the topology of the network.

preprint2020arXiv

Evaluating the Quality of Local Structure Approximation Using Elementary Rule 14

Cellular automata (CA) can be viewed as maps in the space of probability measures. Such maps are normally infinitely-dimensional, and in order to facilitate investigations of their properties, especially in the context of applications, finite-dimensional approximations have been proposed. The most commonly used one is known as the local structure theory, developed by H. Gutowitz et al. in 1987. In spite of the popularity of this approximation in CA research, examples of rigorous evaluations of its accuracy are lacking. In an attempt to fill this gap, we construct a local structure approximation for rule 14, and study its dynamics in a rigorous fashion, without relying on numerical experiments. We then compare the outcome with known exact results.

preprint2020arXiv

Explicit solution of the Cauchy problem for cellular automaton rule 172

Cellular automata (CA) are fully discrete alternatives to partial differential equations (PDE). For PDEs, one often considers the Cauchy problem, or initial value problem: find the solution of the PDE satisfying a given initial condition. For many PDEs of the first order in time, it is possible to find explicit formulae for the solution at the time $t>0$ if the solution is known at $t=0$. Can something similar be achieved for CA? We demonstrate that this is indeed possible in some cases, using elementary CA rule 172 as an example. We derive an explicit expression for the state of a given cell after $n$ iteration of the rule 172, assuming that states of all cells are known at $n=0$. We then show that this expression ("solution of the CA") can be used to obtain an expected value of a given cell after $n$ iterations, provided that the initial condition is drawn from a Bernoulli distribution. This can be done for both finite and infinite lattices, thus providing an interesting test case for investigating finite size effects in CA.

preprint2016arXiv

Mathematical formulae on coins, parts I and II

In the article "The Tale of Two Queens and Two Towering Figures" published in CNJ in 2012 (CNJ vol. 57 No. 5, pp. 304-315), we discussed the contributions of Copernicus and Newton to coin minting and monetary reforms, as well as the commemoration of their achievements on contemporary coins. In the current article we will examine more explicit aspects of the relationship between mathematical sciences and numismatics, namely the presence of mathematical formulae on coins. We will survey examples of coins depicting mathematical formulae, ranging from some universally recognizable ones to more advanced symbolic expressions which are known only to specialists.

preprint2016arXiv

Mathematics on coins, part I: The Tale of Two Queens and Two Towering Figures

In the first article in the series examining mathematics on coins, we discuss two great scientists who are not only featured on many coins, but also contributed to both theory and practice of coinage. The first one is Nicolaus Copernicus, author of the treaty \emph{Monetae cudendae ratio} in which he proposed the law governing competition between money and proposed reform of coinage in the Royal Prussia. The second one is Sir Isaac Newton, warden and Master of the Royal Mint, and overseer of the Great Recoinage of 1696.

preprint2016arXiv

Solving two-dimensional density classification problem with two probabilistic cellular automata

The density classification problem is one of the simplest yet non-trivial computing tasks which seem to be ideally suitable for cellular automata (CA). Unfortunately, there exists no one-dimensional two-state CA which classifies binary strings according to their densities. If, however, in place of simple cells one uses agents which change their behaviour from one rule to another after a fixed number of iterations, the classification can be performed by the traffic rule 184 and the majority rule 232. This two-rule solution cannot be easily generalized to two (or higher) dimensions, because it critically depends on a kinetic phase transition occurring in the rule 184. No rule exhibiting analogous transition is known in two dimensions, most likely because no such rule exists. We propose, therefore, to approach this problem form a slightly different angle, namely by introducing a stochastic component into each of the two rules. If one precedes each iteration of rule 184 by the stochastic "lane changing rule," and each iteration of rule 232 by the stochastic "crowd avoidance" rule, in the limit of infinitely many iterations the classification can be performed correctly with probability 1. This solution can be described either in the language of CA, or using the paradigm of agents which move and proliferate on the 2D lattice, following probabilistic rules.

preprint2015arXiv

An example of degenerate hyperbolicity in a cellular automaton with 3 states

We show that a behaviour analogous to degenerate hyperbolicity can occur in nearest-neighbour cellular automata (CA) with three states. We construct a 3-state rule by "lifting" elementary CA rule 140. Such "lifted" rule is equivalent to rule 140 when arguments are restricted to two symbols, otherwise it behaves as identity. We analyze the structure of multi-step preimages of 0, 1 and 2 under this rule by using minimal finite state machines (FSM), and exploit regularities found in these FSM. This allows to construct explicit expressions for densities of 0s and 1s after $n$ iterations of the rule starting from Bernoulli distribution. When the initial Bernoulli distribution is symmetric, the densities of all three symbols converge to their stationary values in linearly-exponential fashion, similarly as in finite-dimensional dynamical systems with hyperbolic fixed point with degenerate eigenvalues.

preprint2013arXiv

Minimal entropy approximation for cellular automata

We present a method for construction of approximate orbits of measures under the action of cellular automata which is complementary to the local structure theory. The local structure theory is based on the idea of Bayesian extension, that is, construction of a probability measure consistent with given block probabilities and maximizing entropy. If instead of maximizing entropy one minimizes it, one can develop another method for construction of approximate orbits, at the heart of which is the iteration of finitely-dimensional maps, called minimal entropy maps. We present numerical evidence that minimal entropy approximation sometimes spectacularly outperforms the local structure theory in characterizing properties of cellular automata. Density response curve for elementary CA rule 26 is used to illustrate this claim.

preprint2013arXiv

Sequences of preimages in elementary cellular automata

We search for regularities in the sequences of numbers of preimages for elementary cellular automata. For 46 out of 88 "minimal" rules, we find recognizable patterns, usually in the form of second order recurrence equations with constant coefficients. Introducing the concept of asymptotic emulation of CA rules, we then show how the regularities in the sequences of preimage numbers can be used to find rules emulating identity. We also show that the average density of nonzero sites after arbitrary number of steps (starting from disordered configuration) can be computed using the sequences of preimage numbers.

preprint2012arXiv

Classification of two-dimensional binary cellular automata with respect to surjectivity

While the surjectivity of the global map in two-dimensional cellular automata (2D CA) is undecidable in general, in specific cases one can often decide if the rule is surjective or not. We attempt to classify as many 2D CA as possible by using a sequence of tests based on the balance theorem, injectivity of the restriction to finite configurations, as well as permutivity. We introduce the notion of slice permutivity which is shown to imply surjectivity in 2D CA. The tests are applied to 2D binary CA with neighbourhoods consisting of up to five sites, considering all possible contiguous shapes of the neighbourhood. We find that if the size of the neighbourhood is less than five, complete classification of all rules is possible. Among 5-site rules, those with von Neuman neighbourhoods as well as neighbourhoods corresponding to T, V, and Z pentominos can also be completely classified.

preprint2011arXiv

Response Curves and Preimage Sequences of Two-Dimensional Cellular Automata

We consider the problem of finding response curves for a class of binary two-dimensional cellular automata with $L$-shaped neighbourhood. We show that the dependence of the density of ones after an arbitrary number of iterations, on the initial density of ones, can be calculated for a fairly large number of rules by considering preimage sets. We provide several examples and a summary of all known results. We consider a special case of initial density equal to 0.5 for other rules and compute explicitly the density of ones after $n$ iterations of the rule. This analysis includes surjective rules, which in the case of $L$-shaped neighbourhood are all found to be permutive. We conclude with the observation that all rules for which preimage curves can be computed explicitly are either finite or asymptotic emulators of identity or shift.