cuikplotpath.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 "samples.h"

#include <stdlib.h>

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

  unsigned int nx,ny,nz;
  unsigned int i,nv,nvs;  
  boolean *systemVars;

  unsigned int ns;
  double **path;

  Tfilename fparam;
  Tfilename fpath;

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

      /*Read the world from file*/ 
      CS_WD_INIT(&parameters,arg[1],&world);
      nv=CS_WD_GET_SYSTEM_VARS(&systemVars,&world);
          
      if (argc>5)
        {
          CreateFileName(NULL,arg[2],NULL,SOL_EXT,&fpath);
      
          nx=atoi(arg[3]);
          ny=atoi(arg[4]);
          nz=atoi(arg[5]);
        }
      else
        {
          CreateFileName(NULL,arg[1],"_path",SOL_EXT,&fpath);
      
          nx=atoi(arg[2]);
          ny=atoi(arg[3]);
          nz=atoi(arg[4]);
        }

      nvs=0;
      for(i=0;i<nv;i++)
        {
          if (systemVars[i])
              nvs++;
        }
      if ((nx>=nvs)||(ny>=nvs)||(nz>=nvs))
        Error("Projection dimensions are beyond the ambient space dimensionality");

      free(systemVars);

      fprintf(stderr,"Projecting on variables           : ");
      PRINT_VARIABLE_NAME(stderr,CS_WD_GET_SYSTEM_VAR_NAME(nx,&world));fprintf(stderr,",");
      PRINT_VARIABLE_NAME(stderr,CS_WD_GET_SYSTEM_VAR_NAME(ny,&world));fprintf(stderr,",");
      PRINT_VARIABLE_NAME(stderr,CS_WD_GET_SYSTEM_VAR_NAME(nz,&world));fprintf(stderr,"\n");

      if (LoadSamples(&fpath,&nvs,&ns,&path))
        {
          /*If the path was actually loaded*/
          Tfilename fpathplot;
          Tplot3d p3d;

          if (argc>5)
            CreateFileName(NULL,arg[2],NULL,PLOT3D_EXT,&fpathplot);
          else
            CreateFileName(NULL,arg[1],"_path",PLOT3D_EXT,&fpathplot);

          fprintf(stderr,"Generating path plot to           : %s\n",GetFileFullName(&fpathplot));
              
          InitPlot3d(GetFileFullName(&fpathplot),FALSE,argc,arg,&p3d);

          PlotSamples(&parameters,&p3d,nx,ny,nz,ns,path);

          ClosePlot3d(FALSE,0,0,0,&p3d);
          
          DeleteFileName(&fpathplot);
          DeleteSamples(ns,path);
        }
      else
        Error("Can not read the file with the path");

      DeleteParameters(&parameters);
      CS_WD_DELETE(&world);

      DeleteFileName(&fpath);
      DeleteFileName(&fparam);
    }
 else
   {
     fprintf(stderr,"  Wrong number of parameters.\n");
     fprintf(stderr,"  Use:\n");   
     fprintf(stderr,"      cuikplotpath <problem filename> [<path filename>] <xID> <yID> <zID>\n");
     fprintf(stderr,"  where <problem filename> is the description of the problem.\n");
     fprintf(stderr,"        <path filename> is the path to plot. If not given we use <problem filename>_path.sol\n");
     fprintf(stderr,"        <xID> <yID> <zID> are the 3 projection dimensions.\n");
     fprintf(stderr,"          (indices of system variables as those in the _path.sol file from 0)\n");
   }
  return(EXIT_SUCCESS);
}