Approximate inverse computation using Frobenius inner product

Parallel preconditioners are presented for the solution of general linear systems of equations. The computation of these preconditioners is achieved by orthogonal projections related to the Frobenius inner product. So, min_M∈??∥AM?I∥ _F and matrix M₀∈?? corresponding to this minimum (?? being any vectorial subspace of ??_n(?)) are explicitly computed using accumulative formulae in order to reduce computational cost when subspace ?? is extended to another one containing it. Every step, the computation is carried out taking advantage of the previous one, what considerably reduces the amount of work. These general results are illustrated with the subspace of matrices M such that AM is symmetric. The main application is developed for the subspace of matrices with a given sparsity pattern which may be constructed iteratively by augmenting the set of non‐zero entries in each column. Finally, the effectiveness of the sparse preconditioners is illustrated with some numerical experiments. Copyright © 2002 John Wiley & Sons, Ltd.     

Classification problems for shifts on modules over a principal ideal domain

In this paper we study symbolic dynamics over alphabets which are modules over a principal ideal domain, considering topological shifts which are also submodules. Our main result is the classification, up to algebraic and topological conjugacy, of the torsion-free, transitive, finite memory shifts.

     

匀速运动的电偶极子的电场和磁场分布

本文通过对电偶极子电磁势的分析和相对论下电磁场在不同惯性系中的协变分布，计算了作匀速直线运动的电偶极子的电场和磁场的空间分布．     

EXISTENCE OF PERIODIC SOLUTIONS OF SYSTEM OF DIFFERENCE EQUATIONS

This paper studies the nonautonomous nonlinear system of difference equationsΔx(n)=A(n)x(n)+f(n,x(n)),n∈Z,(*) where x(n)∈R~N,A(n)=(a_(ij)(n))N×N is an N×N matrix,with a-(ij)∈C(R,R) for i,j= 1,2,3,...,N,and f=(f_1,f_2,...,f_N)~T∈C(R×R~N,R~N),satisfying A(t+ω)=A(t),f(t+ω,z)=f(t,z) for any t∈R,(t,z)∈R×R~N andωis a positive integer.Sufficient conditions for the existence ofω-periodic solutions to equations (*) are obtained.     

镧系金属一氢化物有效芯势研究

采用Cundari和Stevens等推导的有效芯势对镧系金属一氢化物进行了理论计算，以探讨镧系金属元素与氢的相互作用。结果表明所有镧系金属一氢化物基态时理论上是稳定的，最稳定的是SmH,最不稳定的是DyH;键长计算结果显示，基态时镧系金属一招兵买马花物有独立王国 收缩现象发生；红外振动频率理论计算值与实验结果一致；成键轨道中，金属原子轨道的贡献主要是s轨道和d轨道：从CeH至ErH(GdH)例外）随着外层电子的增加s轨道成分逐渐增大d轨道成分逐渐减小；从TmH和LuH（包括GdH)，成键轨道中金属原子轨道的贡献主要是d轨道，约为90％；约大多数镧系金属一氧化物的成键轨道中金属原子轨道f成分小于1％。     

B_2分子X~3Σ_g~-和A~3Σ_u~-态的势能函数

使用SAC/SAC-CI方法,利用D95(d),6-311g**以及cc-PVTZ等基组,对B2分子的基态(X3Σg-)和第一激发态(A3Σu-)的平衡结构和谐振频率进行了优化计算.通过对3个基组的计算结果的比较,得出了D95(d)基组为3个基组中的最优基组的结论;使用D95(d)基组,利用SAC的GSUM(Group Sum of Operators)方法对基态(X3Σg-),SAC-CI的GSUM方法对激发态(A3Σu-)进行单点能扫描计算,用正规方程组拟合Murrell-Sorbie函数,得到了相应电子态的完整势能函数;从得到的势能函数计算了与基态(X3Σg-)和第一激发态(A3Σu-)相对应的光谱常数(Be,αe,ωe和ωeχe),结果与实验数据吻合.     

On a class of new and practical performance indexes for approximation of fold bifurcations of nonlinear power flow equations

Efficient measurement of the performance index (the distance of a loading parameter from the voltage collapse point) is one of the key problems in power system operations and planning and such an index indicates the severity of a power system with regard to voltage collapse. There exist many interesting methods and ideas to compute this index. However, some successful methods are not yet mathematically justified while other mathematically sound methods are often proposed directly based on the bifurcation theory and they require the initial stationary state to be too close to the unknown turning point to make the underlying methods practical.This paper first gives a survey of several popular methods for estimating the fold bifurcation point including the continuation methods, bifurcation methods and the test function methods (Seydel's direct solution methods, the tangent vector methods and the reduced Jacobian method) and discuss their relative advantages and problems. Test functions are usually based on scaling of the determinant of the Jacobian matrix and it is generally not clear how to determine the behaviour of such functions. As the underlying nonlinear equations are of a particular type, this allows us to do a new analysis of the determinants of the Jacobian and its submatrices in this paper. Following the analysis, we demonstrate how to construct a class of test functions with a predictable analytical behaviour so that a suitable index can be produced. Finally, examples of two test functions from this class are proposed. For several standard IEEE test systems, promising numerical results have been achieved.     

Instability in a delay-PDE model of earthquake occurrence

Summary Ductile deformation prior to brittle fracture in rocks causes fracture to take place with a time delay after the critical stress for fracture is reached, in the presence of an increasing load stress. We discuss the stability of a stochastic model of interactive earthquake occurrence under the influence of time delays resulting from the ductile process. A threshold for oscillatory behavior is found for large enough coupling and discrete time delays. The system is stable if the time delays are spread over a broad time interval, even for large coupling.     

Stabilizing management of fluctuating resources

We consider a resource management problem in which the management objective is to minimize fluctuations in resource stocks. Stabilizing management policies constitute the designing of memoryless state feedback control strategies for a discrete time resource model which contains unknown but bounded fluctuations. We also show that the problem of maximizing sustainable yield in an uncertain fishery can be considered as the problem of stabilizing of the stock level.The paper corresponds to an invited talk at the 14th Symposium on Operations Research, Ulm, September 6–8, 1989.The support by Yrjö Jahnsson Foundation is gratefully acknowledged.This work was supported by NSF and AFOSR under grant ECS 8602524.     

相对论性无自旋粒子在Hartmann势场中运动的精确解

在标量势等于矢量势的条件下,本文获得了具有Hartmann型势的Klein-Gordon方程的精确解.给出了束缚态的精确的能谱方程和归一化的径向波函数,对于散射态,获得了按“k/2π标度”归一化的径向波函数和相移的解析计算公式.讨论了散射振幅的解析性质和波函数、能谱方程以及相移的非相对论近似.