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]

In this study, we have managed to model the energy surface of the aminophenol derivative azomethine molecule mathematically depending on two torsion angles SC1(C9C10C12N14) and SC2(C2C1C6C11). For this purpose, firstly, discrete data obtained from Density Functional Theory calculations have been converted into continuous data with the help of the Fuzzy Inference System. Thus, it is possible to calculate energy values for untested data, which are very costly in terms of time to obtain with other methods/experiments. Then, the continuous and non-smooth surface obtained from the fuzzy inference system and representing the energy values of the molecule has been transformed into a differentiable surface with the help of Bezier surfaces. Thus, an objective function has been obtained in which global optimization methods based on the derivative (or gradient) operator could be used. [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]

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]

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]

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]

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]

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]

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]

A novel approach to zipper fractal interpolation theory for functions of several variables is presented. Multivariate zipper fractal functions are constructed and then perturbed through free choices of base functions, scaling functions, and a binary matrix called signature to obtain their zipper α-fractal versions. In particular, we propose a multivariate Bernstein zipper fractal function and study its coordinate-wise monotonicity which depends on the values of signature. We derive bounds for the graph of a multivariate zipper fractal function by imposing conditions on the scaling factors and the Hölder exponent of the associated germ function and base function. The box dimension result for multivariate Bernstein zipper fractal function is derived. Finally, we study some constrained approximation properties for multivariate zipper Bernstein fractal functions. [ABSTRACT FROM AUTHOR]

In this article, we propose generalized forms of three well-known fractional numerical methods namely Euler, Runge–Kutta 2-step, and Runge–Kutta 4-step, respectively. The new versions we provide of these methods are derived by utilizing a non-uniform grid which is slightly different from previous versions of these algorithms. A new generalized form of the well-known Caputo-type fractional derivative is used to derive the results. All necessary analyses related to the stability, convergence, and error bounds are also provided. The precision of all simulated results is justified by performing multiple numerical experiments, with some meaningful problems solved by implementing the code in Mathematica. Finally, we give a brief discussion on the simulated results which shows that the generalized methods are novel, effective, reliable, and very easy to implement. [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]

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]

The main concern of this paper is to develop a new class of A-stable fourth-order numerical scheme for solving initial value problems. The idea is an evolutionary and heuristic approach; using the Grasshopper optimization, along with the Hermite interpolation for stages, we obtain a class of A-stable methods. Four types of weighting rules are introduced for the current formulation. The fundamental weighting rule (FWR) is the most important rule, which emphasizes on the symmetric and central structure of the method. A systematic strategy is proposed to obtain the FWR based on swarm intelligence and regression. The new techniques are called α I -(e + i)P, where e and i are the number of terminal and internal points, respectively. The numerical experiments demonstrate the reasonable behaviour of the algorithms on several test problems from different applications. Finally, we find that the new formulas are well suited for long time behaviour of the time-history analysis in earthquake engineering. [ABSTRACT FROM AUTHOR]

In this study, the numerical solutions of the modified Fornberg–Whitham (mFW) equation, which describes immigration of the solitary wave and peakon waves with discontinuous first derivative at the peak, have been obtained by the collocation finite element method using quintic trigonometric B-spline bases. Although there are solutions of this equation by semi-analytical and analytical methods in the literature, there are very few studies on the solution of the equation by numerical methods. Any linearization technique has not been used while applying the method. The stability analysis of the applied method is examined by the von-Neumann Fourier series method. To show the performance of the method, we have considered three test problems with nonhomogeneous boundary conditions having analytical solutions. The error norms L2 and L∞ are calculated to demonstrate the accuracy and efficiency of the presented numerical scheme. [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]

Digitized handwritten documents have been used for various purposes, including age detection, a crucial area of research in fields like forensic investigation and medical diagnosis. Automated age recognition is deemed to be a difficult task, due to the great degree of similarity and overlap across people's handwriting. Consequently, the efficiency of the classification system is determined by extracting pertinent features from handwritten documents. This research proposes a set of age-related features suggested by a graphologist to classify handwritten documents into two age groups: youth adult and mature adult. The extracted features are: slant irregularity (SI), pen pressure irregularity (PPI), text lines irregularity (TLI) and the percentage of black and white pixels (PWB). Support Vector Machines (SVM) and Neural Network (NN) classifiers have been used to train, validate and test the proposed approach using two different datasets: the FSHS and the Khatt datasets. When applied to the FSHS dtaset using SVM and NN approaches, the proposed method resulted in a classification accuracy of 71% and 63.5%, respectively. Meanwhile, when applied to the Khatt dataset, our method outperformed state-of-the-art methods with a classification accuracy of 65.2% and 67% utilising SVM and NN classifiers, respectively. These are the best rates available right now in this field. [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]

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]

