Necessary and Sufficient Conditions of Doubly Weighted Hardy-Littlewood-Sobolev Inequality

Using product and convolution theorems on Lorentz spaces, we characterize the sufficient and necessary conditions which ensure the validity of the doubly weighted Hardy-Littlewood-Sobolev inequality. It should be pointed out that we consider whole ranges of p and q, i.e., 0 p ≤∞ and 0 q ≤∞.     

Weak-type endpoint bounds for Riesz means

We use restriction theory to prove optimal weak-type bounds of Bochner-Riesz multipliers and Riesz means of elliptic pseudo-differential operators on compact manifolds, for .

     

A Search for Best Constants in the Hardy-Littlewood Maximal Theorem

Let be the centered maximal operator on the line. Through a numerical search procedure, we have conjectural best constants for the weak-type 1-1 estimate (3/2) and the L^p estimate (the constant B(p,1) such that We prove that these constants are lower bounds for the best constants and discuss the numerical evidence for the conjectures.     

PREDUAL SPACES FOR Q SPACES

To find the predual spaces Pα（R^n） of Qα（R^n） is an important motivation in the study of Q spaces. In this article, wavelet methods are used to solve this problem in a constructive way. First, an wavelet tent atomic characterization of Pα（Rn） is given, then its usual atomic characterization and Poisson extension characterization are given. Finally, the continuity on Pα of Calderon-Zygmund operators is studied, and the result can be also applied to give the Morrey characterization of Pα（Rn）.     

STABILIZED FEM FOR CONVECTION-DIFFUSION PROBLEMS ON LAYER-ADAPTED MESHES

The application of a standard Galerkin finite element method for convection-diffusion problems leads to oscillations in the discrete solution, therefore stabilization seems to be necessary. We discuss several recent stabilization methods, especially its combination with a Galerkin method on layer-adapted meshes. Supercloseness results obtained allow an improvement of the discrete solution using recovery techniques.     

A TWO-SCALE HIGHER-ORDER FINITE ELEMENT DISCRETIZATION FOR SCHROEDINGER EQUATION

In this paper, a two-scale higher-order finite element discretization scheme is proposed and analyzed for a Schroedinger equation on tensor product domains. With the scheme, the solution of the eigenvalue problem on a fine grid can be reduced to an eigenvalue problem on a much coarser grid together with some eigenvalue problems on partially fine grids. It is shown theoretically and numerically that the proposed two-scale higher-order scheme not only significantly reduces the number of degrees of freedom but also produces very accurate approximations.     

FINITE ELEMENTS WITH LOCAL PROJECTION STABILIZATION FOR INCOMPRESSIBLE FLOW PROBLEMS

In this paper we review recent developments in the analysis of finite element methods for incompressible flow problems with local projection stabilization （LPS）. These methods preserve the favourable stability and approximation properties of classical residual-based stabilization （RBS） techniques but avoid the strong coupling of velocity and pressure in the stabilization terms. LPS-methods belong to the class of symmetric stabilization techniques and may be characterized as variational multiscale methods. In this work we summarize the most important a priori estimates of this class of stabilization schemes developed in the past 6 years. We consider the Stokes equations, the Oseen linearization and the NavierStokes equations. Furthermore, we apply it to optimal control problems with linear（ized） flow problems, since the symmetry of the stabilization leads to the nice feature that the operations ＂discretize＂ and ＂optimize＂ commute.     

AN EFFICIENT METHOD FOR MULTIOBJECTIVE OPTIMAL CONTROL AND OPTIMAL CONTROL SUBJECT TO INTEGRAL CONSTRAINTS

We introduce a new and efficient numerical method for multicriterion optimal control and single criterion optimal control under integral constraints. The approach is based on extending the state space to include information on a ＂budget＂ remaining to satisfy each constraint; the augmented Hamilton-Jacobi-Bellman PDE is then solved numerically. The efficiency of our approach hinges on the causality in that PDE, i.e., the monotonicity of characteristic curves in one of the newly added dimensions. A semi-Lagrangian ＂marching＂ method is used to approximate the discontinuous viscosity solution efficiently. We compare this to a recently introduced ＂weighted sum＂ based algorithm for the same problem [25]. We illustrate our method using examples from flight path planning and robotic navigation in the presence of friendly and adversarial observers.     

DIRECT MINIMIZATION FOR CALCULATING INVARIANT SUBSPACES IN DENSITY FUNCTIONAL COMPUTATIONS OF THE ELECTRONIC STRUCTURE

In this article, we analyse three related preconditioned steepest descent algorithms, which are partially popular in Hartree-Fock and Kohn-Sham theory as well as invariant subspace computations, from the viewpoint of minimization of the corresponding functionals, constrained by orthogonality conditions. We exploit the geometry of the admissible manifold, i.e., the invariance with respect to unitary transformations, to reformulate the problem on the Grassmann manifold as the admissible set. We then prove asymptotical linear convergence of the algorithms under the condition that the Hessian of the corresponding Lagrangian is elliptic on the tangent space of the Grassmann manifold at the minimizer.     

RECENT PROGRESS ON SPHERICAL HARMONIC APPROXIMATION MADE BY BNU RESEARCH GROUP --In memory of Professor Sun Yongsheng

As early as in 1990, Professor Sun Yongsheng, suggested his students at Beijing Normal University to consider research problems on the unit sphere. Under his guidance and encouragement his students started the research on spherical harmonic analysis and approximation. In this paper, we incompletely introduce the main achievements in this area obtained by our group and relative researchers during recent 5 years （2001-2005）. The main topics are： convergence of Cesaro summability, a.e. and strong summability of Fourier-Laplace series; smoothness and K-functionals; Kolmogorov and linear widths.     

Cartesian closed categories of FƵ-domains

A subset system Z assigns to each partially ordered set P a certain collection Z（P） of subsets. In this paper, a new kind of subset systems called directable subset systems is introduced. For a directable subset system Z, the concepts of FZ-way-below relation and FZ-domain are introduced. The well-known Scott topology is naturally generalized to the Z-level and the resulting topology is called FZ-Scott topology, and the continuous functions with respect to this topology are characterized by preserving the suprema of directed Z-sets. Then, we mainly consider a generalization of the cartesian closedness of the categories DCPO of directed complete posets, BF of bifinite domains and FS of FS-domains to the Z-level. Corresponding to them, it is proved that, for a suitable subset system Z, the categories FZCPO of Z-complete posets, FSFZ of finitely separated FZ-domains and BFFZ of bifinite FZ-domains are all cartesian closed. Some examples of these categories are given.     

Projections,Birkhoff Orthogonality and Angles in Normed Spaces

Let X be a Minkowski plane, i.e., a real two dimensional normed linear space. We use projections to give a definition of the angle Aq（x, y） between two vectors x and y in X, such that x is Birkhoff orthogonal to y if and only if Aq（x,y）=π/2. Some other properties of this angle are also discussed.     

FINITE ELEMENT APPROXIMATIONS FOR SCHROEDINGER EQUATIONS WITH APPLICATIONS TO ELECTRONIC STRUCTURE COMPUTATIONS

In this paper, both the standard finite element discretization and a two-scale finite element discretization for SchrSdinger equations are studied. The numerical analysis is based on the regularity that is also obtained in this paper for the Schroedinger equations. Very satisfying applications to electronic structure computations are provided, too.     

Spectral gap and convergence rate for discrete-time Markov chains

Let P be a transition matrix which is symmetric with respect to a measure π.The spectral gap of P in L2(π)-space,denoted by gap(P),is defined as the distance between 1 and the rest of the spectrum of P.In this paper,we study the relationship between gap(P) and the convergence rate of Pn.When P is transient,the convergence rate of P n is equal to 1 gap(P).When P is ergodic,we give the explicit upper and lower bounds for the convergence rate of Pn in terms of gap(P).These results are extended to L∞(π)-space.     

K-controllability and Approximate K-controllability of Nonlinear Neutral Systems in Banach Spaces

In this paper, K-controllability and approximate K-controllability of non- linear neutral differential equations in Banach spaces are studied. Sufficient conditions are established for each of these types of controllability. The results are obtained by using Leray-Schauder theory.     

More about the kernel convergence and the ideal convergence

Let I 2N be an ideal and let XI = span{χI ： I ∈ I}, and let pI be the quotient norm of l∞/XI. In this paper, we show first that for each proper ideal I 2N, the ideal convergence deduced by I is equivalent to pI-kernel convergence. In addition, let K = {x＊oχ（·） ： x＊∈ p（e）}, where p（x） = lim supn→∞1/n（∑k=1n｜x（k）｜, and let Iμ = {A N ： μ（A） = 0} for all μ = x＊oχ（·） ∈ K. Then Iμ is a proper ideal. We also show that the ideal convergence deduced by the proper ideal Iμ, the p-kernel convergence and the statistical convergence are also equivalent.     

Visualizing diffusion tensor fields on streamsurfaces with merging ellipsoids and LIC texture

Streamsurfaces in diffusion tensor fields are used to represent structures with pri- marily planar diffusion. So far, however, no effort has been made on the visualization of the anisotropy of diffusion on them, although this information is very important to identify the problematic regions of these structures. We propose two methods to display this anisotropy information. The first one employs a set of merging ellipsoids, which simultaneously character- ize the local tensor details - anisotropy - on them and portray the shape of the streamsurfaces. The weight between the streamsurfaces continuity and the discrete local tensors can be inter- actively adjusted by changing some given parameters. The second one generates a dense LIC （line integral convolution） texture of the two tangent eigenvector fields along the streamsurfaces firstly, and then blends in some color mapping indicating the anisotropy information. For high speed and high quality of texture images, we confine both the generation and the advection of the LIC texture in the image space. Merging ellipsoids method reveals the entire anisotropy information at discrete points by exploiting the geometric attribute of ellipsoids, and thus suits for local and detailed examination of the anisotropy; the texture-based method gives a global representation of the anisotropy on the whole streamsurfaces with texture and color attributes. To reveal the anisotropy information more efficiently, we integrate the two methods and use them at two different levels of details.     

NONLINEAR RANK-ONE MODIFICATION OF THE SYMMETRIC EIGENVALUE PROBLEM

Nonlinear rank-one modification of the symmetric eigenvalue problem arises from eigenvibrations of mechanical structures with elastically attached loads and calculation of the propagation modes in optical fiber. In this paper, we first study the existence and uniqueness of eigenvalues, and then investigate three numerical algorithms, namely Picard iteration, nonlinear Rayleigh quotient iteration and successive linear approximation method （SLAM）. The global convergence of the SLAM is proven under some mild assumptions. Numerical examples illustrate that the SLAM is the most robust method.     

ON SMOOTH LU DECOMPOSITIONS WITH APPLICATIONS TO SOLUTIONS OF NONLINEAR EIGENVALUE PROBLEMS

We study the smooth LU decomposition of a given analytic functional A-matrix A（A） and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A（A）, and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.     

对称熵损失下的一类指数族刻度参数估计

In this paper we investigate the estimator for the rth power of the scale parameter in a class of exponential family under symmetric entropy loss L（θ, δ） = v（θ/δ ＋ δ/θ - 2）. An exact form of the minimum risk equivariant estimator under symmetric entropy loss is given, and the minimaxity of the minimum risk equivariant estimator is proved. The results with regard to admissibility and inadmissibility of a class of linear estimators of the form cT（X） ＋ d are given, where T（X） Gamma（v, θ）.