In this paper, we propose a numerical method for computing Hadamard
finite-part integrals with an integral-power singularity at an endpoint, the
part of the divergent integral which is finite as a limiting procedure. In the
proposed method, we express the desired finite-part integral using a complex
loop integral, and obtain the finite-part integral by evaluating the complex
integral by the trapezoidal rule. Theoretical error estimate and some numerical
examples show the effectiveness of the proposed method.

more |
pdf
| html
None.

mathNAb:
Hidenori Ogata : A numerical method for Hadamard finite-part integrals with an integral power singularity at an endpoint https://t.co/dN121kSNuj https://t.co/bixpuao6xs

None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 0

Unqiue Words: 0

The discretization of surface intrinsic elliptic partial differential
equations (PDEs) pose interesting challenges not seen in flat space. The
discretization of these PDEs typically proceeds by either parametrizing the
surface, triangulating the surface, or embedding the surface in a higher
dimensional flat space. The closest point method (CPM) is an embedding method
that represents surfaces using a function that maps points in the embedding
space to their closest points on the surface. In the CPM, this mapping also
serves as an extension operator that brings surface intrinsic data onto the
embedding space, allowing PDEs to be numerically approximated by standard
methods in a narrow tubular neighborhood of the surface. We focus here on
numerically approximating the positive Helmholtz equation,
$\left(c-\Delta_{\mathcal{S}}\right)u=f,~c\in\mathbb{R}^+$ by the CPM. This
yields a large, sparse, and nonsymmetric system to solve. Herein, we develop
restricted additive Schwarz (RAS) and optimized restricted additive Schwarz
(ORAS)...

more |
pdf
| html
None.

mathNAb:
Ian May, Ronald Haynes, Steven Ruuth : Additive Schwarz solvers and preconditioners for the closest point method https://t.co/dUtQmIh3lP https://t.co/vYXQCicu79

None.

None.

Sample Sizes : None.

Authors: 3

Total Words: 0

Unqiue Words: 0

In this paper, we propose a model and algorithm for sparse representing
Gaussian molecular surface. The original Gaussian molecular surface is
approximated by a relatively small number of radial basis functions (RBFs) with
rotational ellipsoid feature. The sparsity of the RBF representation is
achieved by solving a nonlinear $L_1$ optimization problem. Experimental
results demonstrate that the original Gaussian molecular surface is able to be
represented with good accuracy by much fewer RBFs using our $L_1$ model and
algorithm. The sparse representation of Gaussian molecular surface is useful in
various applications, such as molecular structure alignment, calculating
molecular areas and volumes, and the method in principle can be applied to
sparse representation of general shapes and coarse-grained molecular modeling.

more |
pdf
| html
None.

mathNAb:
Sheng Gui, Minxin Chen, Benzhuo Lu : Sparse Representation of Gaussian Molecular Surface https://t.co/oHQJ78B6bg https://t.co/Gtafzy8yEr

None.

None.

Sample Sizes : None.

Authors: 3

Total Words: 0

Unqiue Words: 0

Many methods for processing scalar and vector valued images, volumes and
other data in the context of inverse problems are based on variational
formulations. Such formulations require appropriate regularization functionals
that model expected properties of the object to reconstruct. Prominent examples
of regularization functionals in a vector-space context are the total variation
(TV) and the Mumford-Shah functional, as well as higher-order schemes such as
total generalized variation models. Driven by applications where the signals or
data live in nonlinear manifolds, there has been quite some interest in
developing analogous methods for nonlinear, manifold-valued data recently. In
this chapter, we consider various variational regularization methods for
manifold-valued data. In particular, we consider TV minimization as well as
higher order models such as total generalized variation (TGV). Also, we discuss
(discrete) Mumford-Shah models and related methods for piecewise constant data.
We develop discrete energies for denoising and...

more |
pdf
| html
mathNAb:
Martin Holler, Andreas Weinmann : Non-smooth variational regularization for processing manifold-valued data https://t.co/BtVaTiPzQc https://t.co/XRrrCLQ4fQ

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 18520

Unqiue Words: 4232

The dynamical formulation of optimal transport, also known as Benamou-Brenier
formulation or Computational Fluid Dynamic formulation, amounts to write the
optimal transport problem as the optimization of a convex functional under a
PDE constraint, and can handle \emph{a priori} a vast class of cost functions
and geometries. Several disretizations of this problem have been proposed,
leading to computations on flat spaces as well as Riemannian manifolds, with
extensions to mean field games and gradient flows in the Wasserstein space.
In this article, we provide a framework which guarantees convergence under
mesh refinement of the solutions of the space-time discretized problems to the
one of the infinite-dimensional one for quadratic optimal transport. The
convergence holds without condition on the ratio between spatial and temporal
step sizes, and can handle arbitrary positive measures as input, while the
underlying space can be a Riemannian manifold. Both the finite volume
discretization proposed by Gladbach, Kopfer and Maas, as...

more |
pdf
| html
None.

