{"title":"Discrete Helmholtz Decompositions of Piecewise Constant and Piecewise Affine Vector and Tensor Fields","authors":"","doi":"10.1007/s10208-024-09642-1","DOIUrl":"https://doi.org/10.1007/s10208-024-09642-1","url":null,"abstract":"<h3>Abstract</h3> <p>Discrete Helmholtz decompositions dissect piecewise polynomial vector fields on simplicial meshes into piecewise gradients and rotations of finite element functions. This paper concisely reviews established results from the literature which all restrict to the lowest-order case of piecewise constants. Its main contribution consists of the generalization of these decompositions to 3D and of novel decompositions for piecewise affine vector fields in terms of Fortin–Soulie functions. While the classical lowest-order decompositions include one conforming and one nonconforming part, the decompositions of piecewise affine vector fields require a nonconforming enrichment in both parts. The presentation covers two and three spatial dimensions as well as generalizations to deviatoric tensor fields in the context of the Stokes equations and symmetric tensor fields for the linear elasticity and fourth-order problems. While the proofs focus on contractible domains, generalizations to multiply connected domains and domains with non-connected boundary are discussed as well.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"46 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140015412","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Maria Alberich-Carramiñana, Jordi Guàrdia, Enric Nart, Adrien Poteaux, Joaquim Roé, Martin Weimann
{"title":"Polynomial Factorization Over Henselian Fields","authors":"Maria Alberich-Carramiñana, Jordi Guàrdia, Enric Nart, Adrien Poteaux, Joaquim Roé, Martin Weimann","doi":"10.1007/s10208-024-09646-x","DOIUrl":"https://doi.org/10.1007/s10208-024-09646-x","url":null,"abstract":"<p>We present an algorithm that, given an irreducible polynomial <i>g</i> over a general valued field (<i>K</i>, <i>v</i>), finds the factorization of <i>g</i> over the Henselianization of <i>K</i> under certain conditions. The analysis leading to the algorithm follows the footsteps of Ore, Mac Lane, Okutsu, Montes, Vaquié and Herrera–Olalla–Mahboub–Spivakovsky, whose work we review in our context. The correctness is based on a key new result (Theorem 4.10), exhibiting relations between generalized Newton polygons and factorization in the context of an arbitrary valuation. This allows us to develop a polynomial factorization algorithm and an irreducibility test that go beyond the classical discrete, rank-one case. These foundational results may find applications for various computational tasks involved in arithmetic of function fields, desingularization of hypersurfaces, multivariate Puiseux series or valuation theory.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"237 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139924244","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Stable Spectral Methods for Time-Dependent Problems and the Preservation of Structure","authors":"Arieh Iserles","doi":"10.1007/s10208-024-09647-w","DOIUrl":"https://doi.org/10.1007/s10208-024-09647-w","url":null,"abstract":"<p>This paper is concerned with orthonormal systems in real intervals, given with zero Dirichlet boundary conditions. More specifically, our interest is in systems with a skew-symmetric differentiation matrix (this excludes orthonormal polynomials). We consider a simple construction of such systems and pursue its ramifications. In general, given any <span>(text {C}^1(a,b))</span> weight function such that <span>(w(a)=w(b)=0)</span>, we can generate an orthonormal system with a skew-symmetric differentiation matrix. Except for the case <span>(a=-infty )</span>, <span>(b=+infty )</span>, only few powers of that matrix are bounded and we establish a connection between properties of the weight function and boundedness. In particular, we examine in detail two weight functions: the Laguerre weight function <span>(x^alpha textrm{e}^{-x})</span> for <span>(x>0)</span> and <span>(alpha >0)</span> and the ultraspherical weight function <span>((1-x^2)^alpha )</span>, <span>(xin (-1,1))</span>, <span>(alpha >0)</span>, and establish their properties. Both weights share a most welcome feature of <i>separability,</i> which allows for fast computation. The quality of approximation is highly sensitive to the choice of <span>(alpha )</span>, and we discuss how to choose optimally this parameter, depending on the number of zero boundary conditions.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"240 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139750190","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On the Existence of Monge Maps for the Gromov–Wasserstein Problem","authors":"Théo Dumont, Théo Lacombe, François-Xavier Vialard","doi":"10.1007/s10208-024-09643-0","DOIUrl":"https://doi.org/10.1007/s10208-024-09643-0","url":null,"abstract":"<p>The Gromov–Wasserstein problem is a non-convex optimization problem over the polytope of transportation plans between two probability measures supported on two spaces, each equipped with a cost function evaluating similarities between points. Akin to the standard optimal transportation problem, it is natural to ask for conditions guaranteeing some structure on the optimizers, for instance, if these are induced by a (Monge) map. We study this question in Euclidean spaces when the cost functions are either given by (i) inner products or (ii) squared distances, two standard choices in the literature. We establish the existence of an optimal map in case (i) and of an optimal 2-map (the union of the graphs of two maps) in case (ii), both under an absolute continuity condition on the source measure. Additionally, in case (ii) and in dimension one, we numerically design situations where optimizers of the Gromov–Wasserstein problem are 2-maps but are not maps. This suggests that our result cannot be improved in general for this cost. Still in dimension one, we additionally establish the optimality of monotone maps under some conditions on the measures, thereby giving insight into why such maps often appear to be optimal in numerical experiments.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"100 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139750239","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Discrete Pseudo-differential Operators and Applications to Numerical Schemes","authors":"Erwan Faou, Benoît Grébert","doi":"10.1007/s10208-024-09645-y","DOIUrl":"https://doi.org/10.1007/s10208-024-09645-y","url":null,"abstract":"<p>We consider a class of discrete operators introduced by O. Chodosh, acting on infinite sequences and mimicking standard properties of pseudo-differential operators. By using a new approach, we extend this class to finite or periodic sequences, allowing a general representation of discrete pseudo-differential operators obtained by finite differences approximations and easily transferred to time discretizations. In particular we can define the notion of order and regularity, and we recover the fundamental property, well known in pseudo-differential calculus, that the commutator of two discrete operators gains one order of regularity. As examples of practical applications, we revisit standard error estimates for the convergence of splitting methods, obtaining in some Hamiltonian cases no loss of derivative in the error estimates, in particular for discretizations of general waves and/or water-waves equations. Moreover, we give an example of preconditioner constructions inspired by normal form analysis to deal with the similar question for more general cases.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"258 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139750263","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Strong Norm Error Bounds for Quasilinear Wave Equations Under Weak CFL-Type Conditions","authors":"Benjamin Dörich","doi":"10.1007/s10208-024-09639-w","DOIUrl":"https://doi.org/10.1007/s10208-024-09639-w","url":null,"abstract":"<p>In the present paper, we consider a class of quasilinear wave equations on a smooth, bounded domain. We discretize it in space with isoparametric finite elements and apply a semi-implicit Euler and midpoint rule as well as the exponential Euler and midpoint rule to obtain four fully discrete schemes. We derive rigorous error bounds of optimal order for the semi-discretization in space and the fully discrete methods in norms which are stronger than the classical <span>(H^1times L^2)</span> energy norm under weak CFL-type conditions. To confirm our theoretical findings, we also present numerical experiments.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"84 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139733608","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Analysis of a Modified Regularity-Preserving Euler Scheme for Parabolic Semilinear SPDEs: Total Variation Error Bounds for the Numerical Approximation of the Invariant Distribution","authors":"Charles-Edouard Bréhier","doi":"10.1007/s10208-024-09644-z","DOIUrl":"https://doi.org/10.1007/s10208-024-09644-z","url":null,"abstract":"<p>We propose a modification of the standard linear implicit Euler integrator for the weak approximation of parabolic semilinear stochastic PDEs driven by additive space-time white noise. This new method can easily be combined with a finite difference method for the spatial discretization. The proposed method is shown to have improved qualitative properties compared with the standard method. First, for any time-step size, the spatial regularity of the solution is preserved, at all times. Second, the proposed method preserves the Gaussian invariant distribution of the infinite dimensional Ornstein–Uhlenbeck process obtained when the nonlinearity is removed, for any time-step size. The weak order of convergence of the proposed method is shown to be equal to 1/2 in a general setting, like for the standard Euler scheme. A stronger weak approximation result is obtained when considering the approximation of a Gibbs invariant distribution, when the nonlinearity is a gradient: one obtains an approximation in total variation distance of order 1/2, which does not hold for the standard method. This is the first result of this type in the literature and this is the major and most original result of this article.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"98 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139710666","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Phaseless Sampling on Square-Root Lattices","authors":"Philipp Grohs, Lukas Liehr","doi":"10.1007/s10208-024-09640-3","DOIUrl":"https://doi.org/10.1007/s10208-024-09640-3","url":null,"abstract":"<p>Due to its appearance in a remarkably wide field of applications, such as audio processing and coherent diffraction imaging, the short-time Fourier transform (STFT) phase retrieval problem has seen a great deal of attention in recent years. A central problem in STFT phase retrieval concerns the question for which window functions <span>(g in {L^2({mathbb R}^d)})</span> and which sampling sets <span>(Lambda subseteq {mathbb R}^{2d})</span> is every <span>(f in {L^2({mathbb R}^d)})</span> uniquely determined (up to a global phase factor) by phaseless samples of the form </p><span>$$begin{aligned} |V_gf(Lambda )| = left{ |V_gf(lambda )|: lambda in Lambda right} , end{aligned}$$</span><p>where <span>(V_gf)</span> denotes the STFT of <i>f</i> with respect to <i>g</i>. The investigation of this question constitutes a key step towards making the problem computationally tractable. However, it deviates from ordinary sampling tasks in a fundamental and subtle manner: recent results demonstrate that uniqueness is unachievable if <span>(Lambda )</span> is a lattice, i.e <span>(Lambda = A{mathbb Z}^{2d}, A in textrm{GL}(2d,{mathbb R}))</span>. Driven by this discretization barrier, the present article centers around the initiation of a novel sampling scheme which allows for unique recovery of any square-integrable function via phaseless STFT-sampling. Specifically, we show that square-root lattices, i.e., sets of the form </p><span>$$begin{aligned} Lambda = A left( sqrt{{mathbb Z}} right) ^{2d}, sqrt{{mathbb Z}} = { pm sqrt{n}: n in {mathbb N}_0 }, end{aligned}$$</span><p>guarantee uniqueness of the STFT phase retrieval problem. The result holds for a large class of window functions, including Gaussians</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"38 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139710668","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Low-Dimensional Invariant Embeddings for Universal Geometric Learning","authors":"Nadav Dym, Steven J. Gortler","doi":"10.1007/s10208-024-09641-2","DOIUrl":"https://doi.org/10.1007/s10208-024-09641-2","url":null,"abstract":"<p>This paper studies separating invariants: mappings on <i>D</i>-dimensional domains which are invariant to an appropriate group action and which separate orbits. The motivation for this study comes from the usefulness of separating invariants in proving universality of equivariant neural network architectures. We observe that in several cases the cardinality of separating invariants proposed in the machine learning literature is much larger than the dimension <i>D</i>. As a result, the theoretical universal constructions based on these separating invariants are unrealistically large. Our goal in this paper is to resolve this issue. We show that when a continuous family of semi-algebraic separating invariants is available, separation can be obtained by randomly selecting <span>(2D+1 )</span> of these invariants. We apply this methodology to obtain an efficient scheme for computing separating invariants for several classical group actions which have been studied in the invariant learning literature. Examples include matrix multiplication actions on point clouds by permutations, rotations, and various other linear groups. Often the requirement of invariant separation is relaxed and only generic separation is required. In this case, we show that only <span>(D+1)</span> invariants are required. More importantly, generic invariants are often significantly easier to compute, as we illustrate by discussing generic and full separation for weighted graphs. Finally we outline an approach for proving that separating invariants can be constructed also when the random parameters have finite precision.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"25 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139710669","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Exotic B-Series and S-Series: Algebraic Structures and Order Conditions for Invariant Measure Sampling","authors":"Eugen Bronasco","doi":"10.1007/s10208-023-09638-3","DOIUrl":"https://doi.org/10.1007/s10208-023-09638-3","url":null,"abstract":"<p>B-Series and generalizations are a powerful tool for the analysis of numerical integrators. An extension named exotic aromatic B-Series was introduced to study the order conditions for sampling the invariant measure of ergodic SDEs. Introducing a new symmetry normalization coefficient, we analyze the algebraic structures related to exotic B-Series and S-Series. Precisely, we prove the relationship between the Grossman–Larson algebras over exotic and grafted forests and the corresponding duals to the Connes–Kreimer coalgebras and use it to study the natural composition laws on exotic S-Series. Applying this algebraic framework to the derivation of order conditions for a class of stochastic Runge–Kutta methods, we present a multiplicative property that ensures some order conditions to be satisfied automatically.</p>","PeriodicalId":55151,"journal":{"name":"Foundations of Computational Mathematics","volume":"6 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139505892","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}