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

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

Stargazers: 2

Subscribers: 1

Subscribers: 1

Forks: 2

Open Issues: 1

Open Issues: 1

None.

Sample Sizes : None.

Authors: 4

Total Words: 12737

Unqiue Words: 3397

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.

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.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

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.

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.

Sample Sizes : None.

Authors: 3

Total Words: 0

Unqiue Words: 0

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

Sample Sizes : None.

Authors: 2

Total Words: 7052

Unqiue Words: 1919

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
mathNAb:
Johannes Wallner : Geometric subdivision and multiscale transforms https://t.co/7m3rD2iI5K https://t.co/tnfkPUkzfy

None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 13133

Unqiue Words: 3129

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} \begin{equation}\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} \end{equation} into a discretized form in a
tighter and more accurate way without affecting the main results of the
article.

more |
pdf
| html
None.

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.

Sample Sizes : None.

Authors: 2

Total Words: 10093

Unqiue Words: 1737

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

Sample Sizes : None.

Authors: 2

Total Words: 13590

Unqiue Words: 2913

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.

mathNAb:
Cecilia Magherini : Weakly regular Sturm-Liouville problems: a corrected spectral matrix method https://t.co/Ykva7iDbO6 https://t.co/gJHwnQJ4WX

None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 6351

Unqiue Words: 1616

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.

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.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

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.

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.

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

Sort results based on if they are interesting or reproducible.

Interesting

Reproducible