### Top 10 Arxiv Papers Today in Numerical Analysis

##### #1. Automatic discrete differentiation and its applications
###### Ai Ishikawa, Takaharu Yaguchi
In this paper, a method for automatically deriving energy-preserving numerical methods for the Euler-Lagrange equation and the Hamilton equation is proposed. The derived energy-preserving scheme is based on the discrete gradient method. In the proposed approach, the discrete gradient, which is a key tool for designing the scheme, is automatically computed by a similar algorithm to the automatic differentiation. Besides, the discrete gradient coincides with the usual gradient if the two arguments required to define the discrete gradient are the same. Hence the proposed method is an extension of the automatic differentiation in the sense that the proposed method derives not only the discrete gradient but also the usual gradient. Due to this feature, both energy-preserving integrators and variational (and hence symplectic) integrators can be implemented in the same programming code simultaneously. This allows users to freely switch between the energy-preserving numerical method and the symplectic numerical method in accordance with...
more | pdf | html
None.
###### Tweets
mathNAb: Ai Ishikawa, Takaharu Yaguchi : Automatic discrete differentiation and its applications https://t.co/Yc99rjEJHC https://t.co/z7pFpZM296
SciFi: Automatic discrete differentiation and its applications. https://t.co/bxQCw1iumH
arxiv_cs_LG: Automatic discrete differentiation and its applications. Ai Ishikawa and Takaharu Yaguchi https://t.co/3nuo2ITeIi
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 7026
Unqiue Words: 1677

##### #2. A Two-stage Classification Method for High-dimensional Data and Point Clouds
###### Xiaohao Cai, Raymond Chan, Xiaoyu Xie, Tieyong Zeng
High-dimensional data classification is a fundamental task in machine learning and imaging science. In this paper, we propose a two-stage multiphase semi-supervised classification method for classifying high-dimensional data and unstructured point clouds. To begin with, a fuzzy classification method such as the standard support vector machine is used to generate a warm initialization. We then apply a two-stage approach named SaT (smoothing and thresholding) to improve the classification. In the first stage, an unconstraint convex variational model is implemented to purify and smooth the initialization, followed by the second stage which is to project the smoothed partition obtained at stage one to a binary partition. These two stages can be repeated, with the latest result as a new initialization, to keep improving the classification quality. We show that the convex model of the smoothing stage has a unique solution and can be solved by a specifically designed primal-dual algorithm whose convergence is guaranteed. We test our...
more | pdf | html
None.
###### Tweets
mathNAb: Xiaohao Cai, Raymond Chan, Xiaoyu Xie, Tieyong Zeng : A Two-stage Classification Method for High-dimensional Data and Point Clouds https://t.co/0zrmfnsrHy https://t.co/OkM3COnjya
arxiv_cscv: A Two-stage Classification Method for High-dimensional Data and Point Clouds https://t.co/UAXeXr2Nrb
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 4
Total Words: 0
Unqiue Words: 0

##### #3. Uniform bounds for invariant subspace perturbations
###### Anil Damle, Yuekai Sun
For a fixed matrix A and perturbation E we develop purely deterministic bounds on how invariant subspaces of A and A+E can differ when measured by a suitable "row-wise" metric rather than via traditional norms such as two or Frobenius. Understanding perturbations of invariant subspaces with respect to such metrics is becoming increasingly important across a wide variety of applications and therefore necessitates new theoretical developments. Under minimal assumptions we develop new bounds on subspace perturbations under the two-to-infinity matrix norm and show in what settings these row-wise differences in the invariant subspaces can be significantly smaller than the two or Frobenius norm differences. We also demonstrate that the constitutive pieces of our bounds are necessary absent additional assumptions and therefore our results provide a natural starting point for further analysis of specific problems.
more | pdf | html
None.
###### Tweets
mathNAb: Anil Damle, Yuekai Sun : Uniform bounds for invariant subspace perturbations https://t.co/EghAgbJH33 https://t.co/PB31CQd13i
StatsPapers: Uniform bounds for invariant subspace perturbations. https://t.co/YoW9CDpswP
###### Github

Numerical experiments for "Row-wise invariant subspace perturbation bounds"

