Multisplitting for regularized least squares with Krylov subspace recycling

The method of multisplitting (MS), implemented as a restricted additive Schwarz type algorithm, is extended for the solution of regularized least squares problems. The presented non‐stationary version of the algorithm uses dynamic updating of the weights applied to the subdomains in reconstituting the global solution. Standard convergence results follow from extensive prior literature on linear MS schemes. Additional convergence results on nonstationary iterations yield convergence conditions for the presented nonstationary MS algorithm. The global iteration uses repeated solves of local problems with changing right hand sides but a fixed system matrix. These problems are solved inexactly using a conjugate gradient least squares algorithm which provides a seed Krylov subspace. Recycling of the seed system Krylov subspace to obtain the solutions of subsequent nearby systems of equations improves the overall efficiency of the MS algorithm, and is apparently novel in this context. The obtained projected solution is not always of sufficient accuracy to satisfy a reasonable inner convergence condition on the local solution. Improvements to accuracy may be achieved by reseeding the solution space either every few steps, or when the successive right hand sides are sufficiently close as measured by a provided tolerance. Restarting and augmenting the solution space are also discussed. Any time a new space is generated it is used for subsequent steps. Numerical simulations validate the use of the recycling algorithm. These numerical experiments use the standard reconstruction of the two dimensional Shepp–Logan phantom, as well as a two dimensional problem from seismic tomography. Copyright © 2011 John Wiley & Sons, Ltd.     

Existence and uniqueness for a nonlinear inverse reaction‐diffusion problem with a nonlinear source in higher dimensions

This paper analyzes the existence and the uniqueness problem for an n‐dimensional nonlinear inverse reaction‐diffusion problem with a nonlinear source. A transformation is used to obtain a new inverse coefficient problem. Then, a parabolic differential operator L_λ is defined to establish the relation between the solution of L_λ = 0 and the new inverse problem. Following this, it is shown that the inverse problem has at least one solution in the class of admissible coefficients. Furthermore, it is proved that this solution is the unique solution of the undertaken inverse problem. A numerical example is given to illustrate ill‐posedness of the inverse problem. Copyright © 2013 John Wiley & Sons, Ltd.     

Reconstructing the past from imprecise knowledge of the present: Effective non‐uniqueness in solving parabolic equations backward in time

Identifying sources of ground water pollution and deblurring astronomical galaxy images are two important applications generating growing interest in the numerical computation of parabolic equations backward in time. However, while backward uniqueness typically prevails in parabolic equations, the precise data needed for the existence of a particular backward solution is seldom available. This paper discusses previously unexplored non‐uniqueness issues, originating from trying to reconstruct a particular solution from imprecise data. Explicit 1D examples of linear and nonlinear parabolic equations are presented, in which there is strong computational evidence for the existence of distinct solutions w^red(x,t) and w^green(x,t), on 0 ≤ t ≤ 1. These solutions have the property that the traces w^red(x,1) and w^green(x,1) at time t = 1 are close enough to be visually indistinguishable, while the corresponding initial values w^red(x,0) and w^green(x,0) are vastly different, well‐behaved, physically plausible functions, with comparable L² norms. This implies effective non‐uniqueness in the recovery of w^red(x,0) from approximate data for w^red(x,1). In all these examples, the Van Cittert iterative procedure is used as a tool to discover unsuspected, valid, additional solutions w^green(x,0). This methodology can generate numerous other examples and indicates that multidimensional problems are likely to be a rich source of striking non‐uniqueness phenomena. Published 2012. This article is a US Government work and is in the public domain in the USA.     

样条最小二乘拟合在估计水塔水流量中的应用

本文采用样条最小二乘法 ,通过对水塔水位的观测数据进行拟合 ,能够估计出任意时刻水塔的水流量及一天内的总用水量 ,同时较好地克服了使用高次多项式时容易导致的正规方程组的病态性问题 .     

Determining the temperature from incomplete boundary data

