Source author record

Nandor Sieben

Nandor Sieben 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

11works
11topics
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

11 published item(s)

preprint2020arXiv

Invariant synchrony subspaces of sets of matrices

A synchrony subspace of R^n is defined by setting certain components of the vectors equal according to an equivalence relation. Synchrony subspaces invariant under a given set of square matrices form a lattice. Applications of these invariant synchrony subspaces include equitable and almost equitable partitions of the vertices of a graph used in many areas of graph theory, balanced and exo-balanced partitions of coupled cell networks, and coset partitions of Cayley graphs. We study the basic properties of invariant synchrony subspaces and provide many examples of the applications. We also present what we call the split and cir algorithm for finding the lattice of invariant synchrony subspaces. Our theory and algorithm is further generalized for non-square matrices. This leads to the notion of tactical decompositions studied for its application in design theory.

preprint2013arXiv

Newton's Method and Symmetry for Semilinear Elliptic PDE on the Cube

We seek discrete approximations to solutions $u:Ω\to R$ of semilinear elliptic partial differential equations of the form $Δu + f_s(u) = 0$, where $f_s$ is a one-parameter family of nonlinear functions and $Ω$ is a domain in $R^d$. The main achievement of this paper is the approximation of solutions to the PDE on the cube $Ω=(0,π)^3 \subseteq R^3$. There are 323 possible isotropy subgroups of functions on the cube, which fall into 99 conjugacy classes. The bifurcations with symmetry in this problem are quite interesting, including many with 3-dimensional critical eigenspaces. Our automated symmetry analysis is necessary with so many isotropy subgroups and bifurcations among them, and it allows our code to follow one branch in each equivalence class that is created at a bifurcation point. Our most complicated result is the complete analysis of a degenerate bifurcation with a 6-dimensional critical eigenspace. This article extends the authors' work in {\it Automated Bifurcation Analysis for Nonlinear Elliptic Partial Difference Equations on Graphs} (Int. J. of Bifurcation and Chaos, 2009), wherein they combined symmetry analysis with modified implementations of the gradient Newton-Galerkin algorithm (GNGA, Neuberger and Swift) to automatically generate bifurcation diagrams and solution graphics for small, discrete problems with large symmetry groups. The code described in the current paper is efficiently implemented in parallel, allowing us to investigate a relatively fine-mesh discretization of the cube. We use the methodology and corresponding library presented in our paper {\it An MPI Implementation of a Self-Submitting Parallel Job Queue} (Int. J. of Parallel Prog., 2012).

preprint2012arXiv

An MPI Implementation of a Self-Submitting Parallel Job Queue

We present a simple and easy to apply methodology for using high-level self-submitting parallel job queues in an MPI environment. Using C++, we implemented a library of functions, MPQueue, both for testing our concepts and for use in real applications. In particular, we have applied our ideas toward solving computational combinatorics problems and for finding bifurcation diagrams of solutions of partial differential equations (PDE). Our method is general and can be applied in many situations without a lot of programming effort. The key idea is that workers themselves can easily submit new jobs to the currently running job queue. Our applications involve complicated data structures, so we employ serialization to allow data to be effortlessly passed between nodes. Using our library, one can solve large problems in parallel without being an expert in MPI. We demonstrate our methodology and the features of the library with several example programs, and give some results from our current PDE research. We show that our techniques are efficient and effective via overhead and scaling experiments.

preprint2011arXiv

Biased Weak Polyform Achievement Games

In a biased weak $(a,b)$ polyform achievement game, the maker and the breaker alternately mark $a,b$ previously unmarked cells on an infinite board, respectively. The maker's goal is to mark a set of cells congruent to a polyform. The breaker tries to prevent the maker from achieving this goal. A winning maker strategy for the $(a,b)$ game can be built from winning strategies for games involving fewer marks for the maker and the breaker. A new type of breaker strategy called the priority strategy is introduced. The winners are determined for all $(a,b)$ pairs for polyiamonds and polyominoes up to size four.

preprint2010arXiv

Automated Bifurcation Analysis for Nonlinear Elliptic Partial Difference Equations on Graphs

We seek solutions $u\in\R^n$ to the semilinear elliptic partial difference equation $-Lu + f_s(u) = 0$, where $L$ is the matrix corresponding to the Laplacian operator on a graph $G$ and $f_s$ is a one-parameter family of nonlinear functions. This article combines the ideas introduced by the authors in two papers: a) {\it Nonlinear Elliptic Partial Difference Equations on Graphs} (J. Experimental Mathematics, 2006), which introduces analytical and numerical techniques for solving such equations, and b) {\it Symmetry and Automated Branch Following for a Semilinear Elliptic PDE on a Fractal Region} wherein we present some of our recent advances concerning symmetry, bifurcation, and automation fo We apply the symmetry analysis found in the SIAM paper to arbitrary graphs in order to obtain better initial guesses for Newton's method, create informative graphics, and be in the underlying variational structure. We use two modified implementations of the gradient Newton-Galerkin algorithm (GNGA, Neuberger and Swift) to follow bifurcation branches in a robust way. By handling difficulties that arise when encountering accidental degeneracies and higher-dimension we can find many solutions of many symmetry types to the discrete nonlinear system. We present a selection of experimental results which demonstrate our algorithm's capability to automatically generate bifurcation diagrams and solution graphics starting with only an edgelis of a graph. We highlight interesting symmetry and variational phenomena.

