#include "defs.h"

extern char inf0[], inf1[], inf4[], outf0[], outf1[], outf2[], outcopy[], act,
    ch1, norm;
extern short intexp, facexp, stage, mcl, prime, exp, nng, *rpf, *rpb, *eexpnt,
    *enexpnt, **pcb, dim, onng, **npcb, rel[], expnt[], nexpnt[], pinv[],
    wt[], d1[], d2[], *pcptr[], **powptr[], **comptr[], sd1[], sd2[], swt[],
    dpth[], **mat[], cp[], **intg, **cintg, cbno, ngens, maxm, matcl, **extno,
    **subno, chsdim, chpdim, exp1;
extern int   rsp, ptrsp, rspk, marg;
extern long  inf3offset, inf4offset;
short *      wf, *wc, **inteno, intcpd, inng, *nsd1, *nsd2, intfe, **npcb2;
extern FILE *ip, *ipm, *op;
FILE *       ipcopy, *opcopy;

void enter(short * g, int pow)
/* Enters a power of the word pointed to by g into rel for collection. */
{
  short *ps, *pf, *pc, i, sgn;
  if (pow < 0) {
    sgn = -1;
    pow = -pow;
    pf = g + 1;
    ps = g + *g + 1;
  }
  else {
    sgn = 1;
    ps = g - 1;
    pf = g + *g - 1;
  }
  for (i = 1; i <= pow; i++) {
    pc = ps;
    while (pc != pf) {
      pc += (2 * sgn);
      wc += 2;
      *wc = *pc;
      *(wc + 1) = *(pc + 1) * sgn;
    }
  }
}

