If hessenberg is called on a matrix over ZZ it calls inv(::BigInt) (or inv(::fmpz) in Nemo. However it is defined incorrectly by Julia for BigInt and not defined at all for fmpz in Nemo.
It's not clear we actually want to define inv for ZZ or not. We've been using div to separate fields from other kinds of rings.
One possibility would be to use divexact instead in hessenberg.
If hessenberg is called on a matrix over ZZ it calls inv(::BigInt) (or inv(::fmpz) in Nemo. However it is defined incorrectly by Julia for BigInt and not defined at all for fmpz in Nemo.
It's not clear we actually want to define inv for ZZ or not. We've been using div to separate fields from other kinds of rings.
One possibility would be to use divexact instead in hessenberg.