reduce_mod(x::ZZMatrix, y::ZZRingElem) versus mod(M::ZZMatrix, p::ZZRingElem)

[fingolfin]

There are two reduce_mod methods with first argument a ZZMatrix. They seem to overlap with several mod and mod! etc. methods (which originally came from Hecke). See also this comment.

What should we do? Redefine one in terms of the other? To me reduce_mod is the oddball, I would have expected it to be called mod. But note that the name reduce_mod is explicitly mentioned in the manual.