BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Robert Schaback

Robert Schaback appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysisphysics.soc-phPopulations and Evolution

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Small Errors Imply Large Instabilities

Numerical Analysts and scientists working in applications often observe that once they improve their techniques to get a better accuracy, some instability creeps in through the back door. This paper shows for a large class of numerical methods that such a Trade-off Principle between error and stability is unavoidable. The setting is confined to recovery of functions from data, but it includes solving differential equations by writing such methods as a recovery of functions under constraints imposed by differential operators and boundary values. It is shown in particular that Kansa's Unsymmetric Collocation Method sacrifices accuracy for stability, when compared to symmetric collocation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modelling Recovered Cases and Death Probabilities for the COVID-19 Outbreak

From March 23rd, the data for the recovered cases of COVID-19 are missing from the standard repository maintained by the Johns Hopkins University in collaboration with the WHO. But since data concerning recovered patients are extremely important for modelling the COVID-19 outbreak, a method for estimating the missing data is provided and tested. As a byproduct, it produces estimates for the probabilities to die $k$ days after confirmation, or to survive after $d$ days.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On COVID-19 Modelling

This contribution analyzes the COVID-19 outbreak by comparably simple mathematical and numerical methods. The final goal is to predict the peak of the epidemic outbreak per country with a reliable technique. This is done by an algorithm motivated by standard SIR models and aligned with the standard data provided by the Johns Hopkins University. To reconstruct data for the unregistered Infected, the algorithm uses current values of the infection fatality rate and a data-driven estimation of a specific form of the recovery rate. All other ingredients are data-driven as well. Various examples of predictions are provided for illustration.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive Approximation of Functions with Discontinuities

One of the basic principles of Approximation Theory is that the quality of approximations increase with the smoothness of the function to be approximated. Functions that are smooth in certain subdomains will have good approximations in those subdomains, and these {\em sub-approximations} can possibly be calculated efficiently in parallel, as long as the subdomains do not overlap. This paper proposes a class of algorithms that first calculate sub-approximations on non-overlapping subdomains, then extend the subdomains as much as possible and finally produce a global solution on the given domain by letting the subdomains fill the whole domain. Consequently, there will be no Gibbs phenomenon along the boundaries of the subdomains. Throughout, the algorithm works for fixed scattered input data of the function itself, not on spectral data, and it does not resample.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Error Analysis of Nodal Meshless Methods

There are many application papers that solve elliptic boundary value problems by meshless methods, and they use various forms of generalized stiffness matrices that approximate derivatives of functions from values at scattered nodes $x_1,\ldots,x_M\in Ω\subset\R^d$. If $u^*$ is the true solution in some Sobolev space $S$ allowing enough smoothness for the problem in question, and if the calculated approximate values at the nodes are denoted by $\tilde u_1,\ldots,\tilde u_M$, the canonical form of error bounds is $$ \max_{1\leq j\leq M}|u^*(x_j)-\tilde u_j|\leq ε\|u^*\|_S $$ where $ε$ depends crucially on the problem and the discretization, but not on the solution. This contribution shows how to calculate such $ε$ {\em numerically and explicitly}, for any sort of discretization of strong problems via nodal values, may the discretization use Moving Least Squares, unsymmetric or symmetric RBF collocation, or localized RBF or polynomial stencils. This allows users to compare different discretizations with respect to error bounds of the above form, without knowing exact solutions, and admitting all possible ways to set up generalized stiffness matrices. The error analysis is proven to be sharp under mild additional assumptions. As a byproduct, it allows to construct worst cases that push discretizations to their limits. All of this is illustrated by numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Fractional Derivatives of Radial Basis Functions

The paper provides the fractional integrals and derivatives of the Rie\-mann-Liouville and Caputo type for the five kinds of radial basis functions (RBFs), including the powers, Gaussian, multiquadric, Matern and thin-plate splines, in one dimension. It allows to use high order numerical methods for solving fractional differential equations. The results are tested by solving two fractional differential equations. The first one is a fractional ODE which is solved by the RBF collocation method and the second one is a fractional PDE which is solved by the method of lines based on the spatial trial spaces spanned by the Lagrange basis associated to the RBFs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimal Stencils in Sobolev Spaces

This paper proves that the approximation of pointwise derivatives of order $s$ of functions in Sobolev space $W_2^m(\R^d)$ by linear combinations of function values cannot have a convergence rate better than $m-s-d/2$, no matter how many nodes are used for approximation and where they are placed. These convergence rates are attained by {\em scalable} approximations that are exact on polynomials of order at least $\lfloor m-d/2\rfloor +1$, proving that the rates are optimal for given $m,\,s,$ and $d$. And, for a fixed node set $X\subset\R^d$, the convergence rate in any Sobolev space $W_2^m(Ω)$ cannot be better than $q-s$ where $q$ is the maximal possible order of polynomial exactness of approximations based on $X$, no matter how large $m$ is. In particular,scalable stencil constructions via polyharmonic kernels are shown to realize the optimal convergence rates, and good approximations of their error in Sobolev space can be calculated via their error in Beppo-Levi spaces. This allows to construct near-optimal stencils in Sobolev spaces stably and efficiently, for use in meshless methods to solve partial differential equations via generalized finite differences (RBF-FD). Numerical examples are included for illustration.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Superconvergence of Kernel-Based Interpolation

