cuikplot3d.c

Go to the documentation of this file.

#include "boolean.h"
#include "plot3d.h"
#include "box.h"
#include "defines.h"
#include "filename.h"
#include "color.h"

#include <stdlib.h> 
#include <stdio.h>
#include <string.h>
#include <math.h>

int main(int argc, char **arg)
{
  Tbox box;
  unsigned int nbox;
  int token;
  Tplot3d plot;
  FILE *PSin; 
  Tinterval *x,*y,*z;
  unsigned int dim_x,dim_y,dim_z;
  double average_x,average_y,average_z;
  double min_x,max_x,min_y,max_y, min_z,max_z;
  double minSize;
  Tfilename fsols,fplot;
  Tcolor plotColor;
  
  if (argc<7)
    {
      fprintf(stderr,"Use:\n");
      fprintf(stderr,"  cuikplot3d <filename> dimx dimy dimz min_size <plotname> \n");
      fprintf(stderr,"    <filename> the input .sol file\n");
      fprintf(stderr,"    dimx dimy dimz -> the three dimensions to be plotted (numbered from 1)\n");
      fprintf(stderr,"    min_size -> minimum size for the boxes\n");
      fprintf(stderr,"    <plotname> the output .cgl file\n");
    }
  else
    {
      CreateFileName(NULL,arg[1],NULL,SOL_EXT,&fsols);

      /*"Parse" the input parameters*/
      
      /*Input file name*/
      PSin=fopen(GetFileFullName(&fsols),"r");
      if (!PSin)
        Error("Input file can not be opened");

      /*Dimension on which to project*/
      dim_x=(unsigned int)atoi(arg[2]);
      dim_y=(unsigned int)atoi(arg[3]);
      dim_z=(unsigned int)atoi(arg[4]);

      /*Check that the dimensions are ok (positive not repeated)*/
      if ((dim_x==0)||(dim_y==0)||(dim_z==0))
        Error("Dimensions are numbered from 1\n");

      if ((dim_x==dim_y)||(dim_x==dim_z)||(dim_y==dim_z))
        Error("Plotting dimensions must be different\n");

      dim_x=dim_x-1;
      dim_y=dim_y-1;
      dim_z=dim_z-1;

      /*Compute the center of mass of the boxes in the selected dimension*/
      average_x=0.0;
      average_y=0.0;
      average_z=0.0;

      nbox=0;
      do {
        token=ReadBox(PSin,&box);
        if (token!=EOF)
          {
            average_x+=IntervalCenter(GetBoxInterval(dim_x,&box));
            average_y+=IntervalCenter(GetBoxInterval(dim_y,&box)); 
            average_z+=IntervalCenter(GetBoxInterval(dim_z,&box)); 
            nbox++;
            DeleteBox(&box);
          }
      } while (token!=EOF);

      if (nbox>0)
        {
          average_x/=(double)nbox;
          average_y/=(double)nbox;
          average_z/=(double)nbox;
        }
      else
        {
          average_x=0.0;
          average_y=0.0;
          average_z=0.0;
        }
          
      /*Rewin the file*/
      fclose(PSin);

      PSin=fopen(GetFileFullName(&fsols),"r");

      average_x=0.0;
      average_y=0.0;
      average_z=0.0;

      /*Get the minSize parameter*/
      minSize=(double)atof(arg[5]);

      /*and now we plot in 3D*/
      CreateFileName(NULL,arg[6],NULL,PLOT3D_EXT,&fplot);
      InitPlot3d(GetFileFullName(&fplot),TRUE,argc,arg,&plot);
      DeleteFileName(&fplot);

      NewColor(DCP3D_R,DCP3D_G,DCP3D_B,&plotColor);
      StartNew3dObject(&plotColor,&plot);
      DeleteColor(&plotColor);

      do {
        token=ReadBox(PSin,&box);
        if (token!=EOF)
          { 
            /*Get the interval of the box in the projection dimentions*/
            x=GetBoxInterval(dim_x,&box);
            y=GetBoxInterval(dim_y,&box);
            z=GetBoxInterval(dim_z,&box);
           
            /*and get the extreme of the intervals*/
            min_x=LowerLimit(x) ;
            max_x=UpperLimit(x) ;
            
            min_y=LowerLimit(y) ;
            max_y=UpperLimit(y) ;
            
            min_z=LowerLimit(z) ;
            max_z=UpperLimit(z) ;

            /*if the box is too small we enlarge it (otherwise we have "invisible" boxes*/
            if (((max_x-min_x)<minSize)&&
                ((max_y-min_y)<minSize)&&
                ((max_z-min_z)<minSize))
              {
                double r;

                r=minSize/2;

                max_x+=r; min_x-=r;
                max_y+=r; min_y-=r;
                max_z+=r; min_z-=r;
              }

            /*and add the box to the plot*/

            PlotBox3d(min_x-average_x,max_x-average_x,
                      min_y-average_y,max_y-average_y,
                      min_z-average_z,max_z-average_z,
                      &plot);

            //if (isolated)
            //  StartNew3dObject(0,0,1,&plot);

            DeleteBox(&box);
          }
      } while (token!=EOF);

      ClosePlot3d(FALSE,average_x,average_y,average_z,&plot);

      DeleteFileName(&fsols);
      fclose(PSin);
    }
  
  return(EXIT_SUCCESS);
}