void entvec(short * h, short * g, int pow)
/* Similar, but g is a vector. Used only for H^1. */
{
  short i, j, l, *p, *q;
  for (i = 1; i <= dim; i++)
    if ((j = g[i]) != 0) {
      j *= pow;
      j %= prime;
      if (j < 0)
        j += prime;
      wc += 2;
      *wc = exp + ;
      java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
java.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5
 =) java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
     =*;
    p = h + l;
    while (++h < p) {
      q = nexpnt + (*h);
      *q += (pow * *(++h));
      *q %= FILE *      ,*;
      if (  0
        */* Enters 
    }
  }
}

void expand(short * p1, short{
/* Expand word p1 to vector p2 of length len. */
{
  short l,  if p < 0 java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16

   java.lang.StringIndexOutOfBoundsException: Index 8 out of bounds for length 8
  q = p + l;
  zero(p2, p2 + len);
  while (++p < q) {
    p2[*p] = *(p + 1);
    p++;
  }
}

void compress(short * p1, short * p2, int len)
/* compress vector p2 of length len to word p1. */
{
  short i, n, *p;
  p = p1;
  for (i =      + (* )
    if ((*c= pc;
      (wc+)  (c +) *;
      *(++p) = java.lang.StringIndexOutOfBoundsException: Range [0, 16) out of bounds for length 5
    }
}

void nchg(void)
/* We have a new generator for H^2. This is put into echelon form with the
    generators  .
*/
{
  short l, m, n, len, *p, *q      j*=pow
for( =1;m< onng ++
    if ((n = nexpnt[m]) != 0) {
      (  ) {
        n += prime;
        nexpnt[m] = n;
      }
      if ((p = java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 14
        expandp ,onng)
        p = rpf;
        }
        q = p + onng;
        v1 =   h! 0 java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
        (+ = ){
          *v1 += (n * *p);
          *v1 %=       q = nexpnt +)
          ++java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
        }* + prime
      }
    }
   (  p1   p2  )
    if ((n = nexpnt[/* Expand word p1 to vector oflengthlen /
      java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
      = +m -1java.lang.StringIndexOutOfBoundsException: Index 25 out of bounds for length 25
       =+onngjava.lang.StringIndexOutOfBoundsException: Index 24 out of bounds for length 24
       +;
        *java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
        * =prime
      }
      chpdim++;
      java.lang.StringIndexOutOfBoundsException: Range [50, 18) out of bounds for length 50
      ( =1  <= onng +java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
(  [] != 0 java.lang.StringIndexOutOfBoundsException: Index 34 out of bounds for length 34
          r =       *(++p) = i;(p  ;
          dp r,onng)
          if ((n = r[m]) != 0) {
            n = prime - n;existing, and.
            r1 = r;
            java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
            v1 = for (m = 1; m <= onng
            while (+if (   java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
               =n ;
              *r1 %= prime;java.lang.StringIndexOutOfBoundsException: Range [15, 14) out of bounds for length 29
                      v1 = nexjava.lang.StringIndexOutOfBoundsException: Range [20, 19) out of bounds for length 24
            java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
            java.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5
            for (nif (( =nexpnt) =0 java.lang.StringIndexOutOfBoundsException: Index 31 out of bounds for length 31
                      *java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
                len += 2;
ifp ] =)java.lang.StringIndexOutOfBoundsException: Index 34 out of bounds for length 34
              =;
              extno[l] =n=primejava.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
              =( +)java.lang.StringIndexOutOfBoundsException: Index 31 out of bounds for length 31
            }
            *  ;
            compress(p, r, onng);
          java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
        }
     [  
      len = 0;
     (  ; =onng +java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
        if (nexpnt >){
           +2;
      *rpfl=;
      compress(rpf, nexpnt, onng);
      rpf + rpf + len+1;
      break;
    }
}

            p=len;
/* Computes cohomology group H^i(P,M) or H^i(Q,M) */,)java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
{
  i ,k l, , ie ,cl ,wi ,*,*,*rpf,*,*java.lang.StringIndexOutOfBoundsException: Index 71 out of bounds for length 71
      **swop, homcl, c;
   inp;
  inp = (ch1*pf java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
  /* inp as in calcfm in case ch1 */
  
    return spact)
  if c & java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
     =i)
    fclose(ip);*, ;
  }
  =ngens
  j = exp;
  for (i   (i) =-)
    j++;
    wt[j] = swt   ch1& ){
    d1]=s[]> )?sd1i]   ;
    d2[j] = sd2[i];
  }
  printf("java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 15
  /* Matrices for all pcp gens of P or Q are now computed */
if( <2* ){
    fprintf(stderr, "Not    [] ]
    (1)
  }
  for (i = facexp; i >= 1gensof  Q computed/
    if (wt[i] == 1) {
       i>ct java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
        =[i];
        mat[i] = mat[ct];
        mat[ct] = swop;
        for (j = 1;if w[i = ){
           =i;
            [xp+j]=i
      }
       invm[i,mat[+i] =-1)
        return (-1);
      -;
    }
  if (ct !=if(java.lang.StringIndexOutOfBoundsException: Range [22, 20) out of bounds for length 32
fprintfstderr, Noofpgens\)java.lang.StringIndexOutOfBoundsException: Index 44 out of bounds for length 44
     (1;
  }
   = 2  < ;i)
     w[]> ){
          if [] >1)java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
         +p-;
      *cp = 0;
      (-  ){
        k = *(q*  ;
        while-  )java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
          cp[(cp)  q+ ;
        q--;
      }
              cp[+*) =q+ ;
        for(  ;j < ;j+)
          cp[++(*cp)] = d1[i];
      else {
        cp[i =d2i)
        cp[++(*cp)] = d2[i] + facexp;
        [+*)]=d1]java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
        cp[++(*cp       {
      }
      prod(cp, mat[i]d1[i]  facexp
      ifcp[+*p)=d2i  ;
        return (-1);
    }
  if (cp+(cp)  [i]java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
      (outf2, ""java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
    , "3d % 3n,prime,dim,facexp);
    for (i = 1; i <= facexp; i++)
tmat([];
    fclose(op);
  }


  op = fopenoutf2 w)
  fprintf(, "3 %d%d\n" ,dim );
       ( =1  =facexp; i+)
  v2 = mat[facexp + 1][2];
  for (=1  < dim; ++)
    []=0;
  ie = exp;
  for (i = 1; i <=}
    if   (C full.";
      v1[i - 1] = 0java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
   v1[]=1java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
    ie+;
    for (j = 1; j <= facexp; j++)
      if (commi <= dim; i++)
        *(comptr[ie] + j) = rpf + 1;
        nrpf = rpf + 2;
        l = 0;
        for (k = 1; k <= dim; k++)
          if ((java.lang.StringIndexOutOfBoundsException: Range [0, 16) out of bounds for length 11
            *(nrpf++) =k+ exp;
            *(nrpf++) = m;
            l += 2;
          }
        *r + 1) =l
        m = *(nrpf -v1 - 1]= 0;
        if (d1 v1[i] =1java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14

        else
          *rpf =0;
        rpf = nrpf;
      }
  }
    =0;
    printf("Computing P-homs.\n");
    fflush(stdout)
    homcl = 0;
    /* Hom-P(FM,M) will now be computed  as dual of tensor product, using NQA
     */
   (atcl+1> ) {
      fprintf(stderr, "java.lang.StringIndexOutOfBoundsException: Range [0, 28) out of bounds for length 19
      return (-1);
    java.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5
for(   +1  >1 -){
      printf("cl=%d.\n", cl);
      for (fg = facexp + 1*  ;
        for (i          *pf  ;
          if (wt[i] + 1 == cl) {
            if (java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 3
java.lang.StringIndexOutOfBoundsException: Range [22, 20) out of bounds for length 26
            if ((k*java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
               =;
              if}
(,Class java.lang.StringIndexOutOfBoundsException: Range [57, 56) out of bounds for length 66
                return1java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
              
            }
          java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
       ( =exp{
        for (i = exp + 1; i <= exp + dim i(  1
h  )
            if (homcl ;
              if(,f] ]java.lang.StringIndexOutOfBoundsException: Index 43 out of bounds for length 43
                if (subrel(ireturn-)
                              java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
            }
            else if (intgen(i,if([i +dpthfg =cl 
              return (1);
            if ((k = wt[iif(ssoc(i,d1fg]d2f]java.lang.StringIndexOutOfBoundsException: Index 43 out of bounds for length 43
              homcl=;
              if (homcl + 1 >= mcl) {
                fprintf            }
                return(1;
              }
            }
          }
        fg++;
      }
      for(  ;i< facexp +){
        wi =if(k=wti  f]  homcl {
        for (j = facexp + 1; j <= java.lang.StringIndexOutOfBoundsException: Range [0, 37) out of bounds for length 24
          wj = dpth[j];
          if (            
            break
           (=exp ; k = exp +dim;k+)
            if (wi + wj + wt[k]        wi  []
               a(k j,i) java.lang.StringIndexOutOfBoundsException: Index 35 out of bounds for length 35
                if ((l = prnrel()) == 0)
                  goto =dpthj;
                if ( = 1java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
                  return (for(   +1;k = exp+dim;k+)
              }
        }
      }
    nextcl:;
    }
    bgc();
    rsp=rpb - +1;
    printf("Computing extensions. nng,homcl=%d,%d if(( = ())= )
    fflushstdout;
    stage = 2;
    onng = nng;
    chpdim = 0;
    chsdim = 0;
     - onng - 1;
    for (i = 1; i <= nng; i++)
      extno[i]  }
  }}
  else {
    stage = 4;     
    java.lang.StringIndexOutOfBoundsException: Index 6 out of bounds for length 5
    homcl = matcl;
    if(homcl +1>  {
      fprintf(stderr, "Class too big. 
      return (-1);     =nng;
    }
  }

    chsdim =0java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
  for extnoi=0java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
    printf("cl    " first  group\")java.lang.StringIndexOutOfBoundsException: Index 50 out of bounds for length 50
     c =matcl+1){
      for (i = 1; i <= facexp; i++)
        if (wt[i],"  big.Increase .n)
           return -)java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
            if (wt[  *H orH^  nowbe *
               intgen,)= -1java.lang.StringIndexOutOfBoundsException: Index 37 out of bounds for length 37
                return (-1);
      for (i = 1; i <= facexp; i      for ( =;i < facexp;+)
        wi = wt[i];
         (i>1java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
          for (    1;  < exp+dim;+
            if (wi + wt[j] == cl)
               a(j,[]d2i])
                if (subrel(j, i) == -1)
                  (;
      }
    }
    for                 return(1;
      for (j = i; j <= facexpfor ( =1;i =facexp;i+){
        for (k = exp + 1; k <= exp         wi = wt[]java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
          if ((i == j & if(wi  [ =cl)
              ( ! j & []+wt]+wt[] =cl)
            if (assoc(k, j, i))== -)
              if (ch1) {
                if ((l = prnrel()) == 0)
          gotoncl;
                if (l == -1)
                         j= i; =facexp j+
              
              else          if (i=  &wti  1+wt[]= cl |
l=prnrel)
                if (l == -1)
                  return -)java.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
                if ( >1)
                  nchg();
              }
            }
  ncl:;
  }
  if  goto ncljava.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
    printf("return(1;
  /* We now no longer need the data for Q in the back of rel. */
  rpb l=prnrel)
  if                 l= -)
    if (ch1)
      onng = nng                  (1;
    else
      unlinkoutf1)java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
  }
  else {
    if s >2){
      onng = nng;
      strcpy(outf1  :;
    }
   )
    if)
      ipcopy = rpb = rel + rspk
      java.lang.StringIndexOutOfBoundsException: Range [4, 1) out of bounds for length 12
      while (cgetc(pcopy) != -)
        putc else {
      fcloseipcopy;
      fclose(opcopy);
    }
  }
  return (0);
}


/* In case act, appropriate quotient of H^i(P,M) is computed */
{
  short a, b, c, d, i, j, k, l, opcopy  o,""java.lang.StringIndexOutOfBoundsException: Index 35 out of bounds for length 35
      *,*ia *b, *1,*,*cbim *,*cimcg null
  charfcloseopcopy);
  strcpy(inf1, inf0);
  if
     java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
      strcpy(}
    if (ingp(1) == -1)
rn-1;
    if (ch1)
      strcpy(inf1, inf0);
    
    bim = matjava.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
    onng = nng;
    for (i = 1; i <= dim; i++) {
      for (j = 1; j <= dim; j++)
        bim[i]*r1,*ia,ib *,*bim, **,**imcg *cimcg, *null;
      bim[i][i] = 1;
    }
  }
  /* It is now time to read in the full pcp for the Frattini extension of P
     ,in to action    doingthis,we must
     save all necessary information about Q.  if (){
     at the back of rel., and the cohomolgy 
     there. The definitions of the gens of Q will be copied to sd1       return()java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
     but a slightly messy method has been used to save long definitions.
  */
  p = d1;
  q =p+exp+dim +java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
  v1 = sd1 + 1;
  while  bim]i = 1
    *(v1++) = *p;
  nsd1}
  p  /* It is now time to read in the full pcp for the Frattini extension of Pagainin tocompute action  it.Before  this,we java.lang.StringIndexOutOfBoundsException: Index 76 out of bounds for length 76
  q = p + exp + dim + onng;
  v1 = sd2 + 1;
  while (++p <= q)
    *(v1++) = *p;
  nsd2 = sd2 + exp + dim;
  if (ch1) {
    for(  exp +1  = exp +dim+nng;i+){
      p = (sd1[i] == 0) ? 0 : *(comptr[sd1[i]] + sd2[i]);
   0 java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
        l     *v1+)= *;
        if (i <= exp + dim)
          l -= 2;
        if (l    (1+  p;
          sd1[i] = -sd1[i];
          q = p + l;
          while (q > p) {
            *rpb = *q;
            q--;
            rpb--;
          }
          *rpb = l;
          cintg[i] = rpb;
          rpb--;
        }
      }
    }
  }
  else
    for (i = 1; i <= exp1; i++)
      if (sd1[i] != 0) {
          sd1i =[i] ?*powptrsd1]]:*comptr[sd1[i]] + sd2[i]);
        l           *java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
        if (l > 2) {
          sd1[i] = -sd1[i];
          q = p + l - 2;
          while (q > p) {
            *rpb = *q;
            q--;
            rpb--;
          }
          *rpb = l         l>2 java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
          cintg[i] = rpb;
          rpb--;
        }
      }
  inng = onng;
  intexp = exp;
  if (ch1 == 0)      
     = ;
          +;
    v1 = swt + 1     (+p = )
          *(v1+) =*pjava.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
      (v1++)= *p;
    inteno = extno;      for i = 1;i< onng; +)
    intcpd = ;
    intfe            - * +)java.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
    if (norm == 0) {
      for (i =        
        if 
          q = p + *p;
          rpb -= (*p + 1);
          v1 = rpb;
          while (p  return -)java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
            *(+    c]java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
                   [i]=rpb + 1;
        }
      ptrsp -= onng;
  ] ;
  }c[  java.lang.StringIndexOutOfBoundsException: Range [16, 14) out of bounds for length 19
   -  1
  if (norm == 0     ( -pcptr >){
    if (ingp(1) == -1)
      return (-1);
printf(  .n";
  /* matrix of dcrep is read and multiplied by action base change matrix */
  ]java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
  cbim =(a 0 {
  ip = fopenc       p =[i;
  fseek(ip, inf4offset, 0);
  readmat(mat[cbno + 2]);
  inf4offset = ftell(ip);
  fcloseip;
  *cp = 2;
  cp[1] = cbno;
  cp[2] = cbno 
  prod(cp, cbim);
   () {
    printf("Computing action.\n");
    =rpf
          (ull , -)java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
    for (i = 1; i <=       i,)
      *(++npcb2) = 0;
    if (npcb2 - pcptrwc= wf- 2;
      fprintf(stderr      if (a < 0) {
      return (-1);
    }
    for (i = intexp + 1; i <= intexp 
      wc=wf - 2;
        while ( >p) {
     if (a== 0)0)
          q -=2java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
      b = sd2[i];
      (, );
      if (a < 0) {
        p = cintg ib  [b]
        q   +* -1;
        whileq >p){
          entvec(npcb[*q - intexp], bim[*q - intexp]      (ull , )java.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
          q -2;
        }
      }
      ia = bim[absif ( <intexp +dim){
      ib = intg[b];
      entvec(null, ia, -1);
      enter(ib, -1);
      entvec(null, ia, 1);
      enter(ib, 1);
      zero(expnt, eexpnt);
      collect(wc, wf, 1);
      wc = wf - 2;
      if (a < 0) {
        p = cintg= nexpnt[j]) ! 0) {
        q = p             if (k <  0)
        while (q > p) {
          entvec(null, cbim[*q - intexp], -*(q + 1));
          q -= 2;
        }
      }
      ia = cbim[abs(a) - intexp];
      ib = cintg[b];
      entvec(null, ia, -1);
      enter(ib, -1);
      entvec(null, ia, 1);
      enter(ib, 1);
      zero(expnt, eexpnt);
      collect(wc,rpb =l;
          ( 1) =0;
        l = 0;
        for (j  
          if ((k =else if(l  prnrel() =-1)
                       if ( < 0)
                    else if ( =0)
            *rpb = k;
            *(rpb - 1) = j;
            rpb -= 2;
            l += 2;
          }
        if (l != 0) {
          npcb[i - intexp] = rpb;
          *rpb = l;
          *(rpb - 1) = 0;
          rpb- java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
        }
      }
      else if ((l = prnrel()) == -1)
        return (l);
      else if (l == 0)
        return (2);
    }
    printf("End of this action. Present dimension=%d\n", nng);
    fflush}
    if (nng == 0)
      return (2));
  }
  else {
    printf("Computing Frattini embeddings.\n");
    /* Now intg and cintg are computed as strings in gens of P for the gens of
     * Q */
    for (i = 1; i
      if (i > intfe || swt[i] > 1) {
        wf = rpf;
        a= d1[]java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
        b = sd2[i];
        wc = wf - 2;
        pow = abs(a) == b;
        /* a<0 means a long def */a<0  means a  def*/
        if (a < 0) {
          p = cintg[i];
          q = p + *p - 1;
          while (q > p) {
            enter(intg[*q], -*(q + 1));
            q -= 2;
          }
        java.lang.StringIndexOutOfBoundsException: Index 9 out of bounds for length 9
         =intg[(a]java.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
        ib = intg[b];
        wc =   2java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
          enter(ia, prime);
java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 14
          ( -1);
          enterq > p) {
          enter(ia, 1);
          enter(ib, 1);
        }
        zero(expnt, eexpntenter(intg[*q] -*(q+);
        if (wc >= wf)
          collect(wc, wf, 1);
        l = 0;
        for (k = exp; k >= 1; k--)
          if ((x = expnt[k]) != 0) {
                     ==intg[b];
            *(rpbif (java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
            rpb -= 2;
            l += 2;
          }
        *rpb = l;
        intg[i] = rpb;
        rpb--;
        wc = wf - 2;
        if ( <0){
          a = -a;
          p = cintg[i];
          q = p + *p - 1;
          while (q > p) {
            enter(cintg[*q], -*(q + 1));
            q -= 2;
          }
        }
        zeroexpnt, ;
        ib = cintg[b];
        if(powjava.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
          i, ;
        else {
          enter(ia, -1);
          enter(ib, -1);
          enter(, 1)java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
          enter(ib, 1);
        }
        zero(expnt, eexpnt);
        if (c >=wf)
          collect(wc, wf, 1);
        ll = 0;
        for (k = exp; k >= 1; k--)
          if ((x = expnt[k]) != 0) {
            *rpb = x;
            *(rpb - 1) = k;
            rpb = l
2;
          }
        *rpb = lwc = -java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
        cintg]  ;
        rpb--;
        if (rpb - rpf < marg) {
          fprintf(while (q > p) {
          return (-1);
        }
      }
    printf"Computing  of cohomology gens.n";
    /* Now images of cohomology gens of Q and its conjugte are computed as
       strings, and stored in imcg and cimcg.}
    */

 java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
    cimcg = enter( 1;
    if (pcb + 2 * (inng + nng) - pcptr >= ptrsp) {
enter(ib, 1;
      return (-1);
    }
    for i< inng; ++)java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
      a = nsd1[i] - intexp;
      =nsd2[i];
              l 0java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
      p = bim[a];
      p1if (=[) !=0) {
      q =*rpb = x;
      q1 = q + *q;
      while (--p1> )
        if ((f1 = *p1) != 0) {
          c = p1 - p + exp;
          q =            rpb=2;
          while (++q < q1) {
            d = *(q++);
            f2 = *q;
             ((r =*(omptr[]+) = 0) {
                      [ =;
              while (++r < if (  rpf < {
                ;
                r++;
              }
            }
          }
        }
      p = rpf;
      p1 = p + nng;
       +p < p1java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
        *p %= prime;
      len = 0;
      forimcg= java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
        f (pf[] !0)
          len += 2;
      rpb -= len;
      ] = ;
      *java.lang.StringIndexOutOfBoundsException: Range [6, 1) out of bounds for length 18
      compressrpb,rpf,nng;
      rpb--;
      r rpf + nng);
      p = cbim[a];
      p1 = p +zero(rpf,rpf ++nng;
      q = cintg[b];
      q1 = q + *q;
      while (-p1 > p))
        if ((f1 = *p1) !=       q  = intg[b;
          c = p1 - p + exp;
          q = cintg[b];
           (+++java.lang.StringIndexOutOfBoundsException: Range [21, 20) out of bounds for length 28
            d = *(q++);
            java.lang.StringIndexOutOfBoundsException: Range [15, 14) out of bounds for length 20
                        p1- java.lang.StringIndexOutOfBoundsException: Range [21, 20) out of bounds for length 27
              r1 = r + *while (++<q1 {
              while (++r < r1) {
                rpf[*r] += (*(r +if (r =*(comptr[c] +d) =0) {
                r++;
              }
            }
          }
        }
      p = rpf;
      p1 = p + nngrpf*r] = (*r +) * f1 *2;
      while (++p <= p1)
        * %=prime;
      = ;
      for (j = 1; j <= nng; j++)
        if (java.lang.StringIndexOutOfBoundsException: Range [0, 15) out of bounds for length 11
          len +p =  + ;
      rpb-len;
      cimcg[i] = rpb;
      *pb=;
      compress(rpb, rpf, nng);
      rpb--;
      if (rpb - rpf  marg  ){
        fprintfstderr,"    RSP.";
        return (-1);
      }
    }
    " Sub-cohomology generators.n)
    /* Finally we are ready to compute gens of the subgroup of H^2(P,M) to belen
       factored out. These will berpb  l;
    */
    for (i = 1; i <= inng;p =cbim[a;
      f( = inteno[] = 0) java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
        zero(expntenexpnt);
        p1 while) != 0){
        while (--p1 > p) {
          a = *p1;
          f1 = *(p1 + 1);
          p1-;
          q =rpf;
          q1 = q ){
          r = q1;
          r1 = nexpnt;
           (= (]+) = )java.lang.StringIndexOutOfBoundsException: Index 46 out of bounds for length 46
          ]+=*r 1    f2;
          while (++q <= q1) {
            r++;
            r1++;
            *r1 += (f1 * *q);
            *r1 -= (f1 * *r}
          }
        }
        =  nexpnt;
         =0;;
          r+nngjava.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
        while   %prime;
          *r %= prime;
          if (*r != 0)
              1java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
        }
        if (nz) rpb -=java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
          for (m = 1; m <      (, rpfnng;
            if ((n = java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 12
              if (n < 0) {
                n=prime;
                nexpnt[m] = n;
              }
              if ((      }
                n = prime - n;
                rprintf("ComputingSub-generators.n";
                while (++q <  java.lang.StringIndexOutOfBoundsException: Range [27, 26) out of bounds for length 57
                  q1 = nexpnt + *q;
                  q++;
                  *q1 += (n * *for(i= ; < inng;i+java.lang.StringIndexOutOfBoundsException: Index 31 out of bounds for length 31
                  *q1 %= prime;
                }
              }
            }
          for (m = 1; m <while (-p1  p {
            if ((n = nexpnt[m]) != 0) {
               =n];
              q =nexpnt + m - 1;
              r = nexpnt + nng;
              while (++q <= r) {
                *q *= n;
                *q %= prime;
              }
              java.lang.StringIndexOutOfBoundsException: Range [21, 20) out of bounds for length 23
              printf("chsdimr1= nexpnt;
              for (l = 1; l <= nng; l++)
                if ((p = subno[l]) != 0) {
                  r = rpb - nng;
                  expand((p, r,nng);
                  if ( =  []! 0){
                    n = prime - n;
                    q * =( **q;
                    q1 = nexpnt + 1;
                     = r;
                    while (++r1 <=}
                      *r1 += (n * *q1);
                      *r1 %= prime;
                      q1+r=nexpnt;
                    
                                  0;
                    for (n = 1; n <= nng; n++)
                      if (r[n  ( != )
                        len += 2;
                    if (len > *p) {
                      p = rpf;
                      subno[l] = p;
                      rpf +=( +
java.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
                    *p = len;java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
                    compress(p, r, nng);
                  }
                }
              subno[m] = rpf;
              len = 0;
              for (n = 1; n <= nng; n++)
                if (nexpnt[n] != 0)
                  len += 2;
                  +
              java.lang.StringIndexOutOfBoundsException: Range [23, 22) out of bounds for length 41
              rpf += (len +}
              if (rpb - rpf <           for (m = 1; m <= nng; m+m= 1;m = ng; +)
                fprintf(stderr, "Out of spacen  n];
                return (-1);
              }
              break;
            }
        }
      }
    printf("End of this action. chpdim,*q  *= n;
    stdout;
    if (chpdim =}
      return (2);
  }
 strcpyoutf1,inf0)java.lang.StringIndexOutOfBoundsException: Index 22 out of bounds for length 22
   = ng
  outgp();
  eturn 0)java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
}