gmHdf5ModelReader.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together 
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND, 
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR 
** A PARTICULAR PURPOSE.
**
************************************************************************/
 
#ifndef _GEMA_HDF5_MODEL_READER_H_
#define _GEMA_HDF5_MODEL_READER_H_
 
#include "gmLog.h"
#include "gmHdf5Model.h"
#include "gmBoundaryCondition.h"
#include "gmCellBoundary.h"
#include "gmDiscontinuitySet.h"
#include "gmNanBoxedValue.h"
 
#include <QStringList>
 
class GmValueAccessor;
class GmCellAccessor;
class GmGaussAccessor;
class GmMeshGaussIndex;
class GmBoundaryConditionAccessor;
class GmPCR;
 
namespace H5 {
  class H5File;
  class DataType;
}
 
class GMC_API_EXPORT GmHdf5ModelReader : public GmHdf5Model
{
public:
 
  GmHdf5ModelReader(QString fileName, const GmLogCategory& logger);
  ~GmHdf5ModelReader();
 
  struct DicInfo
  {
    QString _name;  
    int     _n;     
  };
 
  struct AttrInfo
  {
    QString _name;      
    int     _state;     
    int     _n;         
    int     _d;         
    bool    _functions; 
  };
 
  struct BcInfo
  {
    QString                     _name;      
    GmBoundaryCondition::BcType _type;      
    int                         _n;         
    int                         _d;         
    QVector<AttrInfo>           _attrData;  
  };
 
  struct DiscInfo
  {
    QString _name;  
    int     _n1;    
    int     _n2;    
  };
 
  struct DsInfo
  {
    QString                               _name;     
    GmDiscontinuitySet::DiscontinuityType _type;     
    QVector<DiscInfo>                     _discData; 
  };
 
  struct MeshInfo
  {
    int _d;               
    int _nGeometryNodes;  
    int _nGhostNodes;     
    int _nCells;          
    int _nPSets;          
    int _nCellNodes;      
    int _nGhostCells;     
    int _nInactiveCells;  
 
    QStringList      _functionsDictionary;     
    QStringList      _propertySetsDictionary;  
    QVector<DicInfo> _cellGroupsDictionary;    
    QVector<DicInfo> _nodeSetsDictionary;      
    QVector<DicInfo> _boundaryEdgesDictionary; 
    QVector<DicInfo> _boundaryFacesDictionary; 
 
    QVector<AttrInfo> _nodeInitialData;  
    QVector<AttrInfo> _cellInitialData;  
    QVector<AttrInfo> _gaussInitialData; 
 
    QVector<BcInfo>   _bcList; // Information on each associated boundary condition
    QVector<DsInfo>   _dsList; // Information on each associated discontinuity set
 
    MeshInfo() 
      : _d(0), _nGeometryNodes(0), _nGhostNodes(0), _nCells(0), _nPSets(0), 
        _nCellNodes(0), _nGhostCells(0), _nInactiveCells(0) {}
  };
 
  bool open();
 
  QString fileName() const { return _fileName; }
 
  QStringList     meshList() const;
  const MeshInfo* setCurrentMesh(QString mesh, bool meshInfo, bool bcInfo, bool dsInfo, const GmPCR* pcr);
 
  const MeshInfo* currentMeshInfo() const { return _info; }
 
  bool loadGeometryCoordinates (GmValueAccessor* coordAc, GmValueAccessor* ghostAc) const;
  bool loadGeometryCellTopology(const QStringList& propertySetOrder, const QVector<int>& propertySetSizes, 
                                int* type, int* propertySets, int* nodes, int* cellsPerType) const;
  bool loadAndSetMeshGhostCells   (GmCellMesh* mesh) const;
  bool loadAndSetMeshInactiveCells(GmCellMesh* mesh) const;
  bool loadCellGroupData          (int groupIndex, int* data, int* realNumCells) const;
  bool loadNodeSetData            (int groupIndex, int* data, int* realNumNodes) const;
  bool loadBoundaryData           (int groupIndex, const GmCellMesh* mesh, int* data, int* realNumCells, GmCellBoundaryType type) const;
 
  bool loadNodeData (int attrIndex, GmValueAccessor* ac) const;
  bool loadCellData (int attrIndex, GmCellAccessor*  ac) const;
  bool loadGaussData(int attrIndex, GmGaussAccessor* ac, const QVector<int>& cellTypeList) const;
 
  bool loadBcApplicationPoints(int bcIndex, const GmMesh* mesh, int* data, QVector<int>& skipList) const;
  bool loadBcData             (int bcIndex, int attrIndex, GmValueAccessor* ac, const QVector<int>& skipList) const;
 
  bool loadDsCoordinateData(int dsIndex, int dIndex, double* data) const;
  bool loadDsBorderData    (int dsIndex, int dIndex, const GmCellMesh* mesh, int* data, int* realNumCells) const;
  bool loadDsTriangleData  (int dsIndex, int dIndex, int* data) const;
 
private:
  bool loadSharedDictionariesInfo();
  bool loadMeshInfo();
  bool loadBcInfo();
  bool loadDsInfo();
  bool loadAttributeInfo(GroupType group, QString subGroupName, QVector<AttrInfo>& data, QString excludedTable = "");
 
  QString groupName(GroupType group, QString subGroupName) const { return GmHdf5Model::groupName(_currentMesh, group, subGroupName); }
 
  int checkTableLayout(QString tablePath, int ncol, const H5::DataType& type, bool required, 
                       int* nlines = NULL, int* ncolumns = NULL) const;
 
  bool loadDictionary(DictionaryType type, QStringList& data,      bool required) const;
  bool loadDictionary(DictionaryType type, QVector<DicInfo>& data, bool required) const;
 
  bool loadDictionaryData(int groupIndex, GroupType group, QVector<DicInfo>& dictionary, int* data) const;
  
  bool loadAccessorData(const AttrInfo& attrInfo, GroupType group, QString subGroupName, GmValueAccessor* ac, 
                        const GmMeshGaussIndex* gindex = NULL, const QVector<int>& bcSkipOrCellTypeList = QVector<int>()) const;
 
  bool isFunctionValue(double v) const { return GmNanBoxedIndex(v, false).isIndex(); }
 
  unsigned nanToDindex(double v) const { assert(isFunctionValue(v));  return GmNanBoxedIndex(v, false).toIndex(); }
 
  QString              _fileName;    
  const GmLogCategory& _logger;      
 
  H5::H5File*          _fileh;       
  QString              _currentMesh; 
  MeshInfo*            _info;        
  
  const GmPCR* _pcr;    
  const GmPCR* _defPcr; 
};
 
 
#endif