cuikatan2.c

Go to the documentation of this file.

#include "box.h"

#include "defines.h"
#include "filename.h"

#include <string.h>


int main(int argc, char **arg)
{
  FILE *fileIn,*fileOut;
  Tbox boxIn,boxOut;
  Tinterval iAtan2;
  unsigned int *map,*dimSin,*dimCos,nIn,nOut;
  unsigned int k,l,t,na,nbox;
  int token;
  Tfilename fIn,fOut;
  boolean found;

  if ((argc<4)||((argc%2)==0))
    {
      fprintf(stderr,"Use:\n");
      fprintf(stderr,"  cuikatan2 <input_file> <dim_sin_1> <dim_cos_1> ... <dim_sin_n> <dim_cos_n>  <output_file>\n");
      fprintf(stderr,"    <input_file> the input .sol file\n");
      fprintf(stderr,"    <dim_sin_i> <dim_cos_i> the sin/cos dimensions (numbered from 1)\n");
      fprintf(stderr,"    <output_file> the .sol file to be generated\n");
    }
  else
    {
      CreateFileName(NULL,arg[1],NULL,SOL_EXT,&fIn);
      CreateFileName(NULL,arg[argc-1],NULL,SOL_EXT,&fOut);
  
      fileIn=fopen(GetFileFullName(&fIn),"r");
      if (!fileIn)
        Error("Input file can not be opened");

      fileOut=fopen(GetFileFullName(&fOut),"w");
      if (!fileOut)
        Error("Output file can not be opened");


      na=(argc-3)/2;
      NEW(dimSin,na,unsigned int);
      NEW(dimCos,na,unsigned int);
      
      for(k=0;k<na;k++)
        {
          dimSin[k]=(unsigned int)atoi(arg[2+2*k])-1;
          dimCos[k]=(unsigned int)atoi(arg[3+2*k])-1;
        }

      nbox=0;
      do {
        token=ReadBox(fileIn,&boxIn);
        if (token!=EOF)
          {
            if (nbox==0)
              {
                nIn=GetBoxNIntervals(&boxIn);

                nOut=nIn-na;

                NEW(map,nOut-na,unsigned int);
                l=0;
                for(k=0;k<nIn;k++)
                  {
                    found=FALSE;
                    t=0;
                    while((!found)&&(t<na))
                      {
                        found=((k==dimSin[t])||(k==dimCos[t]));
                        t++;
                      }
                    if (!found)
                      {
                        map[l]=k;
                        l++;
                      }
                  }
              }

            InitBox(nOut,NULL,&boxOut);
            for(k=0;k<nOut-na;k++)
              SetBoxInterval(k,GetBoxInterval(map[k],&boxIn),&boxOut);
            
            for(k=0;k<na;k++)
              {
                IntervalAtan2(GetBoxInterval(dimSin[k],&boxIn),
                              GetBoxInterval(dimCos[k],&boxIn),
                              &iAtan2);
                SetBoxInterval(nOut-na+k,&iAtan2,&boxOut);
                DeleteInterval(&iAtan2);
              }

            PrintBox(fileOut,&boxOut);

            DeleteBox(&boxIn);
            DeleteBox(&boxOut);

            nbox++;
          }     
      } while (token!=EOF);
      

      fclose(fileIn);
      fclose(fileOut);


      if (nbox>0)
        free(map);
      free(dimSin);
      free(dimCos);

      DeleteFileName(&fIn);
      DeleteFileName(&fOut);
      
    }
  return(EXIT_SUCCESS);
}