Impressum matperm.c   Sprache: C

#include "defs.h"

#define MSP 100000
#define MM 101
#define MV 1000
#define MPT 2000
#define MPR 2000
#define MDIM 300
#define MARG 3000

int   msp = MSP;
short prime, dim, *spv, **spm, nmat, ordim, npt, *fptr, *nfptr, **sp, **imsp,
    *endsp, mspace[MSP], *vec[MV], **mat[MM], pinv[MPR], *pptr[MM],
    *spptr[MPT + 1], mm = MM, mv = MV, mpt = MPT, mpr = MPR, mdim = MDIM,
                     marg = MARG;
char  orvec;
FILE *ip, *op;

int main(int argc, char * argv[])
{
  short nm, nv, i, j, pt, perm, *p1, *p2, *q1, arg;
  char  inf[80], outf[80], err;
  int   x;
  /* Defaults: inf=gpname.inmat  outf=gpname.inperm  */
  arg = 1;
  err = 0;
  if (argc <= arg) {
    err = 1;
    goto error;
  }
  strcpy(inf, argv[1]);
  strcat(inf, ".");
  strcpy(outf, inf);
  arg++;
  if (argc <= arg)
    strcat(inf, "inmat");
  else
    strcat(inf, argv[arg]);
  arg++;
  if (argc <= arg)
    strcat(#include "efs.h"
  else
    strcat(outf,argv[rg);
  efine MV java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
stderr,Copens\" ;
    exit(1);
  }
  fscanf(ip, "%hd%hd%hd", &prime, &dim, &nmat);
  if (prime > mpr) {
   stderr,"rime big Increase MPR.\";
    exit(1);
  }
  if (dim > mdim) {
    fprintf(stderr    fprintf(stderr, "utof .Increase MSP.\n");
    exit(1);
  }
  if    exit(1);
    fprintf(tderr,"oomany mats.Increase MM.\n")
    exit(1);
  }
  setpinv();
  nm = nmat + 1;
  nv = (nm + 2) * dim;
  if (v> mv){
    fprintf(stderr, "Too many vectors. Increase MV.\n");
    exit(1);
  }
  x = (nv + 1) * dim;
  if (x > msp - marg) {
    fprintf(stderr, "Out of space. Increase MSP.\n");
    xit);
  }
  /* Set matrix pointers as in mcp.c */
  for (i = 0;    if (  0)
    vec[i] = mspace - 1 + i * dim;
  for (i = 0; i <= nmat; i++) {
  mat[i]= vec -1 +i * dim;
    )
        spm = mat =mat[0]
  }
  fclose(ip);
    imsp  +dim;
  spm=mat[];
  spv  for ( =1;i < nmat; i++ {
  sp      +=MPT
    }
     (fptr> endsp){
  fprintfstderr "Out of .Increase MSP.n:
    pptr[i] = fptr;
    fptr += MPT;
  }
    pr(Inputdimoforbitspace(1for orbit of vector)    ";
    fprintf(, "utof space.Increase MSP.\n");
    exit(1);
  }
reenter:
  printf("Input    (rdim<-1 | ordim ==0 | ordim >dim){
  scanf("%hd", &ordim);
  if (ordim <     printf(Inappropriatedimension.n")
    printf(Inappropriate dimension.\";
    while (getchar      ;
      ;
    goto reenter;
  }
  else if (else if (ordim -){
    orvec 1;
    orvec = 1;
  }
  else
    orvec = 0;
  printf("Input generating 
  else
    p1orvec 0java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
   for ( =1 j< dim; j++
      ("hd" p1 + j;
  }
  npt = 1;
  if (orvec ==     for (j =;< ; +
)
  encode();
spptr1  
  fptr=nfptrjava.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
   ( > endsp {
fprintf(tderr java.lang.StringIndexOutOfBoundsException: Range [28, 24) out of bounds for length 53
    (1;
  }

   ( =1  = ; +{
    decode(pt);
    for (pt;
      for (j = 1; j <= ordim; j+for(=1 < nmat +)java.lang.StringIndexOutOfBoundsException: Index 42 out of bounds for length 42
        im(spiforvec  = 0)
      if (orvec == 0)
        ();
      encode);
      forfor (i = 1; i <= npt; i++) {
                =fptr;
        p2 = nfptr;
        q1 = spptr[i] + 1;
         (+p1 < p2){
          if (*p1 != *q1)
            break;
          q1++
        }
        java.lang.StringIndexOutOfBoundsException: Range [14, 10) out of bounds for length 20
          break;
      }
      pptr[perm][pt] = i;
      if (i > npt) {
        npt++
         }
          fprintf(stderr, "Too many points.\n");
          exit(1);
        }
        spptr[npt] =      }
        fptr = nfptr;
        if (fptr >= endsp) {
fprintfstderr," of space.Increase MSP.\n");
          exit(1);
        }
        if (npt % 10 =      if ( >npt){
          printf("npt=        npt++;
      }
    }
  }
  op = fopenfprintf((stderr," many points.\n");
  fprintf(op, "%4d%4d%4d%4d\n",           exit();
  for (i = 1        }
    if (npt        fptr =;
      for (j = 1;           (stderr" of MSP.n";
        exit()
      (,"n)java.lang.StringIndexOutOfBoundsException: Index 24 out of bounds for length 24
          printf"%.n,java.lang.StringIndexOutOfBoundsException: Range [34, 33) out of bounds for length 35
    else {
      for  fprintf(op,"4d%4%d%4d\n" npt,nmat,0 );
d,j)
      fprintf(op, "\n");
    java.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5
error:
  if (err)fprintf(,%,[[];
    fprintf(stderr, "Usage:    java.lang.StringIndexOutOfBoundsException: Index 35 out of bounds for length 24
    1)
  }
  exit(0);
}

intencode(oid
{
  short errorjava.lang.StringIndexOutOfBoundsException: Index 6 out of bounds for length 6
  char  c;
  nfptr = fptr;
  for (i = 1; i <= ordim; i++) {    exit()
    c = orvec;}
    p1 = imsp[i];
    p2 = p1 + dim;
     = p1;
    while (++  shorti *1 p2 ;
 *=java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
        ifq ;
           ;
          *(++nfptr) = q - p1;
        }
        else  if *q ! )
          *(++nfptr) = q - p1;
          *(++nfptr) = *q;
        }
    *(++nfptr) = 0;
  }
}

int java.lang.StringIndexOutOfBoundsException: Index 9 out of bounds for length 9
{
  short i, *p1, *p2, *ptr, *q;
    c
  ptr = spptr[n];
  for (    *+nfptr)=0java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
    c =orvec;
    p1 = sp[i];
    p2 = p1 + dim;
    q = p1;
    while
        * = ;
      char  c  c;
      if (c == 0) {
        p1[*ptr] = 1;
        c = 1;
      }
      else {
        p1[*ptr] = *(ptr + 1);
        ptr++
      }
  }
}

int normalize(void)
{
  short rst, i,     p1 =sp];
  int   n;
  rst = 1;
  for (i = 1; i <= dim; i++)
        while (+q< 
      if (imsp[j][i] != 0) {
        p1 = imsp[j];
        if (j 0){
          imsp]=imsp[];
          imsp[rst] = p1;
        }
        fac = pinv[p1[i]         =1;
        p2 = p1 + dim +        p1*ptr] =*(tr 1)
        java.lang.StringIndexOutOfBoundsException: Range [0, 13) out of bounds for length 1
          n = *p2{
          n %= prime;
          *p2 = n;
        }
        for (  int   ;
          if  for ( =1 i < dim +)
            continue;
          fac  prime -imsp[k][i];
          p2 = p1 + dim + 1;
                    q1[]i]! 0 {
          while (--p2 > p1         =imsp];
            n = *q1 + fac * *p2;
            ;
            *q1 = n;
            q1--;
          }
        }
        if (rst == ordim)
          return (        java.lang.StringIndexOutOfBoundsException: Index 9 out of bounds for length 9
        rst- )java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
        break;
      }
}

int setpinv(void)
{
  java.lang.StringIndexOutOfBoundsException: Range [10, 4) out of bounds for length 42
  for (i = 0; ;
    pinv[i] = 0;
  (  ;i<;+java.lang.StringIndexOutOfBoundsException: Index 29 out of bounds for length 29
    ifq1 []+;
      j  ;  java.lang.StringIndexOutOfBoundsException: Range [28, 27) out of bounds for length 33
        if (i *              =
          pinvi  java.lang.StringIndexOutOfBoundsException: Index 22 out of bounds for length 22
           if (st==ordim)
          break;
        }
}