Connectivity.h

// -*- tab-width: 2; indent-tabs-mode: nil; coding: utf-8-with-signature -*-
//-----------------------------------------------------------------------------
// Copyright 2000-2022 CEA (www.cea.fr) IFPEN (www.ifpenergiesnouvelles.com)
// See the top-level COPYRIGHT file for details.
// SPDX-License-Identifier: Apache-2.0
//-----------------------------------------------------------------------------
/*---------------------------------------------------------------------------*/
/* Connectivity.h                                              (C) 2000-2015 */
/*                                                                           */
/* Descripteur de connectivité.                                              */
/*---------------------------------------------------------------------------*/
#ifndef ARCANE_CONNECTIVITY_H
#define ARCANE_CONNECTIVITY_H
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#include "arcane/ArcaneTypes.h"
#include "arcane/MeshVariable.h"
#include "arcane/VariableTypes.h"
 
#include <iostream>
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ARCANE_BEGIN_NAMESPACE
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
class ARCANE_CORE_EXPORT Connectivity
{
public:
 
  enum eConnectivityType {
    CT_Null       = 0,
    // From Node
    // CT_NodeToNode = 1<<1,
    CT_NodeToEdge = 1<< 2,
    CT_NodeToFace = 1<< 3, // Utilisé dans UnstructuredMeshUtilities::localIdsFromConnectivity
    CT_NodeToCell = 1<< 4,
    // From Edge
    CT_EdgeToNode = 1<< 5,
    // CT_EdgeToEdge = 1<<6,
    CT_EdgeToFace = 1<< 7,
    CT_EdgeToCell = 1<< 8,
    // From Face
    CT_FaceToNode = 1<< 9,
    CT_FaceToEdge = 1<<10,
    CT_FaceToFace = 1<<11,
    CT_FaceToCell = 1<<12,
    // From Cell
    CT_CellToNode = 1<<13,
    CT_CellToEdge = 1<<14,
    CT_CellToFace = 1<<15,
    // CT_CellToCell = 1<<16,
    CT_DoFToNode     = 1<<17,
    CT_DoFToEdge     = 1<<18,
    CT_DoFToFace     = 1<<19,
    CT_DoFToCell     = 1<<20,
    CT_DoFToDoF      = 1<<21,
    CT_DoFToParticle = 1<<22,
 
    // Existing kind
    CT_HasNode    = 1<<23,
    CT_HasEdge    = 1<<24,
    CT_HasFace    = 1<<25,
    CT_HasCell    = 1<<26,
 
    // Frozen mark
    CT_Frozen     = 1<<27,
 
    // Dimension mark
    CT_Dim1D      = 1<<28,
    CT_Dim2D      = 1<<29,
    CT_Dim3D      = 1<<30,
 
    // Types composés
    // Connectivité par défaut
    CT_Default = CT_NodeToCell + CT_NodeToFace
               + CT_FaceToNode + CT_FaceToCell
               + CT_CellToNode + CT_CellToFace
               + CT_HasNode + CT_HasFace + CT_HasCell,
     
    CT_Default1D = CT_NodeToCell // minimum en 1D
                 + CT_FaceToNode + CT_FaceToCell
                 + CT_CellToNode + CT_CellToFace
                 + CT_HasNode + CT_HasFace + CT_HasCell,
 
    CT_Default2D = CT_NodeToCell // minimum en 2D
                 + CT_FaceToNode + CT_FaceToCell
                 + CT_CellToNode + CT_CellToFace
                 + CT_HasNode + CT_HasFace + CT_HasCell,
 
 
    CT_Default3D = CT_NodeToCell // minimum en 3D
                 + CT_FaceToNode + CT_FaceToCell 
                 + CT_CellToNode + CT_CellToFace
                 + CT_HasNode + CT_HasFace + CT_HasCell,
 
    CT_FullConnectivity2D = CT_NodeToFace + CT_NodeToCell
                          + CT_FaceToNode + CT_FaceToCell
                          + CT_CellToNode + CT_CellToFace
                          + CT_HasNode + CT_HasFace + CT_HasCell,
 
    CT_FullConnectivity3D = CT_NodeToEdge + CT_NodeToFace + CT_NodeToCell
                          + CT_EdgeToNode + CT_EdgeToFace + CT_EdgeToCell
                          + CT_FaceToNode + CT_FaceToEdge + CT_FaceToCell + CT_FaceToFace
                          + CT_CellToNode + CT_CellToEdge + CT_CellToFace +
                          + CT_HasNode + CT_HasEdge + CT_HasFace + CT_HasCell,
 
    CT_EdgeConnectivity = CT_HasEdge
                        + CT_NodeToEdge + CT_FaceToEdge + CT_CellToEdge +
                        + CT_EdgeToNode + CT_EdgeToFace + CT_EdgeToCell,
 
    CT_GraphConnectivity = CT_DoFToNode + CT_DoFToEdge + CT_DoFToFace + CT_DoFToCell + CT_DoFToDoF + CT_DoFToParticle
  };
 
                       
public:
  class Printer {
  public:
    Printer(const Integer connectivity) : m_connectivity(connectivity) { }
    virtual ~Printer() { }
    void print(std::ostream & o) const;
  private:
    const Integer m_connectivity;
  };
 
public:
  Connectivity(VariableScalarInteger connectivity);
  
  virtual ~Connectivity();
 
  bool hasFace() const { return hasConnectivity(CT_HasFace); }
  bool hasEdge() const { return hasConnectivity(CT_HasEdge); }
 
  void enableConnectivity(const Integer c);
  void disableConnectivity(const Integer c);
 
  bool hasConnectivity(const Integer c) const;
  inline static bool hasConnectivity(const Integer connectivity, const Integer c);
 
  bool isFrozen() const;
  void freeze(IMesh * mesh);
  
  static Integer getPrealloc(const Integer connectivity, eItemKind kindA, eItemKind kindB);
 
  static void print(std::ostream & o, const Integer connectivity);
 
  static Integer kindsToConnectivity(eItemKind kindA, eItemKind kindB);
 
private:
  VariableScalarInteger m_connectivity;
 
  void _enableConnectivity(const Integer c);
  void _disableConnectivity(const Integer c);
  static inline bool _hasConnectivity(const Integer connectivity, const Integer c);
 
  static void _checkValid(const Integer c);
  void _checkFrozen() const;
};
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
inline bool Connectivity::
hasConnectivity(const Integer connectivity, const Integer c) 
{ 
  // _checkValid(c);
  return _hasConnectivity(connectivity,c); 
}
 
/*---------------------------------------------------------------------------*/
 
inline bool Connectivity::
_hasConnectivity(const Integer connectivity, const Integer c) 
{ 
  return ((connectivity & c) != 0); 
}
 
/*---------------------------------------------------------------------------*/
 
ARCANE_CORE_EXPORT std::ostream & operator<<(std::ostream & o, const Connectivity::Printer & p);
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ARCANE_END_NAMESPACE
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#endif /* ARCANE_CONNECTIVITY_H */