includeLinks/TopologyInfo.hpp

Go to the documentation of this file.

/* ***************************************************************** 
    MESQUITE -- The Mesh Quality Improvement Toolkit

    Copyright 2004 Lawrence Livermore National Laboratory.  Under 
    the terms of Contract B545069 with the University of Wisconsin -- 
    Madison, Lawrence Livermore National Laboratory retains certain
    rights in this software.

    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
    License as published by the Free Software Foundation; either
    version 2.1 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public License 
    (lgpl.txt) along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
    kraftche@cae.wisc.edu    
   
  ***************************************************************** */

#ifndef MESQUITE_TOPOLOGY_INFO_HPP
#define MESQUITE_TOPOLOGY_INFO_HPP

#include "Mesquite.hpp"

namespace Mesquite
{

class MsqError;

class TopologyInfo
{
  public:
    
    static unsigned dimension( EntityTopology topo )
      { return topo >= MIXED ? 0: instance.dimMap[topo]; }
    
    static unsigned adjacent( EntityTopology topo, unsigned dimension )
      { return (topo >= MIXED  || dimension > 2) ? 0 : instance.adjMap[topo][dimension]; }
    
    static unsigned sides( EntityTopology topo )
      { return (topo >= MIXED || instance.dimMap[topo] < 1) ? 0 : 
          instance.adjMap[topo][instance.dimMap[topo]-1]; }
    
    static unsigned corners( EntityTopology topo )
      { return topo >= MIXED ? 0 : instance.adjMap[topo][0]; }
    
    static unsigned edges( EntityTopology topo )
      { return adjacent(topo, 1); }
      
    static unsigned faces( EntityTopology topo )
      { return adjacent(topo, 2 ); }
      
    static void higher_order( EntityTopology topo, unsigned num_nodes,
                              bool& midedge, bool& midface, bool& midvol,
                              MsqError &err );
    
    static const unsigned* edge_vertices( EntityTopology topo,
                                          unsigned edge_number,
                                          MsqError& err );
    
    static const unsigned* face_vertices( EntityTopology topo,
                                          unsigned face_number,
                                          unsigned& num_vertices_out,
                                          MsqError& err );

    static const unsigned* side_vertices( EntityTopology topo,
                                          unsigned side_dimension, 
                                          unsigned side_number,
                                          unsigned& num_verts_out,
                                          MsqError& err );

     
      
    static void side_number( EntityTopology topo, 
                             unsigned connectivity_length,
                             unsigned node_index,
                             unsigned& side_dimension_out,
                             unsigned& side_number_out,
                             MsqError& err );

  private:
 
    enum {
      MAX_EDGES = 12,
      MAX_FACES = 6,
      MAX_FACE_CONN = 5,
      FIRST_FACE = TRIANGLE,
      LAST_FACE = QUADRILATERAL,
      FIRST_VOL= TETRAHEDRON,
      LAST_VOL = PYRAMID
    };
    
    unsigned dimMap[MIXED];    
    unsigned adjMap[MIXED][3]; 
    unsigned edgeMap[LAST_VOL-FIRST_FACE+1][MAX_EDGES][2] ;
    unsigned faceMap[LAST_VOL-FIRST_VOL+1][MAX_FACES][MAX_FACE_CONN];

    TopologyInfo();
    
    static TopologyInfo instance;
  
};

} //namespace Mesquite

#endif