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.

Parallel Implementations of the Jacobi Linear Algebraic Systems Solve

Parallel Algorithm for Solving Banded Linear Systems

A Parallel Algorithm Based on Galerkin Theory for Block-Tridiagonal Linear Systems

Parallel Algorithm for Solving Nonlinear System of Equations

Research and Implementation of Effective Jacobi Iteration Algorithms on GPU

A Parallel Iterative Method for Solving Periodical Block-Tridiagonal Linear Equations.

Parallel Row Operation Algorithm for Banded Linear Systems

Iterative Solvers for Linear Systems on Parallel Computers

Parallel Algorithms for Solving Semi Linear Elliptic Equations

A Parallel Method With Hybrid Algorithms For Mixed Integer Nonlinear Programming

Parallel Asynchronous Team Algorithms: Convergence and Performance Analysis

Jacobi solver: A fast FPGA-based engine system for Jacobi method

Parallel Block-iterative Methods: Application to the Load Flow Problem

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

Investigation of implement of Jacobi iteration on GPU

Iterative Algorithm for Linear Equations on GPU

Improved Parallel Arnoldi Method for Solving Linear Equations

Parallel Algorithm with Parameters Based on Alternating Direction for Solving Banded Linear Systems

Optimizing Algorithm of Sparse Linear Systems on GPU

J-Parallelio -- automatic parallelization framework for Java virtual machine code

Parallel least squares problems on massively distributed memory computers