Parallelisation of partial differential equations via representation theory

Sheehan Olver
2024-11-12
Abstract:A little utilised but fundamental fact is that if one discretises a partial differential equation using a symmetry-adapted basis corresponding to so-called irreducible representations, the basic building blocks in representational theory, then the resulting linear system can be completely decoupled into smaller independent linear systems. That is, representation theory can be used to trivially parallelise the numerical solution of partial differential equations. This beautiful theory is introduced via a crash course in representation theory aimed at its practical utilisation, its connection with decomposing expansions in polynomials into different symmetry classes, and give examples of solving Schrödinger's equation on simple symmetric geometries like squares and cubes where there is as much as four-fold increase in the number of independent linear systems, each of a significantly smaller dimension than results from standard bases.
Numerical Analysis,Representation Theory
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to achieve the parallel solution of partial differential equations (PDEs) through representation theory. Specifically, the author proposes a method, that is, by using the symmetry - adapted basis corresponding to the irreducible representations of the symmetry group to discretize PDEs, thereby decomposing the originally complex linear system into multiple smaller independent linear subsystems. These subsystems can be solved in parallel, thus significantly improving the computational efficiency. The core idea of the paper is to use the symmetry in representation theory to simplify and accelerate the numerical solution of PDEs. For example, when dealing with the Schrödinger equation on simple symmetric geometries such as squares and cubes, this method can decompose the problem into multiple independent small - scale sub - problems, and the dimension of each sub - problem is much smaller than the result obtained by using standard basis functions. This not only improves the computational efficiency but also makes parallel computing possible. ### Specific Problem Description 1. **Application of Symmetry in PDE Solution**: - The paper first introduces how to decompose functions by understanding the symmetry under symmetric geometries. For example, a one - dimensional function can be decomposed into an even function and an odd function by reflection symmetry. - This decomposition method can be extended to squares and cubes in two - and three - dimensional spaces. Through the symmetry - adapted basis, functions can be decomposed into different symmetry classes, and each symmetry class corresponds to an irreducible representation. 2. **Application of Representation Theory**: - The author elaborates on the basic concepts of representation theory, including reducible representations and irreducible representations, and how to decompose a representation into irreducible representations by an orthogonal matrix. - In particular, the paper discusses how to use the symmetry - adapted basis to decompose PDEs so that they can be solved in parallel. For example, using the symmetry of the dihedral group \(D_4\) on a square can decompose the problem into 6 independent sub - problems; using the symmetry of the octahedral group \(O_h\) on a cube can decompose the problem into 20 independent sub - problems. 3. **Improvement of Numerical Methods**: - The paper also discusses how to apply these theories in actual numerical calculations. For example, using multivariate orthogonal polynomials as basis functions can construct the symmetry - adapted basis while maintaining numerical stability. - The author provides some specific numerical experiment results, showing the advantages of this method in improving computational efficiency. ### Summary In general, the main objective of this paper is to achieve the efficient parallel solution of PDEs through the symmetry - adapted basis in representation theory. This method is not only of great theoretical significance but also shows a significant improvement in computational efficiency in practical applications.