Search

This work is concerned with the interpolation of a function f when using a low number of interpolation points, as required by the finite element method for solving PDEs numerically. The function f is assumed to have a jump or a steep derivative, and our goal is to minimize the oscillations produced by the Gibbs phenomenon while preserving the approximation properties for smoother functions. This is achieved by interpolating the transform f ^ = g ∘ f using Lagrange polynomials, where g is a rational transformation chosen by minimizing a suitable functional depending on the values of f . The mapping g is monotonic and constructed to possess boundary layers that remove the Gibbs phenomenon. No previous knowledge of the location of the jump is required. The extension to functions of several variables is straightforward, of which we provide several examples. Finally, we show how the interpolation fits the finite element method and compare it with known strategies. [ABSTRACT FROM AUTHOR]

Polynomial reconstruction on Cartesian grids is fundamental in many scientific and engineering applications, yet it is still an open problem how to construct for a finite subset K of Z D a lattice T ⊂ K so that multivariate polynomial interpolation on this lattice is unisolvent. In this work, we solve this open problem of poised lattice generation (PLG) via an interdisciplinary research of approximation theory, abstract algebra, and artificial intelligence. Specifically, we focus on the triangular lattices in approximation theory, study group actions of permutations upon triangular lattices, prove an isomorphism between the group of permutations and that of triangular lattices, and dynamically organize the state space of permutations so that a depth-first search of poised lattices has optimal efficiency. Based on this algorithm, we further develop the PLG finite difference method that retains the simplicity of Cartesian grids yet overcomes the disadvantage of legacy finite difference methods in handling irregular geometries. Results of various numerical tests demonstrate the effectiveness of our algorithm and the simplicity, flexibility, efficiency, and fourth-order accuracy of the PLG finite difference method. [ABSTRACT FROM AUTHOR]

We construct some smooth functions defined over a sphere that interpolate large sets of scattered data, using some modified Shepard methods, the least squares thin-plate spline and the inverse multiquadric functions. We illustrate the benefits of our methods in numerical examples for some test functions and some real data applications. [ABSTRACT FROM AUTHOR]

The radial basis function (RBF) approximation is a rapidly developing field of mathematics. In the paper, we are concerned with the measurement of scalar physical quantities at nodes on sphere in the 3D Euclidean space and the spherical RBF interpolation of the data acquired. We employ a multiquadric as the radial basis function and the corresponding trend is a polynomial of degree 2 considered in Cartesian coordinates. Attention is paid to geodesic metrics that define the distance of two points on a sphere. The choice of a particular geodesic metric function is an important part of the construction of interpolation formula. We show the existence of an interpolation formula of the type considered. The approximation formulas of this type can be useful in the interpretation of measurements of various physical quantities. We present an example concerned with the sampling of anisotropy of magnetic susceptibility having extensive applications in geosciences and demonstrate the advantages and drawbacks of the formulas chosen, in particular the strong dependence of interpolation results on condition number of the matrix of the system considered and on round-off errors in general. [ABSTRACT FROM AUTHOR]

In an error estimation of finite element solutions to the Poisson equation, we usually impose the shape regularity assumption on the meshes to be used. In this paper, we show that even if the shape regularity condition is violated, the standard error estimation can be obtained if "bad" elements that violate the shape regularity or maximum angle condition are covered virtually by simplices that satisfy the minimum angle condition. A numerical experiment illustrates the theoretical result. [ABSTRACT FROM AUTHOR]

Standard interpolatory subdivision schemes and their underlying interpolating refinable functions are of interest in CAGD, numerical PDEs, and approximation theory. Generalizing these notions, we introduce and study n s -step interpolatory M -subdivision schemes and their interpolating M -refinable functions with n s ∈ N ∪ { ∞ } and a dilation factor M ∈ N \ { 1 } . We completely characterize C m -convergence and smoothness of n s -step interpolatory subdivision schemes and their interpolating M -refinable functions in terms of their masks. Inspired by n s -step interpolatory stationary subdivision schemes, we further introduce the notion of r-mask quasi-stationary subdivision schemes, and then we characterize their C m -convergence and smoothness properties using only their masks. Moreover, combining n s -step interpolatory subdivision schemes with r-mask quasi-stationary subdivision schemes, we can obtain r n s -step interpolatory subdivision schemes. Examples and construction procedures of convergent n s -step interpolatory M -subdivision schemes are provided to illustrate our results with dilation factors M = 2 , 3 , 4 . In addition, for the dyadic dilation M = 2 and r = 2 , 3 , using r masks with only two-ring stencils, we provide examples of C r -convergent r-step interpolatory r-mask quasi-stationary dyadic subdivision schemes. [ABSTRACT FROM AUTHOR]

A proof that Wronskians of non-zero B-splines are positive is given, using only recursive formulas for B-splines and their derivatives. This could be used to generalize the de Boor–DeVore geometric proof of the Schoenberg–Whitney conditions and total positivity of B-splines to Hermite interpolation. For Wronskians of maximal order with respect to a given degree, positivity follows from a simple formula. [ABSTRACT FROM AUTHOR]

