amroc/converters/HDFToV3/V3Visualizer.h

// -*- C++ -*-

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

#ifndef AMROC_V3VISUALIZER_H
#define AMROC_V3VISUALIZER_H

#include <wsdepend.h>

#include "HDFToV3/V3VisualGrid.h"
#include "HDFToV3/FileEventhandler.h"  

extern "C" {
 void V3SURFACE(int nsurf[][2], int surf[], int scel[][4], char tsurf[][20]);
 void V3GRID(float xyz[][3]);
 void V3SCAL(int *key, float v[]);
 void V3DRAW3D(int *istat);
 void V3EQUIV(int listequiv[][2]);
 void V3VECT (int *key, float v[][3]);
 void V3PEVENTS(int *iwin, int *type, int *xe, int *ye, int *state);
}              
                  

V3VisualGridBase *Grid;
FileEventHandler FileEvents;


class V3Visualizer {
public:
  V3Visualizer() {}
  
  void Show(int step, std::string* solvername, std::string* inputname,
            std::string** filenames, int Nfilenames, V3VisualGridBase* gr) {

    remove(".Visual3.struc");
    Grid = gr;
    std::string sname;
    if (solvername) sname = *solvername;
    else sname = "solver.in";
    std::ifstream infile;
    infile.open(sname.c_str(), std::ios::in);
    if (!infile) {
      std::cerr << sname << " not found. Aborting." << std::endl;
      std::exit(-1);
    } 
    infile.close();

    ControlDevice DisplayCtrl(GetFileControlDevice(sname.c_str(),""));
    Grid->register_at(DisplayCtrl,"");
    DisplayCtrl.update();     
    Grid->SetUpdateName(inputname);
    Grid->update();  

    if (Grid->SetUp(step,filenames,Nfilenames) < 0) {
      std::cerr << "Nothing read in. Aborting." << std::endl;
      std::exit(-1);
    } 
    std::cerr << Grid->Size() << " Grids constructed from all levels."
              << std::endl;  
    FileEvents.SetStep(step);
    
    // The important declarations needed for Visual3;
    
    // Generating standard spec.col if file is missing; 
    std::ofstream outfile;
    int i;
    outfile.precision(4);
    infile.open("spec.col", std::ios::in);
    if (!infile) {
      outfile.open("spec.col", std::ios::out);
      outfile << 256 << "    " << 2 << std::endl;
      for(i=0;i<64;i++){
        outfile << "0.0000    " << i/64.0 << "    1.0000" << std::endl;
      }

      for(i=0;i<64;i++){
        outfile << "0.0000    1.0000    " << 1-i/64.0 << std::endl;
      }

      for(i=0;i<64;i++){
        outfile << i/64.0 << "    1.0000    0.0000" << std::endl;
      }
      for(i=1;i<65;i++){
        outfile << "1.0000    " << 1-i/64.0 << "    0.0000" << std::endl;
      }
      outfile << "0.0000    0.0000    0.0000" << std::endl;
      outfile << "1.0000    1.0000    1.0000" << std::endl;

      outfile.close();
    }
    else{
      infile.close();
    }
    
    
    char cmf[] = "spec.col";
    int len_cmf = std::strlen(cmf); 
    
    char titl[LEN_TKEYS];
    int len_titl = strcpn(titl, Grid->Title());

    int win3d = 1; 
    int num_keys = Grid->NKeys() + (Grid->Equiv() ? 1000 : 0);

    std::cerr << "Limits:" << std::endl;

    int lus = 99;    
    int iopt = 0; // steady grid and steady data   
    int mirror = 0;// no miyrroring  
    int knode = 0;// number of non-block nodes  
    int kcel1 = 0;// no tetrahedras  
    int kcel2 = 0; // no pyramids 
    int kcel3 = 0; // no prisms
    int kcel4 = 0;// hexahedras 
    int knptet = 0;// no poly-tetrahedral strips
    int kptet = 0;// no cells in any polytetrahedra
    int knblock = Grid->Size();// structured blocks
  
    int *blocks;          
    blocks = new int[3*knblock];// structured blocks definitions
    Grid->SetBlocks(blocks);    
    
    if (Grid->Equiv()) 
      Grid->FindNodePairs(); 
    int kequiv = (Grid->Equiv() ? Grid->Kequiv() : 0);
    
    std::cerr << " Found "  << kequiv << " matching node pairs." << std::endl;
    
    std::cerr << "my calculated ksurf: " << Grid->NFacetsOnSurface()
              << std::endl;
    int ksurf = Grid->NFacetsOnSurface(); // number of domain surfaces
    int knsurf = 6; // number of domain surface groups
    
    
    // call Visual3 
    // note string lengths are appended onto the stack as FORTRAN
    // expects them to be there.
    
    V3_INIT(titl, &iopt, cmf, &lus, &win3d, &num_keys, Grid->KeyboardKeys(),
            Grid->TitleKeys(), Grid->FunctionsKeys(), Grid->Mvals(), &mirror,
            &knode, &kequiv, &kcel1, &kcel2, &kcel3, &kcel4, &knptet, &kptet,
            &knblock, blocks, &ksurf, &knsurf, len_titl, len_cmf, LEN_TKEYS);
 
    delete [] blocks; 
  }
};


void V3SURFACE (int nsurf[][2], int surf[], int scel[][4], char tsurf[][20])
{ 
  Grid->SetSurfaces(nsurf, surf, scel, tsurf);
}
  
void V3GRID(float xyz[][3])
{ 
  Grid->SetGrid(xyz);  
}
  
void V3SCAL (int *key, float v[])
{
  Grid->SetScal(key,v);
}

void V3VECT (int *key, float v[][3])
{
  Grid->Vect(key,v);
} 
  
void V3DRAW3D(int* istat)
{
  if (*istat == 0)  
    Grid->Draw3D();            
}

void V3EQUIV(int listequiv[][2])
{  
  Grid->SetNodePairs(listequiv);
} 
        
void V3PEVENTS(int *iwin, int *type, int *xe, int *ye, int *state) {
  FileEvents.ProcessEvent(iwin, type, xe, ye, state);
}

#endif

---