cuikunsimplify.c

Go to the documentation of this file.

 
#include "cuiksystem.h"
#include "parameters.h"
 
#include "defines.h"
#include "error.h"
#include "box_list.h"
#include "filename.h"
 
#include <string.h>
#include <stdlib.h>
#include <time.h>
 
int main(int argc, char **arg)
{
  TCuikSystem cuiksystem;  /* The set of equations */
  Tparameters parameters;  /* Parameters used in the Cuik process */
 
  Tlist solutions;  /* The list with the simplified solutions. */ 
  Titerator it;     /* Iterator on the solution list */
  Tbox boxO;        /* Space where to recompose the un-simplified box. */ 
 
  Tfilename fcuik;
  Tfilename fsol;
  Tfilename fsolsimp;
  Tfilename fparam;
 
  boolean *systemVars;
  
  FILE *fs;
  
  if (argc>2)
    {
      /*Init parameters*/          
      CreateFileName(NULL,arg[1],NULL,PARAM_EXT,&fparam);
      fprintf(stderr,"Reading parameter file: %s\n",GetFileFullName(&fparam));
      InitParametersFromFile(GetFileFullName(&fparam),&parameters);
 
      /*Read the problem from file*/      
      CreateFileName(NULL,arg[1],NULL,CUIK_EXT,&fcuik);
      fprintf(stderr,"Reading equation file : %s\n",GetFileFullName(&fcuik));
      InitCuikSystemFromFile(&parameters,GetFileFullName(&fcuik),&cuiksystem);
      GetCSSystemVars(&systemVars,&cuiksystem);
 
      /*Read the simplified boxes*/
      CreateFileName(NULL,arg[2],NULL,SOL_WITH_DUMMIES_EXT,&fsolsimp);
      fprintf(stderr,"Reading boxes file    : %s\n",GetFileFullName(&fsolsimp));
      if (!ReadListOfBoxes(GetFileFullName(&fsolsimp),&(solutions)))
        Error("Solution file can not be opened");
        
      
      /*Save the un-simplified boxes*/
      CreateFileName(NULL,arg[1],NULL,SOL_EXT,&fsol);
      fprintf(stderr,"Writting boxes file   : %s\n",GetFileFullName(&fsol));
      fs=fopen(GetFileFullName(&fsol),"w");
      if (!fs)
        Error("Could not open the output file for the un-simplified boxes.");
      
 
 
      InitIterator(&it,&solutions);
      First(&it);
      while(!EndOfList(&it))
        {
          RegenerateOriginalBox(&parameters,(Tbox *)GetCurrent(&it),&boxO,&cuiksystem);
          PrintBoxSubset(fs,systemVars,NULL,&boxO);
          DeleteBox(&boxO);
          Advance(&it);
        }
      fclose(fs);
      
      /*Remove the allocated objects*/
      DeleteListOfBoxes(&solutions);
      DeleteParameters(&parameters);
      DeleteCuikSystem(&cuiksystem);
 
      DeleteFileName(&fparam);
      DeleteFileName(&fcuik);
      DeleteFileName(&fsol);
      DeleteFileName(&fsolsimp);
      free(systemVars);
    }
 else
   {
     fprintf(stderr,"  Wrong number of parameters.\n");
     fprintf(stderr,"  Use:\n");   
     fprintf(stderr,"      cuikunsimplify <problem filename>.cuik  <solution filename>.sol\n");
     fprintf(stderr,"  where <problem filename> contains the original kinematic equations\n");
     fprintf(stderr,"        <solution filename> contains the solutions for the simplified system\n");
     fprintf(stderr,"    (the '.cuik' and '.sol' extensions are not required)\n");
   }
  return(EXIT_SUCCESS);
}