### Top 10 Arxiv Papers Today in Numerical Analysis

##### #1. Stability selection enables robust learning of partial differential equations from limited noisy data
###### Suryanarayana Maddu, Bevan L. Cheeseman, Ivo F. Sbalzarini, Christian L. Müller
We present a statistical learning framework for robust identification of partial differential equations from noisy spatiotemporal data. Extending previous sparse regression approaches for inferring PDE models from simulated data, we address key issues that have thus far limited the application of these methods to noisy experimental data, namely their robustness against noise and the need for manual parameter tuning. We address both points by proposing a stability-based model selection scheme to determine the level of regularization required for reproducible recovery of the underlying PDE. This avoids manual parameter tuning and provides a principled way to improve the method's robustness against noise in the data. Our stability selection approach, termed PDE-STRIDE, can be combined with any sparsity-promoting penalized regression model and provides an interpretable criterion for model component importance. We show that in particular the combination of stability selection with the iterative hard-thresholding algorithm from...
more | pdf | html
###### Tweets
arxivml: "Stability selection enables robust learning of partial differential equations from limited noisy data", Suryanaray… https://t.co/OFDpz6MeM4
mathNAb: Suryanarayana Maddu, Bevan L. Cheeseman, Ivo F. Sbalzarini, Christian L. Müller : Stability selection enables robust learning of partial differential equations from limited noisy data https://t.co/xuNa6PpdJH https://t.co/b80ksanPR5
Memoirs: Stability selection enables robust learning of partial differential equations from limited noisy data. https://t.co/sdeNJQ8LQv
MOSAICgroup1: RT @surimk92: @MOSAICgroup1 new paper on PDE learning from limited noisy data is out on arXiv https://t.co/NjZA98jN9R. Exciting work at the…
MOSAICgroup1: RT @surimk92: @MOSAICgroup1 paper https://t.co/NjZA98jN9R on automated PDE learning at a glance. https://t.co/F1hfrKlMgP
microbionaut: RT @surimk92: @MOSAICgroup1 new paper on PDE learning from limited noisy data is out on arXiv https://t.co/NjZA98jN9R. Exciting work at the…
FlatironCCM: RT @surimk92: @MOSAICgroup1 new paper on PDE learning from limited noisy data is out on arXiv https://t.co/NjZA98jN9R. Exciting work at the…
BevanCheeseman: RT @surimk92: @MOSAICgroup1 new paper on PDE learning from limited noisy data is out on arXiv https://t.co/NjZA98jN9R. Exciting work at the…
###### Github

This repository contains codes and data-sets for the PDE inference from limited spatio-temporal data

Repository: PDE-STRIDE
User: SuryanarayanaMK
Language: Jupyter Notebook
Stargazers: 2
Subscribers: 1
Forks: 2
Open Issues: 1
None.
###### Other stats
Sample Sizes : None.
Authors: 4
Total Words: 12737
Unqiue Words: 3397

