Quelle  nqp2.c

  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;
 0jsp;         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
}