/*****************************************************/
/*

cc $INC -O oa.c -o prog $LIB

prog filename

*/
/*****************************************************/
#include <math.h>
#include <stdlib.h>
#include <stdio.h>


#define LINE_LENGTH 500

/*Dynamic Memory Functions*/
int    **imatrix(int n, int m);
int    *ivector(int length);

int  GetAllRuns(int **runs, int m, int *current);

/*****************************************************/
main()
{
int  **runs, Nruns, NRbar, *ivec, *ivec2;
int  i, j, k, N, m, levels=2, strength=3, index=0;
char input[30], output[30], line[LINE_LENGTH];
FILE *IN, *OUT;

m=8;
sprintf(input, "oa8_s3lam2_e6_all.sol");

Nruns=1;
for(i=0; i<m; ++i) Nruns*=2;
ivec=ivector(m);
runs=imatrix(Nruns, m);
Nruns=GetAllRuns(runs, m, ivec);
ivec2=ivector(Nruns);

IN=fopen(input, "r"); 
while(fgets(line, LINE_LENGTH, IN)!=NULL)
{ 
 if(strncmp(line, "cur", 3)==0) 
  {
   N=0;
   index++;
   for(i=0; i<Nruns; ++i)
    {
     fscanf(IN, "%d", &ivec2[i]);
     N+=ivec2[i];
    }

   sprintf(output, "OA.%d.%d.%d.%d.%d.txt", N, m, levels, strength, index);
   OUT=fopen(output, "w");
   for(i=0; i<Nruns; ++i)
    for(j=0; j<ivec2[i]; ++j)
     for(k=0; k<m; ++k)
      {
       fprintf(OUT, "%2d ", runs[i][k]);
       if(k==m-1) fprintf(OUT, "\n");
      }
   fclose(OUT);
  }
}
fclose(IN);

}
/*****************************************************/
int GetAllRuns(int **runs, int m, int *current)
{
int i, j, done=0, Nruns=0;

for(i=0; i<m; i++) current[i]=0;
while(!done)
{
 for(j=m-1; j>=0; j--)
 {
  if(current[j]==2)
   {
    if(j==0) done=1;
    current[j-1]++;
    current[j]=0;
   }
 }
 if(done) break;
 for(j=0; j<m; j++) runs[Nruns][j]=2*current[j]-1;
 Nruns++;
 current[m-1]++;
}

return Nruns;
}
/*****************************************************/
int **imatrix(int n, int m)
{
  int i, **x, *xnow;

if((x = (int**)malloc((unsigned)(n * sizeof(int*))))==NULL){
 printf("\nError: allocation failure 1 in imatrix().\n\n"); exit(-1);}

if((xnow = (int*)malloc((unsigned)(n * m * sizeof(int))))==NULL){
 printf("\nError: allocation failure 2 in imatrix().\n\n"); exit(-1);}

for(i=0; i<n; i++, xnow += m) x[i] = xnow;

 return x;
}
/*****************************************************/
int *ivector(int length)
{
  int i, *x;

if((x = (int*)malloc((unsigned)(length * sizeof(int))))==NULL){
 printf("\nError: allocation failure in ivector().\n\n"); exit(-1);}

 return x;
}
/*****************************************************/