In this study, we investigate a hybrid-type anisotropic weakly over-penalised symmetric interior penalty method for the Poisson equation on convex domains. Compared with the well-known hybrid discontinuous Galerkin methods, our approach is simple and easy to implement. Our primary contributions are the proposal of a new scheme and the demonstration of a proof for the consistency term, which allows us to estimate the anisotropic consistency error. The key idea of the proof is to apply the relation between the Raviart–Thomas finite element space and a discontinuous space. In numerical experiments, we compare the calculation results for standard and anisotropic mesh partitions. [ABSTRACT FROM AUTHOR]

This article explores various techniques for generating fractal-like Bézier curves in both 2D and 3D environments. It delves into methods such as subdivision schemes, Iterated Function System (IFS) theory, perturbation of Bézier curves, and perturbation of Bézier basis functions. The article outlines conditions on subdivision matrices necessary for convergence and demonstrates their use in creating an IFS with an attractor aligned to the convergent point of the subdivision scheme based on specified initial data. Additionally, it discusses conditions for obtaining a one-sided approximation of a given Bézier curve through perturbation. The article also addresses considerations for perturbed Bézier basis functions to construct fractal-like Bézier curves that remain within the convex hull polygon/polyhedron defined by control points. These methods find applications in various fields, including computer graphics, art, and design. [ABSTRACT FROM AUTHOR]

Interpolation-based methods are well-established and effective approaches for the efficient generation of accurate reduced-order surrogate models. Common challenges for such methods are the automatic selection of good or even optimal interpolation points and the appropriate size of the reduced-order model. An approach that addresses the first problem for linear, unstructured systems is the iterative rational Krylov algorithm (IRKA), which computes optimal interpolation points through iterative updates by solving linear eigenvalue problems. However, in the case of preserving internal system structures, optimal interpolation points are unknown, and heuristics based on nonlinear eigenvalue problems result in numbers of potential interpolation points that typically exceed the reasonable size of reduced-order systems. In our work, we propose a projection-based iterative interpolation method inspired by IRKA for generally structured systems to adaptively compute near-optimal interpolation points as well as an appropriate size for the reduced-order system. Additionally, the iterative updates of the interpolation points can be chosen such that the reduced-order model provides an accurate approximation in specified frequency ranges of interest. For such applications, our new approach outperforms the established methods in terms of accuracy and computational effort. We show this in numerical examples with different structures. [ABSTRACT FROM AUTHOR]

In this study, we investigate an anisotropic weakly over-penalised symmetric interior penalty method for the Stokes equation on convex domains. Our approach is a simple discontinuous Galerkin method similar to the Crouzeix–Raviart finite element method. As our primary contribution, we show a new proof for the consistency term, which allows us to obtain an estimate of the anisotropic consistency error. The key idea of the proof is to apply the relation between the Raviart–Thomas finite element space and a discontinuous space. While inf-sup stable schemes of the discontinuous Galerkin method on shape-regular mesh partitions have been widely discussed, our results show that the Stokes element satisfies the inf-sup condition on anisotropic meshes. Furthermore, we provide an error estimate in an energy norm on anisotropic meshes. In numerical experiments, we compare calculation results for standard and anisotropic mesh partitions. [ABSTRACT FROM AUTHOR]

Composition detection and semantic line detection are important research topics in computer vision and play an important auxiliary role in the analysis of image esthetics. However, at present, few researchers have considered the internal relationship between these two related tasks for comprehensive research. In order to solve this problem, we propose a synchronous detection network of composition class and semantic lines based on cross-attention, which can realize the mutual supervision and guidance between composition class detection and semantic line detection, to improve the accuracy of each other’s detection. First, the pre-trained composition detection model and the pre-trained semantic line detection model as two teacher models to provide data labels of composition and semantic line information for the student model. Then, we train a student model with the help of the teacher model. The student model adopts the multi-task learning architecture by combining soft and hard parameter sharing, as we propose. At the same time, we develop a cross-attention module to ensure that both tasks get the help and supervision they need from each other. Experimental results show that our method can draw semantic lines while detecting composition classes, which increases the interpretability of composition class detection. Our composition detection accuracy reaches 92.57%, and for benchmark semantic lines, the accuracy of our AUC_A metric can reach 92.00%. [ABSTRACT FROM AUTHOR]

The focus of this investigation centers on the formulation of a novel spectral method deployed to numerically solve partial integro-differential equations that possess memory features. The approach leverages the hexic shifted Chebyshev polynomials, and addresses the nonlinearity of the computational output using iterative methods. A comprehensive explanation of the methodology is presented, and a thorough convergence analysis is conducted. This technique has the capacity to be effortlessly modified and used to solve a plethora of linear and nonlinear issues, while minimizing computational time. Ultimately, the effectiveness of this novel strategy is demonstrated through the successful resolution of two exemplary problems. The findings of this study suggest that this spectral approach holds significant promise for solving partial integro-differential equations. [ABSTRACT FROM AUTHOR]

