Solution of large‐scale weighted least‐squares problems

A sequence of least‐squares problems of the form min_y∥G^1/2(A^T y?h)∥₂, where G is an n×n positive‐definite diagonal weight matrix, and A an m×n (m?n) sparse matrix with some dense columns; has many applications in linear programming, electrical networks, elliptic boundary value problems, and structural analysis. We suggest low‐rank correction preconditioners for such problems, and a mixed solver (a combination of a direct solver and an iterative solver). The numerical results show that our technique for selecting the low‐rank correction matrix is very effective. Copyright © 2002 John Wiley & Sons, Ltd.     

Coupled intervals in the discrete calculus of variations: necessity and sufficiency

In this work we study nonnegativity and positivity of a discrete quadratic functional with separately varying endpoints. We introduce a notion of an interval coupled with 0, and hence, extend the notion of conjugate interval to 0 from the case of fixed to variable endpoint(s). We show that the nonnegativity of the discrete quadratic functional is equivalent to each of the following conditions: The nonexistence of intervals coupled with 0, the existence of a solution to Riccati matrix equation and its boundary conditions. Natural strengthening of each of these conditions yields a characterization of the positivity of the discrete quadratic functional. Since the quadratic functional under consideration could be a second variation of a discrete calculus of variations problem with varying endpoints, we apply our results to obtain necessary and sufficient optimality conditions for such problems. This paper generalizes our recent work in [R. Hilscher, V. Zeidan, Comput. Math. Appl., to appear], where the right endpoint is fixed.     

Conjugate addition of pyrroles to α,β-unsaturated ketones using copper bromide as a catalyst

Copper bromide was used as a catalyst for the addition of pyrroles to enones. When both the reactants were used in equimolar amounts, mono and dialkylated products were obtained. However, the use of excess enone furnished only dialkylated products. Thus, copper bromide was shown to be an efficient catalyst for the dialkylation of pyrroles.     

An inner product lemma

Given the operator product BA in which both A and B are symmetric positive‐definite operators, for which symmetric positive‐definite operators C is BA symmetric positive‐definite in the C inner product 〈x, y〉C? This question arises naturally in preconditioned iterative solution methods, and will be answered completely here. Copyright © 2004 John Wiley & Sons, Ltd.     

二维椭圆型方程反问题中优化算法的比较

近年来，决定椭圆型方程系数反问题在地磁、地球物理、冶金和生物等实际问题上有着广泛的应用．本文讨论了二维的决定椭圆型方程系数反问题的数值求解方法．由误差平方和最小原则，这个反问题可化为一个变分问题，并进一步离散化为一个最优化问题，其目标函数依赖于要决定的方程系数．本文着重考察非线性共轭梯度法在此最优化问题数值计算中的表现，并与拟牛顿法作为对比．为了提高算法的效率我们适当选择加快收敛速度的预处理矩阵．同时还考察了线搜索方法的不同对优化算法的影响．数值实验的结果表明，非线性共轭梯度法在这类大规模优化问题中相对于拟牛顿法更有效．     

Differentiability of the Value Function of Nonclassical Optimal Growth Models

We consider an optimal growth (multi-sector) model with nonconvex technology. Using the Clarke results on generalized gradients, we prove that the value function has left and right derivatives with respect to the initial capital stock, without requiring supermodularity assumptions.     

Commutators, C0-semigroups and resolvent estimates

We study the existence and the continuity properties of the boundary values on the real axis of the resolvent of a self-adjoint operator H in the framework of the conjugate operator method initiated by Mourre. We allow the conjugate operator A to be the generator of a C₀-semigroup (finer estimates require A to be maximal symmetric) and we consider situations where the first commutator [H,iA] is not comparable to H. The applications include the spectral theory of zero mass quantum field models.     

Optimal shape design for packaging containing heating elements by inverse heat transfer method

The present study is concerned with the determination of the optimal shape for a package containing multiple heating elements. The optimization tool has been developed based on the inverse heat transfer (IHT) approach, incorporating a direct problem solver, a numerical grid generator, a direct-differentiation sensitivity analyzer, and the conjugate gradient method. Shape design that leads to a specified outer surface temperature distribution is predicted by the approach. In this study, the effects of internal heat generation on optimal shapes of the packagings have also been evaluated. Several practical cases with various imbedded heating elements and thermal conditions are studied. Results show that the approach provides an efficient computer-aided design scheme for the shape profile determination.     

一类共轭梯度法的全局收敛性结果

本文证明了在Grippo-Lucidi线搜索下当βk取βk=σ1βPRPk+σ2βnewk,其中σ1≥0,σ2≥0,σ1+σ2＞0,βnewk=gTk(gk-gk-1)/-dTk-1gk-1时一类共轭梯度法的全局收敛性,并给出了此类方法良好的数值效果.