An overlapping additive Schwarz preconditioner for interpolation on the unit sphere with spherical radial basis functions

The problem of interpolation of scattered data on the unit sphere has many applications in geodesy and Earth science in which the sphere is taken as a model for the Earth. Spherical radial basis functions provide a convenient tool for constructing the interpolant. However, the underlying linear systems tend to be ill-conditioned. In this paper, we present an additive Schwarz preconditioner for accelerating the solution process. An estimate for the condition number of the preconditioned system will be discussed. Numerical experiments using MAGSAT satellite data will be presented.     

A fractional partial differential equation based multiscale denoising model for texture image

Traditional integer‐order partial differential equation based image denoising approach can easily lead edge and complex texture detail blur, thus its denoising effect for texture image is always not well. To solve the problem, we propose to implement a fractional partial differential equation (FPDE) based denoising model for texture image by applying a novel mathematical method—fractional calculus to image processing from the view of system evolution. Previous studies show that fractional calculus has some unique properties that it can nonlinearly enhance complex texture detail in digital image processing, which is obvious different with integer‐order differential calculus. The goal of the modeling is to overcome the problems of the existed denoising approaches by utilizing the aforementioned properties of fractional differential calculus. Using classic definition and property of fractional differential calculus, we extend integer‐order steepest descent approach to fractional field to implement fractional steepest descent approach. Then, based on the earlier fractional formulas, a FPDE based multiscale denoising model for texture image is proposed and further analyze optimal parameters value for FPDE based denoising model. The experimental results prove that the ability for preserving high‐frequency edge and complex texture information of the proposed fractional denoising model are obviously superior to traditional integral based algorithms, as for texture detail rich images. Copyright © 2013 John Wiley & Sons, Ltd.     

Neighboring extremals of dynamic optimization problems with path equality constraints

Neighboring extremals of dynamic optimization problems with path equality constraints and with an unknown parameter vector are considered in this paper. With some simplifications, the problem is reduced to solving a linear, time-varying two-point boundary-value problem with integral path equality constraints. A modified backward sweep method is used to solve this problem. Two example problems are solved to illustrate the validity and usefulness of the solution technique. This research was supported in part by the National Aeronautics and Space Administration under NASA Grant No. NCC-2-106. The author is indebted to Professor A. E. Bryson, Jr., Department of Aeronautics and Astronautics, Stanford University, for many stimulating discussions.     

On a Property of Subexponential Distributions

Let F be a distribution function (d.f.) on [0, ) with finite first moment m >0. We define the integrated tail distribution function F ₁ of F by F ₁(t)=m^-1 ₀ ^t (1- F(u))du, t0. In this paper, we obtain sufficient conditions under which implications FSF ₁S and F ₁S FS hold, where S is the class of subexponential distributions.     

Generalized One-Sided Laws for the Larger Observations of a Triangular Array

Consider independent and identically distributed random variables {X,X _nj , 1jn,n1} with density f(x)=px ^–p–1 I(x1), where p>0. We show that there exist unusual generalized Laws of the Iterated Logarithm involving the larger order statistics from our array.     

GLOBAL CLASSICAL SOLUTIONS WITH SMALL INITIAL TOTAL VARIATION FOR QUASILINEAR HYPERBOLIC SYSTEMS

By means of the continuous Glimm functional,a proof is given on the global existence ofclassical solutions to Cauchy problem for general first order quasilinear hyperbolic systems withsmall initial total rariation.     

The effect of a surface energy term on the distribution of phases in an elastic medium with a two‐well elastic potential

We consider the problem of minimizing among functions u:?^d?Ω→?^d, u_∣?Ω=0, and measurable subsets E of Ω. Here f_h⁺, f^? denote quadratic potentials defined on Ω¯×{symmetric d×d matrices}, h is the minimum energy of f_h⁺ and ε(u) is the symmetric gradient of the displacement field u. An equilibrium state û, Ê of J(u,E) is called one‐phase if E=?? or E=Ω, two‐phase otherwise. For two‐phase states, σ∣?E∩Ω∣ measures the effect of the separating surface, and we investigate the way in which the distribution of phases is affected by the choice of the parameters h??, σ>0. Additional results concern the smoothness of two‐phase equilibrium states and the behaviour of inf J(u,E) in the limit σ↓0. Moreover, we discuss the case of additional volume force potentials, and extend the previous results to non‐zero boundary values. Copyright © 2002 John Wiley & Sons, Ltd.     

On the covariance of the level sizes in random recursive trees

In this paper we study the covariance structure of the number of nodes k and l steps away from the root in random recursive trees. We give an analytic expression valid for all k, l and tree sizes N. The fraction of nodes k steps away from the root is a random probability distribution in k. The expression for the covariances allows us to show that the total variation distance between this (random) probability distribution and its mean converges in probability to zero. © 2002 Wiley Periodicals, Inc. Random Struct. Alg., 20: 519–539, 2002     

The Integrated Density of States for Some Random Operators with Nonsign Definite Potentials

We study the integrated density of states of random Anderson-type additive and multiplicative perturbations of deterministic background operators for which the single-site potential does not have a fixed sign. Our main result states that, under a suitable assumption on the regularity of the random variables, the integrated density of states of such random operators is locally Hölder continuous at energies below the bottom of the essential spectrum of the background operator for any nonzero disorder, and at energies in the unperturbed spectral gaps, provided the randomness is sufficiently small. The result is based on a proof of a Wegner estimate with the correct volume dependence. The proof relies upon the L^p-theory of the spectral shift function for p?1 (Comm. Math. Phys.218 (2001), 113-130), and the vector field methods of Klopp (Comm. Math. Phys.167 (1995), 553-569). We discuss the application of this result to Schrödinger operators with random magnetic fields and to band-edge localization.     

Contingent derivatives and regularization for noncoercive inverse problems

Abstract

We study the inverse problem of parameter identification in noncoercive variational problems that commonly appear in applied models. We examine the differentiability of the set-valued parameter-to-solution map using the first-order and the second-order contingent derivatives. We explore the inverse problem using the output least-squares and the modified output least-squares objectives. By regularizing the noncoercive variational problem, we obtain a single-valued regularized parameter-to-solution map and investigate its smoothness and boundedness. We also consider optimization problems using the output least-squares and the modified output least-squares objectives for the regularized variational problem. We give a complete convergence analysis showing that for the output least-squares and the modified output least-squares, the regularized minimization problems approximate the original optimization problems suitably. We also provide the first-order and the second-order adjoint method for the computation of the first-order and the second-order derivatives of the output least-squares objective. We provide discrete formulas for the gradient and the Hessian calculation and present numerical results.