Commuting extensions and cubature formulae

Based on a novel point of view on 1-dimensional Gaussian quadrature, we present a new approach to d-dimensional cubature formulae. It is well known that the nodes of 1-dimensional Gaussian quadrature can be computed as eigenvalues of the so-called Jacobi matrix. The d-dimensional analog is that cubature nodes can be obtained from the eigenvalues of certain mutually commuting matrices. These are obtained by extending (adding rows and columns to) certain noncommuting matrices A₁,...,A_d, related to the coordinate operators x₁,...,x_d, in R^d. We prove a correspondence between cubature formulae and “commuting extensions” of A₁,...,A_d, satisfying a compatibility condition which, in appropriate coordinates, constrains certain blocks in the extended matrices to be zero. Thus, the problem of finding cubature formulae can be transformed to the problem of computing (and then simultaneously diagonalizing) commuting extensions. We give a general discussion of existence and of the expected size of commuting extensions and briefly describe our attempts at computing them.     

Regular simplices, symmetric polynomials and the mean value property

LetP be ann-dimensional regular simplex in ℝⁿ centered at the origin, and let P(k) be thek-skeleton ofP fork = 0, 1,…,n. Then the set of all continuous functions in ℝⁿ satisfying the mean value property with respect to P(k) forms a finite-dimensional linear space of harmonic polynomials. In this paper the function space is explicitly determined by group theoretic and combinatorial arguments for symmetric polynomials.     

New cubature formulae and hyperinterpolation in three variables

A new algebraic cubature formula of degree 2n+1 for the product Chebyshev measure in the d-cube with ≈n ^d /2 ^d−1 nodes is established. The new formula is then applied to polynomial hyperinterpolation of degree n in three variables, in which coefficients of the product Chebyshev orthonormal basis are computed by a fast algorithm based on the 3-dimensional FFT. Moreover, integration of the hyperinterpolant provides a new Clenshaw-Curtis type cubature formula in the 3-cube. Work supported by the National Science Foundation under Grant DMS-0604056, by the “ex-60%” funds of the Universities of Padova and Verona, and by the INdAM-GNCS.     

On locally constructible spheres and balls

Durhuus and Jonsson (1995) introduced the class of “locally constructible” (LC) 3-spheres and showed that there are only exponentially many combinatorial types of simplicial LC 3-spheres. Such upper bounds are crucial for the convergence of models for 3D quantum gravity.     

The distribution of points on the sphere and corresponding cubature formulae

In applications, for instance in optics and astrophysics, thereis a need for high-accuracy integration formulae for functionson the sphere. To construct better formulae than previouslyused, almost equidistantly spaced nodes on the sphere and weightsbelonging to these nodes are required. This problem is closelyrelated to an optimal dispersion problem on the sphere and tothe theories of spherical designs and multivariate Gauss quadratureformulae. We propose a two-stage algorithm to compute optimal point locationson the unit sphere and an appropriate algorithm to calculatethe corresponding weights of the cubature formulae. Points aswell as weights are computed to high accuracy. These algorithmscan be extended to other integration problems. Numerical examplesshow that the constructed formulae yield impressively smallintegration errors of up to 10^-12.     

Weighted spherical semidesigns and cubature formulae for calculating integrals on a sphere

In this paper we study weighted spherical semidesigns, i.e., systems of points on a sphere of a specific type. We propose a new proof of the necessary and sufficient condition for a system of points on a sphere to be a weighted spherical semidesign. This criterion gives new approaches to the construction of cubature formulae for calculating integrals over a sphere with the degree of accuracy of 5 and 9.     

Orthogonal polynomials for exponential weights on , II

Let I=[0,d), where d is finite or infinite. Let W_ρ(x)=x^ρexp(-Q(x)), where and Q is continuous and increasing on I, with limit ∞ at d. We obtain further bounds on the orthonormal polynomials associated with the weight , finer spacing on their zeros, and estimates of their associated fundamental polynomials of Lagrange interpolation. In addition, we obtain weighted Markov and Bernstein inequalities.     

Isotropic covariance matrix polynomials on spheres

This article characterizes the covariance matrix function of a Gaussian or elliptically contoured vector random field that is stationary, isotropic, and mean square continuous on the sphere. By applying the characterization to examine the validity of a matrix function whose entries are polynomials of degrees up to 4, we obtain a necessary and sufficient condition for the polynomial matrix to be an isotropic covariance matrix function on the sphere.     

Orthogonal polynomials on the disc

