-
A method for determining the parameters in a rheological model for viscoelastic materials by minimizing Tikhonov functionals arXiv.cs.NA Pub Date : 2021-02-26 Rebecca Rothermel; Wladimir Panfilenko; Prateek Sharma; Anne Wald; Thomas Schuster; Anne Jung; Stefan Diebels
Mathematical models describing the behavior of viscoelastic materials are often based on evolution equations that measure the change in stress depending on its material parameters such as stiffness, viscosity or relaxation time. In this article, we introduce a Maxwell-based rheological model, define the associated forward operator and the inverse problem in order to determine the number of Maxwell
-
On Rational Krylov and Reduced Basis Methods for Fractional Diffusion arXiv.cs.NA Pub Date : 2021-02-26 Tobias Danczul; Clemens Hofreither
We establish an equivalence between two classes of methods for solving fractional diffusion problems, namely, Reduced Basis Methods (RBM) and Rational Krylov Methods (RKM). In particular, we demonstrate that several recently proposed RBMs for fractional diffusion can be interpreted as RKMs. This changed point of view allows us to give convergence proofs for some methods where none were previously available
-
Multi-fidelity regression using artificial neural networks: efficient approximation of parameter-dependent output quantities arXiv.cs.NA Pub Date : 2021-02-26 Mengwu Guo; Andrea Manzoni; Maurice Amendt; Paolo Conti; Jan S. Hesthaven
Highly accurate numerical or physical experiments are often time-consuming or expensive to obtain. When time or budget restrictions prohibit the generation of additional data, the amount of available samples may be too limited to provide satisfactory model results. Multi-fidelity methods deal with such problems by incorporating information from other sources, which are ideally well-correlated with
-
Data-driven modeling of linear dynamical systems with quadratic output in the AAA framework arXiv.cs.NA Pub Date : 2021-02-25 Ion Victor Gosea; Serkan Gugercin
We extend the AAA (Adaptive-Antoulas-Anderson) algorithm to develop a data-driven modeling framework for linear systems with quadratic output (LQO). Such systems are characterized by two transfer functions: one corresponding to the linear part of the output and another one to the quadratic part. We first establish the joint barycentric representations and the interpolation theory for the two transfer
-
Arbitrarily High-order Maximum Bound Preserving Schemes with Cut-off Postprocessing for Allen-Cahn Equations arXiv.cs.NA Pub Date : 2021-02-26 Jiang Yang; Zhaoming Yuan; Zhi Zhou
We develop and analyze a class of maximum bound preserving schemes for approximately solving Allen--Cahn equations. We apply a $k$th-order single-step scheme in time (where the nonlinear term is linearized by multi-step extrapolation), and a lumped mass finite element method in space with piecewise $r$th-order polynomials and Gauss--Lobatto quadrature. At each time level, a cut-off post-processing
-
Power series expansion neural network arXiv.cs.NA Pub Date : 2021-02-25 Qipin Chen; Wenrui Hao; Juncai He
In this paper, we develop a new neural network family based on power series expansion, which is proved to achieve a better approximation accuracy comparing to existing neural networks. This new set of neural networks can improve the expressive power while preserving comparable computational cost by increasing the degree of the network instead of increasing the depth or width. Numerical results have
-
Variance Reduction via Primal-Dual Accelerated Dual Averaging for Nonsmooth Convex Finite-Sums arXiv.cs.NA Pub Date : 2021-02-26 Chaobing Song; Stephen J. Wright; Jelena Diakonikolas
We study structured nonsmooth convex finite-sum optimization that appears widely in machine learning applications, including support vector machines and least absolute deviation. For the primal-dual formulation of this problem, we propose a novel algorithm called \emph{Variance Reduction via Primal-Dual Accelerated Dual Averaging (\vrpda)}. In the nonsmooth and general convex setting, \vrpda~has the
-
Contrast-independent partially explicit time discretizations for multiscale wave problems arXiv.cs.NA Pub Date : 2021-02-25 Eric T. Chung; Yalchin Efendiev; Wing Tat Leung; Petr N. Vabishchevich
In this work, we design and investigate contrast-independent partially explicit time discretizations for wave equations in heterogeneous high-contrast media. We consider multiscale problems, where the spatial heterogeneities are at subgrid level and are not resolved. In our previous work, we have introduced contrast-independent partially explicit time discretizations and applied to parabolic equations
-
Krylov solvability under perturbations of abstract inverse linear problems arXiv.cs.NA Pub Date : 2021-02-26 Noe Angelo Caruso; Alessandro Michelangeli
When a solution to an abstract inverse linear problem on Hilbert space is approximable by finite linear combinations of vectors from the cyclic subspace associated with the datum and with the linear operator of the problem, the solution is said to be a Krylov solution, i.e., it belongs to the Krylov subspace of the problem. Krylov solvability of the inverse problem allows for solution approximations
-
Sparse Approximations with Interior Point Methods arXiv.cs.NA Pub Date : 2021-02-26 Valentina De Simone; Daniela di Serafino; Jacek Gondzio; Spyridon Pougkakiotis; Marco Viola
Large-scale optimization problems that seek sparse solutions have become ubiquitous. They are routinely solved with various specialized first-order methods. Although such methods are often fast, they usually struggle with not-so-well conditioned problems. In this paper, specialized variants of an interior point-proximal method of multipliers are proposed and analyzed for problems of this class. Computational
-
Approximation of Stochastic Volterra Equations with kernels of completely monotone type arXiv.cs.NA Pub Date : 2021-02-26 Aurélien Alfonsi; Ahmed Kebaier
In this work, we develop a multi-factor approximation for Stochastic Volterra Equations with Lipschitz coefficients and kernels of completely monotone type that may be singular. Our approach consists in truncating and then discretizing the integral defining the kernel, which corresponds to a classical Stochastic Differential Equation. We prove strong convergence results for this approximation. For
-
Derivative-Free Multiobjective Trust Region Descent MethodUsing Radial Basis Function Surrogate Models arXiv.cs.NA Pub Date : 2021-02-26 Manuel Berkemeier; Sebastian Peitz
We present a flexible trust region descend algorithm for unconstrained and convexly constrained multiobjective optimization problems. It is targeted at heterogeneous and expensive problems, i.e., problems that have at least one objective function that is computationally expensive. The method is derivative-free in the sense that neither need derivative information be available for the expensive objectives
-
A Fast Proximal Gradient Method and Convergence Analysis for Dynamic Mean Field Planning arXiv.cs.NA Pub Date : 2021-02-26 Jiajia Yu; Rongjie Lai; Wuchen Li; Stanley Osher
In this paper, we propose an efficient and flexible algorithm to solve dynamic mean-field planning problems based on an accelerated proximal gradient method. Besides an easy-to-implement gradient descent step in this algorithm, a crucial projection step becomes solving an elliptic equation whose solution can be obtained by conventional methods efficiently. By induction on iterations used in the algorithm
-
Data-Free Likelihood-Informed Dimension Reduction of Bayesian Inverse Problems arXiv.cs.NA Pub Date : 2021-02-26 Tiangang Cui; Olivier Zahm
Identifying a low-dimensional informed parameter subspace offers a viable path to alleviating the dimensionality challenge in the sampled-based solution to large-scale Bayesian inverse problems. This paper introduces a novel gradient-based dimension reduction method in which the informed subspace does not depend on the data. This permits an online-offline computational strategy where the expensive
-
Recovery of regular ridge functions on the ball arXiv.cs.NA Pub Date : 2021-02-25 Tatyana Zaitseva; Yuri Malykhin; Konstantin Ryutin
We consider the problem of the uniform (in $L_\infty$) recovery of ridge functions $f(x)=\varphi(\langle a,x\rangle)$, $x\in B_2^n$, using noisy evaluations $y_1\approx f(x^1),\ldots,y_N\approx f(x^N)$. It is known that for classes of functions $\varphi$ of finite smoothness the problem suffers from the curse of dimensionality: in order to provide good accuracy for the recovery it is necessary to make
-
Solitary water wave interactions for the Forced Korteweg-de Vries equation arXiv.cs.NA Pub Date : 2021-02-25 M. V. Flamarion; R. Ribeiro-Jr
The aim of this work is to study solitary water wave interactions between two topographic obstacles for the forced Korteweg-de Vries equation (fKdV). Focusing on the details of the interactions, we identify regimes in which solitary wave interactions maintain two well separated crests and regimes where the number of local maxima varies according to the laws $2\rightarrow 1\rightarrow 2\rightarrow 1\rightarrow
-
Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers arXiv.cs.NA Pub Date : 2021-02-25 P. V. F. Edelmann; L. Horst; J. P. Berberich; R. Andrassy; J. Higl; C. Klingenberg; F. K. Roepke
Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are but tiny perturbations on top of a close balance between gravity and pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can
-
Semi-implicit Hybrid Discrete $\left(\text{H}^T_N\right)$ Approximation of Thermal Radiative Transfer arXiv.cs.NA Pub Date : 2021-02-25 Ryan G. McClarren; James A. Rossmanith; Minwoo Shin
The thermal radiative transfer (TRT) equations form a system that describes the propagation and collisional interactions of photons. Computing accurate and efficient numerical solutions to TRT is challenging for several reasons, the first of which is that TRT is defined on a high-dimensional phase space. In order to reduce the dimensionality, classical approaches such as the P$_N$ (spherical harmonics)
-
Lie Group integrators for mechanical systems arXiv.cs.NA Pub Date : 2021-02-25 Elena Celledoni; Ergys Çokaj; Andrea Leone; Davide Murari; Brynjulf Owren
Since they were introduced in the 1990s, Lie group integrators have become a method of choice in many application areas. These include multibody dynamics, shape analysis, data science, image registration and biophysical simulations. Two important classes of intrinsic Lie group integrators are the Runge--Kutta--Munthe--Kaas methods and the commutator free Lie group integrators. We give a short introduction
-
A continuation method for computing the multilinear Pagerank arXiv.cs.NA Pub Date : 2021-02-25 Alberto Bucci; Federico Poloni
The multilinear Pagerank model [Gleich, Lim and Yu, 2015] is a tensor-based generalization of the Pagerank model. Its computation requires solving a system of polynomial equations that contains a parameter $\alpha \in [0,1)$. For $\alpha \approx 1$, this computation remains a challenging problem, especially since the solution may be non-unique. Extrapolation strategies that start from smaller values
-
ISALT: Inference-based schemes adaptive to large time-stepping for locally Lipschitz ergodic systems arXiv.cs.NA Pub Date : 2021-02-25 Xingjie Li; Fei Lu; Felix X. -F. Ye
Efficient simulation of SDEs is essential in many applications, particularly for ergodic systems that demand efficient simulation of both short-time dynamics and large-time statistics. However, locally Lipschitz SDEs often require special treatments such as implicit schemes with small time-steps to accurately simulate the ergodic measure. We introduce a framework to construct inference-based schemes
-
Approximation of Manifold-valued Functions arXiv.cs.NA Pub Date : 2021-02-24 Ralf Hielscher; Laura Lippert
We consider the approximation of manifold-valued functions by embedding the manifold into a higher dimensional space, applying a vector-valued approximation operator and projecting the resulting vector back to the manifold. It is well known that the approximation error for manifold-valued functions is close to the approximation error for vector-valued functions. This is not true anymore if we consider
-
Decompositions of high-frequency Helmholtz solutions via functional calculus, and application to the finite element method arXiv.cs.NA Pub Date : 2021-02-25 David Lafontaine; Euan A. Spence; Jared Wunsch
Over the last ten years, results from [Melenk-Sauter, 2010], [Melenk-Sauter, 2011], [Esterhazy-Melenk, 2012], and [Melenk-Parsania-Sauter, 2013] decomposing high-frequency Helmholtz solutions into "low"- and "high"-frequency components have had a large impact in the numerical analysis of the Helmholtz equation. These results have been proved for the constant-coefficient Helmholtz equation in either
-
Stein Variational Gradient Descent: many-particle and long-time asymptotics arXiv.cs.NA Pub Date : 2021-02-25 Nikolas Nüsken; D. R. Michiel Renger
Stein variational gradient descent (SVGD) refers to a class of methods for Bayesian inference based on interacting particle systems. In this paper, we consider the originally proposed deterministic dynamics as well as a stochastic variant, each of which represent one of the two main paradigms in Bayesian computational statistics: variational inference and Markov chain Monte Carlo. As it turns out,
-
BGK models for inert mixtures: comparison and applications arXiv.cs.NA Pub Date : 2021-02-25 Sebastiano Boscarino; Seung Yeon Cho; Maria Groppi; Giovanni Russo
Consistent BGK models for inert mixtures are compared, first in their kinetic behavior and then versus the hydrodynamic limits that can be derived in different collision-dominated regimes. The comparison is carried out both analytically and numerically, for the latter using an asymptotic preserving semi-Lagrangian scheme for the BGK models. Application to the plane shock wave in a binary mixture of
-
Survival of the Fittest: Numerical Modeling of Supernova 2014C arXiv.cs.NA Pub Date : 2021-02-24 Felipe Vargas; Fabio De Colle; Daniel Brethauer; Raffaella Margutti; Cristian G. Bernal
Initially classified as a supernova (SN) type Ib, $\sim$ 100 days after the explosion SN\,2014C made a transition to a SN type II, presenting a gradual increase in the H${\alpha}$ emission. This has been interpreted as evidence of interaction between the supernova shock wave and a massive shell previously ejected from the progenitor star. In this paper, we present numerical simulations of the propagation
-
Energy-consistent finite difference schemes for compressible hydrodynamics and magnetohydrodynamics using nonlinear filtering arXiv.cs.NA Pub Date : 2021-02-24 Haruhisa Iijima
In this paper, an energy-consistent finite difference scheme for the compressible hydrodynamic and magnetohydrodynamic (MHD) equations is introduced. For the compressible magnetohydrodynamics, an energy-consistent finite difference formulation is derived using the product rule for the spatial difference. The conservation properties of the internal, kinetic, and magnetic energy equations can be satisfied
-
Neural network guided adjoint computations in dual weighted residual error estimation arXiv.cs.NA Pub Date : 2021-02-24 Julian Roth; Max Schröder; Thomas Wick
In this work, we are concerned with neural network guided goal-oriented a posteriori error estimation and adaptivity using the dual weighted residual method. The primal problem is solved using classical Galerkin finite elements. The adjoint problem is solved in strong form with a feedforward neural network using two or three hidden layers. The main objective of our approach is to explore alternatives
-
A deep neural network approach on solving the linear transport model under diffusive scaling arXiv.cs.NA Pub Date : 2021-02-24 Liu Liu; Tieyong Zeng; Zecheng Zhang
In this work, we propose a learning method for solving the linear transport equation under the diffusive scaling. Due to the multiscale nature of our model equation, the model is challenging to solve by using conventional methods. We employ the physical informed neural network (PINN) framework, a mesh-free learning method that can numerically solve partial differential equations. Compared to conventional
-
Convergence in the maximum norm of ADI-type methods for parabolic problems arXiv.cs.NA Pub Date : 2021-02-24 S. Gonzalez Pinto; D. Hernandez Abreu
Results on unconditional convergence in the Maximum norm for ADI-type methods, such as the Douglas method, applied to the time integration of semilinear parabolic problems are quite difficult to get, mainly when the number of space dimensions $m$ is greater than two. Such a result is obtained here under quite general conditions on the PDE problem in case that time-independent Dirichlet boundary conditions
-
A Provably Componentwise Backward Stable $O(n^2)$ QR Algorithm for the Diagonalization of Colleague Matrices arXiv.cs.NA Pub Date : 2021-02-24 Kirill Serkh; Vladimir Rokhlin
The roots of a monic polynomial expressed in a Chebyshev basis are known to be the eigenvalues of the so-called colleague matrix, which is a Hessenberg matrix that is the sum of a symmetric tridiagonal matrix and a rank-1 matrix. The rootfinding problem is thus reformulated as an eigenproblem, making the computation of the eigenvalues of such matrices a subject of significant practical importance.
-
Multiresolution-based mesh adaptation and error control for lattice Boltzmann methods with applications to hyperbolic conservation laws arXiv.cs.NA Pub Date : 2021-02-24 Thomas BellottiCMAP; Loïc GouarinCMAP; Benjamin GrailleLMO; Marc MassotCMAP
Lattice Boltzmann Methods (LBM) stand out for their simplicity and computational efficiency while offering the possibility of simulating complex phenomena. While they are optimal for Cartesian meshes, adapted meshes have traditionally been a stumbling block since it is difficult to predict the right physics through various levels of meshes. In this work, we design a class of fully adaptive LBM methods
-
Learning optimal multigrid smoothers via neural networks arXiv.cs.NA Pub Date : 2021-02-24 Ru Huang; Ruipeng Li; Yuanzhe Xi
Multigrid methods are one of the most efficient techniques for solving linear systems arising from Partial Differential Equations (PDEs) and graph Laplacians from machine learning applications. One of the key components of multigrid is smoothing, which aims at reducing high-frequency errors on each grid level. However, finding optimal smoothing algorithms is problem-dependent and can impose challenges
-
A learning scheme by sparse grids and Picard approximations for semilinear parabolic PDEs arXiv.cs.NA Pub Date : 2021-02-24 Jean-François Chassagneux; Junchao Chen; Noufel Frikha; Chao Zhou
Relying on the classical connection between Backward Stochastic Differential Equations (BSDEs) and non-linear parabolic partial differential equations (PDEs), we propose a new probabilistic learning scheme for solving high-dimensional semi-linear parabolic PDEs. This scheme is inspired by the approach coming from machine learning and developed using deep neural networks in Han and al. [32]. Our algorithm
-
Operator preconditioning: the simplest case arXiv.cs.NA Pub Date : 2021-02-23 Rob Stevenson; Raymond van Venetië
Using the framework of operator or Calder\'on preconditioning, uniform preconditioners are constructed for elliptic operators discretized with continuous finite (or boundary) elements. The preconditioners are constructed as the composition of an opposite order operator, discretized on the same ansatz space, and two diagonal scaling operators.
-
High order positivity preserving and asymptotic preserving multi-derivative methods arXiv.cs.NA Pub Date : 2021-02-23 Sigal Gottlieb; Zachary J. Grant; Jingwei Hu; Ruiwen Shu
In this work we present multi-derivative implicit-explicit (IMEX) Runge--Kutta schemes. We derive their order conditions up to third order, and show that such methods can preserve positivity (and more generally strong stability) with a time-step restriction independent of the stiff term, under mild assumptions on the operators. We present sufficient conditions under which such methods are positivity
-
Transfer function interpolation remainder formula of rational Krylov subspace methods arXiv.cs.NA Pub Date : 2021-02-23 Yiding Lin
Rational Krylov subspace projection methods are one of successful methods in MOR, mainly because some order derivatives of the approximate and original transfer functions are the same. This is the well known moments matching result. However, the properties of points which are far from the interpolating points are little known. In this paper, we obtain the error's explicit expression which involves
-
Construction of arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes arXiv.cs.NA Pub Date : 2021-02-23 Matthew W. Scroggs; Jørgen S. Dokken; Chris N. Richardson; Garth N. Wells
We develop an approach to generating degree-of-freedom maps for arbitrary order finite element spaces for any cell shape. The approach is based on the composition of permutations and transformations by cell sub-entity. Current approaches to generating degree-of-freedom maps for arbitrary order problems typically rely on a consistent orientation of cell entities that permits the definition of a common
-
Reconstruction, with tunable sparsity levels, of shear-wave velocity profiles from surface wave data arXiv.cs.NA Pub Date : 2021-02-24 Giulio Vignoli; Julien Guillemoteau; Jeniffer Barreto; Matteo Rossi
The analysis of surface wave dispersion curves is a way to infer the vertical distribution of shear-wave velocity. The range of applicability is extremely wide going, for example, from seismological studies to geotechnical characterizations and exploration geophysics. However, the inversion of the dispersion curves is severely ill-posed and only limited efforts have been put into the development of
-
Sparse online variational Bayesian regression arXiv.cs.NA Pub Date : 2021-02-24 Kody J. H. Law; Vitaly Zankin
This work considers variational Bayesian inference as an inexpensive and scalable alternative to a fully Bayesian approach in the context of sparsity-promoting priors. In particular, the priors considered arise from scale mixtures of Normal distributions with a generalized inverse Gaussian mixing distribution. This includes the variational Bayesian LASSO as an inexpensive and scalable alternative to
-
On Unbiased Estimation for Discretized Models arXiv.cs.NA Pub Date : 2021-02-24 Jeremy Heng; Ajay Jasra; Kody J. H. Law; Alexander Tarakanov
In this article, we consider computing expectations w.r.t. probability measures which are subject to discretization error. Examples include partially observed diffusion processes or inverse problems, where one may have to discretize time and/or space, in order to practically work with the probability of interest. Given access only to these discretizations, we consider the construction of unbiased Monte
-
Adaptive two- and three-dimensional multiresolution computations of resistive magnetohydrodynamics arXiv.cs.NA Pub Date : 2021-02-23 Anna Karina Fontes Gomes; Margarete Oliveira Domingues; Odim Mendes; Kai Schneider
Fully adaptive computations of the resistive magnetohydrodynamic (MHD) equations are presented in two and three space dimensions using a finite volume discretization on locally refined dyadic grids. Divergence cleaning is used to control the incompressibility constraint of the magnetic field. For automatic grid adaptation a cell-averaged multiresolution analysis is applied which guarantees the precision
-
Three Ways to Solve Partial Differential Equations with Neural Networks -- A Review arXiv.cs.NA Pub Date : 2021-02-23 Jan Blechschmidt; Oliver G. Ernst
Neural networks are increasingly used to construct numerical solution methods for partial differential equations. In this expository review, we introduce and contrast three important recent approaches attractive in their simplicity and their suitability for high-dimensional problems: physics-informed neural networks, methods based on the Feynman-Kac formula and the Deep BSDE solver. The article is
-
Functional norms, condition numbers and numerical algorithms in algebraic geometry arXiv.cs.NA Pub Date : 2021-02-23 Felipe Cucker; Alperen A. Ergür; Josué Tonelli-Cueto
In numerical linear algebra, a well-established practice is to choose a norm that exploits the structure of the problem at hand in order to optimize accuracy or computational complexity. In numerical polynomial algebra, a single norm (attributed to Weyl) dominates the literature. This article initiates the use of $L_p$ norms for numerical algebraic geometry, with an emphasis on $L_{\infty}$. This classical
-
Deep ReLU Neural Network Approximation for Stochastic Differential Equations with Jumps arXiv.cs.NA Pub Date : 2021-02-23 Lukas Gonon; Christoph Schwab
Deep neural networks (DNNs) with ReLU activation function are proved to be able to express viscosity solutions of linear partial integrodifferental equations (PIDEs) on state spaces of possibly high dimension $d$. Admissible PIDEs comprise Kolmogorov equations for high-dimensional diffusion, advection, and for pure jump L\'{e}vy processes. We prove for such PIDEs arising from a class of jump-diffusions
-
Twelve Ways To Fool The Masses When Giving Parallel-In-Time Results arXiv.cs.NA Pub Date : 2021-02-23 Sebastian Goetschel; Michael Minion; Daniel Ruprecht; Robert Speck
Getting good speedup -- let alone high parallel efficiency -- for parallel-in-time (PinT) integration examples can be frustratingly difficult. The high complexity and large number of parameters in PinT methods can easily (and unintentionally) lead to numerical experiments that overestimate the algorithm's performance. In the tradition of Bailey's article "Twelve ways to fool the masses when giving
-
The nonconforming Trefftz virtual element method: general setting, applications, and dispersion analysis for the Helmholtz equation arXiv.cs.NA Pub Date : 2021-02-23 L. Mascotto; I. Perugia; A. Pichler
We present a survey of the nonconforming Trefftz virtual element method for the Laplace and Helmholtz equations. For the latter, we present a new abstract analysis, based on weaker assumptions on the stabilization, and numerical results on the dispersion analysis, including comparison with the plane wave discontinuous Galerkin method.
-
A preconditioner based on sine transform for two-dimensional Riesz space factional diffusion equations in convex domains arXiv.cs.NA Pub Date : 2021-02-23 Xin Huang; Hai-Wei Sun
In this paper, we develop a fast numerical method for solving the time-dependent Riesz space fractional diffusion equations with a nonlinear source term in the convex domain. An implicit finite difference method is employed to discretize the Riesz space fractional diffusion equations with a penalty term in a rectangular region by the volume-penalization approach. The stability and the convergence of
-
Monotone cubic spline interpolation for functions with a strong gradient arXiv.cs.NA Pub Date : 2021-02-23 Francesc Aràndiga; Antonio Baeza; Dionisio F. Yáñez
Spline interpolation has been used in several applications due to its favorable properties regarding smoothness and accuracy of the interpolant. However, when there exists a discontinuity or a steep gradient in the data, some artifacts can appear due to the Gibbs phenomenon. Also, preservation of data monotonicity is a requirement in some applications, and that property is not automatically verified
-
Explicit high-order generalized-$α$ methods for isogeometric analysis of structural dynamics arXiv.cs.NA Pub Date : 2021-02-23 Pouria Behnoudfar; Gabriele Loli; Alessandro Reali; Giancarlo Sangalli; Victor M. Calo
We propose a new family of high-order explicit generalized-$\alpha$ methods for hyperbolic problems with the feature of dissipation control. Our approach delivers $2k,\, \left(k \in \mathbb{N}\right)$ accuracy order in time by solving $k$ matrix systems explicitly and updating the other $2k$ variables at each time-step. The user can control the numerical dissipation in the discrete spectrum's high-frequency
-
Convergence rates for gradient descent in the training of overparameterized artificial neural networks with biases arXiv.cs.NA Pub Date : 2021-02-23 Arnulf Jentzen; Timo Kröger
In recent years, artificial neural networks have developed into a powerful tool for dealing with a multitude of problems for which classical solution approaches reach their limits. However, it is still unclear why randomly initialized gradient descent optimization algorithms, such as the well-known batch gradient descent, are able to achieve zero training loss in many situations even though the objective
-
A Generalized Eulerian-Lagrangian Discontinuous Galerkin Method for Transport Problems arXiv.cs.NA Pub Date : 2021-02-22 Xue Hong; Jing-Mei Qiu
We propose a generalized Eulerian-Lagrangian (GEL) discontinuous Galerkin (DG) method. The method is a generalization of the Eulerian-Lagrangian (EL) DG method for transport problems proposed in [arXiv preprint arXiv: 2002.02930 (2020)], which tracks solution along approximations to characteristics in the DG framework, allowing extra large time stepping size with stability. The newly proposed GEL DG
-
Solving high-dimensional parabolic PDEs using the tensor train format arXiv.cs.NA Pub Date : 2021-02-23 Lorenz Richter; Leon Sallandt; Nikolas Nüsken
High-dimensional partial differential equations (PDEs) are ubiquitous in economics, science and engineering. However, their numerical treatment poses formidable challenges since traditional grid-based methods tend to be frustrated by the curse of dimensionality. In this paper, we argue that tensor trains provide an appealing approximation framework for parabolic PDEs: the combination of reformulations
-
Higher order phase averaging for highly oscillatory systems arXiv.cs.NA Pub Date : 2021-02-23 Werner Bauer; Colin J. Cotter; Beth Wingate
We introduce a higher order phase averaging method for nonlinear oscillatory systems. Phase averaging is a technique to filter fast motions from the dynamics whilst still accounting for their effect on the slow dynamics. Phase averaging is useful for deriving reduced models that can be solved numerically with more efficiency, since larger timesteps can be taken. Recently, Haut and Wingate (2014) introduced
-
Using a deep neural network to predict the motion of under-resolved triangular rigid bodies in an incompressible flow arXiv.cs.NA Pub Date : 2021-02-23 Henry von Wahl; Thomas Richter
We consider non-spherical rigid body particles in an incompressible fluid in the regime where the particles are too large to assume that they are simply transported with the fluid without back-coupling and where the particles are also too small to make fully resolved direct numerical simulations feasible. Unfitted finite element methods with ghost-penalty stabilisation are well suited to fluid-structure-interaction
-
Actor-Critic Method for High Dimensional Static Hamilton--Jacobi--Bellman Partial Differential Equations based on Neural Networks arXiv.cs.NA Pub Date : 2021-02-22 Mo Zhou; Jiequn Han; Jianfeng Lu
We propose a novel numerical method for high dimensional Hamilton--Jacobi--Bellman (HJB) type elliptic partial differential equations (PDEs). The HJB PDEs, reformulated as optimal control problems, are tackled by the actor-critic framework inspired by reinforcement learning, based on neural network parametrization of the value and control functions. Within the actor-critic framework, we employ a policy
-
On Theoretical and Numerical Aspect of Fractional Differential Equations with Purely Integral Conditions arXiv.cs.NA Pub Date : 2021-02-03 Saadoune Brahimi; Ahcene Merad; Adem Kilicman
In this paper, we are interested in the study of a problem with fractional derivatives having boundary conditions of integral types. The problem represents a Caputo type advection-diffusion equation where the fractional order derivative with respect to time with $1<\alpha <2$. The method of the energy inequalities is used to prove the existence and the uniqueness of solutions of the problem. The finite
-
Convergence error estimates at low regularity for time discretizations of KdV arXiv.cs.NA Pub Date : 2021-02-22 Frédéric Rousset; Katharina Schratz
We consider various filtered time discretizations of the periodic Korteweg--de Vries equation: a filtered exponential integrator, a filtered Lie splitting scheme as well as a filtered resonance based discretisation and establish convergence error estimates at low regularity. Our analysis is based on discrete Bourgain spaces and allows to prove convergence in $L^2$ for rough data $u_{0} \in H^s,$ $s>0$
-
Energy stable arbitrary order ETD-MS method for gradient flows with Lipschitz nonlinearity arXiv.cs.NA Pub Date : 2021-02-22 Wenbin Chen; Shufen Wang; Xiaoming Wang
We present a methodology to construct efficient high-order in time accurate numerical schemes for a class of gradient flows with appropriate Lipschitz continuous nonlinearity. There are several ingredients to the strategy: the exponential time differencing (ETD), the multi-step (MS) methods, the idea of stabilization, and the technique of interpolation. They are synthesized to develop a generic $k^{th}$
-
Singular Euler-Maclaurin expansion on multidimensional lattices arXiv.cs.NA Pub Date : 2021-02-22 Andreas A. Buchheit; Torsten Keßler
We extend the classical Euler-Maclaurin expansion to sums over multidimensional lattices that involve functions with algebraic singularities. This offers a tool for the precise quantification of the effect of microscopic discreteness on macroscopic properties of a system. First, the Euler-Maclaurin summation formula is generalised to lattices in higher dimensions, assuming a sufficiently regular summand
Contents have been reproduced by permission of the publishers.