gmXdmfHdf5.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_XDMF_HDF5_H_
#define _GEMA_XDMF_HDF5_H_
 
#include <QString>
#include <functional>
 
#include "gmValueInfo.h"
#include "gmFileIO.h"
#include "gmFileFilter.h"
#include "gmFileIOUtils.h"
 
class GmLogCategory;
class GmTaskManager;
 
namespace H5 {
  class H5File;
  class H5Object;
  class PredType;
}
 
class GmSimulationData;
class GmMesh;
class GmCellMesh;
class GmElementMesh;
class GmCell;
class GmValueAccessor;
class GmCellAccessor;
class GmGaussAccessor;
class GmDiscontinuitySet;
 
class GmFileFormat;
class GmFileFilter;
 
class UnitConverter;
 
class GmXdmfHdf5
{
    using TransformF = GmFileIOUtils::TransformF;
public:
    struct DsSizeInfo {
        int originalNodeCount;    
        int originalElemCount;    
        int originalExtElemCount; 
        int intersectNodeCount;   
        int intersectElemCount0;  
        int intersectElemCount1;  
        int dataSize;             
    };
 
  GmXdmfHdf5(const GmLogCategory& logger);
  ~GmXdmfHdf5();
 
  bool create(QString fileName, const GmFileFormat* fileFormat, const GmSimulationData* simData,
              bool singlePrecision, int sequenceNumber);
 
  static bool addSequenceInfo(QString fileName, int numSequenceFiles, double startTime, double endTime, QString timeUnit, const GmLogCategory& logger);
 
  bool createGroups(const QStringList& groups);
  bool createGroup (QString group, double time, QString timeUnit);
 
  bool adjustTimeInfo(int sequenceNum, bool multiTime, double startTime, double endTime, QString timeUnit);
 
  bool saveMeshNodeCoord(QString dsName, const GmMesh* mesh, const GmFileFilter* filter, 
                         GmValueAccessor* nodeAc, const double* scaleFactor, const double* offset);
 
  bool saveMeshGaussNodeCoord(QString dsName, const GmElementMesh* mesh, const GmFileFilter* filter, int ruleIndex,
                              GmValueAccessor* nodeAc, const double* scaleFactor, const double* offset);
 
  bool saveMeshTopology(QString dsName, const GmCellMesh* mesh, const GmFileFilter* filter);
 
  bool saveMeshPointCloudTopology(QString dsName, const GmMesh* mesh, const GmFileFilter* filter, int ruleIndex = -1);
 
  bool saveMeshNodeData(QString dsName, const GmMesh* mesh, const GmFileFilter* filter,
                        GmValueAccessor* dataAc, int dimFilter, const QVector<int>& indexer, bool minMaxAttrib);
 
  bool saveMeshNodeIdColumn(QString dsName, const GmMesh* mesh, const GmFileFilter* filter, GmFileIO::ExtraColumns column);
 
  bool saveMeshCellData(QString dsName, const GmCellMesh* mesh, const GmFileFilter* filter,
                        GmCellAccessor* dataAc, int dimFilter, const QVector<int>& indexer, bool minMaxAttrib);
 
  bool saveMeshGaussData(QString dsName, const GmElementMesh* mesh, const GmFileFilter* filter, int ruleIndex,
                         GmGaussAccessor* dataAc, int dimFilter, const QVector<GmFileFilter::Indexer>& preIndex,int maxCellEntries, bool minMaxAttrib);
 
  bool saveDiscontinuitySetData(QString coordDsName, QString cellDsName, const GmDiscontinuitySet* discSet, bool original,
                                const UnitConverter* conv, const double* scaleFactor, const double* offset, 
                                DsSizeInfo& sizeInfo);
 
  bool saveMeshCellExtraColumn(QString dsName, const GmCellMesh* mesh, const GmFileFilter* filter,
                               GmFileIO::ExtraColumns column, int propertySet);
 
  bool saveMeshGaussIdColumn(QString dsName, const GmElementMesh* mesh, const GmFileFilter* filter, int ruleIndex);
 
private:
  
 
  bool saveNodeData(QString dsName, const GmMesh* mesh, const GmFileFilter* filter, GmValueAccessor* dataAc,
                    int dimFilter, const TransformF& transformFunc, const QVector<int>& indexer);
 
 
  template <class T>
  bool fillDiscontinuityPolylineData(const GmDiscontinuitySet* discSet, bool original, const UnitConverter* conv, 
                                     const double* scaleFactor, const double* offset, 
                                     QVector<T>& nodes, QVector<int>& topology, int* nsegments) const;
  template <class T>
  bool fillDiscontinuityTriangleData(const GmDiscontinuitySet* discSet, bool original, const UnitConverter* conv,
      const double* scaleFactor, const double* offset,
      QVector<T>& nodes, QVector<int>& topology0, int* nsegments0,
      QVector<int>& topology1, int* nsegments1) const;
 
  const char* dataBuffer(GmValueAccessor* ac, const GmValueInfo* info) const;
 
  const H5::PredType& hdf5Type(GmStorageType type) const;
 
  void addTableAttributes(H5::H5Object* obj, const GmMesh* mesh,      const char* type, const GmDiscontinuitySet* discSet = NULL);
  void addNodeAttributes (H5::H5Object* obj, const GmMesh* mesh,      const GmFileFilter* filter);
  void addCellAttributes (H5::H5Object* obj, const GmCellMesh* mesh,  const GmFileFilter* filter, int ruleIndex = -1);
  void addAcAttributes   (H5::H5Object* obj, const GmValueInfo* info, QString unit, int dimFilter, const double* scaleFactor = NULL, const double* offset = NULL, const double* min = NULL, const double* max = NULL);
  void addPropertySetAttributes(H5::H5Object* obj, const GmCellMesh* mesh, int propertySet);
  void addDiscSetNodeAttributes(H5::H5Object* obj, const GmDiscontinuitySet* discSet, bool original, const UnitConverter* conv,
                                const double* scaleFactor, const double* offset);
  void addDiscSetTopologyAttributes(H5::H5Object* obj, const GmDiscontinuitySet* discSet, bool original);
 
  const GmLogCategory& _logger;          
  QString              _fileName;        
  H5::H5File*          _fileh;           
  const GmFileFormat*  _fileFormat;      
  bool                 _singlePrecision; 
  GmTaskManager*       _tm;              
};
 
#endif