Multiscale partial differential equations (PDEs), featuring heterogeneous coefficients oscillating across possibly non-separated scales, pose computational challenges for standard numerical techniques. Over the past two decades, a range of specialized methods has emerged that enables the efficient solution of such problems. Two prominent approaches are numerical multiscale methods with problem-adapted

Symbolic summation as an active research topic of symbolic computation provides efficient algorithmic tools for evaluating and simplifying different types of sums arising from mathematics, computer science, physics and other areas. Most of existing algorithms in symbolic summation are mainly applicable to the problem with univariate inputs. A long-term project in symbolic computation is to develop

The orthogonal group synchronization problem, which focuses on recovering orthogonal group elements from their corrupted pairwise measurements, encompasses examples such as high-dimensional Kuramoto model on general signed networks, \(\mathbb {Z}_2\)-synchronization, community detection under stochastic block models, and orthogonal Procrustes problem. The semidefinite relaxation (SDR) has proven its

We introduce a simple initialization of the Maubach bisection routine for adaptive mesh refinement which applies to any conforming initial triangulation and terminates in linear time with respect to the number of initial vertices. We show that Maubach’s routine with this initialization always terminates and generates meshes that preserve shape regularity and satisfy the closure estimate needed for

The bubble transform is a procedure to decompose differential forms, which are piecewise smooth with respect to a given triangulation of the domain, into a sum of local bubbles. In this paper, an improved version of a construction in the setting of the de Rham complex previously proposed by the authors is presented. The major improvement in the decomposition is that unlike the previous results, in

