The commonly used preconditioning technique is incomplete LU（ILU） decomposition

which is not suitable for parallelization due to the inherent model of forward and back substitutions.Therefore

the parallel computation method of power system transient stability based on preconditioned generalized minimal residual algorithm（GMRES） method with hybrid LU（HLU） decomposition is derived.On the basis of the implicit simultaneous solution method

and combined with the block bordered diagonal form

the Jacobian matrix is preconditioned with hybrid LU decomposition by the technique of dividing and ruling

and then the GMRES method is used to solve the linear equations.Using two different scale systems

the test results show that the iterations of HLU decomposition preconditioning method can almost reduce to the half of ILU decomposition preconditioning method

and with graphics processing unit（GPU） parallel technology

the preconditioned GMRES method with hybrid LU decomposition can obtain more than 10 times speedup when the system scale is up to tens of thousands of nodes.Therefore

the proposed method has not only good convergence rate but also high spatial parallelism

which is more suitable for GPU computing.