This article is devoted to introducing new spline quasi-interpolants for the sharp approximation of data in one and several dimensions. The new construction aims to obtain accurate approximations close to singularities of the function from which the data are obtained. The technique relies in an accurate knowledge of the position of the singularity, which can be known or approximated, and that allows for the obtention of accurate approximations of the jumps of the function and its derivatives at this point. With this information, it is possible to compute correction terms for the B-spline bases of the spline that are affected by the presence of singularities. The result is a piecewise smooth reconstruction with full accuracy (meaning by full accuracy, the accuracy of the spline at smooth zones). The numerical experiments presented support the theoretical results obtained. [ABSTRACT FROM AUTHOR]

The Sinc approximation applied to double-exponentially decaying functions is referred to as the DE-Sinc approximation. Because of its high efficiency, this method has been used in various applications. In the Sinc approximation, the mesh size and truncation numbers should be optimally selected to achieve its best performance. However, the standard selection formula has only been "near-optimally" selected because the optimal formula of the mesh size cannot be expressed in terms of elementary functions of truncation numbers. In this study, we propose two improved selection formulas. The first one is based on the concept by an earlier research that resulted in a better selection formula for the double-exponential formula. The formula performs slightly better than the standard one, but is still not optimal. As a second selection formula, we introduce a new parameter to propose truly optimal selection formula. We provide explicit error bounds for both selection formulas. Numerical comparisons show that the first formula gives a better error bound than the standard formula, and the second formula gives a much better error bound than the standard and first formulas. [ABSTRACT FROM AUTHOR]

High-dimensional/high-fidelity nonlinear dynamical systems appear naturally when the goal is to accurately model real-world phenomena. Many physical properties are thereby encoded in the internal differential structure of these resulting large-scale nonlinear systems. The high dimensionality of the dynamics causes computational bottlenecks, especially when these large-scale systems need to be simulated for a variety of situations such as different forcing terms. This motivates model reduction where the goal is to replace the full-order dynamics with accurate reduced-order surrogates. Interpolation-based model reduction has been proven to be an effective tool for the construction of cheap-to-evaluate surrogate models that preserve the internal structure in the case of weak nonlinearities. In this paper, we consider the construction of multivariate interpolants in frequency domain for structured quadratic-bilinear systems. We propose definitions for structured variants of the symmetric subsystem and generalized transfer functions of quadratic-bilinear systems and provide conditions for structure-preserving interpolation by projection. The theoretical results are illustrated using two numerical examples including the simulation of molecular dynamics in crystal structures. [ABSTRACT FROM AUTHOR]

Bivariate Mittag-Leffler (ML) functions are a substantial generalization of the univariate ML functions, which are widely recognized for their significance in fractional calculus. In the present paper, our initial focus is to investigate the fractional calculus properties of the integral and derivative operators with kernels including the Bivariate ML functions. Further, certain fractional Cauchy-type problems including these operators are considered. Also the numerical approximations of the Caputo type derivative operator are investigated. The theoretical results are justified by applications on examples. Furthermore, the theory of applying the same operators with respect to arbitrary monotonic functions is analyzed in this research. [ABSTRACT FROM AUTHOR]

Linear approximation methods suffer from Gibbs oscillations when approximating functions with jumps. Essentially non oscillatory subcell-resolution (ENO-SR) is a local technique avoiding oscillations and with a full order of accuracy, but a loss of regularity of the approximant appears. The goal of this paper is to introduce a new approach having both properties of full accuracy and regularity. In order to obtain it, we propose a three-stage algorithm: first, the data is smoothed by subtracting an appropriate non-smooth data sequence; then a chosen high order linear approximation operator is applied to the smoothed data and finally, an approximation with the proper jump or corner (jump in the first order derivative) discontinuity structure is reinstated by correcting the smooth approximation with the non-smooth element used in the first stage. This new procedure can be applied as subdivision scheme to design curves and surfaces both in point-value and in cell-average contexts. Using the proposed algorithm, we are able to construct approximations with high precision, with high piecewise regularity, and without smearing nor oscillations in the presence of discontinuities. These are desired properties in real applications as computer aided design or car design, among others. [ABSTRACT FROM AUTHOR]

We present a derivative-free separable quadratic modeling and cubic regularization technique for solving smooth unconstrained minimization problems. The derivative-free approach is mainly concerned with building a quadratic model that could be generated by numerical interpolation or using a minimum Frobenius norm approach, when the number of points available does not allow to build a complete quadratic model. This model plays a key role to generate an approximated gradient vector and Hessian matrix of the objective function at every iteration. We add a specialized cubic regularization strategy to minimize the quadratic model at each iteration, that makes use of separability. We discuss convergence results, including worst case complexity, of the proposed schemes to first-order stationary points. Some preliminary numerical results are presented to illustrate the robustness of the specialized separable cubic algorithm. [ABSTRACT FROM AUTHOR]

We investigate the approximation formulas that were proposed by Tanaka & Sugihara (IMA J. Numer. Anal. 39(4):1957–1984, 2019), in weighted Hardy spaces, which are analytic function spaces with certain asymptotic decay. Under the criterion of minimum worst error of n-point approximation formulas, we demonstrate that the formulas are nearly optimal. We also obtain the upper bounds of the approximation errors that coincide with the existing heuristic bounds in asymptotic order by a duality theorem for the minimization problem of potential energy. [ABSTRACT FROM AUTHOR]

