meshtopology.h

利用C

// Copyright (C) 2006 Anders Logg.
// Licensed under the GNU LGPL Version 2.1.
//
// First added:  2006-05-08
// Last changed: 2007-11-30

#ifndef __MESH_TOPOLOGY_H
#define __MESH_TOPOLOGY_H

#include <dolfin/log/dolfin_log.h>
#include <dolfin/common/types.h>
#include <dolfin/common/Array.h>
#include "MeshConnectivity.h"

namespace dolfin
{
  
  /// MeshTopology stores the topology of a mesh, consisting of mesh entities
  /// and connectivity (incidence relations for the mesh entities). Note that
  /// the mesh entities don't need to be stored, only the number of entities
  /// and the connectivity. Any numbering scheme for the mesh entities is
  /// stored separately in a MeshFunction over the entities.
  ///
  /// A mesh entity e may be identified globally as a pair e = (dim, i), where
  /// dim is the topological dimension and i is the index of the entity within
  /// that topological dimension.
  
  class MeshTopology
  {
  public:
    
    /// Create empty mesh topology
    MeshTopology();

    /// Copy constructor
    MeshTopology(const MeshTopology& topology);
    
    /// Destructor
    ~MeshTopology();

    /// Assignment
    const MeshTopology& operator= (const MeshTopology& topology);
 
    /// Return topological dimension
    inline uint dim() const { return _dim; }
    
    /// Return number of entities for given dimension
    inline uint size(uint dim) const
    { dolfin_assert(dim <= _dim); return num_entities[dim]; }

    /// Clear all data
    void clear();

    /// Initialize topology of given maximum dimension
    void init(uint dim);

    /// Set number of entities (size) for given topological dimension
    void init(uint dim, uint size);

    /// Return connectivity for given pair of topological dimensions
    inline MeshConnectivity& operator() (uint d0, uint d1)
    { dolfin_assert(d0 <= _dim && d1 <= _dim); return connectivity[d0][d1]; }

    /// Return connectivity for given pair of topological dimensions
    inline const MeshConnectivity& operator() (uint d0, uint d1) const
    { dolfin_assert(d0 <= _dim && d1 <= _dim); return connectivity[d0][d1]; }

    /// Display data
    void disp() const;

  private:

    friend class MPIMeshCommunicator;

    // Topological dimension
    uint _dim;
  
    // Number of mesh entities for each topological dimension
    uint* num_entities;

    // Connectivity for pairs of topological dimensions
    MeshConnectivity** connectivity;
   
  };

}

#endif