An infeasible-interior-point potential-reduction algorithm for linear programming

n . The method is based on Rockafellar’s proximal point algorithm and a cutting-plane technique. At each step, we use an approximate proximal point p^a(x_k) of x_k to define a v_k∈∂_εkf(p^a(x_k)) with ε_k≤α∥v_k∥, where α is a constant. The method monitors the reduction in the value of ∥v_k∥ to identify when a line search on f should be used. The quasi-Newton step is used to reduce the value of ∥v_k∥. Without the differentiability of f, the method converges globally and the rate of convergence is Q-linear. Superlinear convergence is also discussed to extend the characterization result of Dennis and Moré. Numerical results show the good performance of the method. Received October 3, 1995 / Revised version received August 20, 1998 Published online January 20, 1999     

Wavelet bases of Hermite cubic splines on the interval

In this paper a pair of wavelets are constructed on the basis of Hermite cubic splines. These wavelets are in C¹ and supported on [−1,1]. Moreover, one wavelet is symmetric, and the other is antisymmetric. These spline wavelets are then adapted to the interval [0,1]. The construction of boundary wavelets is remarkably simple. Furthermore, global stability of the wavelet basis is established. The wavelet basis is used to solve the Sturm–Liouville equation with the Dirichlet boundary condition. Numerical examples are provided. The computational results demonstrate the advantage of the wavelet basis. Dedicated to Dr. Charles A. Micchelli on the occasion of his 60th birthday Mathematics subject classifications (2000) 42C40, 41A15, 65L60. Research was supported in part by NSERC Canada under Grants # OGP 121336.     

Piecewise linear wavelets on sierpinski gasket type fractals

Let SG denote the Sierpinski gasket with Hausdorff measure μ of dimensionlog 3/log 2, let PL_k denote the continuous piecewise linear functions with respect to the usual triangulation of SG into 3^k triangles, and let W_k denote the orthogonal complement of PL_k−1 in PL_k. We construct a basis for each W_k, so that the entire collection is a frame for L²(dμ). This wavelet basis is obtained from three wavelet generators by scaling, translation and rotation, and the wavelets are supported either by one corner triangle or a pair of adjacent triangles in the triangulation of level k − 1. Analogous bases are constructed in the von Koch curve, the hexagasket, and the n-dimensional analog of SG.     

Construction of multivariate compactly supported orthonormal wavelets

We propose a constructive method to find compactly supported orthonormal wavelets for any given compactly supported scaling function φ in the multivariate setting. For simplicity, we start with a standard dilation matrix 2I_2×2 in the bivariate setting and show how to construct compactly supported functions ψ₁,. . .,ψ_n with n>3 such that {2^kψ_j(2^kx−ℓ,2^ky−m), k,ℓ,m∈Z, j=1,. . .,n} is an orthonormal basis for L₂(ℝ²). Here, n is dependent on the size of the support of φ. With parallel processes in modern computer, it is possible to use these orthonormal wavelets for applications. Furthermore, the constructive method can be extended to construct compactly supported multi-wavelets for any given compactly supported orthonormal multi-scaling vector. Finally, we mention that the constructions can be generalized to the multivariate setting. Dedicated to Professor Charles A. Micchelli on the occasion of his 60th birthday Mathematics subject classifications (2000) 42C15, 42C30.     

L P boundedness of area function associated with operators on product spaces

Let L be the infinitesimal generator of an analytic semigroup on L2 (Rn) with suitable upper bounds on its heat kernels. Assume that L has a bounded holomorphic functional calculus on L2(Rn). In this paper,we define the Littlewood- Paley g function associated with L on Rn × Rn, denoted by GL(f)(x1, x2), and decomposition, we prove that ‖SL(f)‖p ≈‖GL(f)‖p ≈‖f‖p for 1 ＜ p ＜∞.     

Compactly supported wavelet bases for Sobolev spaces

