/* st05.c Version 1.1 17/10/2007 FREE PRODUCT CJP */
#define MAXSTORE 1000000

#include<stdio.h>
#include<stdlib.h>
#include"blimp.h"

int main(int argc, char **argv)
{
    FILE *f1, *f2, *f3;
    int *store;
    int *gp1,*gp2,*gp3,*b1, *b2;
    int len1,len2,gen1,gen2,rel1,rel2,i,j,n,m;
    int *pt1, *pt2;
    int storeused=0;
    
    store = (int*) malloc(MAXSTORE*sizeof(int));
    gp1 = (int*) malloc(GRPSIZ*sizeof(int));
    gp2 = (int*) malloc(GRPSIZ*sizeof(int));
    gp3 = (int*) malloc(GRPSIZ*sizeof(int));
    b1 = (int*) malloc(GRPSIZ*sizeof(int));
    b2 = (int*) malloc(GRPSIZ*sizeof(int));

    f1 = (FILE*) fopen(argv[1],"r");
    f2 = (FILE*) fopen(argv[2],"r");
    f3 = (FILE*) fopen(argv[3],"w");

    
    if(argc!=4)
    {
        printf("\nUsage: st05 <input1> <input2> <output>\n");
        return 0;

    }
    
    if(f1==NULL)
    {
        printf("\nCould not open file %s\n",argv[1]);
        return 0;
    }

    if(f2==NULL)
    {
        printf("\nCould not open file %s\n",argv[2]);
        return 0;
    }
        
    if(f3==NULL)
    {
        printf("\nCould not open file %s\n",argv[3]);
        return 0;
    }

    // Read file 2, and store the series of blimp 2 in store.
    pt2=(int*) store;
    while(ReadGroup(gp2,f2)!=-1)
    {
        ExtractBlimp(gp2,2,b2);
        len2 = b2[1];
       
       
       
        storeused+=(len2+2);
        if(storeused>MAXSTORE)
        {
            printf("File 2 too large: ran out of store space");
            return 0;
        }
        for(i=0;i<len2+2;i++) // Whole length + blimp name+length
        {
            *(pt2++) = b2[i];
        }
        *(pt2++)=0;
   
    }

    fclose(f2);

    // Read in 1 group from file 1 at a time, * with each blimp2 in store,
    // Then write out incl. changing name.
    n=0;
    m=0;
     
    while(ReadGroup(gp1,f1)+1)
    {
        n++; // # of group in first file
        pt2=(int*) store;
        m=0;
        
        ExtractBlimp(gp1,2,b1);
        len1 = b1[1];
        gen1 = b1[2];
        rel1 = b1[3];
        
        while(*pt2)
        {
            pt1 = b1;
            pt1++; pt1++; pt1++; pt1++;
            m++; // # of group in second file
            pt2++; // move pt2 onto the length
            len2 = *(pt2++);
            gen2 = *(pt2++);
            rel2 = *(pt2++);
          
            // Create new group in gp3
            gp3[0]=1;
            gp3[1]=3;
            gp3[2]=10;
            gp3[3]=n;
            gp3[4]=m;
            gp3[5]=2;
            gp3[6]=len1 + len2 - 2;
            gp3[7]=gen1 + gen2;
            gp3[8]=rel1 + rel2;
            for(i=1;i<=len1-2;i++) // Copy relators from 1 unchanged
            {
                gp3[i+8] = *(pt1++); 
            }
            
            for(j=0;j<len2-2;j++)
            {
                if(*pt2==0) gp3[j+i+8]=0;
                else gp3[j+i+8] = *pt2 + gen1;
                pt2++;
            }
            gp3[i+j+8] = 999;
         
          
            WriteGroup(gp3,f3);
            pt2++;
        }


        
    }
    
    fclose(f1);
    WriteEOF(f3);
    fclose(f3);
}