In 1991, Moore (Nonlinearity 4:199–230, 1991) raised a question about whether hydrodynamics is capable of performing computations. Similarly, in 2016, Tao (J Am Math Soc 29(3):601–674, 2016) asked whether a mechanical system, including a fluid flow, can simulate a universal Turing machine. In this expository article, we review the construction in Cardona et al. (Proc Natl Acad Sci 118(19):e2026818118

We propose a novel non-parametric learning paradigm for the identification of drift and diffusion coefficients of multi-dimensional non-linear stochastic differential equations, which relies upon discrete-time observations of the state. The key idea essentially consists of fitting a RKHS-based approximation of the corresponding Fokker–Planck equation to such observations, yielding theoretical estimates

This article is concerned with the construction and analysis of new time discretizations for the KdV equation on a torus for low-regularity solutions below \(H^1\). New harmonic analysis tools, including averaging approximations to the exponential phase functions and trilinear estimates of the KdV operator, are established for the construction and analysis of time discretizations with higher convergence

This paper introduces and develops the algebraic framework of moment polynomials, which are polynomial expressions in commuting variables and their formal mixed moments. Their positivity and optimization over probability measures supported on semialgebraic sets and subject to moment polynomial constraints is investigated. On the one hand, a positive solution to Hilbert’s 17th problem for pseudo-moments

Sharpness is an almost generic assumption in continuous optimization that bounds the distance from minima by objective function suboptimality. It facilitates the acceleration of first-order methods through restarts. However, sharpness involves problem-specific constants that are typically unknown, and restart schemes typically reduce convergence rates. Moreover, these schemes are challenging to apply

We count singular vector tuples of a system of tensors assigned to the edges of a directed hypergraph. To do so, we study the generalisation of quivers to directed hypergraphs. Assigning vector spaces to the nodes of a hypergraph and multilinear maps to its hyperedges gives a hyperquiver representation. Hyperquiver representations generalise quiver representations (where all hyperedges are edges) and

We introduce an algorithm to decompose matrix representations of the symmetric group over the reals into irreducible representations, which as a by-product also computes the multiplicities of the irreducible representations. The algorithm applied to a d-dimensional representation of \(S_n\) is shown to have a complexity of \({\mathcal {O}}(n^2 d^3)\) operations for determining which irreducible representations

A fundamental challenge in learning an unknown dynamical system is to reduce model uncertainty by making measurements while maintaining safety. In this work, we formulate a mathematical definition of what it means to safely learn a dynamical system by sequentially deciding where to initialize the next trajectory. In our framework, the state of the system is required to stay within a safety region for

While decomposition of one-parameter persistence modules behaves nicely, as demonstrated by the algebraic stability theorem, decomposition of multiparameter modules is known to be unstable in a certain precise sense. Until now, it has not been clear that there is any way to get around this and build a meaningful stability theory for multiparameter module decomposition. We introduce new tools, in particular

In optimization-based approaches to inverse problems and to statistical estimation, it is common to augment criteria that enforce data fidelity with a regularizer that promotes desired structural properties in the solution. The choice of a suitable regularizer is typically driven by a combination of prior domain information and computational considerations. Convex regularizers are attractive computationally

We obtain essentially matching upper and lower bounds for the expected max-sliced 1-Wasserstein distance between a probability measure on a separable Hilbert space and its empirical distribution from n samples. By proving a Banach space version of this result, we also obtain an upper bound, that is sharp up to a log factor, for the expected max-sliced 2-Wasserstein distance between a symmetric probability

We prove optimal convergence rates for eigenvalues and eigenvectors of the graph Laplacian on Poisson point clouds. Our results are valid down to the critical percolation threshold, yielding error estimates for relatively sparse graphs.

In this work, we develop new regularity conditions in nonsmooth analysis that parallel the stratification conditions of Whitney, Kuo, and Verdier. They quantify how subgradients interact with a certain “active manifold” that captures the nonsmooth activity of the function. Based on these new conditions, we show that several subgradient-based methods converge only to local minimizers when applied to

The criticality problem in nuclear engineering asks for the principal eigenpair of a Boltzmann operator describing neutron transport in a reactor core. Being able to reliably design, and control such reactors requires assessing these quantities within quantifiable accuracy tolerances. In this paper, we propose a paradigm that deviates from the common practice of approximately solving the corresponding

We develop Conley’s theory for multivalued maps on finite topological spaces. More precisely, for discrete-time dynamical systems generated by the iteration of a multivalued map which satisfies appropriate regularity conditions, we establish the notions of isolated invariant sets and index pairs, and use them to introduce a well-defined Conley index. In addition, we verify some of its fundamental properties

We present a randomized, inverse-free algorithm for producing an approximate diagonalization of any \(n \times n\) matrix pencil (A, B). The bulk of the algorithm rests on a randomized divide-and-conquer eigensolver for the generalized eigenvalue problem originally proposed by Ballard, Demmel and Dumitriu (Technical Report 2010). We demonstrate that this divide-and-conquer approach can be formulated

Given a domain \(\Omega \subset \mathbb {R}^n\), the de Rham complex of differential forms arises naturally in the study of problems in electromagnetism and fluid mechanics defined on \(\Omega \), and its discretization helps build stable numerical methods for such problems. For constructing such stable methods, one critical requirement is ensuring that the discrete subcomplex is cohomologically equivalent

We analyze constrained and unconstrained minimization problems on patches of tetrahedra sharing a common vertex with discontinuous piecewise polynomial data of degree p. We show that the discrete minimizers in the spaces of piecewise polynomials of degree p conforming in the \(H^1\), \({\varvec{H}}(\textbf{curl})\), or \({\varvec{H}}({\text {div}})\) spaces are as good as the minimizers in these entire

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent

The theory of group classifications has undergone significant changes in the past 25 years. New methods have been introduced, some difficult problems have been solved and group classifications have become widely available through computer algebra systems. This survey describes the state of the art of the group classification problem, its history, its recent advances and some open problems.

We address the noncommutative version of the Edmonds’ problem, which asks to determine the inner rank of a matrix in noncommuting variables. We provide an algorithm for the calculation of this inner rank by relating the problem with the distribution of a basic object in free probability theory, namely operator-valued semicircular elements. We have to solve a matrix-valued quadratic equation, for which

We present a discretization of the dynamic optimal transport problem for which we can obtain the convergence rate for the value of the transport cost to its continuous value when the temporal and spatial stepsize vanish. This convergence result does not require any regularity assumption on the measures, though experiments suggest that the rate is not sharp. Via an analysis of the duality gap we also

A large toolbox of numerical schemes for dispersive equations has been established, based on different discretization techniques such as discretizing the variation-of-constants formula (e.g., exponential integrators) or splitting the full equation into a series of simpler subproblems (e.g., splitting methods). In many situations these classical schemes allow a precise and efficient approximation. This

We consider numerical approximations of ill-posed elliptic problems with conditional stability. The notion of optimal error estimates is defined including both convergence with respect to discretisation and perturbations in data. The rate of convergence is determined by the conditional stability of the underlying continuous problem and the polynomial order of the approximation space. A proof is given