In this paper we investigate compactly supported wavelet bases for Sobolev spaces. Starting with a pair of compactly supported refinable functions φ and in satisfying a very mild condition, we provide a general principle for constructing a wavelet ψ such that the wavelets ψ_jk:=2^j/2ψ(2^j·−k) ( ) form a Riesz basis for . If, in addition, φ lies in the Sobolev space , then the derivatives 2^j/2ψ^(m)(2^j·−k) ( ) also form a Riesz basis for . Consequently, is a stable wavelet basis for the Sobolev space . The pair of φ and are not required to be biorthogonal or semi-orthogonal. In particular, φ and can be a pair of B-splines. The added flexibility on φ and allows us to construct wavelets with relatively small supports.     

Small fractional parts of additive forms

Letf(x)=θ₁ x ₁ ^k +...+θ _s x _s ^k be an additive form with real coefficients, and ∥α∥ = min {|α-u|:uεℤ} denote the distance fromα to the nearest integer. We show that ifθ ₁,…,θ _s , are algebraic ands = 4k then there are integersx ₁,…,x _s , satisfying l ≤x ₁,≤ N and ∥f(x)∥ ≤ N ^E , withE = − 1 + 2/e. Whens = λk, 1 ≤λ ≤ 2k, the exponentE may be replaced byλE/4, and if we drop the condition thatθ ₁,…,θ _s , be algebraic then the result holds for almost all values of θεℝ ^s . Whenk ≥ 6 is small a better exponent is obtained using Heath-Brown’s version of Weyl’s estimate.     

Local means,wavelet bases and wavelet isomorphisms in Besov‐Morrey and Triebel‐Lizorkin‐Morrey spaces

We consider local means with bounded smoothness for Besov‐Morrey and Triebel‐Lizorkin‐Morrey spaces. Based on those we derive characterizations of these spaces in terms of Daubechies, Meyer, Bernstein (spline) and more general r‐regular (father) wavelets, finally in terms of (biorthogonal) wavelets which can serve as molecules and local means, respectively. Hereby both, local means and wavelet decompositions satisfy natural conditions concerning smoothness and cancellation (moment conditions). Moreover, the given representations by wavelets are unique and yield isomorphisms between the considered function spaces and appropriate sequence spaces of wavelet coefficients. These wavelet representations lead to wavelet bases if, and only if, the function spaces coincide with certain classical Besov‐Triebel‐Lizorkin spaces.     

Accuracy of several multidimensional refinable distributions

Compactly supported distributions f₁,..., f_r on ℝ^d are fefinable if each f_i is a finite linear combination of the rescaled and translated distributions f_j(Ax−k), where the translates k are taken along a lattice Γ ⊂ ∝^d and A is a dilation matrix that expansively maps Γ into itself. Refinable distributions satisfy a refinement equation f(x)=Σ_k∈Λ c_k f(Ax−k), where Λ is a finite subset of Γ, the c_k are r×r matrices, and f=(f₁,...,f_r)^T. The accuracy of f is the highest degree p such that all multivariate polynomials q with degree(q)<p are exactly reproduced from linear combinations of translates of f₁,...,f_r along the lattice Γ. We determine the accuracy p from the matrices c_k. Moreover, we determine explicitly the coefficients y_α,i(k) such that x^α=Σ _i=1 ^r Σ_k∈Γy_α,i(k) fi(x+k). These coefficients are multivariate polynomials y_α,i(x) of degree |α| evaluated at lattice points k∈Γ.     

On explicit basis of bivariate spline space

For a subdivision Δ of a region in d-dimensional Euclidean space, we consider computation of dimension and of basis function in spline space S _k ^r (Δ) consisting of all C piecewise polynomial functions over Δ of degree at most k. A computational scheme is presented for computing the dimension and bases of spline space S _k ^r (Δ). This scheme based on the Grobner basis algorithm and the smooth co-factor method for computing multivariate spline. For bivariate splines, explicit basis functions of S _k ^r (Δ) are obtained for any integer k and r when Δ is a cross-cut partition. The Project is partly supported by the Science and Technology New Star Plan of Beijing and Education Committee of Beijing.     

A remark on interpolation by bivariate splines

