Testing the output

> op(L);

> evalm(L&*C&*transpose(L)-F);

Which is clearly the zero matrix on Z/3Z