cuiksilhouette.c

Go to the documentation of this file.

#include "cuiksystem.h"

#include "parameters.h"
#include "filename.h"
#include "defines.h"
#include "error.h"
#include "box_list.h"
#include "random.h"

#include "cp.h"

#include <stdlib.h>

int main(int argc, char **arg)
{
  TCuikSystem cuiksystem;  /* The set of equations */
  Tparameters parameters;  /* Parameters used in the Cuik process */
  Tbox binit;
  Tcp *cp;
  Tlist cps;
  unsigned int ncp;
  Tfilename fparam;
  Tfilename fcuik;
  Tfilename fout;
  FILE *f_out;

  if (argc>1)
    { 
      randomReset();
      
      /*Read the parameters*/
      CreateFileName(NULL,arg[1],NULL,PARAM_EXT,&fparam);
      InitParametersFromFile(GetFileFullName(&fparam),&parameters);

      /*Init the CuikSystem*/
      CreateFileName(NULL,arg[1],NULL,CUIK_EXT,&fcuik);
      InitCuikSystemFromFile(&parameters,GetFileFullName(&fcuik),&cuiksystem);

      /*Start the output file (where the output boxes will be stored)*/
      CreateFileName(NULL,arg[1],"_sil",SOL_EXT,&fout);
      f_out=fopen(GetFileFullName(&fout),"w"); 
      if (!f_out)
        Error("Could not open main output file");

      /*Create a fake critical point with an empty box to bootstratp the process*/
      InitBox(1,NULL,&binit);
      NEW(cp,1,Tcp);
      NewCP(0,&binit,cp);
      DeleteBox(&binit);

      /*Add the first critical point to the list*/
      InitList(sizeof(Tcp *),&cps);
      AddFirstElement(&cp,&cps);
      ncp=0; /*This is a counter of the number of control points processed*/

      /*Process all critical points in the list (when processin new
        critical points are generated) */
      while(!ListEmpty(&cps))
        {
          ExtractFirst(&cp,&cps);

          DealWithCP(f_out,ncp,&parameters,&cuiksystem,&cps,cp);
          ncp++;

          DeleteCP(cp);
          free(cp);
        }

      fclose(f_out);

      DeleteFileName(&fout);
      DeleteFileName(&fcuik);
      DeleteFileName(&fparam);
    }
  else
    {
      fprintf(stderr,"  Wrong number of parameters.\n");
      fprintf(stderr,"  Use:\n");   
      fprintf(stderr,"      cuiksilhouette <problem filename>.cuik \n");
      fprintf(stderr,"  where <problem filename> contains the kinematic equations\n");
      fprintf(stderr,"    (the '.cuik' extension is not required)\n");
    }
  
  return(EXIT_SUCCESS);
}