It is well-known that univariate cubic spline interpolation, if carried out on point sets with fill distance $h$, converges only like ${\cal O}(h^2)$ in $L_2[a,b]$ for functions in $W_2^2[a,b]$ if no additional assumptions are made. But superconvergence up to order $h^4$ occurs if more smoothness is assumed and if certain additional boundary conditions are satisfied. This phenomenon was generalized in 1999 to multivariate interpolation in Reproducing Kernel Hilbert Spaces on domains $Ω\subset R^d$ for continuous positive definite Fourier-transformable shift-invariant kernels on $R^d$. But the sufficient condition for superconvergence given in 1999 still needs further analysis, because the interplay between smoothness and boundary conditions is not clear at all. Furthermore, if only additional smoothness is assumed, superconvergence is numerically observed in the interior of the domain, but without explanation, so far. This paper first generalizes the "improved error bounds" of 1999 by an abstract theory that includes the Aubin-Nitsche trick and the known superconvergence results for univariate polynomial splines. Then the paper analyzes what is behind the sufficient conditions for superconvergence. They split into conditions on smoothness and localization, and these are investigated independently. If sufficient smoothness is present, but no additional localization conditions are assumed, it is proven that superconvergence always occurs in the interior of the domain. If smoothness and localization interact in the kernel-based case on $R^d$, weak and strong boundary conditions in terms of pseudodifferential operators occur. A special section on Mercer expansions is added, because Mercer eigenfunctions always satisfy the sufficient conditions for superconvergence. Numerical examples illustrate the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

All Well--Posed Problems have Uniformly Stable and Convergent Discretizations

This paper considers a large class of linear operator equations, including linear boundary value problems for partial differential equations, and treats them as linear recovery problems for objects from their data. Well-posedness of the problem means that this recovery is continuous. Discretization recovers restricted trial objects from restricted test data, and it is well-posed or stable, if this restricted recovery is continuous. After defining a general framework for these notions, this paper proves that all well-posed linear problems have stable and refinable computational discretizations with a stability that is determined by the well-posedness of the problem and independent of the computational discretization. The solutions of discretized problems converge when enlarging the trial spaces, and the convergence rate is determined by how well the full data of the object solving the full problem can be approximated by the full data of the trial objects. %of the discretization. This allows very simple proofs of convergence rates for generalized finite elements, symmetric and unsymmetric Kansa-type collocation, and other meshfree methods like Meshless Local Petrov-Galerkin techniques. It is also shown that for a fixed trial space, weak formulations have a slightly better convergence rate than strong formulations, but at the expense of numerical integration. Since convergence rates are reduced to those coming from Approximation Theory, and since trial spaces are arbitrary, this also covers various spectral and pseudospectral methods. All of this is illustrated by examples.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A computational tool for comparing all linear PDE solvers -- Optimal methods are meshless

The paper starts out with a computational technique that allows to compare all linear methods for PDE solving that use the same input data. This is done by writing them as linear recovery formulas for solution values as linear combinations of the input data. Calculating the norm of these reproduction formulas on a fixed Sobolev space will then serve as a quality criterion that allows a fair comparison of all linear methods with the same inputs, including finite-element, finite-difference and meshless local Petrov-Galerkin techniques. A number of illustrative examples will be provided. As a byproduct, it turns out that a unique error--optimal method exists. It necessarily outperforms any other competing technique using the same data, e.g. those just mentioned, and it is necessarily meshless, if solutions are written "entirely in terms of nodes" (Belytschko et. al. 1996). On closer inspection, it turns out that it coincides with {\em symmetric meshless collocation} carried out with the kernel of the Hilbert space used for error evaluation, e.g. with the kernel of the Sobolev space used. This technique is around since at least 1998, but its optimality properties went unnoticed, so far.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Direct meshless local Petrov-Galerkin (DMLPG) method: A generalized MLS approximation

The Meshless Local Petrov{Galerkin (MLPG) method is one of the popular meshless methods that has been used very successfully to solve several types of boundary value problems since the late nineties. In this paper, using a generalized moving least squares (GMLS) approximation, a new direct MLPG technique, called DMLPG, is presented. Following the principle of meshless methods to express everything "entirely in terms of nodes", the generalized MLS recovers test functionals directly from values at nodes, without any detour via shape functions. This leads to a cheaper and even more accurate scheme. In particular, the complete absence of shape functions allows numerical integrations in the weak forms of the problem to be done over low{degree polynomials instead of complicated shape functions. Hence, the standard MLS shape function subroutines are not called at all. Numerical examples illustrate the superiority of the new technique over the classical MLPG. On the theoretical side, this paper discusses stability and convergence for the new discretizations that replace those of the standard MLPG. However, it does not treat stability, convergence, or error estimation for the MLPG as a whole. This should be taken from the literature on MLPG.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solving Heat Conduction Problems by the Direct Meshless Local Petrov-Galerkin (DMLPG) method

As an improvement of the Meshless Local Petrov-Galerkin (MLPG), the Direct Meshless Local Petrov-Galerkin (DMLPG) method is applied here to the numerical solution of transient heat conduction problem. The new technique is based on direct recoveries of test functionals (local weak forms) from values at nodes without any detour via classical moving least squares (MLS) shape functions. This leads to an absolutely cheaper scheme where the numerical integrations will be done over low-degree polynomials rather than complicated MLS shape functions. This eliminates the main disadvantage of MLS based methods in comparison with finite element methods (FEM), namely the costs of numerical integration.

0 citations0 reviewsNo code linkedOpen access