659 research outputs found
Bivariate hierarchical Hermite spline quasi--interpolation
Spline quasi-interpolation (QI) is a general and powerful approach for the
construction of low cost and accurate approximations of a given function. In
order to provide an efficient adaptive approximation scheme in the bivariate
setting, we consider quasi-interpolation in hierarchical spline spaces. In
particular, we study and experiment the features of the hierarchical extension
of the tensor-product formulation of the Hermite BS quasi-interpolation scheme.
The convergence properties of this hierarchical operator, suitably defined in
terms of truncated hierarchical B-spline bases, are analyzed. A selection of
numerical examples is presented to compare the performances of the hierarchical
and tensor-product versions of the scheme
Adaptive meshless centres and RBF stencils for Poisson equation
We consider adaptive meshless discretisation of the Dirichlet problem for Poisson equation based on numerical differentiation stencils obtained with the help of radial basis functions. New meshless stencil selection and adaptive refinement algorithms are proposed in 2D. Numerical experiments show that the accuracy of the solution is comparable with, and often better than that achieved by the mesh-based adaptive finite element method
Macro-element interpolation on tensor product meshes
A general theory for obtaining anisotropic interpolation error estimates for
macro-element interpolation is developed revealing general construction
principles. We apply this theory to interpolation operators on a macro type of
biquadratic finite elements on rectangle grids which can be viewed as a
rectangular version of the Powell-Sabin element. This theory also shows
how interpolation on the Bogner-Fox-Schmidt finite element space (or higher
order generalizations) can be analyzed in a unified framework. Moreover we
discuss a modification of Scott-Zhang type giving optimal error estimates under
the regularity required without imposing quasi uniformity on the family of
macro-element meshes used. We introduce and analyze an anisotropic
macro-element interpolation operator, which is the tensor product of
one-dimensional macro interpolation and Lagrange interpolation.
These results are used to approximate the solution of a singularly perturbed
reaction-diffusion problem on a Shishkin mesh that features highly anisotropic
elements. Hereby we obtain an approximation whose normal derivative is
continuous along certain edges of the mesh, enabling a more sophisticated
analysis of a continuous interior penalty method in another paper
Recommended from our members
Smooth parametric surfaces and n-sided patches
The theory of 'geometric continuity' within the subject of CAGD is reviewed. In particular, we are concerned with how parametric surface patches for CAGD can be pieced together to form a smooth Ck surface. The theory is applied to the problem of filling an n-sided hole occurring within a smooth rectangular patch complex. A number of solutions to this problem are surveyed
From discretization to regularization of composite discontinuous functions
Discontinuities between distinct regions, described by different equation sets, cause difficulties for PDE/ODE solvers. We present a new algorithm that eliminates integrator discontinuities through regularizing discontinuities. First, the algorithm determines the optimum switch point between two functions spanning adjacent or overlapping domains. The optimum switch point is determined by searching for a “jump point” that minimizes a discontinuity between adjacent/overlapping functions. Then, discontinuity is resolved using an interpolating polynomial that joins the two discontinuous functions.
This approach eliminates the need for conventional integrators to either discretize and then link discontinuities through generating interpolating polynomials based on state variables or to reinitialize state variables when discontinuities are detected in an ODE/DAE system. In contrast to conventional approaches that handle discontinuities at the state variable level only, the new approach tackles discontinuity at both state variable and the constitutive equations level. Thus, this approach eliminates errors associated with interpolating polynomials generated at a state variable level for discontinuities occurring in the constitutive equations.
Computer memory space requirements for this approach exponentially increase with the dimension of the discontinuous function hence there will be limitations for functions with relatively high dimensions. Memory availability continues to increase with price decreasing so this is not expected to be a major limitation
The Argyris isogeometric space on unstructured multi-patch planar domains
Multi-patch spline parametrizations are used in geometric design and
isogeometric analysis to represent complex domains. We deal with a particular
class of planar multi-patch spline parametrizations called
analysis-suitable (AS-) multi-patch parametrizations (Collin,
Sangalli, Takacs; CAGD, 2016). This class of parametrizations has to satisfy
specific geometric continuity constraints, and is of importance since it allows
to construct, on the multi-patch domain, isogeometric spaces with optimal
approximation properties. It was demonstrated in (Kapl, Sangalli, Takacs; CAD,
2018) that AS- multi-patch parametrizations are suitable for modeling
complex planar multi-patch domains.
In this work, we construct a basis, and an associated dual basis, for a
specific isogeometric spline space over a given AS-
multi-patch parametrization. We call the space the Argyris
isogeometric space, since it is across interfaces and at all
vertices and generalizes the idea of Argyris finite elements to tensor-product
splines. The considered space is a subspace of the entire
isogeometric space , which maintains the reproduction
properties of traces and normal derivatives along the interfaces. Moreover, it
reproduces all derivatives up to second order at the vertices. In contrast to
, the dimension of does not depend on the domain
parametrization, and admits a basis and dual basis which possess
a simple explicit representation and local support.
We conclude the paper with some numerical experiments, which exhibit the
optimal approximation order of the Argyris isogeometric space and
demonstrate the applicability of our approach for isogeometric analysis
Multivariate Splines and Algebraic Geometry
Multivariate splines are effective tools in numerical analysis and approximation theory. Despite an extensive literature on the subject, there remain open questions in finding their dimension, constructing local bases, and determining their approximation power. Much of what is currently known was developed by numerical analysts, using classical methods, in particular the so-called Bernstein-B´ezier techniques. Due to their many interesting structural properties, splines have become of keen interest to researchers in commutative and homological algebra and algebraic geometry. Unfortunately, these communities have not collaborated much. The purpose of the half-size workshop is to intensify the interaction between the different groups by bringing them together. This could lead to essential breakthroughs on several of the above problems
Spline approximation of a random process with singularity
Let a continuous random process defined on be -smooth,
and have an isolated
singularity point at . In addition, let be locally like a -fold
integrated -fractional Brownian motion for all non-singular points. We
consider approximation of by piecewise Hermite interpolation splines with
free knots (i.e., a sampling design, a mesh). The approximation performance
is measured by mean errors (e.g., integrated or maximal quadratic mean errors).
We construct a sequence of sampling designs with asymptotic approximation rate
for the whole interval.Comment: 16 pages, 2 figure typos and references corrected, revised classes
definition, results unchange
- …