In this work, a general class of interpolation and smoothing natural exponential splines with respect to fourth order differential operators with two real parameters is considered. Some sufficient conditions for the associated matrix R to be a diagonally dominant matrix are given. Based on these, fast algorithms for computing the coefficients of this general class of exponential splines are developed. The obtained splines have C 2 continuity and are the minimum solution of the combination of interpolation and smoothing energy integral. The performances of the resulting splines in financial data from the S &P500 index are given. Numerical experiments show that the resulting splines have more freedom to adjust the shape and control the energy of the curves. Cross-validation and generalized cross-validation for determining an appropriate smoothing parameter are also given. [ABSTRACT FROM AUTHOR]

In this paper, we consider the structure-preserving model order reduction problem for multi-input/multi-output bilinear control systems by tangential interpolation. We propose a new type of tangential interpolation problem for structured bilinear systems, for which we develop a new structure-preserving interpolation framework. This new framework extends and generalizes different formulations of tangential interpolation for bilinear systems from the literature and also provides a unifying framework. We then derive explicit conditions on the projection spaces to enforce tangential interpolation in different settings, including conditions for tangential Hermite interpolation. The analysis is illustrated by means of three numerical examples. [ABSTRACT FROM AUTHOR]

The best polynomial approximation and the Chebyshev approximation are both important in numerical analysis. In tradition, the best approximation is regarded as better than the Chebyshev approximation, because it is usually considered in the uniform norm. However, it is not always superior to the latter, as noticed by Trefethen (Math. Today 47:184–188, 2011) for the function f (x) = | x - 0.25 | . Recently, Wang (arxiv, 2106.03456, 2021) observed a similar phenomenon for a function with an algebraic singularity and proved it in theory. In this paper, we find that for functions with logarithmic endpoint singularities, the pointwise errors of the Chebyshev approximation are smaller than those of the best approximations of the same degree, except only at the narrow boundary layer. The pointwise error for the Chebyshev series truncated at the degree n is O (n - κ) ( κ = min { 2 γ + 1 , 2 δ + 1 } ), but increases as O (n - (κ - 1)) at the endpoint singularities. Theorems are given to explain this effect. [ABSTRACT FROM AUTHOR]

Based on Cauchy’s integral formula and conformal maps, this paper presents a new method for constructing barycentric rational interpolation formulae for complex functions, which may contain singularities such as poles, branch cuts, or essential singularities. The resulting interpolations are pole-free, exponentially convergent, and numerically stable, requiring only O (N) operations. Inspired by the logarithm equilibrium potential, we introduce a Möbius transform to concentrate nodes to the vicinity of singularity to get a spectacular improvement on approximation quality. A thorough convergence analysis is provided, alongside numerous numerical examples that illustrate the theoretical results and demonstrate the accuracy and efficiency of the methodology. Meanwhile, the paper also discusses some applications of the method including the numerical solutions of boundary value problems and the zero locations of holomorphic functions. [ABSTRACT FROM AUTHOR]

We introduce interpolation operators with approximation and stability properties suited for parabolic problems in primal and mixed formulations. We derive localized error estimates for tensor product meshes (occurring in classical time-marching schemes) as well as locally in space-time refined meshes. [ABSTRACT FROM AUTHOR]

In that paper, we new study has been carried out on previous studies of one of the most important mathematical models that describe the global economic movement, and that is described as a non-linear fractional financial model of awareness, where the studies are represented at the steps following: One: The schematic of the model is suggested. Two: The disease-free equilibrium point (DFE) and the stability of the equilibrium point are discussed. Three: The stability of the model is fulfilled by drawing the Lyapunov exponents and Poincare map. Fourth: The existence of uniformly stable solutions have discussed. Five: The Caputo is described as the fractional derivative. Six: Fractional optimal control for NFFMA is discussed by clarifying the fractional optimal control through drawing before and after control. Seven: Reduced differential transform method (RDTM) and Sumudu Decomposition Method (SDM) are used to take the resolution of an NFFMA. Finally, we display that SDM and RDTM are highly identical. [ABSTRACT FROM AUTHOR]

Calcium ions (Ca2+) serve as a crucial signaling mechanism in almost all cells. The buffers are proteins that bind free Ca2+ to reduce the cell's Ca2+ concentration. The most studies reported in the past on calcium signaling in various cells have considered the buffer concentration as constant in the cell. However, buffers also diffuse and their concentration varies dynamically in the cells. Almost no work has been reported on interdependent calcium and buffer dynamics in the cells. In the present study, a model is proposed for inter-dependent spatio-temporal dynamics of calcium and buffer by coupling reaction–diffusion equations of Ca2+ and buffer in a hepatocyte cell. Boundary and initial conditions are framed based on the physiological state of the cell. The effect of various parameters viz. inositol 1,4,5-triphosphate receptor (IP3R), diffusion coefficient, SERCA pump and ryanodine receptor (RyR) on spatio-temporal dynamics of calcium and buffer regulating diacylglycerol (DAG) in a normal and obese hepatocyte cell has been studied using finite element simulation. From the results, it is concluded that the dynamics of calcium and buffer impact each other significantly along the spatio-temporal dimensions, thereby affecting the regulation of all the processes including DAG in a hepatocyte cell. The proposed model is more realistic than the existing ones, as the interdependent system dynamics of calcium and buffer have different regulatory impacts as compared to the individual and independent dynamics of these signaling processes in a hepatocyte cell. [ABSTRACT FROM AUTHOR]

