Skip to content

Commit

Permalink
Fix SMTX.InvariantSesquilinearForm for reducible case
Browse files Browse the repository at this point in the history
... by deleting unnecessary 'normalization' (???) code.
  • Loading branch information
fingolfin committed Nov 8, 2024
1 parent c1513db commit 0829bf4
Showing 1 changed file with 13 additions and 13 deletions.
26 changes: 13 additions & 13 deletions lib/meataxe.gi
Original file line number Diff line number Diff line change
Expand Up @@ -3384,19 +3384,19 @@ SMTX.InvariantSesquilinearForm:=function( module )
return fail;
fi;
# Replace iso by a scalar multiple to get iso twisted symmetric
q:=Size(module.field);
r:=RootInt(q,2);
isot:=List( TransposedMat(iso), x -> List(x, y->y^r) );
isot:=iso * isot^-1;
if not IsDiagonalMat(isot) then
Error("Form does not seem to be of the right kind (non-diagonal)!");
fi;
l:=LogFFE(isot[1,1],Z(q));
if l mod (r-1) <> 0 then
Error("Form does not seem to be of the right kind (not (q-1)st root)!");
fi;
iso:=Z(q)^(l/(r-1)) * iso;
iso:=ImmutableMatrix(GF(q), iso);
# q:=Size(module.field);
# r:=RootInt(q,2);
# isot:=List( TransposedMat(iso), x -> List(x, y->y^r) );
# isot:=iso * isot^-1;
# if not IsDiagonalMat(isot) then
# Error("Form does not seem to be of the right kind (non-diagonal)!");
# fi;
# l:=LogFFE(isot[1,1],Z(q));
# if l mod (r-1) <> 0 then
# Error("Form does not seem to be of the right kind (not (q-1)st root)!");
# fi;
# iso:=Z(q)^(l/(r-1)) * iso;
iso:=ImmutableMatrix(module.field, iso);
SMTX.SetInvariantSesquilinearForm(module, iso);
return iso;
end;
Expand Down

0 comments on commit 0829bf4

Please sign in to comment.