FrobeniusNormalForm() does not work with matrix objects as input

[AliaBonnet]

The function FrobeniusNormalForm() returns an error when the input is a matrix object.

gap> A := RandomInvertibleMat(5,GF(5^5));;
gap> FrobeniusNormalForm(A);
[ [ x_1^5+Z(5^5)^2011*x_1^4+Z(5^5)^124*x_1^3+Z(5^5)^393*x_1^2+Z(5^5)^1076\
*x_1+Z(5^5)^2156 ], 
  [ [ Z(5)^0, 0*Z(5), 0*Z(5), 0*Z(5), 0*Z(5) ], 
      [ Z(5^5)^2020, Z(5^5)^1683, Z(5^5)^2703, Z(5^5)^2316, Z(5^5)^1432 ]
        , [ Z(5^5)^633, Z(5^5)^2345, Z(5^5)^109, Z(5^5)^278, Z(5^5)^803 ]
        , [ Z(5^5)^678, Z(5^5)^434, Z(5^5)^1838, Z(5^5)^3009, Z(5^5)^957 
         ], 
      [ Z(5^5)^1637, Z(5^5)^2308, Z(5^5)^1759, Z(5^5)^2789, Z(5^5)^1735 
         ] ], [ 1 ] ]
gap> A := Matrix(GF(5^5),A);
<5x5-matrix over GF(5^5)>
gap> Print(A);
NewMatrix(IsPlistMatrixRep,GF(3125),5,
[ [ Z(5^5)^2020, Z(5^5)^1683, Z(5^5)^2703, Z(5^5)^2316, Z(5^5)^1432 ], 
  [ Z(5^5)^1354, Z(5^5)^118, Z(5^5)^2779, Z(5^5)^2039, Z(5^5)^2628 ], 
  [ Z(5^5)^2878, Z(5^5)^1444, Z(5^5)^809, Z(5^5)^555, Z(5^5)^59 ], 
  [ Z(5^5)^2461, Z(5^5)^1107, Z(5^5)^599, Z(5^5)^683, Z(5^5)^716 ], 
  [ Z(5^5)^528, Z(5^5)^91, Z(5^5)^64, Z(5^5)^2168, Z(5^5)^2972 ] ])
gap> FrobeniusNormalForm(A);
Error, ConvertToVectorRepNC: Vector cannot be written over GF(5) at /home/alia/Downloads/gap-4.15.1/lib/vecmat.gi:1390 called from
ConvertToVectorRep( r, Q_VEC8BIT( v ) 
 ); at /home/alia/Downloads/gap-4.15.1/lib/vec8bit.gi:975 called from
Vector( Unpack( v ) * Unpack( M ), v 
 ) at /home/alia/Downloads/gap-4.15.1/lib/matobj.gi:1195 called from
PROD_LIST_SCL_DEFAULT( list, nonlist 
 ) at /home/alia/Downloads/gap-4.15.1/lib/list.gi:3447 called from
sp[2] * A at /home/alia/Downloads/gap-4.15.1/pkg/nofoma/gap/nofoma.gi:436 called from
MaximalVectorMat( A 
 ) at /home/alia/Downloads/gap-4.15.1/pkg/nofoma/gap/nofoma.gi:635 called from
...  at *stdin*:14
you can 'quit;' to quit to outer loop, or
you can 'return;' to continue
brk> 

[AliaBonnet]

I tried fixing it by converting the input matrix A into a list at the beginning but this is very inefficient and slows down the function significantly. The majority of functions used in FrobeniusNormalform() probably need to be modified to accomodate usage of matrix objects.