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.

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.

Sample Sizes : None.

Authors: 2

Total Words: 7026

Unqiue Words: 1677

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.

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.

Sample Sizes : None.

Authors: 4

Total Words: 0

Unqiue Words: 0

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.

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

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

Stargazers: 0

Subscribers: 1

Subscribers: 1

Forks: 0

Open Issues: 0

Open Issues: 0

None.

Sample Sizes : None.

Authors: 2

Total Words: 10307

Unqiue Words: 2263

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.

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.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

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.

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.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

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.

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.

Sample Sizes : None.

Authors: 5

Total Words: 0

Unqiue Words: 0

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.

mathNAb:
Felipe Cucker, Teresa Krick : On local analysis https://t.co/zySYRg988m https://t.co/I0t52JI7wB

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

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.

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.

Sample Sizes : None.

Authors: 4

Total Words: 0

Unqiue Words: 0

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
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.

Sample Sizes : None.

Authors: 2

Total Words: 9260

Unqiue Words: 1967

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.

mathNAb:
Roman Kozlov : Conservative difference schemes for one-dimensional flows of polytropic gas https://t.co/0wzyvFNvi0 https://t.co/SzhZCo3IJA

None.

None.

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.*

Sort results based on if they are interesting or reproducible.

Interesting

Reproducible