The structure of bilinear cyclic convolution algorithms is explored over finite fields. The algorithms derived are valid for any length not divisible by the field characteristic and are based upon the small length polynomial multiplication algorithms. The multiplicative complexity of these algorithms is small and depends on the field of constants. The linear transformation matrices

(premultiplication), and

(postmultiplication) defining the algorithm have block structures which are related to one another. The rows of

and

and the columns of

are maximal length recurrent sequences. Because of the highly regular structure of

, and

, the algorithms can be very easily designed even for large lengths. The application of these algorithms to the decoding of Reed-Solomon codes is also examined.