We present precise Raviart–Thomas interpolation error estimates on anisotropic meshes. The novel aspect of our theory is the introduction of a new geometric parameter of simplices. It is possible to obtain new anisotropic Raviart–Thoma error estimates using the parameter. We also include corrections to an error in “General theory of interpolation error estimates on anisotropic meshes” (Japan Journal of Industrial and Applied Mathematics, 38 (2021) 163-191), in which Theorem 3 was incorrect. [ABSTRACT FROM AUTHOR]

In this work, a class of interpolation and smoothing exponential splines with respect to a sixth order differential operator with two parameters is constructed. All the square matrices involved in the construction are proved to be tridiagonal symmetric and diagonally dominant, which results in algorithms for computing this class of exponential splines. The obtained splines have C 4 continuity and are the minimum solution of the combination of interpolation and a generalized smoothing energy integral. The performances of the resulting splines in financial data from the S &P500 index and the effect of fitting multi-exponential decay data are given. Numerical experiments show that the resulting splines have more freedom to adjust the shape and control the energy of the curves and perform better than previous methods in fitting multi-exponential decay data. And cross validation and generalized cross validation for determining an appropriate smoothing parameter are also developed. [ABSTRACT FROM AUTHOR]

We propose to extend results on the interpolation theory for scalar functions to the case of differential k-forms. More precisely, we consider the interpolation of fields in P r - Λ k (T) , the finite element spaces of trimmed polynomial k-forms of arbitrary degree r ≥ 1 , from their weights, namely their integrals on k-chains. These integrals have a clear physical interpretation, such as circulations along curves, fluxes across surfaces, densities in volumes, depending on the value of k. In this work, for k = 1 , we rely on the flexibility of the weights with respect to their geometrical support, to study different sets of 1-chains in T for a high order interpolation of differential 1-forms, constructed starting from "good" sets of nodes for a high order multi-variate polynomial representation of scalar fields, namely 0-forms. We analyse the growth of the generalized Lebesgue constant with the degree r and preliminary numerical results for edge elements support the nonuniform choice, in agreement with the well-known nodal case. [ABSTRACT FROM AUTHOR]

In this paper, the solution to a bivariate Appell interpolation problem proposed in a previous work is given. Bounds of the truncation error are considered. Ten new interpolants for real, regular, bivariate functions are constructed. Numerical examples and comparisons with bivariate Bernstein polynomials are considered. [ABSTRACT FROM AUTHOR]

In this paper, we address the problem of constructing C 2 cubic spline functions on a given arbitrary triangulation T . To this end, we endow every triangle of T with a Wang–Shi macro-structure. The C 2 cubic space on such a refined triangulation has a stable dimension and optimal approximation power. Moreover, any spline function in this space can be locally built on each of the macro-triangles independently via Hermite interpolation. We provide a simplex spline basis for the space of C 2 cubics defined on a single macro-triangle which behaves like a Bernstein/B-spline basis over the triangle. The basis functions inherit recurrence relations and differentiation formulas from the simplex spline construction, they form a nonnegative partition of unity, they admit simple conditions for C 2 joins across the edges of neighboring triangles, and they enjoy a Marsden-like identity. Also, there is a single control net to facilitate control and early visualization of a spline function over the macro-triangle. Thanks to these properties, the complex geometry of the Wang–Shi macro-structure is transparent to the user. Stable global bases for the full space of C 2 cubics on the Wang–Shi refined triangulation T are deduced from the local simplex spline basis by extending the concept of minimal determining sets. [ABSTRACT FROM AUTHOR]

This paper deals with the construction of the Algebraic Trigonometric Pythagorean Hodograph (ATPH) cubic-like Hermite interpolant. A characterization of solutions according to the tangents at both ends and a global free shape parameter α is performed. Since this degree of freedom can be used for adjustments, we study how the curve evolves with respect to α. Several examples illustrating the construction process and a simple fitting method to determine the unique ATPH curve passing through a given point are proposed. [ABSTRACT FROM AUTHOR]

An inverse central ordering of the nodes is proposed for the Newton interpolation formula. This ordering may improve the stability for certain distributions of nodes. For equidistant nodes, an upper bound of the conditioning is provided. This bound is close to the bound of the conditioning in the Lagrange interpolation formula, whose conditioning is the lowest. This ordering is related to a pivoting strategy of a matrix elimination procedure called Neville elimination. The results are illustrated with examples. [ABSTRACT FROM AUTHOR]

Learning mappings of data on manifolds is an important topic in contemporary machine learning, with applications in astrophysics, geophysics, statistical physics, medical diagnosis, biochemistry, and 3D object analysis. This paper studies the problem of learning real-valued functions on manifolds through filtered hyperinterpolation of input–output data pairs where the inputs may be sampled deterministically or at random and the outputs may be clean or noisy. Motivated by the problem of handling large data sets, it presents a parallel data processing approach which distributes the data-fitting task among multiple servers and synthesizes the fitted sub-models into a global estimator. We prove quantitative relations between the approximation quality of the learned function over the entire manifold, the type of target function, the number of servers, and the number and type of available samples. We obtain the approximation rates of convergence for distributed and non-distributed approaches. For the non-distributed case, the approximation order is optimal. [ABSTRACT FROM AUTHOR]