mathNAb:
Hugo Lavenant : Unconditional convergence for discretizations of dynamical optimal transport https://t.co/OkCQt1KZeh https://t.co/fdOnlJPfGC

None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 0

Unqiue Words: 0

Reduced order models of nonlinear conservation laws in fluid dynamics do not
typically inherit stability properties of the full order model. We introduce
projection-based hyper-reduced models of nonlinear conservation laws which are
globally conservative and inherit a semi-discrete entropy inequality
independently of the choice of basis and choice of parameters.

more |
pdf
| html
None.

mathNAb:
Jesse Chan : Entropy stable reduced order modeling of nonlinear conservation laws https://t.co/RoCcbgpScX https://t.co/v3m9X9S4P5

None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 0

Unqiue Words: 0

Some fractional Newton methods have been proposed in order to find roots of
nonlinear equations using fractional derivatives. In this paper we introduce a
fractional Newton method with order $\alpha+1$ and compare with another
fractional Newton method with order $2\alpha$. We also introduce a fractional
Traub method with order $2\alpha+1$ and compare with its first step (fractional
Newton method with order $\alpha+1$). Some tests and analysis of the dependence
on the initial estimations are made for each case.

more |
pdf
| html
None.

mathNAb:
Giro Candelario, Alicia Cordero, Juan R. Torregrosa : A fractional Traub method with $(2α+1)$th-order of convergence and its stability https://t.co/xMNtMrqWns https://t.co/0C87lisHiQ

None.

None.

Sample Sizes : None.

Authors: 3

Total Words: 0

Unqiue Words: 0

In this paper, we consider the problem of accelerating the numerical
simulation of time dependent problems by time domain decomposition. The
available algorithms enabling such decompositions present severe efficiency
limitations and are an obstacle for the solution of large scale and high
dimensional problems. Our main contribution is the improvement of the parallel
efficiency of the parareal in time method. The parareal method is based on
combining predictions made by a numerically inexpensive solver (with coarse
physics and/or coarse resolution) with corrections coming from an expensive
solver (with high-fidelity physics and high resolution). At convergence, the
parareal algorithm provides a solution that has the fine solver's high-fidelity
physics and high resolution In the classical version of parareal, the fine
solver has a fixed high accuracy which is the major obstacle to achieve a
competitive parallel efficiency. In this paper, we develop an adaptive variant
of the algorithm that overcomes this obstacle. Thanks to this,...

more |
pdf
| html
None.

colinjcotter:
Potentially significant Parareal paper on the ArXiV today. I haven’t read it yet but abstract says they address efficiency by adapting the fine propagator.
https://t.co/U4hpXZuPOB

mathNAb:
Y. Maday, O. Mula : An Adaptive Parareal Algorithm https://t.co/iOGLY3GnZN https://t.co/c09p6DlsO6

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

Douglas-Rachford Splitting (DRS) methods based on the proximal point
algorithms for the Poisson and Gaussian log-likelihood functions are proposed
for ptychography and phase retrieval. Fixed point analysis shows that the DRS
iterated sequences are always bounded explicitly in terms of the step size and
that the fixed points are linearly stable if and only if the fixed points are
regular solutions. This alleviates two major drawbacks of the classical
Douglas-Rachford (CDR) algorithm: slow convergence when the feasibility problem
is consistent and divergent behavior when the feasibility problem is
inconsistent. Moreover, the fixed point analysis decisively leads to the
selection of an optimal step size which in turn renders the Gaussian DRS method
in a particularly simple form with no tuning parameter (Averaged
Projection-Reflection). When applied to the challenging problem of blind
ptychography, which seeks to recover both the object and the probe
simultaneously, the DRS methods converge geometrically and globally when
properly initialized.

more |
pdf
| html
None.

mathNAb:
A. Fannjiang, Z. Zhang : Fixed Point Analysis of Douglas-Rachford Splitting for Ptychography and Phase Retrieval https://t.co/1CAFbG6V0f https://t.co/XeImOotMDP

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

The so-called haptotaxis equation is a special class of transport equation
that arises from models of biological cell movement along tissue fibers. This
equation has an anisotropic advection-diffusion equation as its macroscopic
limit. An up to second-order accurate asymptotic preserving method is developed
for the haptotaxis equation in space dimension up to three. For this the
micro-macro decomposition proposed by Lemou and Mieussens is generalized in the
context of finite-volume schemes on staggered grids. The spurious modes that
arise from this discretization can be eliminated by combining flux evaluations
from different points in the right way. The velocity space is discretized by an
arbitrary-order linear moment system ($P_N$).

more |
pdf
| html
None.

mathNAb:
Gregor Corbin : Asymptotic preserving $P_N$ methods for haptotaxis equations https://t.co/oECjzfxRDK https://t.co/t38NTTkKec

None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 0

Unqiue Words: 0

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 192,915 papers.*

Sort results based on if they are interesting or reproducible.

Interesting

Reproducible