We study the determining set for bivariate spline spacesS _k ^o on type-1 triangulation of square using B-net techniques. We further construct the interpolation schemes for these spline spaces that are unisolvent for any function f of C^σ.     

Tight frames and geometric properties of wavelet sets

A construction for providing single dyadic orthonormal wavelets in Euclidean space ℝ^d is given. It is called the general neighborhood mapping construction. The fact that one wavelet is sufficient to generate an orthonormal basis for L²(ℝ^d) is the critical issue. The validity of the construction is proved, and the construction is implemented computationally to provide a host of examples illustrating various geometrical properties of such wavelets in the spectral domain. Because of the inherent complexity of these single orthonormal wavelets, the method is applied to the construction of single dyadic tight frame wavelets, and these tight frame wavelets can be surprisingly simple in nature. The structure of the spectral domains of the wavelets arising from the general neighborhood mapping construction raises a basic geometrical question. There is also the question of whether or not the general neighborhood mapping construction gives rise to all single dyadic orthonormal wavelets. Results are proved giving partial answers to both of these questions. Dedicated to Charles A. Micchelli for his 60th birthday Mathematics subject classification (2000) 42C40. John J. Benedetto: Both authors gratefully acknowledge support from ONR Grant N000140210398. The first named author also gratefully acknowledges support from NSF DMS Grant 0139759.     

〈I〉型三角剖分下非张量积连续小波基的构造

多维非张量积小波是近年小波研究领域中的热点问题之一 ,它们与多维张量积小波相比具有更多的优势 .关于高维张量积、非张量积小波 ,目前已有一些很好的工作 (见文[2 ] [3 ] [4 ] ) ,但关于样条小波 ,还有许多问题有待于研究 .本文针对〈I〉型三角剖分下的二维线性元空间 ,讨论其具有紧支集和对称性的半正交样条小波基 .给定 x1 x2 平面上的〈I〉型三角剖分 (图 1 ( a)所示 ) ,记 j=( j1 ,j2 ) ,| j| =j1 + j2 ,πm= { 0≤ |j|≤ mCj1j2 xj11 xj22 ,Cj1,j2 是任意实数 }为次数不超过 m的代数多项式全体 .引入剖分尺度为 1的线性元空间 V0…     

Integral inequalities with weights for maxiamal operators

Let μ be a measure on the upper half-space R ₊ ⁿ⁺¹ , and v a weight onR ⁿ, we give a characterization for the pair (v, μ) such that ∥M(fv)∥_{L Θ (μ}) ⩽ c ∥f∥_{L Θ (μ}), where Φ is an N-function satisfying Δ₂ condition andMf(x,t), is the maximal function onR ₊ ⁿ⁺¹ , which was introduced by Ruiz, F. and Torrea, J.. Supported by NSFC.     

Von Neumann-Jordan Constants of Absolute Normalized Norms on C^n

In this note, we give some estimations of the Von Neumann-Jordan constant C _{N J} (∥·∥_ψ) of Banach space (ℂ ⁿ , ∥·∥_ψ), where ∥·∥_ψ is the absolute normalized norm on ℂ ⁿ given by function ψ. In the case where ψ and φ are comparable, n=2 and C _{N J} (∥·∥_ψ)=1, we obtain a formula of computing C _{N J} (∥·∥_ψ). Our results generalize some results due to Saito and others. Received May 11, 2002, Accepted November 20, 2002 This work is partly supported by NNSF of China (No. 19771056)     

Superstability of the generalized orthogonality equation on restricted domains

Chmielinski has proved in the paper [4] the superstability of the generalized orthogonality equation |〈f(x), f(y)〉| = |〈x,y〉|. In this paper, we will extend the result of Chmielinski by proving a theorem: LetD _n be a suitable subset of ℝⁿ. If a function f:D _n → ℝⁿ satisfies the inequality ∥〈f(x), f(y)〉| |〈x,y〉∥ ≤ φ(x,y) for an appropriate control function φ(x, y) and for allx, y ∈ D _n, thenf satisfies the generalized orthogonality equation for anyx, y ∈ D _n.     

