Applications of the Bernstein-Durrmeyer operators in estimating the norm of Mercer kernel matrices

The paper is related to the norm estimate of Mercer kernel matrices.The lower and upper bound estimates of Rayleigh entropy numbers for some Mercer kernel matrices on[0,1]×[0,1]based on the Bernstein-Durrmeyer operator kernel ale obtained,with which and the approximation property of the Bernstein-Durrmeyer operator the lower and upper bounds of the Rayleigh entropy number and the l2-norm for general Mercer kernel matrices on[0,1]×[0,1]are provided.     

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 THE ORDER OF SUMMABILITY OF THE FOURIER INVERSION FORMULA

In this article we show that the order of the point value, in the sense of Lojasiewicz, of a tempered distribution and the order of summability of the pointwise Fourier inversion formula are closely related. Assuming that the order of the point values and certain order of growth at infinity are given for a tempered distribution, we estimate the order of summability of the Fourier inversion formula. For Fourier series, and in other cases, it is shown that if the distribution has a distributional point value of order k, then its Fourier series is e.v. Cesaro summable to the distributional point value of order k＋1. Conversely, we also show that if the pointwise Fourier inversion formula is e.v. Cesaro summable of order k, then the distribution is the （k ＋ 1）-th derivative of a locally integrable function, and the distribution has a distributional point value of order k ＋ 2. We also establish connections between orders of summability and local behavior for other Fourier inversion problems.     

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

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, θ）.     

THE ALL-ASSOCIATIVITY OF OCTONIONS AND ITS APPLICATIONS

Using an elementary method, we give a new proof of the all-associativity of octonions. As some applications, the known Taylor theorem is improved, and a new definition and new properties of octonionic determinant are also obtained.     

NUMERICAL RECONSTRUCTION OF SMALL PERTURBATIONS IN THE ELECTROMAGNETIC COEFFICIENTS OF A DIELECTRIC MATERIAL

The aim of this paper is to solve numerically the inverse problem of reconstructing small amplitude perturbations in the magnetic permeability of a dielectric material from partial or total dynamic boundary measurements. Our numerical algorithm is based on the resolution of the time-dependent Maxwell equations, an exact controllability method and Fourier inversion for localizing the perturbations. Two-dimensional numerical experiments illustrate the performance of the reconstruction method for different configurations even in the case of limited-view data.     

SIMPLE PROOFS OF THE CAUCHY-SCHWARTZ INEQUALITY AND THE NEGATIVE DISCRIMINANTPROPERTY IN ARCHIMEDEAN ALMOST f-ALGEBRAS

The paper presents simple proofs of the Cauchy-Schwartz inequality and the negative discriminant property in archimedean almost f-algebras^[5], based on a sequence approximation.     

A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF THE CAHN-HILLIARD EQUATION AND THE HELE-SHAW FLOW

This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut ＋ △（ε△Au-ε^-1f（u）） = 0. It is shown that the a posteriori error bounds depends on ε^-1 only in some low polynomial order, instead of exponential order. Using these a posteriori error estimates, we construct at2 adaptive algorithm for computing the solution of the Cahn- Hilliard equation and its sharp interface limit, the Hele-Shaw flow. Numerical experiments are presented to show the robustness and effectiveness of the new error estimators and the proposed adaptive algorithm.     

RECURSIVE KERNELS

This paper is an extension of earlier papers [8, 9] on the ＂native＂ Hilbert spaces of functions on some domain Ωbelong toR^d Rd in which conditionally positive definite kernels are reproducing kernels. Here, the focus is on subspaces of native spaces which are induced via subsets of Ω, and we shall derive a recursive subspace structure of these, leading to recur- sively defined reproducing kernels. As an application, we get a recursive Neville-Aitken- type interpolation process and a recursively defined orthogonal basis for interpolation by translates of kernels.     

ON THE ERROR ESTIMATES OF A NEW OPERATE SPLITTING METHOD FOR THE NAVIER-STOKES EQUATIONS

In this paper, a new operator splitting scheme is introduced for the numerical solution of the incompressible Navier-Stokes equations. Under some mild regularity assumptions on the PDE solution, the stability of the scheme is presented, and error estimates for the velocity and the pressure of the proposed operator splitting scheme are given.     

A NUMERICAL STUDY FOR THE PERFORMANCE OF THE WENO SCHEMES BASED ON DIFFERENT NUMERICAL FLUXES FOR THE SHALLOW WATER EQUATIONS

In this paper we investigate the performance of the weighted essential non-oscillatory （WENO） methods based on different numerical fluxes, with the objective of obtaining better performance for the shallow water equations by choosing suitable numerical fluxes. We consider six numerical fluxes, i.e., Lax-Friedrichs, local Lax-Friedrichs, Engquist-Osher, Harten-Lax-van Leer, HLLC and the first-order centered fluxes, with the WENO finite volume method and TVD Runge-Kutta time discretization for the shallow water equations. The detailed numerical study is performed for both one-dimensional and two-dimensional shallow water equations by addressing the property, and resolution of discontinuities. issues of CPU cost, accuracy, non-oscillatory     

Boundedness of Multilinear Commutators with Rough Kernels on Morrey-Herz Spaces

