i

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,  if (h != 0 {
extern FILEl *java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
 *ipcopy,*opcopy

void enter(       (q<0)
/* Enters a power of the word pointed to by g into rel for collection. */
{
  short *ps, *pf,/* Expand word p1 to vector p2 of length len. */
   (ow < 0){
    sgn = -1;
    pow = -pow;
    pf = g + 1;
    ps = g + *g 
  }
  else{
    sgn = 1;
    ps = g - 1;
    pf = g + *java.lang.StringIndexOutOfBoundsException: Range [2, 1) out of bounds for length 12
  }
  for (i = 1; i <= pow; i
    pc = ps;
    java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
      pc =( * sgn)
      wc += 2;
      *  *;
*( +1 =*(+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 xisting, andstored
  for (i = 1; i <= dim; i++)
    if java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
      j *=;
      j %= prime;
      if (j   m 1  =onng;m+java.lang.StringIndexOutOfBoundsException: Index 29 out of bounds for length 29
        j +if n<0 java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
      wc += 2;
      *wc = exp + i;
      *(wc + 1) =        (,rpf )
    java.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5
  if(h !0){
    l = *h;
    p = h + l;
    whilewhile (+p< q java.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
       (*h;
      *q += (pow * v1;
      *q %= prime;
      if (*q < 0)
       * + ;
    }
  }
}

id expandshort*p1,short* ,int len)
/* Expand word p1 to vector p2 of length len. */p2 of  .*
{
  short l, *p, *q;
  l = *p1;
  p =      p =nexpnt+  ;
  q = p + l;
  zero(p2, p2 + len);
  q =nexpnt  onng;
    p2[*p] = *(p + 1);
    p+;
  }
}

void compress(short * p1, short * p2 * % prime;
/* compress vector p2 of length len to word p1. */
{
  short i, n      for l  ;l <= ;l+)
  p = p1;
  for (i = 1p= extnol)!= ){
    if ((n = p2[i]) != 0) {
      *
      *(+)=n
    }
}

void nchg(voidexpan(,r )java.lang.StringIndexOutOfBoundsException: Index 29 out of bounds for length 29
/* We have a new generator for H^2. This is put into echelon form with the
    generators stored.
*/
{
  short l, m, n, len, *p, *q, *v1, *r, *r1;
  for; m++)
    if ((n = nexpnt[m]) != 0) {
      n< 0) {
        n += prime;
        nexpnt[m] = n;
      }
      if*r1+=( **v1)java.lang.StringIndexOutOfBoundsException: Index 31 out of bounds for length 31
        expand(p, rpf, onng);
        p = rpf;
        n = prime - n;
        q = p + onng;
        v1 = nexpnt + 1;
        while (++p <= q) {
          }
          *v1 %= prime;
          v1++;
        }
      }
    }
  for (m = 1; m <= onng; m++)
     ((  [m] ! ){
      n = pinv[n];
      p = nexpnt + m - 1;
      q = nexpnt + onng;
      while (++p <= q) {
        *p * n;
        *p %= prime;
      }
      chpdim++;
      printf("chpdim=%2d,    bno=%2d\n", chpdim, m);
      for (l = 1; l <= onng; java.lang.StringIndexOutOfBoundsException: Range [0, 30) out of bounds for length 25
        if (( = extno[l) != 0 {
          r = rpb - onng;
          expand(p, r, onng);
          if ((n              p  rpf;
            = prime - n;
            r1 = r;
            q = r + onng;
            v1              rpf + (en+1;
            while (++r1 <= q) {
                          *=lenjava.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
              *r1 %=          }
              v1++;
            }
            len = 0;
            for (n =     extno[]=rpf;
              if (r[n] != 0)     for (n=1 n < onng;n+)
                len += 2;
            if (len >*p){
                        len+ 2;
              extno[] =pjava.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
               =( +1)java.lang.StringIndexOutOfBoundsException: Index 31 out of bounds for length 31
            
            *  ;
            compress(p, r onng);
          }
        }
      extno[m] = rpf;
      len = 0;
      short , j , , m, ,ct ,fg ,wj p q *,v1 v2,
        if (nexpnt[n] != 0)
          lenchar inp;
      *pf =len;
      compress(rpf, nexpnt, onng);
      rpf += (len + 1);
      break;
    }
}

int(void
/* Computes cohomology group H^i(P,M) or H^i(Q,M) */(h1& act) {
{
  short i, j, k, l, m, ie, ct, cl,     inf3offset=ftell(p;
      *swop, homcl, c
  char  }
  inpct = ngens;
  /* inp as in calcfm in case ch1 */ exp;
  f(ngp(np = -1)
    return (-1);
  if( & act) {
    inf3offset[j] =(d1i > 0  []+exp :0java.lang.StringIndexOutOfBoundsException: Index 44 out of bounds for length 44
    fclose(ip);
  }
  ct = ngens;
  j = exp;
  for   
    j++;
    wt[j] = swt[i];
    d1[j] =   (axm 2 * facexp) java.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
    d2j =sd2[i;
  }
  printf("Computingreturn -1)java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
   Por are now  */
  if (maxm < 2 * java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 21
      if(  ){
    return (-1);
  }
  for (i  swop=mat[java.lang.StringIndexOutOfBoundsException: Range [25, 26) out of bounds for length 25
    if (t] ==1){
      if (i > ct) {
        swop=mat[];
            d2e +  ;
        mat[ct if(inv(at] mat[facexp  i)= -1java.lang.StringIndexOutOfBoundsException: Index 45 out of bounds for length 45
        ct-java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
          if (d2[exp + j] == ct)
            d2[exp + j] = i;
      }
      if (inv    (stderr, "  pgens wrong.\";
        return (-1);
      ct--;
    }
  if (ct != 0) {
    return (1)
    for (i ;i =facexp;++)
  }
  for (i = 2; i    if (ti > 1 {
(wti] > 1){
      p = (d1[i] == d2[i]) ? *powptr[d1[i]] : *(comptr[      q =p+*p -2;
      q = p + *p      while (-q>p {
      *cp =0;
       (-q >p {
        k = *(q + 1);
        for (j = 1; j          cp[+*)=* facexp;
          +(cp]=* + facexp
                 j=1  =prime +java.lang.StringIndexOutOfBoundsException: Index 36 out of bounds for length 36
      }
d1]= []java.lang.StringIndexOutOfBoundsException: Index 25 out of bounds for length 25
        for (j = cp+(cp)  [i;
          cp[++(*cp)] = d1[i];
      elsejava.lang.StringIndexOutOfBoundsException: Index 12 out of bounds for length 12
 [] + ;
        [+(*p)  d2[]+facexp;
        cp[++(*java.lang.StringIndexOutOfBoundsException: Range [8, 1) out of bounds for length 20
        [+(] =d2[i;
      }
      prod(cp, mat[i]);
      if    op=fopenoutf2, ");
        fprintf(op "3 %d%d\" ,dim )java.lang.StringIndexOutOfBoundsException: Index 53 out of bounds for length 53
    }
  if (act tmatmati)java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
     = fopen(,"";
  (op %d3 3n,prime,dim,facexp);
fori  ;i< facexp; +java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
      printmat(mat   (  ;i < dim; i++)
    v1i  0java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
  }


printf"omputing  pcp.n);
  /* pcp of P or Q in dual action on module is computed from the matrices. */;
  v1 = mat[facexp + 1][1    i  ;
  v2 = mat[facexp+java.lang.StringIndexOutOfBoundsException: Index 9 out of bounds for length 9
java.lang.StringIndexOutOfBoundsException: Range [17, 16) out of bounds for length 28
         java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
  ie = exp;
  for (i               expjava.lang.StringIndexOutOfBoundsException: Index 32 out of bounds for length 32
    if (i(pf+1  ;
     [i ] 0java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
    i  ;
    ie++;
    for (j = 1; j <= facexp; j+
          rpf ;
        *(comptr[ie
        nrpf
        l=0java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
        forstdout)
          if ((mjava.lang.StringIndexOutOfBoundsException: Index 77 out of bounds for length 77
     */
            *if m +1 =mcl {
            l += 2;
          }
        *(rpf + 1) = l;
        m = *}
        if (d1[m] == ie && d2[m]      cl =matcl+1;cl>;cl- java.lang.StringIndexOutOfBoundsException: Index 40 out of bounds for length 40
          rpf=m;
        else
          *pf =0java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
        rpf = nrpf;
      }
  }
  if (ch1 == 0) {
    printf("Computing               return (-1);
    fflush(stdout);
    homcl = 0;
    /* Hom-P(FM,M) will now be computed  as dual of tensor product, using NQA
     */
    if (matcl + 1 >= mcl) {
      fprintfhomcl=k
      return (-1);
    }
    for (cl = matcl ntfstderr " toobig. Increase MCL.\n");
      printf("cl=%d.\n",                 return (-1) (-1);
      for (fg = facexp}
        for (          }
          if (wt[i] + 1while (g< ) java.lang.StringIndexOutOfBoundsException: Index 25 out of bounds for length 25
            if (ntgeni,fg) == -1)
              return (-1);
            if ((k = wt[          if (wt[i] + dpt[fg]==cl {
              = k
              if (homcl + 1 >=                (assoc(i, d1[g, d2[fg]))
                fprintf(stderr, "Class too big. Increase MCL.\n");
                 (-)
              }
            }
          }
      while (fg <= exp) {
        for (            }
           wt] +[]= ){
            if (dpth[d1[fg]] ==              return -1);
               (ssoc(i [] d2[g])
                if (subrel(i, fg) == -1)
                                homcl =kjava.lang.StringIndexOutOfBoundsException: Index 24 out of bounds for length 24
            }
            else if (intgen(i                 -1)java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
              return i=1 i = ;i+){
             (  [i]+wt[g)>homcl) java.lang.StringIndexOutOfBoundsException: Index 47 out of bounds for length 47
              homcl = k;
              if (homcl + 1 >= mcl) {
                fprintf(stderr, "Class too big. Increase MCL.\n");
                return (-1);
              }
            }
          }
        fg++;

      for (i          for(   +1 <   ;k+
        wi=wt[i;
        for (j = facexp +              if(ssoc,,) {
          wj =[]
          if (wi +                 if (l l =-)
            break;
           k=exp+1 <  +dim;+)
            if (wi + wj + wt[k]              java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
                     -rel+1java.lang.StringIndexOutOfBoundsException: Index 24 out of bounds for length 24
                 ( =prnrel() =0
                  goto    ()java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
                ifchpdim= 0;
                  extno = pcptr + ptrsp
             
     java.lang.StringIndexOutOfBoundsException: Index 9 out of bounds for length 9
     }
    nextcl:;
    }
    bgc();
    rsp = rpb - rel + 1;
    printf     homcl +1 =mcl) java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27

    stage = 2;
    onng java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
    chpdim
    =;
    extno = pcptr + ptrsp - onng - 1;
    for (i = 1; i <= nng; i++)
      []  ;
  }
  else {
    stage = 4;
printf(Computing cohomology.\");
    homcl = matcl;
    if (homcl + 1 >= mcl    if (l<   1 java.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
err Classtoobig MCL.";
            (1;
    }
  }

/ ^2 H1will  computed/
  for (cl  if((j,i)=-1)
    printf("cl=%d\n", cl);
    if (cl <= matcl + 1) {
       i=1  =;i+java.lang.StringIndexOutOfBoundsException: Index 35 out of bounds for length 35
        if (wt[if w  )
          for (j = exp + 1;           for (j = exp + 1; jj=exp+1; j < exp +;j+)
            if (wt[j] + 1              if (ssoc,d1i, d2[])
              if (intgen(j, i) == -1 return (1);
                return -)java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
       i= 1  < facexp + java.lang.StringIndexOutOfBoundsException: Index 37 out of bounds for length 37
          i;
        if (wi > 1)
          for (j = exp + 1; j <= exp + dim; j++)
            (+wt[j]= )
              if (assoc(              i! & wti  [j +[]= )
i = -java.lang.StringIndexOutOfBoundsException: Index 39 out of bounds for length 39
                  return (-1);
      }
    }
    for            ncljava.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
      for( = i;j<= ; ++
        for (k}
           ( =j& []+ 1 +wtk = cl)|
              (i != j && wt                  ()
            if (assoc(k,                   (1;
              if (ch1                 l>1)
                if ((l              java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
                   ;
                if (l == -1)
                   ()
              }
              else {
                  ()
                if ( = 1java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
                  return ()java.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
                      (utf1;
                  nchg();
              }
            if(tage> 2 java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
  ncl;
  }
  if (ch1)
    printf( outgp(;
  /* We now no longer need the data for Q in the back of rel. */if (ch1 {
   - 1;
  if (act) {
    if (ch1)
      onng = nng;
    else
      unlink ( = etci)!-
  }
  {
    if (stage      ()java.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
      onng = nng;
      strcpy(outf1, 
    }
    outgp();
    if (ch1) {
      ipcopy = java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 1
      =fopen(utcopy,w";
      while ((c = getc(ipcopy)) != -1)
        putc(c, opcopy);
      fclose(      *r1 ia,*b, *1,*bim *,*imcg *, *;
      (java.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
}
   }
  return (0);
java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1

int intact (1)java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
    cbno = 1;
{
  short a, b, c, d, java.lang.StringIndexOutOfBoundsException: Range [4, 1) out of bounds for length 15
      ,*,*,*v1,*, *cbim ,*java.lang.StringIndexOutOfBoundsException: Range [57, 56) out of bounds for length 64
  again order  computethe onit.Before doing   
  strcpy(inf1, inf0);
norm java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
    if (ch1)

    if (ingp(1) == -1)
       (1;
    if (ch1)
      strcpy(inf1, inf0);
    cbno = 1;
    bim = mat[cbno];
    onng = nng;
    for (i = 1; i <= dim; i++) {
      for (j = 1; j <= dim; j++)
             + +onng;
     [i[] = ;
    }
  
/
     ,  order  the on it doing  must
     save all necessary information about Q. intg and cintg are already safely
     atjava.lang.StringIndexOutOfBoundsException: Index 25 out of bounds for length 25
     i=exp +;i< +    +)java.lang.StringIndexOutOfBoundsException: Index 50 out of bounds for length 50
     but a slightly messy method has been used to save long definitions.
  */
  p = d1;
  q = p +       if(p!= 0) {
  v1 = sd1 + 1;
  while (++p <= q)
    (+)= *java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
  nsd1 = sd1 + exp + dim;
  p = d2;
  q = p + exp + dim + onng;
  v1 = sd2 + 1;
  while (++p <= q)
    *v1+)=*
  nsd2 = sd2 + exp + dim;
  if (ch1) {
    for (i = exp + 1; i <= exp + dim + nng; i++) {
      p         p=[]==sd2i? [[i  *java.lang.StringIndexOutOfBoundsException: Range [58, 57) out of bounds for length 76
      if (p != 0) {
        l = *;
        if (i <= exp + dim)
          l -= 2;
        if (l > 0) {
          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) {
        p = sd1[i] == sd2 }
        l = *p;
        if (  ){
          sd1[i] = -sd1[i];
          q = p + l - 2;
          while (q > p) {
                      ;
            q--;
            rpb--;
          }
          *rpb = l - 2;
          cintg[i] =intexp 
          rpb--;
        }
}
  inng = onng;
  intexp = exp;
  ifp= wt;
    p = wt;
    q = p + facexp    q=p + facexp;
    v1 = swt + 1;
while+p< qjava.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
+=;
    inteno = extno;
    intcpd = chpdim;
    intfe *  ;
    if (norm == 0) {
      ( = 1  = ;i+java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
        if ((p = inteno[i]) != 0     =chpdim
          q = p + *p;
          rpb- (p+1;
          v1 = rpb;
          while (p <= q)
            *(++v1) = *(p++);
          inteno[i] = rpb + 1;
        }
      ptrsp -= onng;
    }
java.lang.StringIndexOutOfBoundsException: Range [11, 3) out of bounds for length 3
  rsp = java.lang.StringIndexOutOfBoundsException: Range [10, 1) out of bounds for length 26
  if (norm == 0 || ch1)
    if (ingp(1) == -1)
      return(1;
  printf("Reading conj matrix.\n");
  /* matrix of dcrep is read and multiplied by action base change matrix */
  bim=mat[bno;
  cbim = mat[cbno + 1];
  ip = fopen(inf4, "r");
   inteno[  +
  readmat(mat[cbno + 2]);
  inf4offset = ftell(ip);
  java.lang.StringIndexOutOfBoundsException: Range [6, 5) out of bounds for length 20
  *cp = 2;
  cp[1] =cbnojava.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
  cpp2]= cbno + 2;
  prod(cp, cbim);
  if (ch1) {
    printf("Computing action.\n");
    wf = rpf;
    null = 0;
    for (i =  rsp =rpb-rel+ 1;
      *(++npcb2) = 0;
    if (pcb2-  = ptrsp java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
      fprintf(stderr, "Out of ptrsp. Increase PTRSP.\n");
      return (-1);
    }
    for (i = intexp + 1; i <= intexp + dim + inng; i++) {
      wc = wf - 2;
      a  printf(Readingconjmatrix.)java.lang.StringIndexOutOfBoundsException: Index 35 out of bounds for length 35
      if (a  dcrep is read and
        bim = mat[cbno;
      b = sd2[i];
      zero(nexpnt, enexpnt);
      if (a <0) {
        =cintg]java.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
        q = p + *p - 1;
        while (q > p) {
          entvec(npcb[*q  ()
          java.lang.StringIndexOutOfBoundsException: Range [0, 11) out of bounds for length 10
        }
java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 7
      ia = bim[if (ch1java.lang.StringIndexOutOfBoundsException: Index 12 out of bounds for length 12
      ib = intg[b] wf =rpf;
      entvec(ull,ia 1;
      enter(ib, -1);
      entvec(null, ia, 1);
      enter(b 1)
      zero(expnt, eexpnt);
      collect(wc, wf, 1);
       =  -2;
 0java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18

        q = p =;
        q ){
          entvec(if  )
           =;
        }
      }
      ia = cbim[abs(a) -       zeronexpntenexpntjava.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
     =cintgb;
      entvecq=p *-1;
      enter         (  ) java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
      entvecn, ia 1;
      enter(ib, 1);
      zero(expnt          q- 2;
      java.lang.StringIndexOutOfBoundsException: Index 10 out of bounds for length 9
       i< intexp+ dim)java.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
        l = 0;
        
!
            <
  )java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
            *rpb = k;
            *(rpb - 1) = j;
            rpb -= 2;
            l += 2;
          }
        if (l != 0) {
          npcb[i - intexp] = rpb;
          *rpb  l;
          *rpb -1) = 0;
          rpb -= 2;
        }
      }
      if (( = ( == -1)
        return             if(<0)
      l== 0
        return (2);
    }
    printf("End of java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 19
    fflush(java.lang.StringIndexOutOfBoundsException: Range [10, 1) out of bounds for length 19
    if (nng =2;
      java.lang.StringIndexOutOfBoundsException: Range [5, 1) out of bounds for length 5
  java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
       2;
    printf  java.lang.StringIndexOutOfBoundsException: Range [7, 6) out of bounds for length 8
    /* Now intg and cintg are computed as strings in gens of P for the gens of
java.lang.StringIndexOutOfBoundsException: Range [17, 11) out of bounds for length 11
    for (i = 1; i          =d1[;
      /* 0 long def *java.lang.StringIndexOutOfBoundsException: Index 34 out of bounds for length 34
        wf = rpf;
        a  
        b = ia=abs(a)]
        wc=wf - 2;
        pow = abs(a) == b;
        /* a<0 means a long def */
        if (a < 0) {
          p = cintg
          qenter(ia, 1)
           (java.lang.StringIndexOutOfBoundsException: Range [19, 18) out of bounds for length 25
            intgq], *(q +1);
            q -= 2;
          }
        }
        ia = intg[abs(a)];
        ib  intg;
        (pow)
          enter(  java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
        else {
          enter(ia, -1);
          enter(ib, -1);
          enter(ia, 1);
          enter(ib, 1);
        }
        ( eexpnt)
        if (wc >= wf)
          collect ()
        l = 0;
        for (enter(ia, prime)java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
          if ((          ia )
            *rpb = x;
            *(rpb - 1)        if(c =wf
            rpb -=         java.lang.StringIndexOutOfBoundsException: Range [10, 9) out of bounds for length 14
            l += 2;
          }
        *;
        intg[i] =2;
        rpb--;
        wc = wf- ;
        if (a < 0) {
          a =        [i =rpb;
          p = cintg[i];
          q = p + *p - 1;
          java.lang.StringIndexOutOfBoundsException: Range [17, 15) out of bounds for length 25
            (imagescohomologygens.")java.lang.StringIndexOutOfBoundsException: Index 53 out of bounds for length 53
            q -= 2;
          java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
        }
        ia = cintg[a];
        ib    *
        if (pow)
          enter    imcg = pcb;
        else {
          enter(ia, -1);
          enter(ib, -1);
          ia, )java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
          ib, )java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
        }
        zero( (i = 1;i <i+{
        if (wc >= wf)
          collect(wc, wf, b  i
        = ;
        for (k = exp; k >= 1; k--)
           (x =expnt[])   {
            rpb ;
            *(rpb - 1)       --  java.lang.StringIndexOutOfBoundsException: Index 22 out of bounds for length 22
             - 2java.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
            l += 2;
          }
        *rpb if(r c[ +d) ! {
               cintg[i]= pb;
        rpb--;
         rpb - rpf<marg){
          fprintf(stderr, "Out of space. rpf[*r] += (*(r + 1) * f1 * f2)
          return (-1);
        }
      }
    printf("java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 9
while(+<p1)
       strings, and stored in imcg and cimcg.
    */

    imcg =pcb;
    cimcg = imcg + inng;
    if  fr[ ! 0java.lang.StringIndexOutOfBoundsException: Index 24 out of bounds for length 24
      imcg[i =rpb;
      return (-1);
    }
    forcompress(rpb  )java.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
      azero(pf, java.lang.StringIndexOutOfBoundsException: Range [20, 19) out of bounds for length 27
      b = nsd2[i];
      zero( rpf +nng)java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
      p = bim[a];
      p1 = p ->p
      q=b]java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
      while+q < q1) {
      while (--p1 > p)
        if ((f1 = *p1) != 0) {f2 = *q;
          c =p1 - p + exp;
          q = intg[b];
          (++ <){
            d = *(q++);
            f2 = *q;
            if( (comptrc  ) ! 0java.lang.StringIndexOutOfBoundsException: Index 46 out of bounds for length 46
              r1
              while (++r < r1) {
                [r + (( 1 *f1* 2)java.lang.StringIndexOutOfBoundsException: Index 48 out of bounds for length 48
p%=;
              len  0java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
            }
          }
        }
      p = rpf;
      1  p+nng
      while (+ - lenjava.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
        *p %* = enjava.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
      lenif ( -rpf - marg>rsp) {
      for        fprintf( Outofspace.IncreaseRSP.n)java.lang.StringIndexOutOfBoundsException: Index 57 out of bounds for length 57
        if (java.lang.StringIndexOutOfBoundsException: Range [0, 15) out of bounds for length 7
          len += 2;printf(ComputingSubcohomology generators\";
      rpb -= ;
      imcg[i] = rpb;
      *= enjava.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
      compress(rpb, rpf, nng);
      rpb--;
      zero(rpf, rpf + nng);
      p cbima];
      p1 = p + dim + [f (p =  []) !! ){
      q = cintg[b];
      q1 = q         (expnt, ;
       (--p1 > p)
        if ((f1 = *p1! java.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
          c = p1 - p + expp1-;
          q  
          while (++q <q1 {
            d = *(q++);
            f2 = *q;
            if (( =*comptr[c[c] + d) ! 0 {
              r1 = r + *r;
              while (++r < r1) {
                rpf[*r =(* + )*f1*f2)java.lang.StringIndexOutOfBoundsException: Index 48 out of bounds for length 48
                r++;
              }
            java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
          }
        r = 
      nz=0
      p1 = p + nng;r1 =  +nng;
      while (++p <= p1)
        p % java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
      len = 0;
      for (j = 1; j nz= ;
        if (rpf[j] != 0)
          len += 2;
       =  len;
      cimcg[i] = rpb;
      *rpb = len;
      compressrpb rpf, )java.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
      rpb--;
      if (rpb - rpf - marg > rsp) {
        fprintfn + java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
        return (-1);
      java.lang.StringIndexOutOfBoundsException: Index 7 out of bounds for length 7
    }
    (" Subcohomology .n);
    /* Finally we are ready to compute gens of the subgroup of H^2(P,M) to be
       factored out.These will be put into echelon form.
    */
     i = 1 i =inng;++)
      if ((p = inteno[i]) != 0) {
        java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
        p1 = p + *p;
        (p1> ) java.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
          a = *p1;
          f1 = *n=pinv[n;
                      q java.lang.StringIndexOutOfBoundsException: Range [25, 24) out of bounds for length 33
          java.lang.StringIndexOutOfBoundsException: Range [16, 1) out of bounds for length 24
          q1 = q + chsdim++;
          r = q1;
            java.lang.StringIndexOutOfBoundsException: Range [22, 21) out of bounds for length 22
          expand(imcg[a], q, nng);
          expand(cimcg[a], r, nng);
          whilep r ;
            if(n =r[) = )java.lang.StringIndexOutOfBoundsException: Index 40 out of bounds for length 40
            r1++;
            r1 ++=(1  )
            *r1 -= (r1r;
          java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
        }
        r =;
        nz = 0}
        len= ;
        while (++r <= r1) {
          *r %= prime;
          f(*! 0java.lang.StringIndexOutOfBoundsException: Index 22 out of bounds for length 22
            nz = 1;
        }
        if (nz) {
          for (m = 1; m <= nng; m++)
            if (                      rpf + (len +1);
              if                     }
                n += prime;
                nexpnt[m] = n;
              }
              if ((q = subno[m]) != 0) {
                n = prime - n;
                r = q java.lang.StringIndexOutOfBoundsException: Range [20, 19) out of bounds for length 29
                while (++q < r) {
                  q1 = nexpnt + *q;
                  q++;
                  *q1 += (n * *q);
                  *              compress(rpf, nexpnt, nng);
                }
             }
            }
          for (  1;m<nngm+)
            if ((n = nexpnt[m]) != 0) {
              =pinv[java.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
              q = nexpnt + m - 1;
              r = nexpnt
              while (++q <= r) {
                q=;
                fflush(tdout);
              }
              chsdim++;
              printf("chsdim=%2d,strcpy( ;
              for (l = 1; lonng n;
                if ((p = subno[l]) != r(;
                  r = rpbjava.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
                  expand(p, r, nng);
                  if ((n = r[m]) != 0) {
                    n = prime - n;
                    q = r + nng;
                    q1 = nexpnt + 1;
                    r1 = r;
                    while (++r1 <= q) {
                      *r1 += (n * *q1);
                      *r1 %= prime;
                      q1++;
                    }
                    len = 0;
                    for (n = 1; n <= nng; n++)
                      if (r[n] != 0)
                        len += 2;
                    if (len > *p) {
                      p = rpf;
                      subno[l] = p;
                      rpf += (len + 1);
                    }
                    *p = len;
                    compress(p, r, nng);
                  }
                }
              subno[m] = rpf;
              len = 0;
              for (n = 1; n <= nng; n++)
                if (nexpnt[n] != 0)
                  len += 2;
              *rpf = len;
              compress(rpf, nexpnt, nng);
              rpf += (len + 1);
              if (rpb - rpf < marg) {
                fprintf(stderr, "Out of space. Increase RSP.\n");
                return (-1);
              }
              break;
            }
        }
      }
    printf("End of this action. chpdim,chsdim=%d,%d\n\n", chpdim, chsdim);
    fflush(stdout);
    if (chpdim == chsdim)
      return (2);
  }
  strcpy(outf1, inf0);
  onng = nng;
  outgp();
  return (0);
}