gap&gt; h:=HermitianPolarSpace(2, 7^2);
H(2, 7^2)
gap&gt; pg := AmbientSpace( h );
ProjectiveSpace(2, 49)
gap&gt; lines := Lines( pg );
&lt;lines of ProjectiveSpace(2, 49)&gt;
gap&gt; curve := AsList( Points( h ) );;
gap&gt; Size(curve);
344
gap&gt; Collected( List(lines, t -&gt; Number(curve, c-&gt; c in t)));
[ [ 1, 344 ], [ 8, 2107 ] ]
gap&gt; time;
26412
