plot3d.c

Go to the documentation of this file.

#include "plot3d.h"
#include "error.h"

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

/* Documented but not included in the doxygen documentation because it
   screws up the IRI fomat. */
#define CYLINDER_STRING "{INST geom { CBEZ223 2 1 0 1 1 0 0 1 0 2 1 1 1 1 1 0 1 1 2 0 1 1 0 1 0 0 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 } transforms{ = TLIST 0.5 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0.5 0 0 0 0 0 1 0 0 -1 0 0 0 0 0 1 0.5 0 0 0 0 -1 0 0 0 0 -1 0 0 0 0 1 0.5 0 0 0 0 0 -1 0 0 1 0 0 0 0 0 1}}"

void InitPlot3d(char *name,boolean axes,int argc, char **arg,Tplot3d *p)
{
  char command[200];
  unsigned int i;

  if (name!=NULL)
    {
      p->f=fopen(name,"w");

      if (!p->f)
        Error("Can not create output file in InitPlot3d");
    }
  else
    p->f=stdout;

  //start the ooglfile
  if (argc>0)
    {
      fprintf(p->f,"#\n# Geomview plot generated with the following command line:\n#      ");
      for(i=0;i<(unsigned int)argc;i++)
        fprintf(p->f,"%s ",arg[i]);
      fprintf(p->f,"\n#\n");
    }

  fprintf(p->f,"(progn ");
  fprintf(p->f,"(normalization World none)\n");  
  fprintf(p->f,"(bbox-draw World no)\n");

  /*do not show the main pannel of geomview*/
  /*fprintf(p->f,"(ui-panel geomview off)\n");*/

  /*show the tools pannel of geomview*/
  /*fprintf(p->f,"(ui-panel tools on)\n");*/

  fprintf(p->f,"(backcolor World 1.0 1.0 1.0)\n"); /*background color de fons*/

  fprintf(p->f,"(dice World 15.0)\n"); /*granularity of curbed objects*/
  fprintf(p->f," )\n");


  if (name==NULL)
    {
      NEW(p->fileName,7,char);
      strcpy(p->fileName,"stdout");
    }
  else
    {
      i=strlen(name)+1;
      NEW(p->fileName,i,char);
      strcpy(p->fileName,name);
      i=0;
      while (p->fileName[i]!=0)
        {
          if ((p->fileName[i]=='/')||(p->fileName[i]=='.'))
            p->fileName[i]='_';
          i++;
        }
    }

  if (axes)
    {
      //axes (we assign them identifier 0)
      fprintf(p->f,"(geometry sol_%s_0 ",p->fileName);
      if (p->f!=stdout) fclose(p->f);

      if (p->f==stdout)
        sprintf(command,"cat %s/share/axes.list",_CUIK_SUITE_MAIN_DIR);
      else
        sprintf(command,"cat %s/share/axes.list >> %s",_CUIK_SUITE_MAIN_DIR,name);
      
      system(command);

      if (p->f!=stdout) p->f=fopen(name,"a");
      fprintf(p->f," )\n");
    }

  p->nobj=0; /* Last object added to the scene */

  NewColor(DCP3D_R,DCP3D_G,DCP3D_G,&(p->color));
  
  p->inObject=FALSE;
}

void Start3dBlock(Tplot3d *p)
{
  fprintf(p->f,"(progn \n");
}

void Close3dBlock(Tplot3d *p)
{
  fprintf(p->f,")\n");
  fflush(p->f);
}

void  Delete3dObject(unsigned int nobj,Tplot3d *p)
{
  fprintf(p->f,"(delete sol_%s_%u )\n",p->fileName,nobj);
}

unsigned int StartNew3dObject(Tcolor *c,Tplot3d *p)
{
  Close3dObject(p);

  p->inObject=TRUE;
  CopyColor(&(p->color),c);

  p->nobj++; /* We increas the object identifier (empty objects also use an identifier although they are never created) */

  fprintf(p->f,"(geometry sol_%s_%u {LIST \n",p->fileName,p->nobj); 

  return(p->nobj); /*Return the identifier that will be assigned to the new object*/
}

void Close3dObject(Tplot3d *p)
{
  if (p->inObject)
    {
      /*close the previous appearance and object, if any*/
      fprintf(p->f,"})\n"); 
      fprintf(p->f,"(merge-ap sol_%s_%u {shading flat})\n",p->fileName,p->nobj);
      
      p->inObject=FALSE;
      SetColor3dObject(p->nobj,&(p->color),p);
    }
}

void SetColor3dObject(unsigned int nobj,Tcolor *c,Tplot3d *p)
{
  if (!p->inObject)
    {
      fprintf(p->f,"(merge-ap sol_%s_%u {material {kd 0.75 ",p->fileName,nobj);
      fprintf(p->f," ambient ");PrintColor(p->f,c);
      fprintf(p->f," *diffuse ");PrintColor(p->f,c);
      fprintf(p->f," edgecolor ");PrintColor(p->f,c);
      fprintf(p->f,"}})\n");
    }
}

void Delay3dObject(double t,Tplot3d *p)
{
  if (t!=0.0)
    fprintf(p->f,"(sleep-for %f)\n",t);
}

void Move3dObject(unsigned int nobj,THTransform *t,Tplot3d *p)
{ 
  fprintf(p->f,"(xform-set sol_%s_%u {",p->fileName,nobj);
  HTransformPrintT(p->f,t);
  fprintf(p->f,"})\n");
}

