More on a New Concept of Entropy and Information

An alternative notion of entropy called CRE is proposed in [Ra1] Rao et al. (IEEE Trans. Inf. Theory 50, 2004). This preserves many of the properties of Shannon Entropy and possesses mathematical properties, which we hope will be of use in statistical estimates. In this article, we develop some more mathematical properties of CRE, show its relation to the L log L class, and characterize among others the Weibull distribution.     

Algebraic analysis of multigrid algorithms

We study the convergence rate of multilevel algorithms from an algebraic point of view. This requires a detailed analysis of the constant in the strengthened Cauchy–Schwarz inequality between the coarse‐grid space and a so‐called complementary space. This complementary space may be spanned by standard hierarchical basis functions, prewavelets or generalized prewavelets. Using generalized prewavelets, we are able to derive a constant in the strengthened Cauchy–Schwarz inequality which is less than 0.31 for the L² and H¹ bilinear form. This implies a convergence rate less than 0.15. So, we are able to prove fast multilevel convergence. Furthermore, we obtain robust estimations of the convergence rate for a large class of anisotropic ellipic equations, even for some that are not H¹‐elliptic. Copyright © 1999 John Wiley & Sons, Ltd.     

Kernels and best approximations related to the system of ultraspherical polynomials

We study the uniformly bounded orthonormal system of functions

where is the normalized system of ultraspherical polynomials. We investigate some approximation properties of the system and we show that these properties are similar to one's of the trigonometric system. First, we obtain estimates of L^p-norms of the kernels of the system . These estimates enable us to prove Nikol'skiı˘-type inequalities for -polynomials. Next, we prove directly that is a basis in each , where w is an arbitrary A_p-weight function. Finally, we apply these results to get sharp inequalities for the best -approximations in L^q in terms of the best -approximations in . For the trigonometric system such inequalities have been already known.     

Hecke Eigenfunctions of Quantized Cat Maps Modulo Prime Powers

This paper continues the work done in Olofsson [Commun Math Phys 286(3):1051–1072, 2009] about the supremum norm of eigenfunctions of desymmetrized quantized cat maps. N will denote the inverse of Planck’s constant and we will see that the arithmetic properties of N play an important role. We prove the sharp estimate ||ψ||_∞ = O(N ^1/4) for all normalized eigenfunctions and all N outside of a small exceptional set. We are also able to calculate the value of the supremum norms for most of the so called newforms. For a given N = p ⁿ , with n > 2, the newforms can be divided in two parts (leaving out a small number of them in some cases), the first half all have supremum norm about ${2/\sqrt{1\pm 1/p}}This paper continues the work done in Olofsson [Commun Math Phys 286(3):1051–1072, 2009] about the supremum norm of eigenfunctions of desymmetrized quantized cat maps. N will denote the inverse of Planck’s constant and we will see that the arithmetic properties of N play an important role. We prove the sharp estimate ||ψ||_∞ = O(N ^1/4) for all normalized eigenfunctions and all N outside of a small exceptional set. We are also able to calculate the value of the supremum norms for most of the so called newforms. For a given N = p ⁿ , with n > 2, the newforms can be divided in two parts (leaving out a small number of them in some cases), the first half all have supremum norm about 2/?{1±1/p}{2/\sqrt{1\pm 1/p}} and the supremum norm of the newforms in the second half have at most three different values, all of the order N ^1/6. The only dependence of A is that the normalization factor is different if A has eigenvectors modulo p or not. We also calculate the joint value distribution of the absolute value of n different newforms.     

Function spaces of varying smoothness I