We consider the axisymmetric formulation of the equilibrium problem for a hot plasma in a tokamak. We adopt a non-overlapping mortar element approach, that couples C 0 piece-wise linear Lagrange finite elements in a region that does not contain the plasma and C 1 piece-wise cubic reduced Hsieh–Clough–Tocher finite elements elsewhere, to approximate the magnetic flux field on a triangular mesh of the poloidal tokamak section. The inclusion of ferromagnetic parts is simplified by assuming that they fit within the axisymmetric modeling and a new formulation of the Newton algorithm for the problem solution is stated, both in the static and quasi-static evolution cases. [ABSTRACT FROM AUTHOR]

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training. [ABSTRACT FROM AUTHOR]

In this paper, a kind of bivariate Bernoulli-type multiquadric quasi-interpolation operator is studied by combining the known multiquadric quasi-interpolation operator with the generalized Taylor polynomial as the expansion in the bivariate Bernoulli polynomials. Some error bounds and convergence rates of the combined operators are studied. A selection of numerical examples is presented to compare the performances of the obtained scheme. Furthermore, our method can be applied to time-dependent differential equations. Its advantage is that the algorithm is very simple and easy to implement. [ABSTRACT FROM AUTHOR]

We propose AAA rational approximation as a method for interpolating or approximating smooth functions from equispaced samples. Although it is always better to approximate from large numbers of samples if they are available, whether equispaced or not, this method often performs impressively even when the sampling grid is coarse. In most cases it gives more accurate approximations than other methods. We support this claim with a review and discussion of nine classes of existing methods in the light of general properties of approximation theory as well as the “impossibility theorem” for equispaced approximation. We make careful use of numerical experiments, which are summarized in a sequence of nine figures. Among our new contributions is the observation, summarized in Fig. 7, that methods such as polynomial least-squares and Fourier extension may be either exponentially accurate and exponentially unstable, or less accurate and stable, depending on implementation. [ABSTRACT FROM AUTHOR]

We consider new degrees of freedom for higher order differential forms on cubical meshes. The approach is inspired by the idea of Rapetti and Bossavit to define higher order Whitney forms and their degrees of freedom using small simplices. We show that higher order differential forms on cubical meshes can be defined analogously using small cubes and prove that these small cubes yield unisolvent degrees of freedom. Importantly, this approach is compatible with discrete exterior calculus and expands the framework to cover higher order methods on cubical meshes, complementing the earlier strategy based on simplices. [ABSTRACT FROM AUTHOR]

We describe a general approach for constructing a broad class of operators approximating high-dimensional curves based on geometric Hermite data. The geometric Hermite data consists of point samples and their associated tangent vectors of unit length. Extending the classical Hermite interpolation of functions, this geometric Hermite problem has become popular in recent years and has ignited a series of solutions in the 2D plane and 3D space. Here, we present a method for approximating curves, which is valid in any dimension. A basic building block of our approach is a Hermite average — a notion introduced in this paper. We provide an example of such an average and show, via an illustrative interpolating subdivision scheme, how the limits of the subdivision scheme inherit geometric properties of the average. Finally, we prove the convergence of this subdivision scheme, whose limit interpolates the geometric Hermite data and approximates the sampled curve. We conclude the paper with various numerical examples that elucidate the advantages of our approach. [ABSTRACT FROM AUTHOR]

Hermite interpolation is classically used to reconstruct smooth data when the function and its first order derivatives are available at certain nodes. If first order derivatives are not available, it is easy to set a system of equations imposing some regularity conditions at the data nodes in order to obtain them. This process leads to the construction of a Hermite spline. The problem of the described Hermite splines is that the accuracy is lost if the data contains singularities. The consequence is the appearance of oscillations, if there is a jump discontinuity in the function, that globally affects the accuracy of the spline, or the smearing of singularities, if the discontinuities are in the derivatives of the function. This paper is devoted to the construction and analysis of a new technique that allows for the computation of accurate first order derivatives of a function close to singularities using a Hermite spline. The idea is to correct the system of equations of the spline in order to attain the desired accuracy even close to the singularities. Once we have computed the first order derivatives with enough accuracy, a correction term is added to the Hermite spline in the intervals that contain a singularity. The aim is to reconstruct piecewise smooth functions with O (h 4) accuracy even close to the singularities. The process of adaption will require some knowledge about the position of the singularity and the jumps of the function and some of its derivatives at the singularity. The whole process can be used as a post-processing, where a correction term is added to the classical cubic Hermite spline. Proofs for the accuracy and regularity of the corrected spline and its derivatives are given. We also analyse the mechanism that eliminates the Gibbs phenomenon close to jump discontinuities in the function. The numerical experiments presented confirm the theoretical results obtained. [ABSTRACT FROM AUTHOR]