An iterative procedure for determining temperature fields from Cauchy data given on a part of the boundary is presented. At each iteration step, a series of mixed well‐posed boundary value problems are solved for the heat operator and its adjoint. A convergence proof of this method in a weighted L²‐space is included, as well as a stopping criteria for the case of noisy data. Moreover, a solvability result in a weighted Sobolev space for a parabolic initial boundary value problem of second order with mixed boundary conditions is presented. Regularity of the solution is proved. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Characterization of l-cysteine capped CdTe quantum dots and application to test Cu(II) deficiency in biological samples from critically ill patients

The catalytic activity of copper ion gives, from the physiological point of view, a central role in many biological processes. Variations in the composition and location of cellular copper have been addressed given their physiological and pathological consequences. In this paper l-cysteine capped CdTe quantum dots is used for the fluorimetric determination of Cu(II) in biological samples from healthy individuals and patients admitted to the Intensive Care Units (ICU). An acceptable homogeneity in the CdTe QDs size has been obtained with an average value of 3 nm. No significant alterations in the spectral properties were observed for 2 months when stored in vacutainers at 6 °C and a concentration of approximately 2 μM. Data from oxidative stress markers such superoxide dismutase, total antioxidant capacity and DNA damage can be correlated with a Cu(II) deficiency for the ICU patients as measured by flame-atomic absorption spectroscopy (FAAS) and inductively coupled plasma source mass spectrometry (ICP-MS). Aqueous solutions 0.3 μM of l-cysteine capped CdTe QDs in MOPS buffer (6 mM, pH 7.4) used at 21 °C in the range 15–60 min after preparation of the sample for the measurements of fluorescence gives contents in Cu(II) for erythrocytes in good agreement with those obtained in FAAS and ICP-MS but the comparative ease of use makes the fluorimetric technique more suitable than the other two techniques for routine analysis.     

Adaptive iterative regularization schemes for two‐dimensional integral‐algebraic systems

The main idea of this paper is to utilize the adaptive iterative schemes based on regularization techniques for moderately ill‐posed problems that are obtained by a system of linear two‐dimensional Volterra integral equations with a singular matrix in the leading part. These problems may arise in the modeling of certain heat conduction processes as well as in the dynamic simulation packages such as compressible flow through a plant piping network. Owing to the ill‐posed nature of the first kind Volterra equation that appears in the system, we will focus on the two families of regularization algorithms, ie, the Landweber and Lavrentiev type methods, where we treat both the exact and perturbed data. Our aim is to work directly with the original Volterra equations without any kind of reduction. Two fast iterative algorithms with reasonable computational complexity are developed. Numerical experiments on a few test problems are used to illustrate the validity and efficiency of the proposed iterative methods in comparison with the classical regularization methods.     

Local well‐posedness and ill‐posedness for the fractal Burgers equation in homogeneous Sobolev spaces

In this paper, we consider local well‐posedness and ill‐posedness questions for the fractal Burgers equation. First, we obtain the well‐posedness result in the critical Sobolev space. We also present an unconditional uniqueness result. Second, we show the ill‐posedness from the point of view of the Picard iterations in the supercritical space. Copyright © 2008 John Wiley & Sons, Ltd.     

关于不可压、无粘流体的Euler方程初值问题的适定性(Ⅱ)

以分层理论为基础,讨论了不可压、无粘流体的Euler方程的形式可解性,并给出了各类不适定初值问题存在形式解的条件与计算方法。并讨论了超曲面上和超平面上初值问题的适定性并给出了存在不唯一解的例证。     

关于不可压、无粘流体的Euler方程初值问题的适定性(Ⅰ)

以分层理论为基础,讨论了Euler方程不适定的初值问题以及不适定问题的形式可解性,并给出了某些不适定初值问题存在形式解的条件与计算方法。特别讨论了R4中的超平面{t=0}上初值问题的适定性并给出了存在不唯一解的例证。