I am currently planning to write a python package for matrix operations on GF(q) (finite fields).
The best implementation is part of sage software which depends on so many other packages and can't be used independently.
(alpha version released for more information see take a look at the package package homepage)