### Top 10 Arxiv Papers Today in Numerical Analysis

##### #1. A numerical method for Hadamard finite-part integrals with an integral power singularity at an endpoint
###### Hidenori Ogata
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 0
Unqiue Words: 0

##### #2. Additive Schwarz solvers and preconditioners for the closest point method
###### Ian May, Ronald Haynes, Steven Ruuth
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 0
Unqiue Words: 0

##### #3. Sparse Representation of Gaussian Molecular Surface
###### Sheng Gui, Minxin Chen, Benzhuo Lu
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.
###### Tweets
mathNAb: Sheng Gui, Minxin Chen, Benzhuo Lu : Sparse Representation of Gaussian Molecular Surface https://t.co/oHQJ78B6bg https://t.co/Gtafzy8yEr
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 0
Unqiue Words: 0

##### #4. Non-smooth variational regularization for processing manifold-valued data
###### Martin Holler, Andreas Weinmann
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
###### Tweets
mathNAb: Martin Holler, Andreas Weinmann : Non-smooth variational regularization for processing manifold-valued data https://t.co/BtVaTiPzQc https://t.co/XRrrCLQ4fQ
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 18520
Unqiue Words: 4232

##### #5. Unconditional convergence for discretizations of dynamical optimal transport
###### Hugo Lavenant
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.
###### Tweets
mathNAb: Hugo Lavenant : Unconditional convergence for discretizations of dynamical optimal transport https://t.co/OkCQt1KZeh https://t.co/fdOnlJPfGC
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 0
Unqiue Words: 0

##### #6. Entropy stable reduced order modeling of nonlinear conservation laws
###### Jesse Chan
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.
###### Tweets
mathNAb: Jesse Chan : Entropy stable reduced order modeling of nonlinear conservation laws https://t.co/RoCcbgpScX https://t.co/v3m9X9S4P5
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 0
Unqiue Words: 0

##### #7. A fractional Traub method with $(2α+1)$th-order of convergence and its stability
###### Giro Candelario, Alicia Cordero, Juan R. Torregrosa
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 0
Unqiue Words: 0

##### #8. An Adaptive Parareal Algorithm
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #9. Fixed Point Analysis of Douglas-Rachford Splitting for Ptychography and Phase Retrieval
###### A. Fannjiang, Z. Zhang
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #10. Asymptotic preserving $P_N$ methods for haptotaxis equations
###### Gregor Corbin
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.
###### Tweets
mathNAb: Gregor Corbin : Asymptotic preserving $P_N$ methods for haptotaxis equations https://t.co/oECjzfxRDK https://t.co/t38NTTkKec
None.
None.
###### Other stats
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.

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