Similarity of involutive matrices over a local ring of characteristic $2^k$

An algorithm is constructed for finding the canonical form according to a similarity relation and the transformation matrix for an arbitrary involutory matrix over a finite commutative local ring of characteristic $2^k$.