##### #2. Discontinuous Galerkin methods for short pulse type equations via hodograph transformations
###### Qian Zhang, Yinhua Xia
In the present paper, we consider the discontinuous Galerkin (DG) methods for solving short pulse (SP) type equations. The short pulse equation has been shown to be completely integrable, which admits the loop-soliton, cuspon-soliton solutions as well as smooth-soliton solutions. Through hodograph transformations, these nonclassical solutions can be profiled as the smooth solutions of the coupled dispersionless (CD) system or the sine-Gordon equation. Thus, DG methods can be developed for the CD system or the sine-Gordon equation to simulate the loop-soliton or cuspon-soliton solutions of the SP equation. The conservativeness or dissipation of the Hamiltonian or momentum for the semi-discrete DG schemes can be proved. Also we modify the above DG schemes and obtain an integration DG scheme. Theoretically the a-priori error estimates have been provided for the momentum conserved DG scheme and the integration DG scheme. We also propose the DG scheme and the integration DG scheme for the sine-Gordon equation, in case the SP equation...
more | pdf | html
None.
###### Tweets
mathNAb: Qian Zhang, Yinhua Xia : Discontinuous Galerkin methods for short pulse type equations via hodograph transformations https://t.co/Uq960k4H2u https://t.co/6QCaTx8Ugy
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #3. Extremum Global Sensitivity Analysis with Least Squares Polynomials and their Ridges
###### Chun Yui Wong, Pranay Seshadri, Geoffrey Parks
Global sensitivity analysis is a powerful set of ideas and computational heuristics for understanding the importance and interplay between different parameters in a data-set. Such a data-set is characterized by a set of vector-valued input parameters and a set of scalar-valued output quantities of interest, where we typically assume that the inputs are independent and information on their joint density can be obtained. Alternatively, if the inputs are correlated, one requires information on the marginals and their correlations. In either case, if the output quantities of interest are smooth and continuous, polynomial least squares approximations can be used to extract Sobol' indices. In this paper, we build on these previously well-known ideas by studying two different aspects of this paradigm. First, we study whether sensitivity indices can be computed efficiently if one leverages a polynomial ridge approximation---a polynomial least squares fit over a subspace. We discuss a recipe that utilizes this special dependence structure...
more | pdf | html
None.
###### Tweets
mathNAb: Chun Yui Wong, Pranay Seshadri, Geoffrey Parks : Extremum Global Sensitivity Analysis with Least Squares Polynomials and their Ridges https://t.co/H2oloLE8WW https://t.co/IURKY24lEv
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 0
Unqiue Words: 0

##### #4. On a convergent DSA preconditioned source iteration for a DGFEM method for radiative transfer
###### Olena Palii, Matthias Schlottbom
We consider the numerical approximation of the radiative transfer equation using discontinuous angular and continuous spatial approximations for the even parts of the solution. The even-parity equations are solved using a diffusion synthetic accelerated source iteration. We provide a convergence analysis for the infinite-dimensional iteration as well as for its discretized counterpart. The diffusion correction is computed by a subspace correction, which leads to a convergence behavior that is robust with respect to the discretization. The proven theoretical contraction rate deteriorates for scattering dominated problems. We show numerically that the preconditioned iteration is in practice robust in the diffusion limit. Moreover, computations for the lattice problem indicate that the presented discretization does not suffer from the ray effect. The theoretical methodology is presented for plane-parallel geometries with isotropic scattering, but the approach and proofs generalize to multi-dimensional problems and more general...
more | pdf | html
###### Tweets
mathNAb: Olena Palii, Matthias Schlottbom : On a convergent DSA preconditioned source iteration for a DGFEM method for radiative transfer https://t.co/alfPBiJJ1X https://t.co/WxAFZwj4do
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 7052
Unqiue Words: 1919

##### #5. Geometric subdivision and multiscale transforms
###### Johannes Wallner
Any procedure applied to data, and any quantity derived from data, is required to respect the nature and symmetries of the data. This axiom applies to refinement procedures and multiresolution transforms as well as to more basic operations like averages. This chapter discusses different kinds of geometric structures like metric spaces, Riemannian manifolds, and groups, and in what way we can make elementary operations geometrically meaningful. A nice example of this is the Riemannian metric naturally associated with the space of positive definite matrices and the intrinsic operations on positive definite matrices derived from it. We disucss averages first and then proceed to refinement operations (subdivision) and multiscale transforms. In particular, we report on the current knowledge as regards convergence and smoothness.
more | pdf | html
###### Tweets
mathNAb: Johannes Wallner : Geometric subdivision and multiscale transforms https://t.co/7m3rD2iI5K https://t.co/tnfkPUkzfy
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 13133
Unqiue Words: 3129

##### #6. Corrections on A numerical method for solving nonlinear Volterra--Fredholm integral equations
###### Ngo Thanh Binh, Khuat Van Ninh
Some corrections are made in our article, which was published in Appl. Anal. Optim. Vol. 3 (2019), No. 1, 103--127. These corrections are intended to transform the equation \eqref{eq:1.1} $$\label{eq:1.1} x(t) + \int\limits_a^t {K_1(t,s,x(s)) ds} + \int\limits_a^b {K_2(t,s,x(s)) ds} = g(t),\;\,a \le t \le b \tag{1.1}$$ into a discretized form in a tighter and more accurate way without affecting the main results of the article.
more | pdf | html
None.
###### Tweets
mathNAb: Ngo Thanh Binh, Khuat Van Ninh : Corrections on A numerical method for solving nonlinear Volterra--Fredholm integral equations https://t.co/f0OlEPQNKJ https://t.co/vYd6uJZwlI
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 10093
Unqiue Words: 1737