We propose a fast solver for the variable-order (VO) time-fractional diffusion equation. Due to the impact of the time-dependent VO function, the resulting coefficient matrix of the large linear system assembling discrete equations of all time levels is a block lower triangular matrix without the block Toeplitz structure. Here, we approximate the off-diagonal blocks by low-rank matrices based on the polynomial interpolation, which can be constructed in O (M log 2 M) operations with the same number storage requirement, where M is the number of time steps. Furthermore, a divide-and-conquer method is developed to fast solve the approximated linear system. The proposed solver can be implemented in O (N M log 2 M) complexity with N being the degree of freedom in space. The accuracy of approximation is theoretically studied, and the stability and convergence of the proposed fast method are also investigated. Numerical experiments are carried out to exemplify the accuracy and efficiency of the proposed method. [ABSTRACT FROM AUTHOR]

Panel-based, kernel-split quadrature is currently one of the most efficient methods available for accurate evaluation of singular and nearly singular layer potentials in two dimensions. However, it can fail completely for the layer potentials belonging to the modified Helmholtz, modified biharmonic, and modified Stokes equations. These equations depend on a parameter, denoted α, and kernel-split quadrature loses its accuracy rapidly when this parameter grows beyond a certain threshold. This paper describes an algorithm that remedies this problem, using per-target adaptive sampling of the source geometry. The refinement is carried out through recursive bisection, with a carefully selected rule set. This maintains accuracy for a wide range of the parameter α, at an increased cost that scales as log α . Using this algorithm allows kernel-split quadrature to be both accurate and efficient for a much wider range of problems than previously possible. [ABSTRACT FROM AUTHOR]

Classical rational interpolation usually enjoys better approximation properties than polynomial interpolation because it avoids wild oscillations and exhibits exponential convergence for approximating analytic functions. We develop a rational interpolation operator, which not only preserves the advantage of classical rational interpolation, but also has a finite Lebesgue constant. In particular, it is convergent for approximating any continuous function, and the convergence rate of the interpolants approximating a function is obtained using the modulus of continuity. [ABSTRACT FROM AUTHOR]

In this paper we analyze a greedy procedure to approximate a linear functional defined in a reproducing kernel Hilbert space by nodal values. This procedure computes a quadrature rule which can be applied to general functionals. For a large class of functionals, that includes integration functionals and other interesting cases, but does not include differentiation, we prove convergence results for the approximation by means of quasi-uniform and greedy points which generalize in various ways several known results. A perturbation analysis of the weights and node computation is also discussed. Beyond the theoretical investigations, we demonstrate numerically that our algorithm is effective in treating various integration densities, and that it is even very competitive when compared to existing methods for Uncertainty Quantification. [ABSTRACT FROM AUTHOR]

In this work, we construct a class of Hermite-type interpolation basis functions based on the sixth-order ordinary differential equation S (6) (t) - τ 4 S (2) (t) = 0 . Using them, we propose a kind of C 2 tension interpolation splines with a local tension parameter τ i . For C 2 interpolation, the given interpolant has O (h 2) convergence. Some applications of the C 2 tension interpolation splines on the construction of interest rate term structure in Chinese financial market are given. Moreover, a kind of generalized non-uniform B-splines of the space spanned by span 1 , t , ... , t n - 4 , sin (τ t) , cos (τ t) , sinh (τ t) , cosh (τ t) is constructed. [ABSTRACT FROM AUTHOR]

The algebraic characterization of dual univariate interpolating subdivision schemes is investigated. Specifically, we provide a constructive approach for finding dual univariate interpolating subdivision schemes based on the solutions of certain associated polynomial equations. The proposed approach also makes it possible to identify conditions for the existence of the sought schemes. [ABSTRACT FROM AUTHOR]

We analyze a time-delay Caputo-type fractional mathematical model containing the infection rate of Beddington–DeAngelis functional response to study the structure of a vector-borne plant epidemic. We prove the unique global solution existence for the given delay mathematical model by using fixed point results. We use the Adams–Bashforth–Moulton P-C algorithm for solving the given dynamical model. We give a number of graphical interpretations of the proposed solution. A number of novel results are demonstrated from the given practical and theoretical observations. By using 3-D plots we observe the variations in the flatness of our plots when the fractional order varies. The role of time delay on the proposed plant disease dynamics and the effects of infection rate in the population of susceptible and infectious classes are investigated. The main motivation of this research study is examining the dynamics of the vector-borne epidemic in the sense of fractional derivatives under memory effects. This study is an example of how the fractional derivatives are useful in plant epidemiology. The application of Caputo derivative with equal dimensionality includes the memory in the model, which is the main novelty of this study. [ABSTRACT FROM AUTHOR]