We present a new image scaling method both for downscaling and upscaling, running with any scale factor or desired size. The resized image is achieved by sampling a bivariate polynomial which globally interpolates the data at the new scale. The method's particularities lay in both the sampling model and the interpolation polynomial we use. Rather than classical uniform grids, we consider an unusual sampling system based on Chebyshev zeros of the first kind. Such optimal distribution of nodes permits to consider near-best interpolation polynomials defined by a filter of de la Vallée-Poussin type. The action ray of this filter provides an additional parameter that can be suitably regulated to improve the approximation. The method has been tested on a significant number of different image datasets. The results are evaluated in qualitative and quantitative terms and compared with other available competitive methods. The perceived quality of the resulting scaled images is such that important details are preserved, and the appearance of artifacts is low. Competitive quality measurement values, good visual quality, limited computational effort, and moderate memory demand make the method suitable for real-world applications. [ABSTRACT FROM AUTHOR]

In this paper, based on the Tikhonov regularized approximation formulae derived by An and Wu, we focus on their error estimates, the selection of the regularization parameter and the numerical stability. First, we quantify certain terms in the previous L 2 error bound and the uniform error bound. From these two error bounds, we know that the selection of the regularization parameter is important. Then, we employ two strategies, the balancing principle and Brezinski–Rodriguez–Seatzu estimators, to select the regularization parameter. Some numerical experiments are given to illustrate that both parameter choice strategies can select suitable regularization parameters. The two parameter choice strategies are also compared by testing some oscillatory functions. Finally, we study the numerical stability of Tikhonov regularized barycentric interpolation formula and Tikhonov regularized modified Lagrange interpolation formula. Both interpolation formulae are shown to be forward stable, and the Tikhonov regularized modified Lagrange interpolation formula is also backward stable. We give numerical examples to illustrate the theoretical results. [ABSTRACT FROM AUTHOR]

Finite elements with polynomial basis functions on the simplex with a symmetric distribution of nodes should have a unique polynomial representation. Unisolvence not only requires that the number of nodes equals the number of independent polynomials spanning a polynomial space of a given degree, but also that the Vandermonde matrix controlling their mapping to the Lagrange interpolating polynomials can be inverted. Here, a necessary condition for unisolvence is presented for polynomial spaces that have non-decreasing degrees when going from the edges and the various faces to the interior of the simplex. It leads to a proof of a conjecture on a necessary condition for unisolvence, requiring the node pattern to be the same as that of the regular simplex. [ABSTRACT FROM AUTHOR]

We propose an efficient mesh adaptive method for the numerical solution of time-dependent partial differential equations considered in the fixed space-time cylinder Ω × (0 , T) . We employ the space-time discontinuous Galerkin method which enables us to use different meshes at different time levels in a natural way. The mesh adaptive algorithm is based on control of the interpolation error in the L ∞ (0 , T ; L q (Ω)) -norm. The goal is to construct a sequence of conforming triangular meshes in such a way that the interpolation error bound is under a given tolerance and the number of degrees of freedom is minimal. The resulting grids consist of anisotropic mesh elements with varying polynomial approximation degrees with respect to space. We present a theoretical framework of this approach as well as several numerical examples demonstrating the accuracy, efficiency, and applicability of the method. [ABSTRACT FROM AUTHOR]

Numerical analysts and scientists working in applications often observe that once they improve their techniques to get a better accuracy, some instability of the evaluation 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 evaluation stability is unavoidable. It is an instance of a no free lunch theorem. Here, evaluation is the mathematical map that takes input data to output data. This is independent from the numerical routine that calculates the output. Therefore, evaluation stability is different from computational stability. 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. The trade-off principle bounds the product of two terms from below. The first is related to errors, and the second turns out to be related to evaluation instability. Under certain conditions satisfied for splines and kernel-based interpolation, both can be minimized. Then the lower bound is attained, and the error term is the inverse of the instability term. As a byproduct, it is shown that Kansa’s Unsymmetric Collocation Method sacrifices accuracy for improved evaluation stability, when compared to symmetric collocation. [ABSTRACT FROM AUTHOR]

In this work, our motivation is to design an impressive new numerical approximation on non-uniform grid points for the Caputo fractional derivative in time 0 C D t α with the order α ∈ (1,2). An adaptive high-order stable implicit difference scheme is developed for the time-fractional diffusion wave equations (TFDWEs) by using estimation of order O (N t α - 5) for the Caputo derivative in the time domain on non-uniform mesh and well-known second-order central difference approximation for estimating the spatial derivative on a uniform mesh. The designed algorithm allows one to build adaptive nature where the scheme is adjusted according to the behaviour of α in order to keep the numerical errors very small and converge to the solution very fast as compared to the previously investigated scheme. We rigorously analyze the local truncation errors, unconditional stability of the proposed method, and its convergence of (5 − α)-th order in time and second-order in space for all values of α ∈ (1,2). A reduced order technique is implemented by using moving mesh refinement and assemble with the derived scheme in order to improve the temporal accuracy at several starting time levels. Furthermore, the numerical stability of the derived adaptive scheme is verified by imposing random external noises. Some numerical tests are given to show that the numerical results are consistent with the theoretical results. [ABSTRACT FROM AUTHOR]

