testTimingOverGL := function(p,a)
	time X:=getInsolubleMatrixGroupsOfDimThreeOverGL(p,a);
	return p, a, #&cat(X);
end function;

S:= [[5,1],[5,4],[5,5],[5,10],[5,12],[17,6],[97,6],[131,10]];

for s in S do
	testTimingOverGL(s[1],s[2]);
end for;