Repository: rowwise-perturbation
User: asdamle
Language: MATLAB
Stargazers: 0
Subscribers: 1
Forks: 0
Open Issues: 0
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 10307
Unqiue Words: 2263

##### #4. Guaranteed a posteriori error bounds for low rank tensor approximate solutions
###### Sergey Dolgov, Tomáš Vejchodský
We propose guaranteed and fully computable upper bound on the energy norm of the error in low rank Tensor Train (TT) approximate solutions of (possibly) high dimensional reaction-diffusion problems. The error bound is obtained from Euler--Lagrange equations for a complementary flux reconstruction problem, which are solved in the low rank TT representation using the block Alternating Linear Scheme. This bound is guaranteed to be above the energy norm of the total error, including the discretization error, the tensor approximation error, and the error in the solver of linear algebraic equations. Numerical examples with the Poisson equation and the Schr\"odinger equation with the Henon-Heiles potential in up to 40 dimensions are presented to illustrate the efficiency of this approach.
more | pdf | html
None.
###### Tweets
mathNAb: Sergey Dolgov, Tomáš Vejchodský : Guaranteed a posteriori error bounds for low rank tensor approximate solutions https://t.co/bozqWBxa56 https://t.co/tQWsSwh8qT
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #5. Fast algorithm for computing nonlocal operators with finite interaction distance
###### Xiaochuan Tian, Bjorn Engquist
Developments of nonlocal operators for modeling processes that traditionally have been described by local differential operators have been increasingly active during the last few years. One example is peridynamics for brittle materials and another is nonstandard diffusion including the use of fractional derivatives. A major obstacle for application of these methods is the high computational cost from the numerical implementation of the nonlocal operators. It is natural to consider fast methods of fast multipole or hierarchical matrix type to overcome this challenge. Unfortunately the relevant kernels do not satisfy the standard necessary conditions. In this work a new class of fast algorithms is developed and analyzed, which is some cases reduces the computational complexity of applying nonlocal operators to essentially the same order of magnitude as the complexity of standard local numerical methods.
more | pdf | html
None.
###### Tweets
mathNAb: Xiaochuan Tian, Bjorn Engquist : Fast algorithm for computing nonlocal operators with finite interaction distance https://t.co/DyUkd88fCz https://t.co/i258iRzUDn
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #6. Lagrangian discretization of crowd motion and linear diffusion
###### Quentin Merigot, Hugo Leclerc, Quentin Mérigot, Filippo Santambrogio, Federico Stra
We study a model of crowd motion following a gradient vector field, with possibly additional interaction terms such as attraction/repulsion, and we present a numerical scheme for its solution through a Lagrangian discretization. The density constraint of the resulting particles is enforced by means of a partial optimal transport problem at each time step. We prove the convergence of the discrete measures to a solution of the continuous PDE describing the crowd motion in dimension one. In a second part, we show how a similar approach can be used to construct a Lagrangian discretization of a linear advection-diffusion equation, interpreted as a gradient flow in Wasserstein space. We provide also a numerical implementation in 2D to demonstrate the feasibility of the computations.
more | pdf | html
None.
###### Tweets
mathNAb: Quentin Merigot (LMO), Hugo Leclerc (LMO), Quentin Mérigot (MGMI), Filippo Santambrogio (LMO), Federico Stra (LMO) : Lagrangian discretization of crowd motion and linear diffusion https://t.co/bJSupPCNEr https://t.co/kzJNMLnHJU
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 5
Total Words: 0
Unqiue Words: 0