In this paper, a new Hermite interpolatory subdivision scheme for curve interpolation is introduced. The scheme is constructed from the Rational [3/2] Bernstein Bezier polynomial. We call it the [3/2]-scheme. The limit function of the [3/2]-scheme interpolates both the function values and their derivatives. The proposed scheme has three shape parameters w 0 , w 1 and w 2 . It is shown that if w 1 = w 0 + w 2 2 , then the [3/2]-scheme reproduces linear polynomial and is C 1 provided w 0 and w 2 lie in a region of convergence. The scheme also satisfies the shape preserving properties, i.e., monotonicity and convexity. We also compare the [3/2]-scheme with other existing schemes like the [2/2]-scheme and the Merrien scheme introduced recently. An error analysis shows that the [3/2]-scheme is better than the [2/2]-scheme and the Merrien scheme. Further, it is observed that in case w 0 = w 1 = w 2 , the [3/2]-scheme reduces to the Merrien scheme. [ABSTRACT FROM AUTHOR]

Quadrature formulas (QFs) based on radial basis functions (RBFs) have become an essential tool for multivariate numerical integration of scattered data. Although numerous works have been published on RBF-QFs, their stability theory can still be considered as underdeveloped. Here, we strive to pave the way towards a more mature stability theory for global and function-independent RBF-QFs. In particular, we prove stability of these for compactly supported RBFs under certain conditions on the shape parameter and the data points. As an alternative to changing the shape parameter, we demonstrate how the least-squares approach can be used to construct stable RBF-QFs by allowing the number of data points used for numerical integration to be larger than the number of centers used to generate the RBF approximation space. Moreover, it is shown that asymptotic stability of many global RBF-QFs is independent of polynomial terms, which are often included in RBF approximations. While our findings provide some novel conditions for stability of global RBF-QFs, the present work also demonstrates that there are still many gaps to fill in future investigations. [ABSTRACT FROM AUTHOR]

Data-dependent greedy algorithms in kernel spaces are known to provide fast converging interpolants, while being extremely easy to implement and efficient to run. Despite this experimental evidence, no detailed theory has yet been presented. This situation is unsatisfactory, especially when compared to the case of the data-independent P-greedy algorithm, for which optimal convergence rates are available, despite its performances being usually inferior to the ones of target data-dependent algorithms. In this work, we fill this gap by first defining a new scale of greedy algorithms for interpolation that comprises all the existing ones in a unique analysis, where the degree of dependency of the selection criterion on the functional data is quantified by a real parameter. We then prove new convergence rates where this degree is taken into account, and we show that, possibly up to a logarithmic factor, target data-dependent selection strategies provide faster convergence. In particular, for the first time we obtain convergence rates for target data adaptive interpolation that are faster than the ones given by uniform points, without the need of any special assumption on the target function. These results are made possible by refining an earlier analysis of greedy algorithms in general Hilbert spaces. The rates are confirmed by a number of numerical examples. [ABSTRACT FROM AUTHOR]

We construct new families of direct serendipity and direct mixed finite elements on general planar, strictly convex polygons that are H1 and H(div) conforming, respectively, and possess optimal order of accuracy for any order. They have a minimal number of degrees of freedom subject to the conformity and accuracy constraints. The name arises because the shape functions are defined directly on the physical elements, i.e., without using a mapping from a reference element. The finite element shape functions are defined to be the full spaces of scalar or vector polynomials plus a space of supplemental functions. The direct serendipity elements are the precursors of the direct mixed elements in a de Rham complex. The convergence properties of the finite elements are shown under a regularity assumption on the shapes of the polygons in the mesh, as well as some mild restrictions on the choices one can make in the construction of the supplemental functions. Numerical experiments on various meshes exhibit the performance of these new families of finite elements. [ABSTRACT FROM AUTHOR]

When one solves differential equations by a spectral method, it is often convenient to shift from Chebyshev polynomials Tn(x) with coefficients an to modified basis functions that incorporate the boundary conditions. For homogeneous Dirichlet boundary conditions, u(±1) = 0, popular choices include the "Chebyshev difference basis" ςn(x) ≡ Tn+2(x)−Tn(x) with coefficients here denoted by bn and the "quadratic factor basis" ϱn(x) ≡ (1 − x2)Tn(x) with coefficients cn. If u(x) is weakly singular at the boundary, then the coefficients an decrease proportionally to O (A (n) / n κ) for some positive constant κ, where A(n) is a logarithm or a constant. We prove that the Chebyshev difference coefficients bn decrease more slowly by a factor of 1/n while the quadratic factor coefficients cn decrease more slowly still as O (A (n) / n κ − 2 ) . The error for the unconstrained Chebyshev series, truncated at degree n = N, is O (| A (N) | / N κ) in the interior, but is worse by one power of N in narrow boundary layers near each of the endpoints. Despite having nearly identical error norms in interpolation, the error in the Chebyshev basis is concentrated in boundary layers near both endpoints, whereas the error in the quadratic factor and difference basis sets is nearly uniformly oscillating over the entire interval in x. Meanwhile, for Chebyshev polynomials, the values of their derivatives at the endpoints are O (N 2) , but only O (N) for the difference basis. Furthermore, we give the asymptotic coefficients and rigorous error estimates of the approximations in these three bases, solved by the least squares method. We also find an interesting fact that on the face of it, the aliasing error is regarded as a bad thing; actually, the error norm associated with the downward curving spectral coefficients decreases even faster than the error norm of infinite truncation. But the premise is under the same basis, and when involving different bases, it may not be established yet. [ABSTRACT FROM AUTHOR]