This paper deals with function spaces of varying smoothness. It is a modified version of corresponding parts of [8]. Corresponding spaces of positive smoothness s (x) will be considered in part II. We define the spaces B_p (?ⁿ ), where the function ??: x ? s (x) is negative and determines the smoothness pointwise. First we prove basic properties and then we use different wavelet decompositions to get information about the local smoothness behavior. The main results are characterizations of the spaces B_p (?ⁿ ) by weighted sequence space norms of the wavelet coefficients. These assertions are used to prove an interesting connection to the so‐called two‐microlocal spaces C^{s,s ′} (x⁰). (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Estimates of the deviations from the exact solutions for variational inequalities describing the stationary flow of certain viscous incompressible fluids

This paper is concerned with computable and guaranteed upper bounds of the difference between exact solutions of variational inequalities arising in the theory of viscous fluids and arbitrary approximations in the corresponding energy space. Such estimates (also called error majorants of functional type) have been derived for the considered class of nonlinear boundary‐value problems in (Math. Meth. Appl. Sci. 2006; 29:2225–2244) with the help of variational methods based on duality theory from convex analysis. In the present paper, it is shown that error majorants can be derived in a different way by certain transformations of the variational inequalities that define generalized solutions. The error bounds derived by this techniques for the velocity function differ from those obtained by the variational method. These estimates involve only global constants coming from Korn‐ and Friedrichs‐type inequalities, which are not difficult to evaluate in case of Dirichlet boundary conditions. For the case of mixed boundary conditions, we also derive another form of the estimate that contains only one constant coming from the following assertion: the L² norm of a vector‐valued function from H¹(Ω) in the factor space generated by the equivalence with respect to rigid motions is bounded by the L² norm of the symmetric part of the gradient tensor. As for some ‘simple’ domains such as squares or cubes, the constants in this inequality can be found analytically (or numerically), we obtain a unified form of an error majorant for any domain that admits a decomposition into such subdomains. Copyright © 2009 John Wiley & Sons, Ltd.     

Statistically Efficient Thinning of a Markov Chain Sampler

It is common to subsample Markov chain output to reduce the storage burden. Geyer shows that discarding k ? 1 out of every k observations will not improve statistical efficiency, as quantified through variance in a given computational budget. That observation is often taken to mean that thinning Markov chain Monte Carlo (MCMC) output cannot improve statistical efficiency. Here, we suppose that it costs one unit of time to advance a Markov chain and then θ > 0 units of time to compute a sampled quantity of interest. For a thinned process, that cost θ is incurred less often, so it can be advanced through more stages. Here, we provide examples to show that thinning will improve statistical efficiency if θ is large and the sample autocorrelations decay slowly enough. If the lag ? ? 1 autocorrelations of a scalar measurement satisfy ρ_? > ρ_{? + 1} > 0, then there is always a θ < ∞ at which thinning becomes more efficient for averages of that scalar. Many sample autocorrelation functions resemble first order AR(1) processes with ρ_? = ρ^|?| for some ? 1 < ρ < 1. For an AR(1) process, it is possible to compute the most efficient subsampling frequency k. The optimal k grows rapidly as ρ increases toward 1. The resulting efficiency gain depends primarily on θ, not ρ. Taking k = 1 (no thinning) is optimal when ρ ? 0. For ρ > 0, it is optimal if and only if θ ? (1 ? ρ)²/(2ρ). This efficiency gain never exceeds 1 + θ. This article also gives efficiency bounds for autocorrelations bounded between those of two AR(1) processes. Supplementary materials for this article are available online.     

Large Mass Global Solutions for a Class of L 1-Critical Nonlocal Aggregation Equations and Parabolic-Elliptic Patlak-Keller-Segel Models

We consider a class of L ¹ critical nonlocal aggregation equations with linear or nonlinear porous media-type diffusion which are characterized by a long-range interaction potential that decays faster than the Newtonian potential at infinity. The fast decay breaks the L ¹ scaling symmetry and we prove that all ‘sufficiently spread out’ initial data, even with supercritical mass, results in global, decaying solutions. In particular, we produce decaying solutions with arbitrary mass in cases for which finite time blow up solutions or non-decaying solutions are also known to exist for sufficiently large mass. This is in contrast to the classical parabolic-elliptic PKS for which essentially all solutions with supercritical mass blow up in finite time. The results with linear diffusion are proved using properties of the Fokker-Planck semi-group whereas the results with nonlinear diffusion are proved using a more interesting bootstrap argument coupling the entropy-entropy dissipation methods of the porous media equation together with higher L ^p estimates similar to those used in small-data and local theory for PKS-type equations.     

（L^p,L^q） Estimates for Multilinear Oscillatory Singular Integralswith Smooth Phases

In this paper,the authors establish the weighted (L^p,L^q) estimates for a class of multilinear oscillatory singular integrals with smooth phases.Certain endpoint estimates are also considered.     

Weak Dispersive Estimates for Schrödinger Equations with Long Range Potentials

We prove local smoothing estimates for the Schrödinger initial value problem with data in the energy space L ²(? ^d ), d ≥ 2 and a general class of potentials. In the repulsive setting we have to assume just a power like decay (1 + |x|)^?γ for some γ > 0. Also attractive perturbations are considered. The estimates hold for all time and as a consequence a weak dispersion of the solution is obtained. The proofs are based on similar estimates for the corresponding stationary Helmholtz equation and Kato H-smooth theory.     

On Geodesic <Emphasis Type="Italic">E</Emphasis>-Convex Sets,Geodesic <Emphasis Type="Italic">E</Emphasis>-Convex Functions and <Emphasis Type="Italic">E</Emphasis>-Epigraphs

In this paper, we introduce a new class of sets and a new class of functions called geodesic E-convex sets and geodesic E-convex functions on a Riemannian manifold. The concept of E-quasiconvex functions on R ⁿ is extended to geodesic E-quasiconvex functions on Riemannian manifold and some of its properties are investigated. Afterwards, we generalize the notion of epigraph called E-epigraph and discuss a characterization of geodesic E-convex functions in terms of its E-epigraph. Some properties of geodesic E-convex sets are also studied.     

On a decay rate for 1D-viscous compressible barotropic fluid equations

The Navier-Stokes equations of a compressible barotropic fluid in 1D with zero velocity boundary conditions are considered. We study the case of large initial data in H ¹ as well as the mass force such that the stationary density is positive. The uniform lower bound for the density is proved. By constructing suitable Lyapunov functionals, decay rate estimates in L ²-norm and H ¹-norm are given. The decay rate is exponential if so the decay rate of the nonstationary part of the mass force is. The results are proved in the Eulerian coordinates for a wide class of increasing state functions including with any γ > 0 as well as functions of arbitrarily fast growth. We also extend the results for equations of a multicomponent compressible barotropic mixture (in the absence of chemical reactions). Received December 20, 2000; accepted February 27, 2001.     

Variational problems for maps of bounded variation with values inS 1

The main goal of this paper is to characterize the weak limits of sequences of smooth maps from a Riemannian manifold intoS ¹. This is achieved in terms of Cartesian currents. Applications to the existence of minimizers of area type functionals in the class of maps with values inS ¹ satisfying Dirchlet and homological conditions are then discussed. The so called dipole problem is solved, too.This work has been partially supported by the Ministero dell'Universitá e della Ricerca Scientifica, by C.N.R. contract n. 91.01343.CT01, and by the European Research Project GADGET. It was partially carried out while the first and the third author were visiting the Mathematisches Institut der Universität Bonn under the support respectively of the Alexander von Humboldt-Stiftung and of the SBF 256.     

Opportunity-based age replacement with a one-cycle criterion

In this paper age replacement (AR) and opportunity-based age replacement (OAR) for a unit are considered, based on a one-cycle criterion, both for a known and unknown lifetime distribution. In the literature, AR and OAR strategies are mostly based on a renewal criterion, but in particular when the lifetime distribution is not known and data of the process are used to update the lifetime distribution, the renewal criterion is less appropriate and the one-cycle criterion becomes an attractive alternative. Conditions are determined for the existence of an optimal replacement age T^* in an AR model and optimal threshold age T_opp^* in an OAR model, using a one-cycle criterion and a known lifetime distribution. In the optimal threshold age T_opp^*, the corresponding minimal expected costs per unit time are equal to the expected costs per unit time in an AR model. It is also shown that for a lifetime distribution with increasing hazard rate, the optimal threshold age is smaller than the optimal replacement age. For unknown lifetime distribution, AR and OAR strategies are considered within a nonparametric predictive inferential (NPI) framework. The relationship between the NPI-based expected costs per unit time in an OAR model and those in an AR model is investigated. A small simulation study is presented to illustrate this NPI approach.     

Uniform a priori estimates for a class of horizontal minimal equations

In the product space we obtain uniform a priori C ⁰ horizontal length estimates, uniform a priori C ¹ boundary gradient estimates, as well as uniform modulus of continuity, for a class of horizontal minimal equations. In two independent variables, we derive an uniform global a priori C ¹ estimates and we infer an existence result.     

Complexity and T ‐invariant of Abelian and Milnor groups,and complexity of 3‐manifolds

We investigate the notion of complexity for finitely presented groups and the related notion of complexity for three‐dimensional manifolds. We give two‐sided estimates on the complexity of all the Milnor groups (the finite groups with free action on S³), as well as for all finite Abelian groups. The ideas developed in the process also allow to construct two‐sided bounds for the values of the so‐called T ‐invariant (introduced by Delzant) for the above groups, and to estimate from below the value of T ‐invariant for an arbitrary finitely presented group. Using the results of this paper and of previous ones, we then describe an infinite collection of Seifert threemanifolds for which we can asymptotically determine the complexity in an exact fashion up to linear functions. We also provide similar estimates for the complexity of several infinite families of Milnor groups. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The structure of harmonizable isotropic random fields

The class of harmonizable processes and fields are a natural extension of the class of stationary processes and fields. Random fields admit an additional property called isotropy. The classical spectral and covariance representations for stationary isotropic random fields are extended to the harmonizable isotropic case. A classification of these fields is obtained based upon the smoothness properties of their covariances. In contrast to the stationary case, it is also shown that there exist non-trivial harmonizable isotropic fields which satisfy the Laplace operator in the L ²-sense     

Asymptotic expansions andL ∞-error estimates for mixed finite element methods for second order elliptic problems

Summary Asymptotic expansions for mixed finite element approximations of the second order elliptic problem are derived and Richardson extrapolation can be applied to increase the accuracy of the approximations. A new procedure, which is called the error corrected method, is presented as a further application of the asymptotic error expansion for the first order BDM approximation of the scalar field. The key point in deriving the asymptotic expansions for the error is an establishment ofL ¹-error estimates for mixed finite element approximations for the regularized Green's functions. As another application of theL ¹-error estimates for the regularized Green's functions, we shall present maximum norm error estimates for mixed finite element methods for second order elliptic problems.     

Compatible wavelet-Galerkin method to solve stokes interior problem with circular boundary

This article addresses Neumann boundary value interior problem of Stokes equations with circular boundary. By using natural boundary element method, the Stokes interior problem is reduced into an equivalent natural integral equation with a hyper-singular kernel, which is viewed as Hadamard finite part. Based on trigonometric wavelet functions, the compatible wavelet space is constructed so that it can serve as Galerkin trial function space. In proposed compatible wavelet-Galerkin method, the simple and accurate computational formulae of the entries in stiffness matrix are obtained by singularity removal technique. It is also proved that the stiffness matrix is almost a block diagonal matrix, and its diagonal sub-blocks all are both symmetric and circulant submatrices. These good properties indicate that a 2 ^J+3 × 2 ^J+3 stiffness matrix can be determined only by its 2 ^J + 3J + 1 entries. It greatly decreases the computational complexity. Some error estimates for the compatible wavelet-Galerkin projection solutions are established. Finally, several numerical examples are given to demonstrate the validity of the proposed approach.