##### #7. On local analysis
###### Felipe Cucker, Teresa Krick
We extend to Gaussian distributions a result providing smoothed analysis estimates for condition numbers given as relativized distances to illposedness. We also introduce a notion of local analysis meant to capture the behavior of these condition numbers around a point.
more | pdf | html
None.
###### Tweets
mathNAb: Felipe Cucker, Teresa Krick : On local analysis https://t.co/zySYRg988m https://t.co/I0t52JI7wB
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #8. Local time-stepping for adaptive multiresolution using natural extension of Runge--Kutta methods
###### Müller Moreira Lopes, Margarete Oliveira Domingues, Kai Schneider, Odim Mendes
A space-time fully adaptive multiresolution method for evolutionary non-linear partial differential equations is presented introducing an improved local time-stepping method. The space discretisation is based on classical finite volumes, endowed with cell average multiresolution analysis for triggering the dynamical grid adaptation. The explicit time scheme features a natural extension of Runge--Kutta methods which allow local time-stepping while guaranteeing accuracy. The use of a compact Runge--Kutta formulation permits further memory reduction. The precision and computational efficiency of the scheme regarding CPU time and memory compression are assessed for problems in one, two and three space dimensions. As application Burgers equation, reaction-diffusion equations and the compressible Euler equations are considered. The numerical results illustrate the efficiency and superiority of the proposed local time-stepping method with respect to the reference computations.
more | pdf | html
None.
###### Tweets
mathNAb: Müller Moreira Lopes, Margarete Oliveira Domingues, Kai Schneider, Odim Mendes : Local time-stepping for adaptive multiresolution using natural extension of Runge--Kutta methods https://t.co/t9De7mvXKW https://t.co/p1LTktNBqA
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 4
Total Words: 0
Unqiue Words: 0

##### #9. Asymptotically exact a posteriori error estimates of eigenvalues by the Crouzeix-Raviart element and enriched Crouzeix-Raviart element
###### Jun Hu, Limin Ma
Two asymptotically exact a posteriori error estimates are proposed for eigenvalues by the nonconforming Crouzeix--Raviart and enriched Crouzeix-- Raviart elements. The main challenge in the design of such error estimators comes from the nonconformity of the finite element spaces used. Such nonconformity causes two difficulties, the first one is the construction of high accuracy gradient recovery algorithms, the second one is a computable high accuracy approximation of a consistency error term. The first difficulty was solved for both nonconforming elements in a previous paper. Two methods are proposed to solve the second difficulty in the present paper. In particular, this allows the use of high accuracy gradient recovery techniques. Further, a post-processing algorithm is designed by utilizing asymptotically exact a posteriori error estimators to construct the weights of a combination of two approximate eigenvalues. This algorithm requires to solve only one eigenvalue problem and admits high accuracy eigenvalue approximations...
more | pdf | html
###### Tweets
mathNAb: Jun Hu, Limin Ma : Asymptotically exact a posteriori error estimates of eigenvalues by the Crouzeix-Raviart element and enriched Crouzeix-Raviart element https://t.co/dPxvzZ85b4 https://t.co/8jKRCApWHe
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 9260
Unqiue Words: 1967

##### #10. Conservative difference schemes for one-dimensional flows of polytropic gas
###### Roman Kozlov
The paper considers one-dimensional flows of polytropic (calorically ideal) gas. These flows include three cases of gas dynamics: plain one-dimensional flows (one-dimensional space), radially symmetric flows in two-dimensional space and spherically symmetric flows in three-dimensional space. Starting with the difference schemes which have conservation laws of mass and energy (as well as conservation of momentum and the center of mass motion for the plain one-dimensional flows), we find difference schemes which also have additional conservation laws for the special values of the adiabatic exponent $\gamma = 1 + 1 /d$, where $d$ is the space dimension.
more | pdf | html
None.
###### Tweets
mathNAb: Roman Kozlov : Conservative difference schemes for one-dimensional flows of polytropic gas https://t.co/0wzyvFNvi0 https://t.co/SzhZCo3IJA
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 4010
Unqiue Words: 1071

Assert is a website where the best academic papers on arXiv (computer science, math, physics), bioRxiv (biology), BITSS (reproducibility), EarthArXiv (earth science), engrXiv (engineering), LawArXiv (law), PsyArXiv (psychology), SocArXiv (social science), and SportRxiv (sport research) bubble to the top each day.

Papers are scored (in real-time) based on how verifiable they are (as determined by their Github repos) and how interesting they are (based on Twitter).

To see top papers, follow us on twitter @assertpub_ (arXiv), @assert_pub (bioRxiv), and @assertpub_dev (everything else).

To see beautiful figures extracted from papers, follow us on Instagram.

Tracking 129,963 papers.

###### Search
Sort results based on if they are interesting or reproducible.
Interesting
Reproducible
Online
###### Stats
Tracking 129,963 papers.