---

ExodusMeshReader.h

Go to the documentation of this file.

#ifndef EXODUSMESHREADER_H
#define EXODUSMESHREADER_H

#include "SundanceDefs.h"

#include "MeshReader.h"

#include "Expr.h"


namespace Sundance
{

  using namespace TSF;
  using std::string;

  using std::ostream;

  /** \ingroup UserLevelGeometry
   * Reads a 2D or 3D mesh stored in the Exodus file format.
   */

  class ExodusMeshReader : public MeshReaderBase
    {
    public:
      /** Construct an ExodusMeshReader for reading from a given file */
      ExodusMeshReader(const string& filename)
        : filename_(filename),
        hasMesh_(false),
        hasNodalNames_(false),
        hasElementalNames_(false),
        nodalVarNames_(1),
        elementalVarNames_(1)
        {;}

      /** read a mesh from the reader's file */
      virtual Mesh getMesh() const ;

      /** read in a list of Nodal Field by index */
      Expr getNodalField(const TSFArray<int>& attributeIndices, int timeStep = 1);
      /** read in a single Nodal Field by index */
      Expr getNodalField(const int attributeIndex, int timeStep = 1)
        {
          TSFArray<int> attributeIndices(1);
          attributeIndices[0] = attributeIndex;
          return getNodalField(attributeIndices, timeStep);
        }

      /** read in a list of elemental Fields by index */
      Expr getElementalField(const TSFArray<int>& attributeIndices, int timeStep = 1);
      /** read in a single elemental Field by index */
      Expr getElementalField(const int attributeIndex, int timeStep = 1)
        {
          TSFArray<int> attributeIndices(1);
          attributeIndices[0] = attributeIndex;
          return getElementalField(attributeIndices, timeStep);
        }

      /** read in a list of nodal Fields by name */
      Expr getNodalField(const TSFArray<string>& attributeNames, int timeStep = 1);
      /** read in a single nodal Field by name */
      Expr getNodalField(const string attributeName, int timeStep = 1)
        {
          TSFArray<string> attributeNames(1);
          attributeNames[0] = attributeName;
          return getNodalField(attributeNames, timeStep);
        }

      /** read in a list of elemental Fields by name */
      Expr getElementalField(const TSFArray<string>& attributeNames, int timeStep = 1);
      /** read in a single elemental Field by name */
      Expr getElementalField(const string attributeName, int timeStep = 1)
        {
          TSFArray<string> attributeNames(1);
          attributeNames[0] = attributeName;
          return getElementalField(attributeNames, timeStep);
        }

      TSFArray<string>& getNodalVarNames();
      TSFArray<string>& getElementalVarNames();
    private:

      TSFArray<TSFArray<double> > readNodalField(const TSFArray<int>& attributeIndices, int timeStep = 1);

      TSFArray<TSFArray<double> > readElementalField(const TSFArray<int>& attributeIndices, int timeStep = 1);


      static const int CPUWORDSIZE;
      static const int IOWORDSIZE;

      string filename_;
      mutable Mesh mesh_;

      mutable bool hasMesh_;
      bool hasNodalNames_;
      bool hasElementalNames_;
      TSFArray<string> nodalVarNames_;
      TSFArray<string> elementalVarNames_;

    };


}
#endif

---

Contact:

Kevin Long (krlong@ca.sandia.gov)


Documentation generated by