On Parallel Implementation of Dynamic Equations for Solving Eigenvalue Problems of Complex Matrices

黄祖兰,蒋耀林,陈明敏,林小拉

DOI: https://doi.org/10.1142/9789812792037_0033

2002-01-01

Abstract:This paper provides parallel implementation of dynamic equations for solving complex eigenvalue problems, which are prevalent in engineering fields. The principle of the approach is that the eigenvalue problems is first transformed into the associated problems of solving large systems of nonlinear dynamic equations by optimization techniques and then waveform relaxation in circuit simulation is used to compute the resulting cubic dynamic equations.Because the stable state solutions of the nonlinear dynamic equations correspond to the solutions of the eigenvalue problems, often these dynamic equations were mapped into a neural network, then the transient response of the neural network is obtained by hardware or software implementation of the neural network. This paper introduces the waveform relaxation method to numerically compute the transient solutions of the dynamic equations in order to avoid the complicated neural network circuit design. The waveform relaxation method is a "divide and conquer" approach, that is, at a differential equation level it decouples a large dynamic system into some small dynamic subsystems in the time domain and these subsystems are then solved independently from each other in modern parallel processor environments. The proposed algorithm can efficiently determine all eigenvalues and their eigenvectors for a general complex matrix. It is the first time that the waveform relaxation method is used to compute the cubic dynamic equations. The algorithm is successfully implemented on the IBM RS/6000 SuperPOWER2 system, which is a message passing multi processor. The computer language is ANSI C. The parallel MPI software is applied to message passing between different processors. Authors run the parallel program on different numbers of processors, and acquire the speedup corresponding to the numbers of processors in use respectively. As illustrated by the experimental results on general complex matrices with dimensions 3000×3000, the algorithm of the paper is highly parallelized. The speedup of the parallel algorithm has a linear characteristic property. And the communication overhead of parallel processing is rather low. In view of the linear speedup and the low communication overhead, the parallel implementation reported here is effective for transient computations of large dynamic systems in practical applications.