RICHARDSON EXTRAPOLATION AND DEFECT CORRECTION OF FINITE ELEMENT METHODS FOR OPTIMAL CONTROL PROBLEMS

Asymptotic error expansions in H^1-norm for the bilinear finite element approximation to a class of optimal control problems are derived for rectangular meshes. With the rectan- gular meshes, the Richardson extrapolation of two different schemes and an interpolation defect correction can be applied. The higher order numerical approximations are used to generate a posteriori error estimators for the finite element approximation.     

ON BLOCK PRECONDITIONERS FOR PDE-CONSTRAINED OPTIMIZATION PROBLEMS

Recently, Bal proposed a block-counter-diagonal and a block-counter-triangular precon- ditioning matrices to precondition the GMRES method for solving the structured system of linear equations arising from the Galerkin finite-element discretizations of the distributed control problems in （Computing 91 （2011） 379-395）. He analyzed the spectral properties and derived explicit expressions of the eigenvalues and eigenvectors of the preconditioned matrices. By applying the special structures and properties of the eigenvector matrices of the preconditioned matrices, we derive upper bounds for the 2-norm condition numbers of the eigenvector matrices and give asymptotic convergence factors of the preconditioned GMRES methods with the block-counter-diagonal and the block-counter-triangular pre- conditioners. Experimental results show that the convergence analyses match well with the numerical results.     

HIGHLY OSCILLATORY DIFFUSION-TYPE EQUATIONS

We explore new asymptotic-numeric solvers for partial differential equations with highly oscillatory forcing terms. Such methods represent the solution as an asymptotic series, whose terms can be evaluated by solving non-oscillatory problems and they guarantee high accuracy at a low computational cost. We consider two forms of oscillatory forcing terms, namely when the oscillation is in time or in space： each lends itself to different treatment. Numerical examples highlight the salient features of the new approach.     

ON NEWTON-HSS METHODS FOR SYSTEMS OF NONLINEAR EQUATIONS WITH POSITIVE-DEFINITE JACOBIAN MATRICES

The Hermitian and skew-Hermitian splitting （HSS） method is an unconditionally convergent iteration method for solving large sparse non-Hermitian positive definite system of linear equations. By making use of the HSS iteration as the inner solver for the Newton method, we establish a class of Newton-HSS methods for solving large sparse systems of nonlinear equations with positive definite Jacobian matrices at the solution points. For this class of inexact Newton methods, two types of local convergence theorems are proved under proper conditions, and numerical results are given to examine their feasibility and effectiveness. In addition, the advantages of the Newton-HSS methods over the Newton-USOR, the Newton-GMRES and the Newton-GCG methods are shown through solving systems of nonlinear equations arising from the finite difference discretization of a two-dimensional convection-diffusion equation perturbed by a nonlinear term. The numerical implemen- tations also show that as preconditioners for the Newton-GMRES and the Newton-GCG methods the HSS iteration outperforms the USOR iteration in both computing time and iteration step.     

NORM CALCULATIONS OF A CLASS OF INTEGRAL OPERATORS

In this paper we calculate the norm of a special class of integral operators acting on LP （C^n, dvs）, where dvs is the Gaussian measure on C^n.     

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.     

BOUNDEDNESS OF COMMUTATORS OF GENERALIZED FRACTIONAL INTEGRAL OPERATORS ON WEIGHTED HERZ-TYPE HARDY SPACES

In this paper, we study the boundedness of higher order commutators of gen- eralized fractional integral operators on weighted Lp spaces and Herz-type Hardy spaces.     

SOME RESULTS CONCERNING GROWTH OF POLYNOMIALS

In this paper, we establish inequalities for polynomials with restricted zeros, which in particular yields interesting generalizations of some Zygmund type inequalities for polynomial.     

Strong Convergence of an Implicit Iteration Process for a Finite Family of Asymptotically Ф-pseudocontractive Mappings

Strong convergence theorems for approximation of common fixed points of asymptotically Ф-quasi-pseudocontractive mappings and asymptotically C-strictly- pseudocontractive mappings are proved in real Banach spaces by using a new composite implicit iteration scheme with errors. The results presented in this paper extend and improve the main results of Sun, Gu and Osilike published on J. Math. Anal. Appl.     

ON APPROXIMATION PROPERTIES OF NON-CONVOLUTION TYPE NONLINEAR INTEGRAL OPERATORS-Dedicated to Professor Paul L. Butzer on his 80th Birthday

