cuikatlas.c

Go to the documentation of this file.

#include "world.h"
#include "parameters.h"

#include "defines.h"
#include "error.h"
#include "filename.h"
#include "atlas.h"
#include "random.h"
#include "geom.h"
#include "samples.h"

#include <stdlib.h>
#include <time.h>

int main(int argc, char **arg)
{
  TAtlasBase world;        /* The set of mechanism and obstacles.  */
  Tparameters parameters;  /* Parameters used in the Cuik process. */

  Tfilename fparam;
  Tfilename fatlas;
  
  double *s1;              /* One point on the manifold (values only for the system variables). */

  Tatlas atlas;
  Tstatistics st;

  unsigned int ri;

  if (argc>1)
    { 
      /*Init parameters*/
      CreateFileName(NULL,arg[1],NULL,PARAM_EXT,&fparam);
      fprintf(stderr,"Reading parameters from           : %s\n",GetFileFullName(&fparam));
      InitParametersFromFile(GetFileFullName(&fparam),&parameters);

      /*Read the world/cuik from file*/
      CS_WD_INIT(&parameters,arg[1],&world);
 
      /* Read start sample */
      ReadOneSample(&parameters,arg[1],CS_WD_GET_NUM_SYSTEM_VARS(&world),&s1);

      /* Random seed initialization */
      ri=(unsigned int)time(NULL);
      //ri=1405333052;
      randomSet(ri);
     
      fprintf(stderr,"Random seed                       : %u\n",ri);
      
      /* The statistics here are just used to compute execution time.
         In some cases, we use parallel atlas construction. Then
         CPU time accumulates the time for all the used CPUs, which
         is significantly larger than the real execution time. 
         Wall clock time has the problem of also accumulating the
         time of concurrent processes. */
      #ifdef _OPENMP
        InitStatistics(2,0,&st); /* just to force wall clock time */
      #else
        InitStatistics(0,0,&st); /* use real cpu time */
      #endif

      /* Define the atlas */
      BuildAtlasFromPoint(&parameters,s1,FALSE,&world,&atlas);

      fprintf(stderr,"Atlas compleated.\n");
      fprintf(stderr,"Elapsed time %f\n", GetElapsedTime(&st));
      
      DeleteStatistics(&st);

      /* Save the results */
      CreateFileName(NULL,arg[1],NULL,ATLAS_EXT,&fatlas);
      fprintf(stderr,"Writing atlas to                  : %s\n",GetFileFullName(&fatlas));
      SaveAtlas(&parameters,&fatlas,&atlas);
      DeleteAtlas(&atlas);
      DeleteFileName(&fatlas);

      /* Release memory */
      free(s1);

      DeleteParameters(&parameters);

      CS_WD_DELETE(&world);

      DeleteFileName(&fparam);
    }
 else
   {
     fprintf(stderr,"  Wrong number of parameters.\n");
     fprintf(stderr,"  Use:\n");   
     fprintf(stderr,"      cuikatlas <problem filename>.%s \n",CS_WD_EXT);
     fprintf(stderr,"  where <problem filename> the equations/world description\n");
     fprintf(stderr,"    (the '.%s' extension is not required)\n",CS_WD_EXT);
   }
  return(EXIT_SUCCESS);
}