We present precise anisotropic interpolation error estimates for smooth functions using a new geometric parameter and derive inverse inequalities on anisotropic meshes. In our theory, the interpolation error is bounded in terms of the diameter of a simplex and the geometric parameter. Imposing additional assumptions makes it possible to obtain anisotropic error estimates. This paper also includes corrections to an error in Theorem 2 of our previous paper, "General theory of interpolation error estimates on anisotropic meshes" (Japan Journal of Industrial and Applied Mathematics, 38 (2021) 163–191). [ABSTRACT FROM AUTHOR]

We present a precise anisotropic interpolation error estimate for the Morley finite element method (FEM) and apply it to fourth-order elliptic equations. We do not impose the shape-regularity mesh condition in the analysis. Anisotropic meshes can be used for this purpose. The main contributions of this study include providing a new proof of the term consistency. This enables us to obtain an anisotropic consistency error estimate. The core idea of the proof involves using the relationship between the Raviart-Thomas and Morley finite-element spaces. Our results indicate optimal convergence rates and imply that the modified Morley FEM may be effective for errors. [ABSTRACT FROM AUTHOR]

Under appropriate technical assumptions, the simple-loop theory allows to derive various types of asymptotic expansions for the eigenvalues of Toeplitz matrices generated by a function f. Independently and under the milder hypothesis that f is even and monotone over [0,π], matrix-less algorithms have been developed for the fast eigenvalue computation of large Toeplitz matrices, within a linear complexity in the matrix order: behind the high efficiency of such algorithms there are the expansions predicted by the simple-loop theory, combined with the extrapolation idea. Here we focus our attention on a change of variable, followed by the asymptotic expansion of the new variable, and we adapt the matrix-less algorithm to the considered new setting. Numerical experiments show a higher precision (till machine precision) and the same linear computation cost, when compared with the matrix-less procedures already presented in the relevant literature. Among the advantages, we concisely mention the following: (a) when the coefficients of the simple-loop function are analytically known, the algorithm computes them perfectly; (b) while the proposed algorithm is better or at worst comparable to the previous ones for computing the inner eigenvalues, it is vastly better for the computation of the extreme eigenvalues; a mild deterioration in the quality of the numerical experiments is observed when dense Toeplitz matrices are considered, having generating function of low smoothness and not satisfying the simple-loop assumptions. [ABSTRACT FROM AUTHOR]

In neuroscience, the distribution of a decision time is modelled by means of a one-dimensional Fokker–Planck equation with time-dependent boundaries and space-time-dependent drift. Efficient approximation of the solution to this equation is required, e.g., for model evaluation and parameter fitting. However, the prescribed boundary conditions lead to a strong singularity and thus to slow convergence of numerical approximations. In this article we demonstrate that the solution can be related to the solution of a parabolic PDE on a rectangular space-time domain with homogeneous initial and boundary conditions by transformation and subtraction of a known function. We verify that the solution of the new PDE is indeed more regular than the solution of the original PDE and proceed to discretize the new PDE using a space-time minimal residual method. We also demonstrate that the solution depends analytically on the parameters determining the boundaries as well as the drift. This justifies the use of a sparse tensor product interpolation method to approximate the PDE solution for various parameter ranges. The predicted convergence rates of the minimal residual method and that of the interpolation method are supported by numerical simulations. [ABSTRACT FROM AUTHOR]

The barycentric forms of polynomial and rational interpolation have recently gained popularity, because they can be computed with simple, efficient, and numerically stable algorithms. In this paper, we show more generally that the evaluation of any function that can be expressed as r (x) = ∑ i = 0 n a i (x) f i / ∑ j = 0 m b j (x) in terms of data values f i and some functions a i and b j for i = 0 , ... , n and j = 0 , ⋯ , m with a simple algorithm that first sums up the terms in the numerator and the denominator, followed by a final division, is forward and backward stable under certain assumptions. This result includes the two barycentric forms of rational interpolation as special cases. Our analysis further reveals that the stability of the second barycentric form depends on the Lebesgue constant associated with the interpolation nodes, which typically grows with n, whereas the stability of the first barycentric form depends on a similar, but different quantity, that can be bounded in terms of the mesh ratio, regardless of n. We support our theoretical results with numerical experiments. [ABSTRACT FROM AUTHOR]