amroc/converters/HDFToV3/FileEventhandler.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2002 Ralf Deiterding, Matthias Kraatz
// Brandenburgische Universitaet Cottbus
//
// Copyright (C) 2003-2007 California Institute of Technology
// Ralf Deiterding, ralf@cacr.caltech.edu

#ifndef AMROC_FILEEVENTHANDLER_H
#define AMROC_FILEEVENTHANDLER_H

#include "HDFToV3/Eventhandler.h"

#include <iostream>
#include <fstream>
#include <sstream>
#include <cstdlib>   
#include <cstdio>   
#include <cmath>
#include <cctype>


class FileEventHandler : public VisualEventHandler {
public:
  int step;
  FileEventHandler() : VisualEventHandler() {}

  void SetUpQueue() { 

    std::ifstream infile;
    char charBuffer[256];
    std::string arg,fileID,stepstr;
    int endOfBlock=0;
    std::stringstream temp;
    
    
    /* connect to the X server */
    Display *display = XOpenDisplay (NULL);      
    if (display == NULL) {
      std::fprintf (stderr, "cannot connect to server\n");
      std::exit (EXIT_FAILURE);
    }    
    
    infile.open("display.in", std::ios::in);
    
    if(SeekArg(infile,"SCRIPTBEGIN")){
      SeekArg(infile,"{");
      std::cerr << "Processing general script commands..." << std::endl; 
      endOfBlock=0;
      do{
        arg=GetNextArg(infile,charBuffer);

        std::cerr << arg << std::endl;
        if(charBuffer[0]=='#')
          SeekEol(infile);
        ProcessArg(infile,arg,display);
        if(arg == "}" || arg == "END")
          endOfBlock=1;
      }while(!endOfBlock);
      std::cerr << "Done with general script." << std::endl;
    }

    fileID="FILEID_";
    temp << step;
    temp >> stepstr;
    fileID += stepstr;
    std::cerr << ">" << fileID << "<" << std::endl;
    if(SeekArg(infile,fileID)){
      SeekArg(infile,"{");
      std::cerr << "Processing script commands for step " << step << "..."
                << std::endl; 
      endOfBlock=0;
      do{
        arg=GetNextArg(infile,charBuffer);
        if(charBuffer[0]=='#')
          SeekEol(infile);
        ProcessArg(infile,arg,display);
        if(arg == "}" || arg == "END")
          endOfBlock=1;
      }while(!endOfBlock);
      std::cerr << "Done with script for step " << step << "..." << std::endl; 
    }
    infile.close();
  };
  void SetStep(int stepID){step=stepID;};


private:  
  std::string GetNextArg(std::ifstream& infile,char* charBuffer){
    int i=0;
    std::string str;
    char c;
    do{
      c=infile.get();
    } while((!infile.eof()) && (c<=' '|| c==',' || c==';'));
    charBuffer[i++] = std::toupper(c);
    do{
      c=infile.get();
      charBuffer[i++] = std::toupper(c);      
    } while ((!infile.eof()) && c>' ' && c!=',' && c!=';');
    if(infile.eof())
      return std::string("END");
    charBuffer[i-1]='\0';
    str=charBuffer;
    return str;
  }
  
  int SeekArg(std::ifstream& infile,std::string str){
    char charBuffer[256];
    int exists=0,endOfFile=0;
    std::string cmd;
    do{
      cmd=GetNextArg(infile,charBuffer);
      if(cmd == str)
        exists=1;
      if(cmd == "END" || infile.eof())
        endOfFile=1;      
    }while(!exists && !endOfFile);
    if(!endOfFile)
      return 1;
    else
      return 0;
    
  }
  void SeekEol(std::ifstream& infile){
    char c;
    do{
      c=infile.get();
    } while((!infile.eof()) && c!='\n');
  }
    