preprint2010arXiv

Bounds on the Rubbling and Optimal Rubbling Numbers of Graphs

A pebbling move on a graph removes two pebbles at a vertex and adds one pebble at an adjacent vertex. Rubbling is a version of pebbling where an additional move is allowed. In this new move, one pebble each is removed at vertices $v$ and $w$ adjacent to a vertex $u$, and an extra pebble is added at vertex $u$. A vertex is reachable from a pebble distribution if it is possible to move a pebble to that vertex using rubbling moves. The rubbling number is the smallest number $m$ needed to guarantee that any vertex is reachable from any pebble distribution of $m$ pebbles. The optimal rubbling number is the smallest number $m$ needed to guarantee a pebble distribution of $m$ pebbles from which any vertex is reachable. We give bounds for rubbling and optimal rubbling numbers. In particular, we find an upper bound for the rubbling number of $n$-vertex, diameter $d$ graphs, and estimates for the maximum rubbling number of diameter 2 graphs. We also give a sharp upper bound for the optimal rubbling number, and sharp upper and lower bounds in terms of the diameter.

preprint2010arXiv

Computing Eigenfunctions on the Koch Snowflake: A New Grid and Symmetry

In this paper we numerically solve the eigenvalue problem $Δu + λu = 0$ on the fractal region defined by the Koch Snowflake, with zero-Dirichlet or zero-Neumann boundary conditions. The Laplacian with boundary conditions is approximated by a large symmetric matrix. The eigenvalues and eigenvectors of this matrix are computed by ARPACK. We impose the boundary conditions in a way that gives improved accuracy over the previous computations of Lapidus, Neuberger, Renka & Griffith. We extrapolate the results for grid spacing $h$ to the limit $h \rightarrow 0$ in order to estimate eigenvalues of the Laplacian and compare our results to those of Lapdus et al. We analyze the symmetry of the region to explain the multiplicity-two eigenvalues, and present a canonical choice of the two eigenfunctions that span each two-dimensional eigenspace.

preprint2010arXiv

Critical branching random walk in an IID environment

Using a high performance computer cluster, we run simulations regarding an open problem about d-dimensional critical branching random walks in a random IID environment The environment is given by the rule that at every site independently, with probability p>0, there is a cookie, completely suppressing the branching of any particle located there. Abstract. The simulations suggest self averaging: the asymptotic survival probability in n steps is the same in the annealed and the quenched case; it is \frac{2}{qn}, where q:=1-p. This particular asymptotics indicates a non-trivial phenomenon: the tail of the survival probability (both in the annealed and the quenched case) is the same as in the case of non-spatial unit time critical branching, where the branching rule is modified: branching only takes place with probability q for every particle at every iteration.

preprint2010arXiv

Rectangular Polyomino Set Weak (1,2)-achievement Games

In a polyomino set (1,2)-achievement game the maker and the breaker alternately mark one and two previously unmarked cells respectively. The maker's goal is to mark a set of cells congruent to one of a given set of polyominoes. The breaker tries to prevent the maker from achieving his goal. The teams of polyominoes for which the maker has a winning strategy is determined up to size 4. In set achievement games, it is natural to study infinitely large polyominoes. This enables the construction of super winners that characterize all winning teams up to a certain size.

preprint2010arXiv

Symmetry and Automated Branch Following for a Semilinear Elliptic PDE on a Fractal Region

We apply the Gradient-Newton-Galerkin-Algorithm (GNGA) of Neuberger & Swift to find solutions to a semilinear elliptic Dirichlet problem on the region whose boundary is the Koch snowflake. In a recent paper, we described an accurate and efficient method for generating a basis of eigenfunctions of the Laplacian on this region. In that work, we used the symmetry of the snowflake region to analyze and post-process the basis, rendering it suitable for input to the GNGA. The GNGA uses Newton's method on the eigenfunction expansion coefficients to find solutions to the semilinear problem. This article introduces the bifurcation digraph, an extension of the lattice of isotropy subgroups. For our example, the bifurcation digraph shows the 23 possible symmetry types of solutions to the PDE and the 59 generic symmetry-breaking bifurcations among these symmetry types. Our numerical code uses continuation methods, and follows branches created at symmetry-breaking bifurcations, so the human user does not need to supply initial guesses for Newton's method. Starting from the known trivial solution, the code automatically finds at least one solution with each of the symmetry types that we predict can exist. Such computationally intensive investigations necessitated the writing of automated branch following code, whereby symmetry information was used to reduce the number of computations per GNGA execution and to make intelligent branch following decisions at bifurcation points.