root/trunk/matml/src/ternary/geomview.c

/***************************************
  $Header$

  This file handles the Geomview interaction.
  ***************************************/


#include "gibbs.h" /*+ Ternary prototypes, typedefs, etc. +*/
#include <unistd.h>  /*+ For execlp() +*/
#include <stdio.h>   /*+ For FILE, fdopen(), etc. +*/


/*++++++++++++++++++++++++++++++++++++++
  This function launches a new Geomview process and returns a FILE pointer and
  optionally the Geomview version string.

  int GeomviewBegin It returns zero or an error code.

  FILE **geompipe Address of the FILE * object where GeomviewBegin will you can
  later fprintf commands to Geomview.  Your code should declare a FILE * and
  then pass its address to GeomviewBegin.

  char *version Character array of at least 100 characters where GeomviewBegin
  will put the Geomview version string (ignored if NULL).
  ++++++++++++++++++++++++++++++++++++++*/

int GeomviewBegin (FILE **geompipe, char *version)
{
  int ierr, gv_pid, gv_pipe[2], to_gv_pipe[2];

  pipe (gv_pipe); /* from geomview stdout */
  pipe (to_gv_pipe); /* to geomview stdin */
  gv_pid = fork ();

  if(gv_pid==0) /* child */
    {
      close (0);
      dup (to_gv_pipe[0]);
      close (to_gv_pipe[0]);
      close (to_gv_pipe[1]);
      close (1);
      dup (gv_pipe[1]);
      close (gv_pipe[0]);
      close (gv_pipe[1]);
      /* signal(SIGINT,SIG_IGN); */
      execlp (GEOMVIEW,GEOMVIEW,"-c","(interest (pick world))","-",NULL);
      perror (GEOMVIEW); /* only error gets here */
    }

  /* Program execution resumes here */
  close (gv_pipe[1]);
  close (to_gv_pipe[0]);
  *geompipe = fdopen (to_gv_pipe[1], "w"); /* hooked to stdin of geomview */

  if (version)
    {
      fprintf (*geompipe, "(echo (geomview-version) \"\n\")\n");
      fflush (*geompipe);
      read (gv_pipe[0], version, 99);
    }

  return 0;
}


/*++++++++++++++++++++++++++++++++++++++
  This sends a bunch of triangles to Geomview in COFF format.

  int GeomviewDisplayTriangleCOFF It returns zero or an error code.

  FILE *geompipe FILE pointer hooked to Geomview stdin.

  char *name Full name of this Geomview object.

  char *label Short label of this Geomview object.

  char *appearance Appearance tag for this Geomview object.

  int npoints Number of points.

  visual_point *points Point coordinates and colors.

  int ntriangles Number of triangles.

  int *vertices Vertex indices for each triangle.
++++++++++++++++++++++++++++++++++++++*/

int GeomviewDisplayTriangleCOFF
(FILE *geompipe, char *name, char *label, char *appearance,
 int npoints, visual_point *points, int ntriangles, int *vertices)
{
  int i;

  fprintf (geompipe, "(geometry \"%s\" { : %s })\n", name, label);
  fprintf (geompipe, "(read geometry { define %s \n", label);
  fprintf (geompipe, "appearance { %s }\nCOFF\n", appearance);
  fprintf (geompipe, "%d %d 0\n", npoints, ntriangles);

  for (i=0; i<npoints; i++)
    fprintf (geompipe, "%g %g %g %g %g %g %g\n", points[i].x, points[i].y,
             points[i].z, points[i].red, points[i].green, points[i].blue,
             points[i].alpha);

  for (i=0; i<ntriangles; i++)
    fprintf (geompipe, "3 %d %d %d\n", vertices[3*i], vertices[3*i+1],
             vertices[3*i+2]);

  fprintf (geompipe, "})\n");
  fflush (geompipe);

  return 0;
}


/*++++++++++++++++++++++++++++++++++++++
  This sends a bunch of facet objects to Geomview in COFF format.

  int GeomviewDisplayFacets It returns zero or an error code.

  FILE *geompipe FILE pointer hooked to Geomview stdin.

  char *name Full name of this Geomview object.

  char *label Short label of this Geomview object.

  char *appearance Appearance tag for this Geomview object.

  int npoints Number of points.

  visual_point *points Point coordinates and colors.

  int numfacets Number of facet objects to display.

  hull_facet *facets Array of facet objects.

  facet_type type Facet type to display, or ALL_TYPES to display all of them.
  ++++++++++++++++++++++++++++++++++++++*/

