grd.c

#include "defs.h"

# define PSP   1000000
# define SPACE 3200000
# define SVSP   50000
# define  NPT   100001
# define  MP     500
# define  MB      80
# define  MREL    1024
# define DNWDS    MREL/32
/* SPACE is space available for coset enumeration.
   MREL is max no of relations allowed.
   We will be using a boolean array done on the relator set, so DNWDS is the
   max no. of 32 bit words we need for the relator set.
*/

char inf[80],outf[80],outfg[80],firstnew,gap;
/* Default inf=gpname.outperm
   outf=inf.rel anyway.
*/
int imsp[SPACE],*done[DNWDS],*imcos[MP];
int    perm[PSP],sv[SVSP],cp[5*NPT],orb[NPT+1],
       base[MB],lorb[MB],pno[MP/2],*pptr[MP],*svptr[MB],
       inv[MP],rno[MB],mp=MP,mrel=MREL,
       dnwds=DNWDS,mpt=NPT,mb=MB;
int    psp=PSP,space=SPACE,svsp=SVSP;

int 
main (int argc, char *argv[])
{ int   arg;  char c,err;
  err=0; firstnew=0; arg=1; gap=0;
  if (argc<=arg) {err=1; goto error;}
  while (argv[arg][0]=='-')
  { c=argv[arg][1];
    if (c=='g') gap=1;
    else if (c=='f') firstnew=1;
    else {err=1; goto error;}
    arg++; if (argc<=arg) {err=1; goto error;}
  }
  strcpy(inf,argv[arg]); strcat(inf,".");
  if (gap) {strcpy(outfg,inf); strcat(outfg,"relg");}
  arg++; if (argc<=arg) strcat(inf,"outperm"); else strcat(inf,argv[arg]);
  strcpy(outf,inf); strcat(outf,".rel");
  if (grprog()== -1) exit(1);
  error:  if (err)
  { fprintf(stderr,"Usage:   grrun [-f] [-g] gpname [inf]\n");
    exit(1);
  }
  exit(0);
}