rocin_block.h

Go to the documentation of this file.

/* *******************************************************************
 * Rocstar Simulation Suite                                          *
 * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
 *                                                                   *
 * Illinois Rocstar LLC                                              *
 * Champaign, IL                                                     *
 * www.illinoisrocstar.com                                           *
 * sales@illinoisrocstar.com                                         *
 *                                                                   *
 * License: See LICENSE file in top level of distribution package or *
 * http://opensource.org/licenses/NCSA                               *
 *********************************************************************/
/* *******************************************************************
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
 * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
 * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
 * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
 * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
 *********************************************************************/
// $Id: rocin_block.h,v 1.7 2008/12/06 08:43:19 mtcampbe Exp $

#ifndef _ROCIN_BLOCK_H_
#define _ROCIN_BLOCK_H_

#include <map>

struct VarInfo_HDF4 {
  VarInfo_HDF4(const std::string &name, char pos, COM_Type dType, 
               const std::string& units, int nc, int32 i, int nitems, 
               int ng, bool is_null)
  : m_name(name), m_position(pos), m_dataType(dType), m_units(units), 
    m_indices(nc, i), m_nitems( nitems), m_ng(ng), m_is_null(nc,is_null)
  {}

  std::string m_name;            
  char m_position;              
  COM_Type m_dataType;          
  std::string m_units;          
  std::vector<int32> m_indices; 
  int m_nitems;                 
  int m_ng;                     
  std::vector<bool>  m_is_null; 
};

struct GridInfo_HDF4 {
   inline GridInfo_HDF4(int32* size, int ng)
   : m_name(":st0:"),
     m_numElements(size[2] - 1
                   * std::max(int32(1), (size[1] - 1))
                   * std::max(int32(1), (size[0] - 1))),
     m_numGhostElements(ng), m_index(FAIL)
   { m_name[3] += (size[1] > 1 ? (size[0] > 1 ? 3 : 2) : 1);
     m_size[0] = size[0]; m_size[1] = size[1]; m_size[2] = size[2]; }

   inline GridInfo_HDF4(std::string name, int ne, int ng, int32 i)
   : m_name(name), m_numElements(ne), m_numGhostElements(ng), m_index(i)
   { m_size[0] = 0; m_size[1] = 0; m_size[2] = 0; }

   int32 m_size[3];        
   std::string m_name;     
   int m_numElements;      
   int m_numGhostElements; 
   int32 m_index;          
};

struct Block_HDF4 {
   inline Block_HDF4(const std::string& file, const std::string& geomFile,
                     int32 indices[], int paneId, const std::string& time,
                     const std::string& units, int numNodes, int ghostNodes)
   : m_file(file), m_geomFile(geomFile), m_paneId(paneId), time_level(time),
     m_units(units), m_numNodes(numNodes), m_numGhostNodes(ghostNodes)
   { m_indices[0] = indices[0]; m_indices[1] = indices[1];
     m_indices[2] = indices[2]; }

   std::string m_file;     
   std::string m_geomFile; 
   int32 m_indices[3];     
   int m_paneId;           
   std::string time_level; 
   std::string m_units;    
   int m_numNodes;         
   int m_numGhostNodes;    
   std::vector<GridInfo_HDF4> m_gridInfo; 
   std::vector<VarInfo_HDF4>  m_variables; 
};
typedef std::multimap<std::string, Block_HDF4*> BlockMM_HDF4;

#ifdef USE_CGNS

struct VarInfo_CGNS {
  VarInfo_CGNS(const std::string &name, char pos, COM_Type dType, 
               const std::string& units, int nc, int i, int nitems, 
               int ng, bool is_null)
  : m_name(name), m_position(pos), m_dataType(dType), m_units(units), 
    m_indices(nc, i), m_nitems(nitems), m_ng(ng), m_is_null(nc,is_null)
  {}

  std::string m_name;          
  char m_position;             
  COM_Type m_dataType;         
  std::string m_units;         
  std::vector<int> m_indices;  
  int m_nitems;                
  int m_ng;                    
  std::vector<bool> m_is_null; 
};

struct GridInfo_CGNS {
   inline GridInfo_CGNS(int* size, int cellDim, int ng)
   : m_name(":st0:"), m_numElements(size[cellDim]),
     m_numGhostElements(ng)
   { m_name[3] += cellDim; m_size[0] = size[0]; m_size[1] = m_size[2] = 1;
     if (cellDim > 1) { m_size[1] = size[1]; m_numElements *= size[cellDim+1];
       if (cellDim > 2) { m_size[2] = size[2]; m_numElements *= size[5]; }
   } }

   inline GridInfo_CGNS(std::string name, int ne, int ng)
   : m_name(name), m_numElements(ne), m_numGhostElements(ng)
   { m_size[0] = 0; m_size[1] = 0; m_size[2] = 0; }

   int m_size[3];          
   std::string m_name;     
   int m_numElements;      
   int m_numGhostElements; 
   // We don't really need to save the Element index, since we use all of
   // them in the Zone.
};

struct Block_CGNS {
   inline Block_CGNS(const std::string& file, int B, int Z, int G, int paneId,
                     const std::string& time, const std::string& units,
                     int numNodes, int ghostNodes)
   : m_file(file), m_B(B), m_Z(Z), m_G(G), m_W(0), m_P(0), m_C(0), m_E(0),
     m_N(0), m_paneId(paneId), time_level(time), m_units(units),
     m_numNodes(numNodes), m_numGhostNodes(ghostNodes)
   {}

   std::string m_file;     
   int m_B;                
   int m_Z;                
   int m_G;                
   int m_W;                
   int m_P;                
   int m_C;                
   int m_E;                
   int m_N;                
   int m_paneId;           
   std::string time_level; 
   std::string m_units;    
   int m_numNodes;         
   int m_numGhostNodes;    
   std::vector<GridInfo_CGNS> m_gridInfo; 
   std::vector<VarInfo_CGNS>  m_variables; 
};
typedef std::multimap<std::string, Block_CGNS*> BlockMM_CGNS;

#endif // USE_CGNS

#endif