int GeomviewDisplayFacets
(FILE *geompipe, char *name, char *label, char *appearance,
 int npoints, visual_point *points, int numfacets, hull_facet *facets,
 facet_type type)
{
  int i, n=0;

  // Count number of relevant facets
  if (type == ALL_TYPES)
    n = numfacets;
  else
    for (i=0; i<numfacets; i++)
      if (facets [i].type == type)
        n++;

  fprintf (geompipe, "(geometry \"%s\" { : %s })\n", name, label);
  fprintf (geompipe, "(read geometry { define %s \n", label);
  fprintf (geompipe, "appearance { %s }\nCOFF\n", appearance);
  fprintf (geompipe, "%d %d 0\n", npoints, n);

  for (i=0; i<npoints; i++)
    fprintf (geompipe, "%g %g %g %g %g %g %g\n", points[i].x, points[i].y,
             points[i].z, points[i].red, points[i].green, points[i].blue,
             points[i].alpha);

  for (i=0; i<numfacets; i++)
    if (type == ALL_TYPES || type == facets[i].type)
      fprintf (geompipe, "3 %d %d %d\n", facets[i].vertex[0], facets[i].vertex[1], facets[i].vertex[2]);

  fprintf (geompipe, "})\n");
  fflush (geompipe);

  return 0;
}


/*++++++++++++++++++++++++++++++++++++++
  Displays a phase boundary.  Note this pipes all of the points back to
  geomview again.

  int GeomviewDisplayPhaseBoundary It returns zero or an error code.

  FILE *geompipe FILE pointer hooked to Geomview stdin.

  char *name Full name of this Geomview object.

  char *label Short label of this Geomview object.

  char *appearance Appearance tag for this Geomview object.

  int npoints Number of points.

  visual_point *points Point coordinates and colors.

  phase_boundary theboundary Phase boundary object to display, or NULL to
  display nothing.
  ++++++++++++++++++++++++++++++++++++++*/

int GeomviewDisplayPhaseBoundary
(FILE *geompipe, char *name, char *label, char *appearance,
 int npoints, visual_point *points, phase_boundary *theboundary)
{
  int i, n_edges=0, c;

  if (theboundary)
    {
      n_edges = theboundary->n_edges;
      c = theboundary->compos;
    }
  else
    npoints = 0;

  fprintf (geompipe, "(geometry \"%s\" { : %s })\n", name, label);
  fprintf (geompipe, "(read geometry { define %s \n", label);
  fprintf (geompipe, "appearance { %s }\nCOFF\n", appearance);
  fprintf (geompipe, "%d %d 0\n", npoints, n_edges);

#ifdef DEBUG
  printf ("(geometry \"%s\" { : %s })\n", name, label);
  printf ("(read geometry { define %s \n", label);
  printf ("appearance { %s }\nCOFF\n", appearance);
  printf ("%d %d 0\n", npoints, n_edges);
#endif

  for (i=0; i<npoints; i++)
    {
      fprintf (geompipe, "%g %g %g %g %g %g %g\n", points[i].x, points[i].y,
               points[i].z, points[i].red, points[i].green, points[i].blue,
               points[i].alpha);
#ifdef DEBUG
      printf ("%g %g %g %g %g %g %g\n", points[i].x, points[i].y,
              points[i].z, points[i].red, points[i].green, points[i].blue,
              points[i].alpha);
#endif
    }
  for (i=0; i<n_edges; i++)
    {
      fprintf (geompipe, "3 %d %d %d\n", theboundary->edges [c*i],
               (c>1) ? theboundary->edges [c*i+1] : theboundary->edges [c*i],
               (c>2) ? theboundary->edges [c*i+2] : theboundary->edges [c*i]);
#ifdef DEBUG
      printf ("3 %d %d %d\n", theboundary->edges [c*i],
              (c>1) ? theboundary->edges [c*i+1] : theboundary->edges [c*i],
              (c>2) ? theboundary->edges [c*i+2] : theboundary->edges [c*i]);
#endif
    }
  fprintf (geompipe, "})\n");
#ifdef DEBUG
  printf ("})\n");
#endif
  fflush (geompipe);

  return 0;
}


/*++++++++++++++++++++++++++++++++++++++
  GeomviewEnd exits geomview, closes its FILE, and sets its FILE to NULL so
  nobody tries to use it again.

  int GeomviewEnd It returns zero or an error code.

  FILE **geompipe Address of the FILE * object originally opened by
  GeomviewBegin().
  ++++++++++++++++++++++++++++++++++++++*/

int GeomviewEnd (FILE **geompipe)
{
  fprintf (*geompipe, "(exit)");
  fclose (*geompipe);
  *geompipe = NULL;
  return 0;
}