# File: CompRegion100.gap
# Description: GAP script for computing fundamental
# regions for the action of the unit group of a 
# quaternion algebra over the Rationals with discriminants
# up to 100.
# Author: Assaf Wool
# email: assafwool@yahoo.com
# Date: July 2007

Read("QuatAlg.gap");
DiscList:=[];
logfile:=OutputTextFile("Grp100.txt", false);
OutputLogTo(logfile);

for D in [1..100] do
  fac:=FactorsInt(D);
  len:=Length(fac);
  if RemInt(len,2)=0 and len=Length(Set(fac)) then
    Add(DiscList,D);
  fi;
od;
for D in DiscList do
  Print("D= ",D,"\n");
  ok:=createQuatObj(D);
  if ok=1 then 
    Print("algebra:\nB(",-QuatGlob.u,",",QuatGlob.v,") type= ",QuatGlob.type);
    if IsBound(QuatGlob.aux) then Print(" aux= ",QuatGlob.aux);
    fi;
    Print("\ngen= ",QuatGlob.gen, "\nimap= ",QuatGlob.imap,"\nrel= ", QuatGlob.rel,"\ngrel= ",QuatGlob.grel);
    Print("\n----------------------\n");
  fi;
od;
OutputLogTo();
CloseStream(logfile);