Top 10 Arxiv Papers Today in Numerical Analysis


2.075 Mikeys
#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
Figures
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
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 2
Total Words: 7026
Unqiue Words: 1677

2.04 Mikeys
#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
Figures
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
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 4
Total Words: 0
Unqiue Words: 0

2.016 Mikeys
#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
Figures
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
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 2
Total Words: 10307
Unqiue Words: 2263

2.011 Mikeys
#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
Figures
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
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

2.011 Mikeys
#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
Figures
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
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

2.011 Mikeys
#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
Figures
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
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 5
Total Words: 0
Unqiue Words: 0

2.011 Mikeys
#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
Figures
None.
Tweets
mathNAb: Felipe Cucker, Teresa Krick : On local analysis https://t.co/zySYRg988m https://t.co/I0t52JI7wB
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

2.011 Mikeys
#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
Figures
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
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 4
Total Words: 0
Unqiue Words: 0

2.011 Mikeys
#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
Figures
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
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 2
Total Words: 9260
Unqiue Words: 1967

2.011 Mikeys
#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
Figures
None.
Tweets
mathNAb: Roman Kozlov : Conservative difference schemes for one-dimensional flows of polytropic gas https://t.co/0wzyvFNvi0 https://t.co/SzhZCo3IJA
Github
None.
Youtube
None.
Other stats
Sample Sizes : None.
Authors: 1
Total Words: 4010
Unqiue Words: 1071

About

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
Categories
All
Astrophysics
Cosmology and Nongalactic Astrophysics
Earth and Planetary Astrophysics
Astrophysics of Galaxies
High Energy Astrophysical Phenomena
Instrumentation and Methods for Astrophysics
Solar and Stellar Astrophysics
Condensed Matter
Disordered Systems and Neural Networks
Mesoscale and Nanoscale Physics
Materials Science
Other Condensed Matter
Quantum Gases
Soft Condensed Matter
Statistical Mechanics
Strongly Correlated Electrons
Superconductivity
Computer Science
Artificial Intelligence
Hardware Architecture
Computational Complexity
Computational Engineering, Finance, and Science
Computational Geometry
Computation and Language
Cryptography and Security
Computer Vision and Pattern Recognition
Computers and Society
Databases
Distributed, Parallel, and Cluster Computing
Digital Libraries
Discrete Mathematics
Data Structures and Algorithms
Emerging Technologies
Formal Languages and Automata Theory
General Literature
Graphics
Computer Science and Game Theory
Human-Computer Interaction
Information Retrieval
Information Theory
Machine Learning
Logic in Computer Science
Multiagent Systems
Multimedia
Mathematical Software
Numerical Analysis
Neural and Evolutionary Computing
Networking and Internet Architecture
Other Computer Science
Operating Systems
Performance
Programming Languages
Robotics
Symbolic Computation
Sound
Software Engineering
Social and Information Networks
Systems and Control
Economics
Econometrics
General Economics
Theoretical Economics
Electrical Engineering and Systems Science
Audio and Speech Processing
Image and Video Processing
Signal Processing
General Relativity and Quantum Cosmology
General Relativity and Quantum Cosmology
High Energy Physics - Experiment
High Energy Physics - Experiment
High Energy Physics - Lattice
High Energy Physics - Lattice
High Energy Physics - Phenomenology
High Energy Physics - Phenomenology
High Energy Physics - Theory
High Energy Physics - Theory
Mathematics
Commutative Algebra
Algebraic Geometry
Analysis of PDEs
Algebraic Topology
Classical Analysis and ODEs
Combinatorics
Category Theory
Complex Variables
Differential Geometry
Dynamical Systems
Functional Analysis
General Mathematics
General Topology
Group Theory
Geometric Topology
History and Overview
Information Theory
K-Theory and Homology
Logic
Metric Geometry
Mathematical Physics
Numerical Analysis
Number Theory
Operator Algebras
Optimization and Control
Probability
Quantum Algebra
Rings and Algebras
Representation Theory
Symplectic Geometry
Spectral Theory
Statistics Theory
Mathematical Physics
Mathematical Physics
Nonlinear Sciences
Adaptation and Self-Organizing Systems
Chaotic Dynamics
Cellular Automata and Lattice Gases
Pattern Formation and Solitons
Exactly Solvable and Integrable Systems
Nuclear Experiment
Nuclear Experiment
Nuclear Theory
Nuclear Theory
Physics
Accelerator Physics
Atmospheric and Oceanic Physics
Applied Physics
Atomic and Molecular Clusters
Atomic Physics
Biological Physics
Chemical Physics
Classical Physics
Computational Physics
Data Analysis, Statistics and Probability
Physics Education
Fluid Dynamics
General Physics
Geophysics
History and Philosophy of Physics
Instrumentation and Detectors
Medical Physics
Optics
Plasma Physics
Popular Physics
Physics and Society
Space Physics
Quantitative Biology
Biomolecules
Cell Behavior
Genomics
Molecular Networks
Neurons and Cognition
Other Quantitative Biology
Populations and Evolution
Quantitative Methods
Subcellular Processes
Tissues and Organs
Quantitative Finance
Computational Finance
Economics
General Finance
Mathematical Finance
Portfolio Management
Pricing of Securities
Risk Management
Statistical Finance
Trading and Market Microstructure
Quantum Physics
Quantum Physics
Statistics
Applications
Computation
Methodology
Machine Learning
Other Statistics
Statistics Theory
Feedback
Online
Stats
Tracking 129,963 papers.