Related papers: High-degree compression functions on alternative m…
The elliptic interface problems with discontinuous and high-contrast coefficients appear in many applications and often lead to huge condition numbers of the corresponding linear systems. Thus, it is highly desired to construct high order…
In this note we obtain effective lower bounds for the canonical heights of non-torsion points on $E(\mathbb{Q})$ by making use of suitable elliptic curve ideal class pairings $$\Psi_{E,-D}: E(\mathbb{Q})\times E_{-D}(\mathbb{Q})\mapsto…
In recent years, a new class of mixed finite elements -- compatible-strain mixed finite elements (CSMFEs) -- has emerged that uses the differential complex of nonlinear elasticity. Their excellent performance in benchmark problems, such as…
This paper develops the high-order entropy stable (ES) finite difference schemes for multi-dimensional compressible Euler equations with the van der Waals equation of state (EOS) on adaptive moving meshes. Semi-discrete schemes are first…
In this paper we establish new simple local geometric criteria for discrete entropic curvature introduced in [47] that are powerful enough to capture many geometric properties of complex models arising in mathematical physics. These results…
Symmetry is a key property of numerical methods. The geometric properties of symmetric schemes make them an attractive option for integrating Hamiltonian systems, whilst their ability to exactly recover the initial condition without the…
The article is devoted to the classical problems about the relationships between elliptic functions and hyperelliptic functions of genus 2. It contains new results, as well as a derivation from them of well-known results on these issues.…
We propose a new method for unconstrained optimization of a smooth and strongly convex function, which attains the optimal rate of convergence of Nesterov's accelerated gradient descent. The new algorithm has a simple geometric…
We consider a scalar function depending on a numerical solution of an initial value problem, and its second-derivative (Hessian) matrix for the initial value. The need to extract the information of the Hessian or to solve a linear system…
We consider the problem of maximizing submodular functions; while this problem is known to be NP-hard, several numerically efficient local search techniques with approximation guarantees are available. In this paper, we propose a novel…
We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems…
Standard numerical algorithms like the fast multipole method or $\mathcal{H}$-matrix schemes rely on low-rank approximations of the underlying kernel function. For high-frequency problems, the ranks grow rapidly as the mesh is refined, and…
Trust region and cubic regularization methods have demonstrated good performance in small scale non-convex optimization, showing the ability to escape from saddle points. Each iteration of these methods involves computation of gradient,…
The complexity of the elliptic curve method of factorization (ECM) is proven under the celebrated conjecture of existence of smooth numbers in short intervals. In this work we tackle a different version of ECM which is actually much more…
This work is related to PHG (Parallel Hierarchical Grid). PHG is a toolbox for developing parallel adaptive finite element programs, which is under active development at the State Key Laboratory of Scientific and Engineering Computing. The…
In this work, we establish higher-order div-curl type estimates in the sense of Coifman, Lions, Meyer & Semmes, in a local setting for elliptic homogeneous linear differential operators with smooth coefficients acting on localizable Hardy…
We develop \emph{geometric optimisation} on the manifold of Hermitian positive definite (HPD) matrices. In particular, we consider optimising two types of cost functions: (i) geodesically convex (g-convex); and (ii) log-nonexpansive (LN).…
We consider the Euler type integral associated to the configuration space of points on an elliptic curve, which is an analogue of the hypergeometric function associated to the configuration space of points on a projective line. We calculate…
Though completely integrable Camassa-Holm (CH) equation and Degasperis-Procesi (DP) equation are cast in the same peakon family, they possess the second- and third-order Lax operators, respectively. From the viewpoint of algebro-geometrical…
The numerical approximation of solutions of parametric or stochastic hyperbolic PDEs is still a serious challenge. Because of shock singularities, most methods from the elliptic and parabolic regime, such as reduced basis methods, POD or…