We consider the space P_n of orthogonal polynomials of degree n on the unit disc for a general radially symmetric weight function. We show that there exists a single orthogonal polynomial whose rotations through the angles , j=0,1,…,n forms an orthonormal basis for P_n, and compute all such polynomials explicitly. This generalises the orthonormal basis of Logan and Shepp for the Legendre polynomials on the disc.Furthermore, such a polynomial reflects the rotational symmetry of the weight in a deeper way: its rotations under other subgroups of the group of rotations forms a tight frame for P_n, with a continuous version also holding. Along the way, we show that other frame decompositions with natural symmetries exist, and consider a number of structural properties of P_n including the form of the monomial orthogonal polynomials, and whether or not P_n contains ridge functions.     

A class of bivariate orthogonal polynomials and cubature formula

Summary. The existence of Gaussian cubature for a given measure depends on whether the corresponding multivariate orthogonal polynomials have enough common zeros. We examine a class of orthogonal polynomials of two variables generated from that of one variable. Received February 9, 1993 / Revised version received January 18, 1994     

Explicit formulae for sums of products of Bernoulli polynomials,including poly-Bernoulli polynomials

We study sums of products of Bernoulli polynomials, including poly-Bernoulli polynomials. As a main result, for any positive integer $m$ , explicit expressions of sums of $m$ products are given. This result extends that of the first author, as well as the famous Euler formula about sums of two products of Bernoulli numbers.     

Note on fast evaluation of best L2 approximation polynomials on simplices

We show how best the L₂ approximation polynomial to a given square integrable function on a simplex can be computed very effectively.     

Derivative polynomials and closed-form higher derivative formulae

In a recent paper, Adamchik [1] expressed in a closed-form symbolic derivatives of four functions belonging to the class of functions whose derivatives are polynomials in terms of the same functions. In this sequel, simple closed-form higher derivative formulae which involve the Carlitz-Scoville higher order tangent and secant numbers are derived for eight trigonometric and hyperbolic functions.     

Isospectral potentials and conformally equivalent isospectral metrics on spheres,balls and Lie groups

We construct pairs of conformally equivalent isospectral Riemannian metrics ?1g and ?2g on spheres Sⁿ and balls Bⁿ⁺¹ for certain dimensions n, the smallest of which is n=7, and on certain compact simple Lie groups. In the case of Lie groups, the metric g is left-invariant. In the case of spheres and balls, the metric g not the standard metric but may be chosen arbitrarily close to the standard one. For the same manifolds (M, g) we also show that the functions ?₁ and ?₂ are isospectral potentials for the Schrödinger operator ?₂\gD + \gf. To our knowledge, these are the first examples of isospectral potentials and of isospectral conformally equivalent metrics on simply connected closed manifolds.     

Orthogonal polynomials,Padé approximations and A-stability

This paper consists of a brief survey of implicit Runge-Kutta methods based on Gauss-Legendre and closely related quadrature formulas, together with a new proof of a classical result. The classical result concerns the A-stability of those implicit Runge-Kutta methods whose stability functions are Padé approximations to the exponential function and lie on the diagonal and the first two subdiagonals of the Padé table.     

Orthogonal projections of cubes and regular simplices into a plane

We show that for every \({k\ge 2}\) and \({n\ge k}\), there is an \({n}\)-dimensional unit cube in \({\mathbb{R}^n}\) which is mapped to a regular \({2k}\)-gon by an orthogonal projection in \({\mathbb{R}^n}\) onto a \({2}\)-dimensional subspace. Moreover, by increasing dimension \({n}\), arbitrary large regular \({2k}\)-gon can be obtained in such a way. On the other hand, for every \({m\ge 3}\) and \({n\ge m-1}\), there is an \({n}\)-dimensional regular simplex of unit edge in \({\mathbb{R}^n}\) which is mapped to a regular \({m}\)-gon by an orthogonal projection onto a plane. Moreover, contrary to the cube case, arbitrary small regular \({m}\)-gon can be obtained in such a way, by increasing dimension \({n}\).     

On the Lipschitz classification of normed spaces, unit balls, and spheres

For every normed space , we note its closed unit ball and unit sphere by and , respectively. Let and be normed spaces such that is Lipschitz homeomorphic to , and is Lipschitz homeomorphic to .

We prove that the following are equivalent:

1. is Lipschitz homeomorphic to .

2. is Lipschitz homeomorphic to .

3. is Lipschitz homeomorphic to .

This result holds also in the uniform category, except (2 or 3) 1 which is known to be false.