HDF5OrthoReader Class Reference

#include <hdf5orthoreader.h>

Inheritance diagram for HDF5OrthoReader:

[legend]

Collaboration diagram for HDF5OrthoReader:

[legend]

List of all members.

Public Member Functions
	HDF5OrthoReader (bool debuginfo=false)
OrthoMesh &	getMesh (string strMeshName)
void	readScalars (VecDouble &vec, std::string dataSetName)
hid_t	getFile ()
	~HDF5OrthoReader ()
OrthoMesh *	loadOrthoMesh (std::string strMeshName)
void	iterateTree (hid_t id, HDF5Operator &op)
void	list ()
Protected Member Functions
HDF5OrthoInfo	readOrthoMeshInfo (hid_t triaGroup)
Private Types
typedef __gnu_cxx::hash_map < string, OrthoMesh *, StringHashFun, struct eqstr >	MeshHeap
Private Attributes
MeshHeap	_meshCash

---

Detailed Description

HDF5OrthoRead

Definition at line 19 of file hdf5orthoreader.h.

---

Member Typedef Documentation

typedef __gnu_cxx::hash_map<string,OrthoMesh*,StringHashFun,struct eqstr> HDF5OrthoReader::MeshHeap [private]

Definition at line 22 of file hdf5orthoreader.h.

---

Constructor & Destructor Documentation

HDF5OrthoReader::HDF5OrthoReader

(

bool

debuginfo = false

)

Definition at line 4 of file hdf5orthoreader.cpp.

{
  m_file =-1;
  if (!debuginfo)
  {
    H5Eset_auto(H5E_DEFAULT,NULL,stderr);
    H5Eset_auto(H5E_DEFAULT,NULL,stdout);
  }

}

HDF5OrthoReader::~HDF5OrthoReader

(

)

[inline]

Definition at line 32 of file hdf5orthoreader.h.

{};

---

Member Function Documentation

hid_t HDF5OrthoReader::getFile

(

)

Definition at line 45 of file hdf5orthoreader.cpp.

{
  assert(m_file > 0);
  return m_file; 
}

OrthoMesh & HDF5OrthoReader::getMesh

(

string

strMeshName

)

Definition at line 17 of file hdf5orthoreader.cpp.

{
  if (_meshCash.find(strMeshName) != _meshCash.end())
  {
    return *(_meshCash[strMeshName]);
  }
  else
  {
    _meshCash[strMeshName]=loadOrthoMesh(strMeshName);
    return *(_meshCash[strMeshName]);
  }
}

void HDF5OrthoReader::iterateTree	(	hid_t	id,
		HDF5Operator &	op
	)

Definition at line 72 of file hdf5orthoreader.cpp.

{

  switch(H5Iget_type(id))
  {
  case H5I_GROUP:
    {
      char strElem[STR_SIZE];
      hsize_t nElems;

      op.processGroup(id);

      //Iterate among all attributes
      unsigned nAtts = H5Aget_num_attrs(id);
      for (unsigned i=0;i<nAtts;i++)
      {
        hid_t attId = H5Aopen_idx(id,i);
        iterateTree(attId,op);
        H5Aclose(attId);
      }
      H5Gget_num_objs(id,&nElems);
      for(unsigned i=0;i<nElems;i++)
      {
        H5Gget_objname_by_idx(id,i,strElem,STR_SIZE);
        hsize_t type =H5Gget_objtype_by_idx(id,i);
        switch(type)
        {
        case H5G_LINK:
          {
            printf("Ignoring link %s\n",strElem);
            break;
          }
        case H5G_GROUP:
          {
            hid_t idChild=H5Gopen(id,strElem,H5P_DEFAULT);
            iterateTree(idChild,op);
            closeId(idChild);
            break;
          }
        case H5G_DATASET:
          {
            hid_t idChild=H5Dopen(id,strElem,H5P_DEFAULT);
            iterateTree(idChild,op);
            closeId(idChild);
            break;
          }
        case H5G_TYPE:
          {
            printf("Ignoring type %s\n",strElem);
            break;
          }
        }
      }
      break;
    }
  case H5I_DATASET:
    {
      op.processDataSet(id);
      //Iterate among all attributes
      unsigned nAtts = H5Aget_num_attrs(id);
      for (unsigned i=0;i<nAtts;i++)
      {
        hid_t attId = H5Aopen_idx(id,i);
        iterateTree(attId,op);
        H5Aclose(attId);
      }
      break;
    }
  case H5I_ATTR:
    {
      op.processAtt(id);
      break;
    }
  default:
    {
      throw new Exception("HDF5OrthoReader::iterateTree() Expect HDF5 object to be a dataset, a group or attribute");
    }
  }
}

void HDF5OrthoReader::list

(

)

Definition at line 52 of file hdf5orthoreader.cpp.

{
  hid_t triaG = openGroup("Triangulations");
  hid_t test = H5Gopen1(triaG,"./tria0");
  char name2[1000];
  H5Iget_name(test,name2,1000);
  printf("Test %s\n",name2);
  hsize_t numObjs;
  H5Gget_num_objs(triaG,&numObjs);
  for (unsigned i =0;i<numObjs;i++)
  {
    char name[1000];
    H5Gget_objname_by_idx(triaG,i,name,1000);
    printf("%s\n",name);
  }
  
}

OrthoMesh * HDF5OrthoReader::loadOrthoMesh

(

std::string

strMeshName

)

Definition at line 31 of file hdf5orthoreader.cpp.

{
  char path[500];
  sprintf(path,"Triangulations/%s",strMeshName.c_str());
  hid_t Gtria = H5Gopen1(getFile(),path);
  if (Gtria < 0)
    throw new Exception("HDF5OrthoReader::loadOrthoMesh Could not find mesh \"%s\"",path);
  HDF5OrthoInfo info;
  info.readOrthoMeshInfo(Gtria);
  OrthoMesh* pMesh = new OrthoMesh(info.P0,info.P1,info.nElems[0],info.nElems[1],info.nElems[2],info.wells);
  return pMesh;
}

HDF5OrthoInfo HDF5OrthoReader::readOrthoMeshInfo

(

hid_t

triaGroup

)

[protected]

void HDF5OrthoReader::readScalars	(	VecDouble &	vec,
		std::string	dataSetName
	)

Definition at line 183 of file hdf5orthoreader.cpp.

{
  hid_t dataset=H5Dopen(getFile(),dataSetName.c_str(),H5P_DEFAULT);
  if (dataset < 0)
    throw new Exception("Could not read dataset \"%s\"",dataSetName.c_str());

  hid_t dataspace=H5Dget_space(dataset);
  unsigned rank=H5Sget_simple_extent_ndims(dataspace);
  hsize_t dims[rank],maxdims[rank];
  H5Sget_simple_extent_dims(dataspace,dims,maxdims);
  unsigned nElems=dims[0];
  for (unsigned i=1;i<rank;i++)
  {
    nElems*=dims[i];
  }
  if (vec.size() != nElems)
    vec.reinit(nElems);
  //READ THE DATA
  H5Dread(dataset,H5T_NATIVE_DOUBLE,H5S_ALL,H5S_ALL,H5P_DEFAULT,vec.begin());
  H5Sclose(dataspace);
  H5Dclose(dataset);
   
}

---

Member Data Documentation

MeshHeap HDF5OrthoReader::_meshCash [private]

Definition at line 23 of file hdf5orthoreader.h.

---

The documentation for this class was generated from the following files:

• hdf5orthoreader.h
• hdf5orthoreader.cpp