---

View.C

Go to the documentation of this file.

/*
 * Copyright (C) 2003 Eric Bohm
 *
 *This program is free software; you can redistribute it and/or modify
 *it under the terms of the GNU General Public License as published by
 *the Free Software Foundation version 2 of the License.
 * 
 *This program is distributed in the hope that it will be useful, but
 *WITHOUT ANY WARRANTY; without even the implied warranty of
 *MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 *General Public License for more details.
 * 
 *You should have received a copy of the GNU General Public License along
 *with this program; if not, write to the Free Software Foundation, Inc.,
 *675 Mass Ave, Cambridge, MA 02139, USA.
 *
 */
using namespace std;

#include <GL/glu.h>
#include <GL/glut.h>
#include <stdio.h>
#include <vector>
#include <algorithm>
#include "View.h"



bool View::loadfile(const char *filename)
{
  unsigned int time=12;
  unsigned int chunk=12;
  char base[80];
  char stuff[80];
  char ext[80];
  memset(base,'\0',80);
  // strip path crap to eliminate .
  char *slash=rindex(filename,'/');
  if(slash>=filename)
    {
      strcpy(base,slash+1);
      slash=rindex(filename,'.');
      strcpy(ext,slash+1);
    }
  else
    {
     strcpy(base,filename);
     slash=rindex(filename,'.');
     strcpy(ext,slash+1);
    }
  sscanf(base,"%[^0123456789.]",stuff);
  chunk=atoi(base+strlen(stuff));
  time=atoi(ext);
  map <unsigned int, RGB>::iterator chiter=chunks.find(chunk);
  map <unsigned int, Timestep>::iterator timeiter=timesteps.find(time);
  if(timeiter==timesteps.end())
    {
      timesteps[time]=Timestep();
    }
  else
    {
    }
  if(chiter==chunks.end())
    { 
      GLfloat tetcolor[3];
      tetcolor[0] = drand48() / 2.0 + 0.5;
      tetcolor[1] = drand48() / 2.0 + 0.5;
      tetcolor[2] = drand48() / 2.0 + 0.5;
      chunks[chunk]=RGB(tetcolor);
      if(timesteps[time].loadfile(filename,tetcolor,time,chunk))
        {
          filenum++;
        }
      else
        return false;
    }
  else
    {
      if(timesteps[time].loadfile(filename,chiter->second.rgb,time,chunk))
        {
          filenum++;
        }
      else
        return false;
    }
  if(time>uptodate_step.size())
    {
      uptodate_step.resize(time);
      uptodate_step[time]=false;
    }
  return true;
}




void View::display_step(unsigned int step)
{
  map <unsigned int, Timestep>::iterator stepiter=timesteps.find(step);
  if(stepiter!=timesteps.end())
    {
      currstep=step;
      now=stepiter;
      if(glIsList(step) && uptodate_step[step])
        { // if it exists and is up to date just redisplay
          glCallList(step);
        }
      else
        { // draw it
          uptodate_step[step]=true;
          now->second.draw_polygons(ScreenWidth,ScreenHeight,colormap,aspect_ratio,scale,xtri,ytri,ztri,xtran,ytran,ztran,
                             fillpoly,showlabel,showlength,opaque,pointmode,solo_chunk);
        }
    }
  else
    {
      fprintf(stderr,"Error no such timestep %u\n",step);
    }
}

void View::recenter()
{
  // In a better world we would use for_each here.
  // But that doesn't work nicely when you need to modify the object.
  map < unsigned int, Timestep>::iterator titer=timesteps.begin();
  while(titer!=timesteps.end())
    {
      titer->second.set_random_color(points);
      titer->second.calc_and_move();
      titer++;
    }
}


void View::make_random_color_set()
{
  //simple solution here.  cat all the point vectors together.
  // sort the collection.
  // unique_copy the collection
  // assign a random color to each point
  vector <Coordinate> collection;
  map < unsigned int, Timestep>::iterator titer=timesteps.begin();
  while(titer!=timesteps.end())
    {
      vector <Coordinate> tstamp=titer->second.get_points();
      collection.insert(collection.end(),tstamp.begin(),tstamp.end());
      titer++;
    }
  sort(collection.begin(),collection.end());
  insert_iterator< vector <Coordinate> > ins(points,points.begin());
  //  unique(collection.begin(),collection.end());
  //  points=collection;
  unique_copy(collection.begin(),collection.end(),ins);
  //  for_each(collection.begin(),collection.end(),mem_fun_ref(&Coordinate::edump));

  for(unsigned int i=0;i<points.size();i++)
    {
      points[i].set_rgb(drand48() / 2.0 + 0.5,drand48() / 2.0 + 0.5,drand48() / 2.0 + 0.5);
    }

}

---