In the present paper we state some approximation theorems concerning point- wise convergence and its rate for a class of non-convolution type nonlinear integral opera- tors of the form：Tλ（f;x）=B∫AKλ（t,x,f（t））dr,x∈〈a,b〉λλA.In particular, we obtain the pointwise convergence and its rate at some characteristic points x0 off as （x,λ） → （x0, λ0） in LI 〈A,B 〉, where 〈 a,b 〉 and 〈A,B 〉 are is an arbitrary intervals in R, A is a non-empty set of indices with a topology and X0 an accumulation point of A in this topology. The results of the present paper generalize several ones obtained previously in the papers [191-[23]     

ALTERNATELY LINEARIZED IMPLICIT ITERATION METHODS FOR SOLVING QUADRATIC MATRIX EQUATIONS

A numerical solution of the quadratic matrix equations associated with a nonsingular M-matrix by using the alternately linearized implicit iteration method is considered. An iteration method for computing a nonsingular M-matrix solution of the quadratic matrix equations is developed, and its corresponding theory is given. Some numerical examples are provided to show the efficiency of the new method.     

A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS

We propose an effective stopping criterion for higher-order fast sweeping schemes for static Hamilton-Jacobi equations based on ratios of three consecutive iterations. To design the new stopping criterion we analyze the convergence of the first-order Lax-Friedrichs sweeping scheme by using the theory of nonlinear iteration. In addition, we propose a fifth-order Weighted PowerENO sweeping scheme for static Hamilton-Jacobi equations with convex Hamiltonians and present numerical examples that validate the effectiveness of the new stopping criterion.     

A New Kind of Iteration Method for Finding Approximate Periodic Solutions to Ordinary Differential Equations

In this paper, a new kind of iteration technique for solving nonlinear ordinary differential equations is described and used to give approximate periodic solutions for some well-known nonlinear problems. The most interesting features of the proposed methods are its extreme simplicity and concise forms of iteration formula for a wide range of nonlinear problems.     

Weighted Estimates for the Maximal Commutator of Singular Integral Operator on Spaces of Homogeneous Type

Weighted estimates with general weights are established for the maximal operator associated with the commutator generated by singular integral operator and BMO function on spaces of homogeneous type, where the associated kernel satisfies the HSlder condition on the first variable and some condition which is fairly weaker than the Holder condition on the second variable.     

Liouville type theorems for a system of integral equations on upper half space

In this paper,we consider the following system of integral equations on upper half space {u(x) = ∫Rn + (1/|x-y|n-α-1/|-y|n-α) λ1up1(y) + μ1vp2(y) + β1up3(y)vp4(y) dy;v(x) = ∫Rn + (1/|x-y|n-α-1/|-y|n-α)(λ2uq1(y) + μ2vq2(y) + β2uq3(y)vq4(y) dy,where Rn + = {x =(x1,x2,...,xn) ∈ Rn|xn 0}, =(x1,x2,...,xn-1,-xn) is the reflection of the point x about the hyperplane xn= 0,0 α n,λi,μi,βi≥ 0(i = 1,2) are constants,pi≥ 0 and qi≥ 0(i = 1,2,3,4).We prove the nonexistence of positive solutions to the above system with critical and subcritical exponents via moving sphere method.     

A Method for Solving Fredholm Integral Equations of the First Kind Based on Chebyshev Wavelets

In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] are used and are utilized as a basis in Galerkin method to approximate the solution of integral equations. Then, in some examples the mentioned wavelets are compared with each other.     

ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS

The main purpose of this work is to provide a novel numerical approach for the Volterra integral equations based on a spectral approach. A Legendre-collocation method is proposed to solve the Volterra integral equations of the second kind. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical results confirm the theoretical prediction of the exponential rate of convergence. The result in this work seems to be the first successful spectral approach （with theoretical justification） for the Volterra type equations.     

Realization of numbers as the degrees of maps between manifolds

We determine the set of degrees between some classes of oriented closed （n - 1）-connected 2n-manifolds by using the arithmetic theory of quadratic forms.     

Bifurcations of limit cycles from a quintic Hamiltonian system with a heteroclinic cycle

In this paper,we consider Li′enard systems of the form dx/dt=y,dy/dt=x+bx3-x5+ε(α+βx2+γx4)y,where b∈R,0|ε|1,(α,β,γ)∈D∈R3 and D is bounded.We prove that for |b|1(b0) the least upper bound of the number of isolated zeros of the related Abelian integrals I(h)=∮Γh(α+βx2+γx4)ydx is 2(counting the multiplicity) and this upper bound is a sharp one.     

NON-STANDARD COMMUTATORS FOR ROUGH OSCILLATORY SINGULAR INTEGRALS

In this paper, the author studies a class of non-standard commutators with higher order remainders for oscillatory singular integral operators with phases more general than polynomials. For 1 〈 p 〈 ∞, the L^p-boundedness of such operators are obtained provided that their kernels belong to the spaces L^q（s^n-1） for some q 〉 1.