Source author record

James W. Swift

James W. Swift 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

7works
9topics
3close 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

7 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.

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

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

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.

preprint1996arXiv

Stability of periodic solutions in series arrays of Josephson junctions with internal Capacitance

A mystery surrounds the stability properties of the splay-phase periodic solutions to a series array of N Josephson junction oscillators. Contrary to what one would expect from dynamical systems theory, the splay state appears to be neutrally stable for a wide range of system parameters. It has been explained why the splay state must be neutrally stable when the Stewart-McCumber parameter beta is zero. In this paper we complete the explanation of the apparent neutral stability; we show that the splay state is typically hyperbolic -- either asymptotically stable or unstable -- when beta > 0. We conclude that there is only a single unit Floquet multiplier, based on accurate and systematic computations of the Floquet multipliers for beta ranging from 0 to 10. However, N-2 multipliers are extremely close to 1 for beta larger than about 1. In addition, two more Floquet multipliers approach 1 as beta becomes large. We visualize the global dynamics responsible for these nearly degenerate multipliers, and then estimate them accurately by a multiple time-scale analysis. For N=4 junctions the analysis also predicts that the system converges toward either the in-phase state, the splay state, or two clusters of two oscillators, depending on the parameters.