  void ProcessArg(std::ifstream& infile,std::string cmd,Display *display){
    int opt,i,j,k;
    char charBuffer[256];
    long dummy=0,iv[12],winid[8],keybin[32];     
    char blank=' ';
    float rdummy,rv[12],matA[9],matB[9];
    float Pi=4*std::atan(float(1.0)),phi[3];
    
    opt=0;
    V3_GETSTATE(&opt,winid,rv,&blank,1);
    opt = -1;
    V3_GETSTATE(&opt,keybin,rv,&blank,1);
    if(cmd == "SETVIEW"){
      for(i=0;i<12;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> rv[i]; 
      }
      EventQueue().push(SingleEvent(1, 0, 0, 0, 
                                    new SetState(1, &dummy, 1, rv, 12, &blank, 1)));
    }  
    
    if(cmd == "SETCUT"){
      for(i=0;i<9;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> rv[i]; 
      }
      EventQueue().push(SingleEvent(1, 0, 0, 0, 
                                    new SetState(2, &dummy, 1, rv, 12, &blank, 1)));
    }  

    if(cmd == "SCALARSTATE"){
      std::cerr << "Setting ScalarState" << std::endl;
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
      arg >> iv[0];
      for(i=0;i<2;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> rv[i] ;
      }
      
      
      EventQueue().push(SingleEvent(1, 0, 0, 0, 
                                    new SetState(5, iv, 1, rv, 2, &blank, 1)));
    }
    
    if(cmd == "CELLSCALARSTATE"){
      std::cerr << "Setting CellScalarState" << std::endl;
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
      arg >> iv[0];
      for(i=0;i<2;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> rv[i] ;
      }
      EventQueue().push(SingleEvent(1, 0, 0, 0, 
                                    new SetState(9, iv, 1, rv, 2, &blank, 1)));
    }
    
    if(cmd == "THRESHSTATE"){
      for(i=0;i<2;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> iv[i]; 
      }
      for(i=0;i<2;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> rv[i]; 
      }
      EventQueue().push(SingleEvent(1, 0, 0, 0, 
                                    new SetState(7, iv, 2, rv, 2, &blank, 1)));
    }
    
    if(cmd == "ISOSURFACE"){
      EventQueue().push(SingleEvent(1, winid[2], 2, keybin[6]));
      opt=5;
      
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
      arg >> iv[0];
      V3_GETSTATE(&opt,iv,rv,&blank,1);
      
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg1(std::string(charBuffer,sizeof(charBuffer)-1));
      arg1 >> rv[2];
      rdummy=(rv[2]-rv[0])/(rv[1]-rv[0])*436.0+37.5;
      dummy=(int)rdummy;
      std::cerr << "Isosurface   " << rv[0] << "    " << rv[2] << "    "
                << rv[1] << "        " << dummy << std::endl;

      EventQueue().push(SingleEvent(1,winid[4],4,3,dummy,50));
    }
    
    if(cmd == "SURFSTATE"){
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
      arg >> iv[0];
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg2(std::string(charBuffer,sizeof(charBuffer)-1));
      arg2 >> iv[1];
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg3(std::string(charBuffer,sizeof(charBuffer)-1));
      arg3 >> iv[2];
      iv[1]+=2*iv[2];
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg4(std::string(charBuffer,sizeof(charBuffer)-1));
      arg4 >> iv[2];
      iv[1]+=4*iv[2];
      cmd=GetNextArg(infile,charBuffer);
      std::istringstream arg5(std::string(charBuffer,sizeof(charBuffer)-1));
      arg5 >> iv[2];
      iv[1]+=8*iv[2];
      EventQueue().push(SingleEvent(1, 0, 0, 0, 
                                    new SetState(10, iv, 2, &rdummy, 1, &blank, 1)));
    }
    
    if(cmd == "TRANSLATION"){
      
      
      opt=1;
      V3_GETSTATE(&opt,iv,rv,&blank,0);
      
      for(i=0;i<3;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> rdummy;
        rv[i*4+3]+=rdummy;
      }
      EventQueue().push(SingleEvent(1, 0, 0, 0, 
                                    new SetState(1, &dummy, 1, rv, 12, &blank, 1)));
    }
    
    if(cmd == "ROTATION"){
      opt=1;
      V3_GETSTATE(&opt,iv,rv,&blank,1);
      for(i=0;i<3;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> phi[i];
        phi[i]*=Pi/180;
      }
      
      matA[0]=std::cos(phi[1])*std::cos(phi[2]);
      matA[1]=std::cos(phi[1])*std::sin(phi[2]);
      matA[2]=std::sin(phi[1]);
      matA[3]=-std::sin(phi[0])*std::sin(phi[1])*std::cos(phi[2])
              -std::cos(phi[0])*std::sin(phi[2]);
      matA[4]=-std::sin(phi[0])*std::sin(phi[1])*std::sin(phi[2])
              +std::cos(phi[0])*std::cos(phi[2]);
      matA[5]=std::sin(phi[0])*std::cos(phi[1]);
      matA[6]=-std::cos(phi[0])*std::sin(phi[1])*std::cos(phi[2])
              +std::sin(phi[0])*std::sin(phi[2]);
      matA[7]=-std::cos(phi[0])*std::sin(phi[1])*std::sin(phi[2])
              -std::sin(phi[0])*std::cos(phi[2]);
      matA[8]=std::cos(phi[0])*std::cos(phi[1]);
      for(i=0;i<3;i++){
        for(j=0;j<3;j++){
          matB[i*3+j]=0;
          for(k=0;k<3;k++){
            matB[i*3+j]+=matA[i*3+k]*rv[k*4+j];
          }
        }
      }
      for(i=0;i<3;i++){
        for(j=0;j<3;j++){
          rv[i*4+j]=matB[i*3+j];
        }
      }
      EventQueue().push(SingleEvent(1, 0, 0, 0, 
                                    new SetState(1, &dummy, 1, rv, 12, &blank, 1)));
      
    }
    if(cmd == "2DVIEW"){

      opt=2; 
      for(i=0;i<9;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> rv[i]; 
      }
      EventQueue().push(SingleEvent(1, 0, 0, 0,  
                                      new SetState(2, iv, 1, rv, 9, &blank, 1))); 
      
      EventQueue().push(SingleEvent(1, winid[2], 2, keybin[2] ));
      EventQueue().push(SingleEvent(1, winid[2], 2, keybin[2] ));
    }

    if(cmd == "EDGEOUTLINEOFF"){
      opt=14;
      V3_GETSTATE(&opt,iv,rv,&blank,0);
      iv[0]=0;
      EventQueue().push(SingleEvent(1, 0, 0, 0,  
                                      new SetState(14, iv, 2, rv, 7, &blank, 1))); 
    }
    if(cmd == "EDGEOUTLINEON"){
      opt=14;
      V3_GETSTATE(&opt,iv,rv,&blank,0);
      iv[0]=1;
      EventQueue().push(SingleEvent(1, 0, 0, 0,  
                                      new SetState(14, iv, 2, rv, 7, &blank, 1))); 
    }

    if(cmd == "TRANSPARENTPLANEOFF"){
      opt=14;
      V3_GETSTATE(&opt,iv,rv,&blank,0);
      iv[1]=0;
      EventQueue().push(SingleEvent(1, 0, 0, 0,  
                                      new SetState(14, iv, 2, rv, 7, &blank, 1))); 
    }
    if(cmd == "TRANSPARENTPLANEON"){
      opt=14;
      V3_GETSTATE(&opt,iv,rv,&blank,0);
      iv[1]=1;
      EventQueue().push(SingleEvent(1, 0, 0, 0,  
                                      new SetState(14, iv, 2, rv, 7, &blank, 1))); 
    }
    if(cmd == "EDGESOFFTRANSPARENTON"){
      opt=14;
      V3_GETSTATE(&opt,iv,rv,&blank,0);
      iv[0]=0;
      iv[1]=0;
      EventQueue().push(SingleEvent(1, 0, 0, 0,  
                                      new SetState(14, iv, 2, rv, 7, &blank, 1))); 
    }

    if(cmd == "HARDCOPY"){
      EventQueue().push(SingleEvent(5, winid[2], 2, '~'));
    }
    if(cmd == "2DHARDCOPY"){
      EventQueue().push(SingleEvent(5, winid[1], 2, '~'));  
    }
    if(cmd == "ESCAPE"){
      EventQueue().push(SingleEvent(5, winid[2], 2, keybin[19]));
    }
    if(cmd == "DISPLAYSIZE"){
      for(i=0;i<2;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> iv[i];
      }
      XResizeWindow(display,winid[iv[0]],iv[1],iv[1]);
      XFlush (display);  
    }
    if(cmd == "MOVEWINDOW"){
      for(i=0;i<3;i++){
        cmd=GetNextArg(infile,charBuffer);
        std::istringstream arg(std::string(charBuffer,sizeof(charBuffer)-1));
        arg >> iv[i];
      }
      XMoveWindow(display,winid[iv[0]],iv[1],iv[2]);
      XFlush (display);  
    }
    
    if(cmd == "SETLIMITS"){
      EventQueue().push(SingleEvent(1, winid[4], 2, 'f'));      
    }
    if(cmd == "SETCONTOURS"){
      EventQueue().push(SingleEvent(1, winid[4], 2, 'c'));   
    }
    if(cmd == "PLANARCUTOFF"){
      EventQueue().push(SingleEvent(1, winid[2], 2, keybin[1]));     
    }
    if(cmd == "PLANARCUTON"){
      EventQueue().push(SingleEvent(1, winid[2], 2, keybin[2]));     
      EventQueue().push(SingleEvent(1, winid[2], 2, keybin[2]));     
    }
    if(cmd == "CUTPLANE"){
      EventQueue().push(SingleEvent(1, winid[2], 2, keybin[3]));  
    }
    if(cmd == "CONTOURTOGGLE"){
      EventQueue().push(SingleEvent(1, winid[1], 2, keybin[7]));
      
    }
    if(cmd == "RENDERTOGGLE"){
      EventQueue().push(SingleEvent(1, winid[1], 2, keybin[8]));
    }
    if(cmd == ""){}
  }
  
};

#endif

---