CartesianGrid.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
//-----------------------------------------------------------------------------
/*---------------------------------------------------------------------------*/
/* CartesianGrid.h                                             (C) 2000-2021 */
/*                                                                           */
/* Grille cartésienne avec noeuds, faces et mailles.                         */
/*---------------------------------------------------------------------------*/
#ifndef ARCANE_CARTESIANMESH_CARTESIANGRID_H
#define ARCANE_CARTESIANMESH_CARTESIANGRID_H
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#include "arcane/cartesianmesh/v2/CartesianNumbering.h"
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
namespace Arcane::CartesianMesh::V2
{
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
template <typename IdType>
class CartesianGrid
{
 public:
  using IdType3 = IdType[3];
 
  using CartesianNumberingType = CartesianNumbering<IdType>;
 
  using CartesianNumberingType3 = CartesianNumberingType[3];
 
 public:
  //  param[in] dimension  Dimension de la grille cartésienne (au plus 3)
  CartesianGrid(const IdType3& ncells_dir, Integer dimension)
  : m_dimension(dimension)
  {
 
    for (Integer d(0); d < m_dimension; ++d) {
      m_ncells_dir[d] = ncells_dir[d];
      m_nnodes_dir[d] = m_ncells_dir[d] + 1;
    }
    m_cart_num_cell.initNumbering(m_ncells_dir, dimension);
 
    for (Integer d(m_dimension); d < 3; ++d) {
      m_ncells_dir[d] = 1;
      m_nnodes_dir[d] = 1;
    }
    m_cart_num_node.initNumbering(m_nnodes_dir, dimension);
 
    // Pour les faces, on va numéroter les faces selon X, puis selon Y et selon Z
    // On distingue les faces par leurs orientations (dnorm)
    IdType nfaces_norm_first = 0; // Premier numero de face selon la normale dnorm
    for (Integer dnorm(0); dnorm < m_dimension; ++dnorm) {
 
      // Les directions orthogonales à la normale
      Integer d1 = (dnorm + 1) % 3;
      Integer d2 = (dnorm + 2) % 3;
 
      // Dans la direction de la normale, on a m_ncells_dir[dnorm]+1 faces
      // et dans les directions m_ncells_dir[d] faces
      m_nfaces_norm_dir[dnorm][dnorm] = m_ncells_dir[dnorm] + 1;
      m_nfaces_norm_dir[dnorm][d1] = m_ncells_dir[d1];
      m_nfaces_norm_dir[dnorm][d2] = m_ncells_dir[d2];
 
      m_cart_num_face[dnorm].initNumbering(m_nfaces_norm_dir[dnorm], dimension, nfaces_norm_first);
      nfaces_norm_first += m_cart_num_face[dnorm].nbItem();
    }
 
    for (Integer dnorm(m_dimension); dnorm < 3; ++dnorm) {
      m_nfaces_norm_dir[dnorm][0] = 1;
      m_nfaces_norm_dir[dnorm][1] = 1;
      m_nfaces_norm_dir[dnorm][2] = 1;
    }
  }
 
  const CartesianNumberingType& cartNumCell() const
  {
    return m_cart_num_cell;
  }
 
  const CartesianNumberingType& cartNumNode() const
  {
    return m_cart_num_node;
  }
 
  const CartesianNumberingType& cartNumFace(Integer dir) const
  {
    ARCANE_ASSERT(dir < m_dimension, ("La direction doit être strictement inférieure à la dimension"));
    return m_cart_num_face[dir];
  }
 
  const CartesianNumberingType3& cartNumFace3() const
  {
    return m_cart_num_face;
  }
 
  CartesianNumberingType* cartNumCellPtr()
  {
    return &m_cart_num_cell;
  }
 
  CartesianNumberingType* cartNumNodePtr()
  {
    return &m_cart_num_node;
  }
 
  CartesianNumberingType* cartNumFacePtr(Integer dir)
  {
    ARCANE_ASSERT(dir < m_dimension, ("La direction doit être strictement inférieure à la dimension"));
    return &(m_cart_num_face[dir]);
  }
 
  CartesianNumberingType3* cartNumFace3Ptr()
  {
    return &m_cart_num_face;
  }
 
  Integer dimension() const
  {
    return m_dimension;
  }
 
 protected:
  IdType3 m_ncells_dir = { 1, 1, 1 }; // Nb de mailles par direction
  IdType3 m_nnodes_dir = { 1, 1, 1 }; // Nb de noeuds par direction
  IdType3 m_nfaces_norm_dir[3]; 
 
  Integer m_dimension = 0;
 
  CartesianNumberingType m_cart_num_cell;
  CartesianNumberingType m_cart_num_node;
  CartesianNumberingType3 m_cart_num_face;
};
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
} // namespace Arcane::CartesianMesh::V2
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#endif