Parallel implementation of the semi-conjugate residual method for systems of linear algebraic equations

A preconditioned semi-conjugate residual algorithm for solving systems of algebraic equations with a nonsymmetric square matrix is presented. Most of the calculations in the algorithm is accounted for vector operations whose number grows quadratically with the number of stored directing vectors. A parallel implementation of the algorithm is studied experimentally, the Eisenstat modification being chosen as the preconditioning. An algorithm for parallelizing matrix operations is presented for special problems.