gap> mat := IdentityMat(4,GF(3));
[ [ Z(3)^0, 0*Z(3), 0*Z(3), 0*Z(3) ], [ 0*Z(3), Z(3)^0, 0*Z(3), 0*Z(3) ],
[ 0*Z(3), 0*Z(3), Z(3)^0, 0*Z(3) ], [ 0*Z(3), 0*Z(3), 0*Z(3), Z(3)^0 ] ]
gap> phi := PolarityOfProjectiveSpace(mat,GF(3));
<polarity of PG(3, GF(3)), >
gap> points := AbsolutePoints(phi);
<All elements of Q+(3, 3)>
gap> List(points);
[ <a point in Q+(3, 3)>, <a point in Q+(3, 3)>, <a point in Q+(3, 3)>,
<a point in Q+(3, 3)>, <a point in Q+(3, 3)>, <a point in Q+(3, 3)>,
<a point in Q+(3, 3)>, <a point in Q+(3, 3)>, <a point in Q+(3, 3)>,
<a point in Q+(3, 3)>, <a point in Q+(3, 3)>, <a point in Q+(3, 3)>,
<a point in Q+(3, 3)>, <a point in Q+(3, 3)>, <a point in Q+(3, 3)>,
<a point in Q+(3, 3)> ]