Related papers: Generalized Yee methods: Scalable symplectic finit…
We present a new approach to the problem of proving global stability, based on symplectic geometry and with a focus on systems with several conserved quantities. We also provide a proof of instability for integrable systems whose momentum…
Generalized additive models (GAMs) are a widely used class of models of interest to statisticians as they provide a flexible way to design interpretable models of data beyond linear models. We here propose a scalable and well-calibrated…
The real symplectic Stiefel manifold is the manifold of symplectic bases of symplectic subspaces of a fixed dimension. It features in a large variety of applications in physics and engineering. In this work, we study this manifold with the…
Neural operators (NOs) struggle with high-contrast multiscale partial differential equations (PDEs), where fine-scale heterogeneities cause large errors. To address this, we use the Generalized Multiscale Finite Element Method (GMsFEM) that…
Reduced basis methods are popular for approximately solving large and complex systems of differential equations. However, conventional reduced basis methods do not generally preserve conservation laws and symmetries of the full order model.…
The proximal Galerkin finite element method is a high-order, low-iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of point-wise bound constraints in infinite-dimensional function spaces.…
In this paper, we apply the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) to first solving a nonlinear poroelasticity problem. The arising system consists of a nonlinear pressure equation and a…
In this paper, we investigate a class of approximate Gaussian processes (GP) obtained by taking a linear combination of compactly supported basis functions with the basis coefficients endowed with a dependent Gaussian prior distribution.…
A class of Hamiltonian stochastic differential equations with multiplicative L\'{e}vy noise in the sense of Marcus, and the construction and numerical implementation methods of symplectic Euler scheme, are considered. A general symplectic…
Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative…
We study the symplectic structure and dynamics of Yang-Mills theory in the presence of a boundary. We introduce a decomposition of the fields on a Cauchy slice such that the symplectic form splits cleanly into independent bulk and edge…
This paper presents new methods for analyzing and evaluating generalized plans that can solve broad classes of related planning problems. Although synthesis and learning of generalized plans has been a longstanding goal in AI, it remains…
Many equilibrated flux recovery methods for finite element solutions rely on ad hoc or method-specific techniques, limiting their generalizability and efficiency. In this work, we introduce the Equilibrated Averaging Residual Method (EARM),…
Krylov subspace recycling is a powerful tool for solving long series of large, sparse linear systems that change slowly. In PDE constrained shape optimization, these appear naturally, as hundreds or more optimization steps are needed with…
Maxwell's equations describe the evolution of electromagnetic fields, together with constraints on the divergence of the magnetic and electric flux densities. These constraints correspond to fundamental physical laws: the nonexistence of…
This paper focuses on the numerical approximation of the linearized shallow water equations using hybridizable discontinuous Galerkin (HDG) methods, leveraging the Hamiltonian structure of the evolution system. First, we propose an…
This paper introduces a new weak Galerkin (WG) finite element method for second order elliptic equations on polytopal meshes. This method, called WG-FEM, is designed by using a discrete weak gradient operator applied to discontinuous…
The Energy Conserving Semi-Implicit Method (ECSIM) introduced by Lapenta (2017) has many advantageous properties compared to the classical semi-implicit and explicit PIC methods. Most importantly, energy conservation eliminates the growth…
We describe a new finite element method (FEM) to construct continuous equilibrium distribution functions of stellar systems. The method is a generalization of Schwarzschild's orbit superposition method from the space of discrete functions…
Finite differences, finite elements, and their generalizations are widely used for solving partial differential equations, and their high-order variants have respective advantages and disadvantages. Traditionally, these methods are treated…