void PlotBox3d(double min_x,double max_x,
               double min_y,double max_y,
               double min_z,double max_z,
               Tplot3d *p)
{
  if (!p->inObject)
    Error("Adding geometry to a non-existing object");

  fprintf(p->f,"{OFF 8 6 1\n");
 
  /*vertices*/
  fprintf(p->f,"%f %f %f\n",min_x,min_y,min_z);
  fprintf(p->f,"%f %f %f\n",min_x,max_y,min_z);
  fprintf(p->f,"%f %f %f\n",max_x,max_y,min_z);
  fprintf(p->f,"%f %f %f\n",max_x,min_y,min_z);

  fprintf(p->f,"%f %f %f\n",min_x,min_y,max_z);
  fprintf(p->f,"%f %f %f\n",min_x,max_y,max_z);
  fprintf(p->f,"%f %f %f\n",max_x,max_y,max_z);
  fprintf(p->f,"%f %f %f\n",max_x,min_y,max_z);

  /*faces*/
  fprintf(p->f,"4 3 2 1 0\n");
  fprintf(p->f,"4 4 5 6 7\n");
  fprintf(p->f,"4 0 1 5 4\n"); 
  fprintf(p->f,"4 2 3 7 6\n"); 
  fprintf(p->f,"4 3 0 4 7\n");
  fprintf(p->f,"4 1 2 6 5}\n"); 

  fflush(p->f);
}

void PlotTriangle3d(double x1,double y1,double z1,
                    double x2,double y2,double z2,
                    double x3,double y3,double z3,
                    Tplot3d *p)
{
   if (!p->inObject)
    Error("Adding geometry to a non-existing object");

  /*vertices*/
  fprintf(p->f,"{OFF 3 1 1\n");
  fprintf(p->f,"%f %f %f\n",x1,y1,z1);
  fprintf(p->f,"%f %f %f\n",x2,y2,z2);
  fprintf(p->f,"%f %f %f\n",x3,y3,z3);
  
  /*faces*/
  fprintf(p->f,"3 0 1 2}\n"); 

  fflush(p->f);
}

void Plot3dObject(unsigned int nv,unsigned int nf,unsigned int ne,
                  double **v,
                  unsigned int *nvf,
                  unsigned int **fv,Tplot3d *p)
{
  unsigned int i,j;

  if (!p->inObject)
    Error("Adding geometry to a non-existing object");

  fprintf(p->f,"OFF\n");
  fprintf(p->f,"%u %u %u\n",nv,nf,ne);
  fprintf(p->f,"\n");
  for(i=0;i<nv;i++)
    fprintf(p->f,"%f %f %f\n",v[i][0],v[i][1],v[i][2]);

  fprintf(p->f,"\n");
  for(i=0;i<nf;i++)
    {
      fprintf(p->f,"%u   ",nvf[i]);
      for(j=0;j<nvf[i];j++)
        fprintf(p->f,"%u ",fv[i][j]);

      fprintf(p->f,"\n");
    }
  fprintf(p->f,"\n");
}

void PlotSphere(double r,double x,double y,double z,Tplot3d *p)
{
  if (!p->inObject)
    Error("Adding geometry to a non-existing object");
 
  fprintf(p->f,"{SPHERE %f %f %f %f}\n",r,x,y,z);
}

void PlotCylinder(double r,double *p1,double *p2,Tplot3d *p)
{
  THTransform f;

  if (!p->inObject)
    Error("Adding geometry to a non-existing object");

  HTransformX2Vect(r,r,p1,p2,&f);

  fprintf(p->f,"{INST\n");
  fprintf(p->f,"  geom       %s \n",CYLINDER_STRING);
  fprintf(p->f,"  transform ");
  HTransformPrintT(p->f,&f);
  fprintf(p->f,"}\n");
  
  HTransformDelete(&f);
}

void PlotVect3d(unsigned int n,
                double *x,double *y,double *z,
                Tplot3d *p)
{
  unsigned int i;

  if (!p->inObject)
    Error("Adding geometry to a non-existing object");
 
  fprintf(p->f,"{VECT 1 %u 1 \n   %u \n   1\n",n,n);

  for(i=0;i<n;i++)
    fprintf(p->f,"   %f %f %f\n",x[i],y[i],z[i]);

  PrintColor(p->f,&(p->color));
  fprintf(p->f," 0}\n"); /*alpha component is 0*/

  fflush(p->f);
  
}

void Take3dSnapshot(char *name,Tplot3d *p)
{
  /*fprintf(p->f,"(snapshot Camera \"| convert - %s.gif\")\n",name);*/
  fprintf(p->f,"(snapshot Camera \"%s.ppm\")\n",name);
}

void ClosePlot3d(boolean quit,double average_x,double average_y,double average_z,Tplot3d *p)
{
  unsigned int i;

  Close3dObject(p);

  /*translate the objects*/
  if ((average_x!=0.0)&&(average_y!=0.0)&&(average_z!=0.0))
    {
      /* Object 0 (reserved for the global axes, if any) are not plotted */
      for(i=1;i<=p->nobj;i++)
        fprintf(p->f,"(xform-set sol_%s_%u {1 0 0 0 0 1 0 0 0 0 1 0 %f %f %f 1})\n",p->fileName,i,average_x,average_y,average_z);
    }

  if (quit)
    fprintf(p->f,"(quit)\n");

  if (p->f!=stdout) fclose(p->f);
  
  free(p->fileName);
}