Error analysis of efficient evaluation algorithms for tensor product surfaces

Backward stability of the Casteljau algorithm and two more efficient algorithms for polynomial tensor product surfaces with interest in CAGD is shown. The conditioning of the corresponding bases are compared. These algorithms are also compared with the corresponding Horner algorithm and their higher accuracy is shown. A running error analysis of the algorithms is also carried out providing algorithms which calculate “a posteriori” sharp error bounds simultaneously to the evaluation of the surface without increasing significantly the computational cost.     

Fast evaluation of nonlinear functionals of tensor product wavelet expansions

For a nonlinear functional f, and a function u from the span of a set of tensor product interpolets, it is shown how to compute the interpolant of f (u) from the span of this set of tensor product interpolets in linear complexity, assuming that the index set has a certain multiple tree structure. Applications are found in the field of (adaptive) tensor product solution methods for semilinear operator equations by collocation methods, or after transformations between the interpolet and (bi-) orthogonal wavelet bases, by Galerkin methods.     

Twisted tensor product codes

We present two families of constacyclic linear codes with large automorphism groups. The codes are obtained from the twisted tensor product construction.      

Locally tensor product functions

We present in this paper a family of functions which are tensor product functions in subdomains, while not having the usual drawback of functions which are tensor product functions in the whole domain. With these functions we can add more points in some region without adding points on lines parallel to the axes. These functions are linear combinations of tensor product polynomial B-splines, and the knots of different B-splines are less connected together than with usual polynomial B-splines. Approximation of functions, or data, with such functions gives satisfactory results, as shown by numerical experimentation. AMS subject classification 41A15, 41A63, 65Dxx     

An equivariant tensor product on Mackey functors

For all subgroups H of a cyclic p-group G we define norm functors that build a G-Mackey functor from an H-Mackey functor. We give an explicit construction of these functors in terms of generators and relations based solely on the intrinsic, algebraic properties of Mackey functors and Tambara functors. We use these norm functors to define a monoidal structure on the category of Mackey functors where Tambara functors are the commutative ring objects.     

Analysis of tensor product multigrid

We consider anisotropic second order elliptic boundary value problems in two dimensions, for which the anisotropy is exactly aligned with the coordinate axes. This includes cases where the operator features a singular perturbation in one coordinate direction, whereas its restriction to the other direction remains neatly elliptic. Most prominently, such a situation arises when polar coordinates are introduced.The common multigrid approach to such problems relies on line relaxation in the direction of the singular perturbation combined with semi-coarsening in the other direction. Taking the idea from classical Fourier analysis of multigrid, we employ eigenspace techniques to separate the coordinate directions. Thus, convergence of the multigrid method can be examined by looking at one-dimensional operators only. In a tensor product Galerkin setting, this makes it possible to confirm that the convergence rates of the multigrid V-cycle are bounded independently of the number of grid levels involved. In addition, the estimates reveal that convergence is also robust with respect to a singular perturbation in one coordinate direction.Finally, we supply numerical evidence that the algorithm performs satisfactorily in settings more general than those covered by the proof.     

The tensor product of f-algebras

In this paper we prove that the Fremlin tensor product of two f-algebras can be endowed with an f-algebra structure and satisfies an appropriate universal property. In particular, the Riesz tensor product of C(X) and C(Y), where X and Y are topological spaces, is an f-subalgebra of C(X × Y).     

Simplicial join via tensor product

Let be a field, and M and N two finitely generated graded modules over standard graded -algebras A and B, respectively. We will study generalized, sequentially, almost, and approximately Cohen–Macaulay as well as clean, and pretty clean properties of the -module through the corresponding properties of M and N. The behavior of these properties with respect to the simplicial join of two simplicial (multi)complexes will be revealed as corollaries.     

Lagrange interpolation on subgrids of tensor product grids

This note shows that a wide class of algebraically motivated constructions for Lagrange interpolation polynomials always yields a tensor product interpolation space as long as the nodes form a tensor product grid or a lower subset thereof.

     

Differential geometry of tensor product immersions

In this article we obtain the best possible estimates of the type number of tensor product immersions and investigate tensor product immersions with lowest possible type. Several classification theorems in this respect are then proved.     

Nonlinear tensor product approximation of functions

We are interested in approximation of a multivariate function $f(x_1,\dots ,x_d)$ by linear combinations of products $u^1\left(x_1\right)\cdots u^d\left(x_d\right)$ of univariate functions $u^i\left(x_i\right)$, $i=1,\dots ,d$. In the case $d=2$ it is the classical problem of bilinear approximation. In the case of approximation in the $L_2$ space the bilinear approximation problem is closely related to the problem of singular value decomposition (also called Schmidt expansion) of the corresponding integral operator with the kernel $f(x_1,x_2)$. There are known results on the rate of decay of errors of best bilinear approximation in $L_p$ under different smoothness assumptions on $f$. The problem of multilinear approximation (nonlinear tensor product approximation) in the case $d\ge 3$ is more difficult and much less studied than the bilinear approximation problem. We will present results on best multilinear approximation in $L_p$ under mixed smoothness assumption on $f$.     

Syzygies and tensor product of modules

We give an application of the New Intersection Theorem and prove the following: let $R$ be a local complete intersection ring of codimension $c$ and let $M$ and $N$ be nonzero finitely generated $R$ -modules. Assume $n$ is a nonnegative integer and that the tensor product $M\otimes _{R}N$ is an $(n+c)$ th syzygy of some finitely generated $R$ -module. If ${{\mathrm{Tor}}}^{R}_{>0}(M,N)=0$ , then both $M$ and $N$ are $n$ th syzygies of some finitely generated $R$ -modules.     

Families of twisted tensor product codes

Using geometric properties of the variety ${\mathcal V_{r,t}}$ , the image under the Grassmannian map of a Desarguesian (t ? 1)-spread of PG(rt ? 1, q), we introduce error correcting codes related to the twisted tensor product construction, producing several families of constacyclic codes. We determine the precise parameters of these codes and characterise the words of minimum weight.     

Ergodic limits of tensor product semigroups

Let e^tSande^?tT be (C₀)-semigroups on a Banach space X. Their tensor product $\text{L}$(t) is defined by $\text{L}$(t)A = e^tSAe^tT (A?B(X)) and has the generator Δ formally of the form ΔA = SA ? AT. Under the assumption that {$\text{L}$(t); t ? 0} is bounded, we investigate the Abel limit and the Cesàro limit of $\text{L}$(t)A at ∞. If $\text{gW}{}_{\mathrm{s}}\text{[Ω}{}_{\mathrm{u}}\text{]}$ denotes the set of operators A for which the Abel limit P_s(A) [resp. P_u(A)] exists in the strong [resp. uniform] operator topology, then

$\text{N(Δ)⊕}\text{R(Δ)}\text{= ω}{}_{\mathrm{u}}\text{? ω}{}_{\mathrm{s}}\text{? N(Δ) +}\text{R(Δ)}$

and the limit defines a projection P_s[P_u] from $\text{Ω}{}_{\mathrm{s}}$ [resp. $\text{Ω}{}_{\mathrm{u}}$] onto N(Δ) with N(Δ) with $\text{R(Δ)}\text{= N(P}{}_{\mathrm{u}}\text{) ? N(P}{}_{\mathrm{u}}\text{) ?}\text{R(Δ)}$. If, in addition, S and T are Hilbert space normal operators such that _gq(S) ∩ _gq(T) ≠ φ, then $\text{Ω}{}_{\mathrm{u}}$ contains all compact operators.