A general construction of nonseparable multivariate orthonormal wavelet bases

The construction of nonseparable and compactly supported orthonormal wavelet bases of L ²(R ⁿ ); n ≥ 2, is still a challenging and an open research problem. In this paper, we provide a special method for the construction of such wavelet bases. The wavelets constructed by this method are dyadic wavelets. Also, we show that our proposed method can be adapted for an eventual construction of multidimensional orthogonal multiwavelet matrix masks, candidates for generating multidimensional multiwavelet bases.      

Some simple Haar-type wavelets in higher dimensions

An orthonormal wavelet system in ℝ^d, d ∈ ℕ, is a countable collection of functions {ψ _j,k ^ℓ }, j ∈ ℤ, k ∈ ℤ^d, ℓ = 1,..., L, of the form that is an orthonormal basis for L² (ℝ^d), where a ∈ GL_d (ℝ) is an expanding matrix. The first such system to be discovered (almost 100 years ago) is the Haar system for which L = d = 1, ψ¹(x) = ψ(x) = κ[0,1/2)^(x) − κ_[l/2,1) ^(x), a = 2. It is a natural problem to extend these systems to higher dimensions. A simple solution is found by taking appropriate products Φ(x₁, x₂, ..., x_d) = φ₁ (x₁)φ₂(x₂) ... φ_d(x_d) of functions of one variable. The obtained wavelet system is not always convenient for applications. It is desirable to find “nonseparable” examples. One encounters certain difficulties, however, when one tries to construct such MRA wavelet systems. For example, if a = ( _1-1 ^{1 1} ) is the quincunx dilation matrix, it is well-known (see, e.g., [5]) that one can construct nonseparable Haar-type scaling functions which are characteristic functions of rather complicated fractal-like compact sets. In this work we shall construct considerably simpler Haar-type wavelets if we use the ideas arising from “composite dilation” wavelets. These were developed in [7] and involve dilations by matrices that are products of the form a^jb, j ∈ ℤ, where a ∈ GL_d(ℝ) has some “expanding” property and b belongs to a group of matrices in GL_d(ℝ) having |det b| = 1.     

On estimation and detection of a function of infinitely many variables

We observe an unknown function of infinitely many variables f = f(t), t = (t₁, ..., t_n, ... ) ∈, [0, 1]^∞, in the Gaussian white noise of level ε > 0. We suppose that in each variable there exists a 1-periodical σ-smooth extension of the function f(t) to IR ^∞. Taking a quantity σ > 0 and a positive sequence a = {a_k}, we consider the set that consists of functions f such that . We consider the cases a_k = k^α and a_k = exp(λk), α > 0, λ > 0. We would like to estimate a function f ∈ or to test the null hypothesis H₀: f = 0 against the alternatives f ∈ , where the set consists of functions f ∈ such that ∥f∥₂ ≥ r. In the estimation problem, we obtain the asymptotics (as ε → 0) of the minimax quadratic risk. In the detection problem, we study the sharp asymptotics of minimax separation rates f _ɛ ^* that provide distiguishability in the problems. Bibliography: 12 titles. __________ Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 328, 2005, pp. 91–113.     

A nodal basis of Cμ-rational spline functions on triangulations

Let ρ be a triangulation of a polygonal domain D⊂R² with vertices V={v_i:l≤i≤N_v} and RS^k(D, ρ)={u∈C^k(D): ≠ T∈ρ, u/_T is a rational function}. The purpose of this paper is to study the existence and construction of C^μ-rational spline functions on any triangulation ρ for CAGD. The Hermite problem H^μ(V,U)={find u∈U: D^αu(v_i)=D^αf(v_i),|α|≤μ} is solved by the generalized wedge function method in rational spline function family, i.e. U=RS^μ. this solution needs only the knowledge of partial derivatives of order≤μ at v_i. The explicit repesentations of all C^μ-GWF(generalized wedge functions)and the interpolating operator with degree of precision at least 2μ+1 for any triangulation are given.