This paper is concerning the commutators generated by the multilinear singular integral with rough kernels and BMO functions. The boundedness of the multilinear commutators Tb→ （f→） is established on the Morrey-Herz space by using the John-Nirenberg inequality.     

SUPERCONVERGENCE ANALYSIS OF FINITE ELEMENT METHODS FOR OPTIMAL CONTROL PROBLEMS OF THE STATIONARY BENARD TYPE

In this paper, we consider the finite element approximation of the distributed optimal control problems of the stationary Benard type under the pointwise control constraint. The states and the co-states are approximated by polynomial functions of lowest-order mixed finite element space or piecewise linear functions and the control is approximated by piecewise constant functions. We give the superconvergence analysis for the control; it is proved that the approximation has a second-order rate of convergence. We further give the superconvergence analysis for the states and the co-states. Then we derive error estimates in L^∞-norm and optimal error estimates in L^2-norm.     

ON THE QUASI-RANDOM CHOICE METHOD FOR THE LIOUVILLE EQUATION OF GEOMETRICAL OPTICS WITH DISCONTINUOUS WAVE SPEED

We study the quasi-random choice method （QRCM） for the Liouville equation of ge- ometrical optics with discontinuous locM wave speed. This equation arises in the phase space computation of high frequency waves through interfaces, where waves undergo partial transmissions and reflections. The numerical challenges include interface, contact discon- tinuities, and measure-valued solutions. The so-called QRCM is a random choice method based on quasi-random sampling （a deterministic alternative to random sampling）. The method not only is viscosity-free but also provides faster convergence rate. Therefore, it is appealing for the prob!em under study which is indeed a Hamiltonian flow. Our analy- sis and computational results show that the QRCM 1） is almost first-order accurate even with the aforementioned discontinuities; 2） gives sharp resolutions for all discontinuities encountered in the problem; and 3） for measure-valued solutions, does not need the level set decomposition for finite difference/volume methods with numerical viscosities.     

SHARP MAXIMAL FUNCTION ESTIMATES FOR MULTILINEAR SINGULAR INTEGRALS WITH NON-SMOOTH KERNELS

In this article we obtain weighted norm estimates for multilinear singular integrals with non-smooth kernels and the boundedness of certain multilinear commutators by making use of a sharp maximal function.     

ERROR ESTIMATES FOR THE RECURSIVE LINEARIZATION OF INVERSE MEDIUM PROBLEMS

This paper is devoted to the mathematical analysis of a general recursive linearization algorithm for solving inverse medium problems with multi-frequency measurements. Under some reasonable assumptions, it is shown that the algorithm is convergent with error estimates. The work is motivated by our effort to analyze recent significant numerical results for solving inverse medium problems. Based on the uncertainty principle, the recursive linearization allows the nonlinear inverse problems to be reduced to a set of linear problems and be solved recursively in a proper order according to the measurements. As an application, the convergence of the recursive linearization algorithm [Chen, Inverse Problems 13（1997）, pp.253-282] is established for solving the acoustic inverse scattering problem.     

THE OPTIMAL CONVERGENCE ORDER OF THE DISCONTINUOUS FINITE ELEMENT METHODS FOR FIRST ORDER HYPERBOLIC SYSTEMS

In this paper, a discontinuous finite element method for the positive and symmetric, first-order hyperbolic systems （steady and nonsteady state） is constructed and analyzed by using linear triangle elements, and the O（h^2）-order optimal error estimates are derived under the assumption of strongly regular triangulation and the Ha-regularity for the exact solutions. The convergence analysis is based on some superclose estimates of the interpolation approximation. Finally, we discuss the Maxwell equations in a two-dimensional domain, and numerical experiments are given to validate the theoretical results.     

A NATURAL GRADIENT ALGORITHM FOR THE SOLUTION OF LYAPUNOV EQUATIONS BASED ON THE GEODESIC DISTANCE

A new framework based on the curved Riemannian manifold is proposed to calculate the numerical solution of the Lyapunov matrix equation by using a natural gradient descent algorithm and taking the geodesic distance as the objective function. Moreover, a gradient descent algorithm based on the classical Euclidean distance is provided to compare with this natural gradient descent algorithm. Furthermore, the behaviors of two proposed algorithms and the conventional modified conjugate gradient algorithm are compared and demonstrated by two simulation examples. By comparison, it is shown that the convergence speed of the natural gradient descent algorithm is faster than both of the gradient descent algorithm and the conventional modified conjugate gradient algorithm in solving the Lyapunov equation.     

A NON-TRIVIAL PRODUCT OF FILTRATION s＋6 IN THE STABLE HOMOTOPY GROUPS OF SPHERES

By a method improving that of [1], the authors prove the existence of a non-trivial product of filtration, s ＋ 6, in the stable homotopy groups of sphere, πt-6S, which is represented up to non-zero scalar by β^-s＋2ho（hmbn-1 -hnbm-1） ∈ ExtA^s＋6,t＋s（Zp, Zp） in the Adams spectral sequence, where p ≥ 7, n ≥ m ＋ 2 ≥ 5, q = 2（p- 1）, 0 ≤ s 〈 p - 2, t= （s ＋ 2 ＋ （s ＋ 2）p ＋ p^m ＋ p^n）q. The advantage of this method is to extend the range of s without much complicated argument as in [1].     

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.