/* st08.c CJP Version 1.01 17/10/07 */

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

int main(int argc, char **argv)
{
    int i,j,k;
    FILE *f1,*f2,*f3;
    int *gp,*b,*inv,*rel,*rel2;
    int good,changed,len,current,test;
    int *pt,*ptc,*ptt,*pt2;
    NTY n;
    NTY *nt;

    nt = &n;
    ntycons(nt);
    
    if(argc!=4)
    {
        printf("Usage: st08 <input> <NoTetrafound> <TetraFound>\n");
        return 0;
    }
        
    f1 = (FILE*) fopen(argv[1],"r");
    if(f1==NULL)
    {
        printf("Error opening file %s",argv[1]);
        return 0;
    }
    
    f2 = (FILE*) fopen(argv[2],"w");
    if(f2==NULL)
    {
        printf("Error opening file %s",argv[2]);
        return 0;
    }
    
    f3 = (FILE*) fopen(argv[3],"w");
    if(f3==NULL) 
    {
        printf("Error opening file %s",argv[3]);
        return 0;
    }
    
    good=1;
 
    gp = (int*) malloc(sizeof(int)*GRPSIZ);
    if(gp==NULL) good=0;
    b = (int*) malloc(sizeof(int)*GRPSIZ);
    if(b==NULL) good=0;
    inv = (int*) malloc(sizeof(int)*GRPSIZ);
    if(inv==NULL) good=0;
    rel = (int*) malloc(sizeof(int)*GRPSIZ);
    if(rel==NULL) good=0;
    rel2 = (int*) malloc(sizeof(int)*GRPSIZ);
    if(rel2==NULL) good=0;
    
    if(good==0)
    {
        printf("Error allocating memory\n");
        return 0;
    }
  
    while(1)
    {
        changed=0;

        i = ReadGroup(gp,f1);
        if(i==-1) break;

        ExtractBlimp(gp,3,inv);
        k=inv[1];
        for(i=0;i<k;i++)
        {
            inv[i] = inv[i+2]; 
        }
        ntystart(nt,inv[0],inv);
    
        // Insert all triangles into NTY
        ExtractBlimp(gp,6,b);
        pt = &b[3];
        for(i=0;i<b[2];i++)
        {
            len=0;
            while((*pt)!=0)
            {
                rel[len]=(*pt);
                len++;
                pt++;
            }
            rel[len]=0;
            pt++;

            if(len==3)
            {
                ntyinsert(nt,rel,3);
            }
        }
 
        // Search for Tetrahedra
        current = 1;
        test = ntynexnts(nt,current,-1);
        while(1)
        {
            current = ntynexnts(nt,current,-1);
            if(current==2) break;

            test = 1;
            ptc = ntyntptr(nt,current);
             
            while(test!=2)
            {
                test = ntynexnts(nt,test,-1);
                if(test==2) break;

                ptt = ntyntptr(nt,test);                
            
                if(ptc[0]!=ptt[0]) continue;
               
                rel[0] = ptc[2];
                rel[1] = inv[ptt[2]];
                for(i=1;i<inv[0];i++)
                {
                    rel[2] = i;
                    
                    j=ntyfind(nt,rel,3);
                    if(j==0) continue;
                    rel2[0] = i;
                    rel2[1] = inv[ptc[1]];
                    rel2[2] = ptt[1];
                    if(ntyinsert(nt,rel2,len)!=0)
                    {

                        changed=1;
                    }
                }
               

            }
        }
        // Write out
        if(changed==0) WriteGroup(gp,f2);
        if(changed==1)
        {
            b[0] = 6;
            b[1] = 1;
            b[2] = 0;
            pt = &b[3];
            current=1;
            while(1)
            {
                current = ntynexrel(nt,current);
                if(current==2) break;
                pt2 = ntyntptr(nt,ntyrelnt(nt,current));
                i = ntyntlen(nt,ntyrelnt(nt,current));
                b[2]++;
                
                for(j=0;j<i;j++)
                {
                    (*pt) = pt2[j];
                    pt++;
                    b[1]++;
                }
                (*pt) = 0;
                pt++;
                b[1]++;
            }
            DeleteBlimp(gp,6);
            InsertBlimp(gp,b);
            WriteGroup(gp,f3);        
        }
        fflush(f2);
        fflush(f3);
    }
    WriteEOF(f2);
    WriteEOF(f3);
    fclose(f1);
    fclose(f2);
    fclose(f3);
    
}