---

MeshData.h

Go to the documentation of this file.

#ifndef MESHDATA_H
#define MESHDATA_H

#include "SundanceDefs.h"

#include <string>
#include "TSFArray.h"
#include "Cell.h"
#include "TSFRBTree.h"
#include "Point.h"
#include "TSFIntHashtable.h"
#include "MPIComm.h"



namespace Sundance
{

  using namespace TSF;
  using std::string;

  using std::ostream;

  /** \ingroup LowLevelGeometry

  MeshData contains the data tables defining a mesh. It lives
  inside a mesh, stored by ref-counted pointer, so you can make lots of
  copies of the mesh object which all refer back to the same
  data. MeshData has no public accessors. Classes Mesh and CellBase are friends and
  access the mesh data directly.

  */



  class MeshData
    {
    public:
      /** */
      MeshData(int spatialDim, const MPIComm& comm);
      MeshData(int spatialDim, int rank, int nProc);
      ~MeshData(){;}

      void diagnostics() const ;

      int lookupLabel(const string& label);
      int lookupTypeCode(const string& typeName) const ;

      void textWrite(ostream& os) const ;

      friend class Mesh;
      friend class CellBase;
      friend class ZeroCell;
      friend class MaximalCell;
      friend class IntermediateCell;
      friend class Cell;
      friend class BoundaryCellSet;
    private:
      MeshData(const MeshData& /*other*/){TSFError::raise("mesh data copy ctor");}

      // multiprocessor information
      MPIComm comm_;
      int procID_;
      int numProcs_;

      // mesh identification.
      // meshID_ is the ID number for this mesh.
      // topID_ is the ID number for the next mesh created.
      static int topID_;
      int meshID_;

      // spatial dimension of this mesh.
      int spatialDim_;

      TSFArray<int> localCellCount_;

      // a flat array of every point in the local processor's submesh.
      TSFArray<Point> points_;
      // list of global indices
      TSFArray<int> pointGlobalIndex_;
      // list of which processor owns each point
      TSFArray<int> pointOwnerProcID_;
      // list of the processors that contact, but don't own, each point.
      TSFArray<TSFArray<int> > pointSharerProcID_;
      // map global point index to local point index
      TSFIntHashtable pointIndexMap_;

      // pointToZeroCellMap lets us find a zero cell given a point index.
      // This will help us locate a cell's facets based given its node numbers.
      TSFArray<int> pointToZeroCellMap_;
      // table of cells, indexed by (dimension, index).
      TSFArray<TSFArray<Cell> > cells_;
      // map global cell index to local cell index
      TSFArray<TSFIntHashtable> cellIndexMap_;

      // processor IDs for cells
      TSFArray<TSFArray<int> > ownerProcID_;

      // faceting information for intermediate cells
      TSFArray<TSFArray<char> > myFacetIndex_;

      // labeling
      TSFArray<string> labels_;
      TSFRBTree<string, short int> labelToIndexMap_;
      TSFArray<TSFArray<short int> > labelIndices_;

      // identification of boudnary cells
      TSFArray<bool> isBoundary_;

      // Auxiliary connectivity information, used for maximal cells only (which
      // is why we keep it here instead of in all cells).

      // for each maximal cell, list the other maxiaml cells that share a node.
      // maximalNeighbors_[cellIndex][neighborIndex]
      TSFArray<TSFArray<int> > maximalNeighbors_;
      // for each maximal cell, list the number of nodes shared by the neighbors
      // neighborWeights_[cellIndex][neighborIndex]
      TSFArray<TSFArray<int> > neighborWeights_;

      static string emptyLabel_;
    };

}
#endif

---

Contact:

Kevin Long (krlong@ca.sandia.gov)


Documentation generated by