A regularization algorithm for matrices of bilinear and sesquilinear forms

We give an algorithm that uses only unitary transformations and for each square complex matrix constructs a *congruent matrix that is a direct sum of a nonsingular matrix and singular Jordan blocks.