##### #7. Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations
###### Kailiang Wu, Chi-Wang Shu
This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not possess an entropy pair. To address this issue, a symmetrizable RMHD system, which admits a convex entropy pair, is proposed by adding a source term into the equations. Arbitrarily high-order accurate entropy stable finite difference schemes are then developed on Cartesian meshes based on the symmetrizable RMHD system. The crucial ingredients of these schemes include (i) affordable explicit entropy conservative fluxes which are technically derived through carefully selected parameter variables, (ii) a special high-order discretization of the source term in the symmetrizable RMHD system, and (iii) suitable high-order dissipative operators based on essentially non-oscillatory reconstruction to ensure the entropy stability. Several benchmark numerical tests demonstrate the accuracy and robustness of the...
more | pdf | html
###### Tweets
mathNAb: Kailiang Wu, Chi-Wang Shu : Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations https://t.co/pa8yn9cTpB https://t.co/F9TeZ3AiOX
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 13590
Unqiue Words: 2913

##### #8. Weakly regular Sturm-Liouville problems: a corrected spectral matrix method
###### Cecilia Magherini
In this paper, we consider weakly regular Sturm-Liouville eigenproblems with unbounded potential at both endpoints of the domain. We propose a Galerkin spectral matrix method for its solution and we study the error in the eigenvalue approximations it provides. The result of the convergence analysis is then used to derive a low-cost and very effective formula for the computation of corrected numerical eigenvalues. Finally, we present and discuss the results of several numerical experiments which confirm the validity of the approach.
more | pdf | html
None.
###### Tweets
mathNAb: Cecilia Magherini : Weakly regular Sturm-Liouville problems: a corrected spectral matrix method https://t.co/Ykva7iDbO6 https://t.co/gJHwnQJ4WX
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 6351
Unqiue Words: 1616

##### #9. Richardson extrapolation for the discrete iterated modified projection solution
###### Gobinda Rakshit, Rekha P. Kulkarni
Approximate solutions of Urysohn integral equations using projection methods involve integrals which need to be evaluated using a numerical quadrature formula. It gives rise to the discrete versions of the projection methods. For $r \geq 1,$ a space of piece-wise polynomials of degree $\leq r - 1$ with respect to an uniform partition is chosen to be the approximating space and the projection is chosen to be the interpolatory projection at $r$ Gauss points. Asymptotic expansion for the iterated modified projection solution is available in literature. In this paper, we obtain an asymptotic expansion for the discrete iterated modified projection solution and use Richardson extrapolation to improve the order of convergence. Our results indicate a choice of a numerical quadrature which preserves the order of convergence in the continuous case.
more | pdf | html
None.
###### Tweets
mathNAb: Gobinda Rakshit, Rekha P. Kulkarni : Richardson extrapolation for the discrete iterated modified projection solution https://t.co/jd2DfMszT9 https://t.co/o5J1bGgsBP
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #10. A filtered Boris algorithm for charged-particle dynamics in a strong magnetic field
###### Ernst Hairer, Christian Lubich, Bin Wang
A modification of the standard Boris algorithm, called filtered Boris algorithm, is proposed for the numerical integration of the equations of motion of charged particles in a strong non-uniform magnetic field in the asymptotic scaling known as maximal ordering. With an appropriate choice of filters, second-order error bounds in the position and in the parallel velocity, and first-order error bounds in the normal velocity are obtained with respect to the scaling parameter. The proof compares the modulated Fourier expansions of the exact and the numerical solutions. Numerical experiments illustrate the error behaviour of the filtered Boris algorithm.
more | pdf | html
None.
###### Tweets
mathNAb: Ernst Hairer, Christian Lubich, Bin Wang : A filtered Boris algorithm for charged-particle dynamics in a strong magnetic field https://t.co/cFegt3MIDG https://t.co/Ivk8ojlO3K
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 3
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 160,434 papers.

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