ItemGroupImpl.cc

// -*- tab-width: 2; indent-tabs-mode: nil; coding: utf-8-with-signature -*-
//-----------------------------------------------------------------------------
// Copyright 2000-2024 CEA (www.cea.fr) IFPEN (www.ifpenergiesnouvelles.com)
// See the top-level COPYRIGHT file for details.
// SPDX-License-Identifier: Apache-2.0
//-----------------------------------------------------------------------------
/*---------------------------------------------------------------------------*/
/* ItemGroupImpl.cc                                            (C) 2000-2024 */
/*                                                                           */
/* Implémentation d'un groupe d'entités de maillage.                         */
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#include "arcane/core/ItemGroupImpl.h"
 
#include "arcane/utils/String.h"
#include "arcane/utils/ITraceMng.h"
#include "arcane/utils/ArgumentException.h"
 
#include "arcane/core/ItemGroupObserver.h"
#include "arcane/core/IItemFamily.h"
#include "arcane/core/ItemGroup.h"
#include "arcane/core/IMesh.h"
#include "arcane/core/MeshPartInfo.h"
#include "arcane/core/MeshKind.h"
#include "arcane/core/ItemPrinter.h"
#include "arcane/core/IItemOperationByBasicType.h"
#include "arcane/core/ItemGroupComputeFunctor.h"
#include "arcane/core/IVariableSynchronizer.h"
#include "arcane/core/ParallelMngUtils.h"
#include "arcane/core/internal/ItemGroupInternal.h"
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
namespace Arcane
{
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*!
 * \internal
 * \brief Classe d'un groupe nul.
 */
class ItemGroupImplNull
: public ItemGroupImpl
{
 public:
 
  ItemGroupImplNull() : ItemGroupImpl() {}
  virtual ~ItemGroupImplNull() {} //!< Libére les ressources
 
 public:
 
  //! Retourne le nom du groupe
  const String& name() const { return m_name; }
  const String& fullName() const { return m_name; }
 
 public:
 
 public:
 
  virtual void convert(NodeGroup& g) { g = NodeGroup(); }
  virtual void convert(EdgeGroup& g) { g = EdgeGroup(); }
  virtual void convert(FaceGroup& g) { g = FaceGroup(); }
  virtual void convert(CellGroup& g) { g = CellGroup(); }
 
 public:
 
 private:
 
  String m_name;
};
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::shared_null= 0;
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
class ItemGroupImplItemGroupComputeFunctor
: public ItemGroupComputeFunctor
{
 private:
 
  typedef void (ItemGroupImpl::*FuncPtr)();
 
 public:
 
  ItemGroupImplItemGroupComputeFunctor(ItemGroupImpl * parent, FuncPtr funcPtr)
  : ItemGroupComputeFunctor()
  , m_parent(parent)
  , m_function(funcPtr) { }
 
 public:
 
  void executeFunctor() override
  {
    (m_parent->*m_function)();
  }
 
 private:
 
  ItemGroupImpl * m_parent;
  FuncPtr m_function;
};
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
checkSharedNull()
{
  // Normalement ce test n'est vrai que si on a une instance globale
  // de 'ItemGroup' ce qui est déconseillé. Sinon, _buildSharedNull() a été
  // automatiquement appelé lors de l'initialisation (dans arcaneInitialize()).
  if (!shared_null)
    _buildSharedNull();
  return shared_null;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl::
ItemGroupImpl(IItemFamily* family,const String& name)
: m_p (new ItemGroupInternal(family,name))
{
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl::
ItemGroupImpl(IItemFamily* family,ItemGroupImpl* parent,const String& name)
: m_p(new ItemGroupInternal(family,parent,name))
{
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl::
ItemGroupImpl()
: m_p (new ItemGroupInternal())
{
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl::
~ItemGroupImpl()
{
  delete m_p;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
const String& ItemGroupImpl::
name() const
{
  return m_p->name();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
const String& ItemGroupImpl::
fullName() const
{
  return m_p->fullName();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
Integer ItemGroupImpl::
size() const
{
  return m_p->itemsLocalId().size();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
empty() const
{
  return m_p->itemsLocalId().empty();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
Int32ConstArrayView ItemGroupImpl::
itemsLocalId() const
{
  return m_p->itemsLocalId();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
beginTransaction()
{
  if (m_p->m_transaction_mode)
    ARCANE_FATAL("Transaction mode already started");
  m_p->m_transaction_mode = true;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
endTransaction()
{
  if (!m_p->m_transaction_mode)
    ARCANE_FATAL("Transaction mode not started");
  m_p->m_transaction_mode = false;
  if (m_p->m_need_recompute) {
    m_p->m_need_recompute = false;
    m_p->m_need_invalidate_on_recompute = false;
    m_p->notifyInvalidateObservers();
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
Int32Array& ItemGroupImpl::
unguardedItemsLocalId(const bool self_invalidate)
{
  ITraceMng* trace = m_p->m_mesh->traceMng();
  trace->debug(Trace::Medium) << "ItemGroupImpl::unguardedItemsLocalId on group " << name()
                              << " with self_invalidate=" << self_invalidate;
 
  if (m_p->m_compute_functor && !m_p->m_transaction_mode)
    ARCANE_FATAL("Direct access for computed group in only available during a transaction");
 
  _forceInvalidate(self_invalidate);
  return m_p->mutableItemsLocalId();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
parent() const
{
  return m_p->m_parent;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
IMesh* ItemGroupImpl::
mesh() const
{
  return m_p->mesh();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
IItemFamily* ItemGroupImpl::
itemFamily() const
{
  return m_p->m_item_family;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
null() const
{
  return m_p->null();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
isOwn() const
{
  return m_p->m_is_own;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
setOwn(bool v)
{
  const bool is_own = m_p->m_is_own;
  if (is_own == v)
    return;
  if (!is_own) {
    if (m_p->m_own_group)
      ARCANE_THROW(NotSupportedException,"Setting Own with 'Own' sub-group already defined");
  }
  else {
    // On a le droit de remettre setOwn() à 'false' pour le groupe de toutes
    // les entités. Cela est nécessaire en reprise si le nombre de parties
    // du maillage est différent du IParallelMng associé à la famille
    if (!isAllItems())
      ARCANE_THROW(NotSupportedException,"Un-setting Own on a own group");
  }
  m_p->m_is_own = v;
  // (HP) TODO: Faut il notifier des observers ?
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
eItemKind ItemGroupImpl::
itemKind() const
{
  return m_p->kind();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
ownGroup()
{
  ItemGroupImpl* ii = m_p->m_own_group;
  // Le flag est déjà positionné dans le ItemGroupInternal::_init ou ItemGroupImpl::setOwn
  if (!ii) {
    if (m_p->m_is_own){
      ii = this;
      m_p->m_own_group = ii;
    } else {
      ii = createSubGroup("Own",m_p->m_item_family,new OwnItemGroupComputeFunctor());
      m_p->m_own_group = ii;
      ii->setOwn(true);
    }
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
ghostGroup()
{
  ItemGroupImpl* ii = m_p->m_ghost_group;
  if (!ii) {
    ii = createSubGroup("Ghost",m_p->m_item_family,new GhostItemGroupComputeFunctor());
    m_p->m_ghost_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
interfaceGroup()
{
  if (itemKind()!=IK_Face)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_interface_group;
  if (!ii) {
    ii = createSubGroup("Interface",m_p->m_item_family,new InterfaceItemGroupComputeFunctor());
    m_p->m_interface_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
nodeGroup()
{
  if (itemKind()==IK_Node)
    return this;
  ItemGroupImpl* ii = m_p->m_node_group;
  if (!ii){
    ii = createSubGroup("Nodes",m_p->m_mesh->nodeFamily(),new ItemItemGroupComputeFunctor<Node>());
    m_p->m_node_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
edgeGroup()
{
  if (itemKind()==IK_Edge)
    return this;
  ItemGroupImpl* ii = m_p->m_edge_group;
  if (!ii){
    ii = createSubGroup("Edges",m_p->m_mesh->edgeFamily(),new ItemItemGroupComputeFunctor<Edge>());
    m_p->m_edge_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
faceGroup()
{
  if (itemKind()==IK_Face)
    return this;
  ItemGroupImpl* ii = m_p->m_face_group;
  if (!ii){
    ii = createSubGroup("Faces",m_p->m_mesh->faceFamily(),new ItemItemGroupComputeFunctor<Face>());
    m_p->m_face_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
cellGroup()
{
  if (itemKind()==IK_Cell)
    return this;
  ItemGroupImpl* ii = m_p->m_cell_group;
  if (!ii){
    ii = createSubGroup("Cells",m_p->m_mesh->cellFamily(),new ItemItemGroupComputeFunctor<Cell>());
    m_p->m_cell_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
innerFaceGroup()
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_inner_face_group;
  if (!ii){
    ii = createSubGroup("InnerFaces",m_p->m_mesh->faceFamily(),
                        new InnerFaceItemGroupComputeFunctor());
    m_p->m_inner_face_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
outerFaceGroup()
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_outer_face_group;
  if (!ii){
    ii = createSubGroup("OuterFaces",m_p->m_mesh->faceFamily(),
                        new OuterFaceItemGroupComputeFunctor());
    m_p->m_outer_face_group = ii;
  }
  return ii;
}
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
//! AMR
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
activeCellGroup()
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_active_cell_group;
 
  if (!ii){
    ii = createSubGroup("ActiveCells",m_p->m_mesh->cellFamily(),
                        new ActiveCellGroupComputeFunctor());
    m_p->m_active_cell_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
ownActiveCellGroup()
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_own_active_cell_group;
  // Le flag est déjà positionné dans le ItemGroupInternal::_init ou ItemGroupImpl::setOwn
  if (!ii) {
    ii = createSubGroup("OwnActiveCells",m_p->m_mesh->cellFamily(),
                        new OwnActiveCellGroupComputeFunctor());
    m_p->m_own_active_cell_group = ii;
    ii->setOwn(true);
  }
  return ii;
}
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
levelCellGroup(const Integer& level)
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_level_cell_group[level];
  if (!ii){
    ii = createSubGroup(String::format("LevelCells{0}",level),
                        m_p->m_mesh->cellFamily(),
                        new LevelCellGroupComputeFunctor(level));
    m_p->m_level_cell_group[level] = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
ownLevelCellGroup(const Integer& level)
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_own_level_cell_group[level];
  // Le flag est déjà positionné dans le ItemGroupInternal::_init ou ItemGroupImpl::setOwn
  if (!ii) {
    ii = createSubGroup(String::format("OwnLevelCells{0}",level),
                        m_p->m_mesh->cellFamily(),
                        new OwnLevelCellGroupComputeFunctor(level));
    m_p->m_own_level_cell_group[level] = ii;
    ii->setOwn(true);
  }
  return ii;
}
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
activeFaceGroup()
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_active_face_group;
  if (!ii){
    ii = createSubGroup("ActiveFaces",m_p->m_mesh->faceFamily(),
                        new ActiveFaceItemGroupComputeFunctor());
    m_p->m_active_face_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
ownActiveFaceGroup()
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_own_active_face_group;
  if (!ii){
    ii = createSubGroup("OwnActiveFaces",m_p->m_mesh->faceFamily(),
                        new OwnActiveFaceItemGroupComputeFunctor());
    m_p->m_own_active_face_group = ii;
  }
  return ii;
}
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
innerActiveFaceGroup()
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_inner_active_face_group;
  if (!ii){
    ii = createSubGroup("InnerActiveFaces",m_p->m_mesh->faceFamily(),
                        new InnerActiveFaceItemGroupComputeFunctor());
    m_p->m_inner_active_face_group = ii;
  }
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
outerActiveFaceGroup()
{
  if (itemKind()!=IK_Cell)
    return checkSharedNull();
  ItemGroupImpl* ii = m_p->m_outer_active_face_group;
  if (!ii){
    ii = createSubGroup("OuterActiveFaces",m_p->m_mesh->faceFamily(),
                        new OuterActiveFaceItemGroupComputeFunctor());
    m_p->m_outer_active_face_group = ii;
  }
  return ii;
}
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
createSubGroup(const String& suffix, IItemFamily* family, ItemGroupComputeFunctor * functor)
{
  String sub_name = name() + "_" + suffix;
  auto finder = m_p->m_sub_groups.find(sub_name);
  if (finder != m_p->m_sub_groups.end()) {
    ARCANE_FATAL("Cannot create already existing sub-group ({0}) in group ({1})",
                 suffix, name());
  }
  ItemGroup ig = family->createGroup(sub_name,ItemGroup(this));
  ItemGroupImpl* ii = ig.internal();
  ii->setComputeFunctor(functor);
  functor->setGroup(ii);
  // Observer par défaut : le sous groupe n'est pas intéressé par les infos détaillées de transition
  attachObserver(ii,newItemGroupObserverT(ii,
                                          &ItemGroupImpl::_executeInvalidate));
  m_p->m_sub_groups[sub_name] = ii;
  ii->invalidate(false);
  return ii;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImpl* ItemGroupImpl::
findSubGroup(const String& suffix)
{
  String sub_name = name() + "_" + suffix;
  auto finder = m_p->m_sub_groups.find(sub_name);
  if (finder != m_p->m_sub_groups.end()) {
    return finder->second.get();
  }
  else {
    // ou bien erreur ?
    return checkSharedNull();
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
changeIds(Int32ConstArrayView old_to_new_ids)
{
  ITraceMng* trace = m_p->m_mesh->traceMng();
  if (m_p->m_compute_functor) {
    trace->debug(Trace::High) << "ItemGroupImpl::changeIds on "
                              << name() << " : skip computed group";
    return;
  }
 
  // ItemGroupImpl ne fait d'habitude pas le checkNeedUpdate lui meme, plutot le ItemGroup
  checkNeedUpdate();
  if (isAllItems()) {
    // Int32ArrayView items_lid = m_p->itemsLocalId();
    // for( Integer i=0, is=items_lid.size(); i<is; ++i ) {
    //   Integer lid = items_lid[i];
    //   ARCANE_ASSERT((lid = old_to_new_ids[lid]),("Non uniform order on allItems"));;
    // }
    // m_p->notifyCompactObservers(&old_to_new_ids);
    m_p->notifyInvalidateObservers();
    return;
  }
 
  Int32ArrayView items_lid = m_p->itemsLocalId();
  for( Integer i=0, is=items_lid.size(); i<is; ++i ){
    Integer old_id = items_lid[i];
    items_lid[i] = old_to_new_ids[old_id];
  }
 
  m_p->updateTimestamp();
  // Pour l'instant, il ne faut pas trier les entités des variables
  // partielles car les valeurs correspondantes des variables ne sont
  // pas mises à jour (il faut implémenter changeGroupIds pour cela)
  // NOTE: est-ce utile de le faire ?
  // NOTE: faut-il le faire pour les autre aussi ? A priori, cela n'est
  // utile que pour garantir le meme resultat parallele/sequentiel
  if (m_p->m_observer_need_info) {
    m_p->notifyCompactObservers(&old_to_new_ids);
  } else {
    // Pas besoin d'infos, on peut changer arbitrairement leur ordre
    // TODO: #warning "(HP) Connexion de ce triage d'item avec la famille ?"
    std::sort(std::begin(items_lid),std::end(items_lid));
    m_p->notifyCompactObservers(nullptr);
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
invalidate(bool force_recompute)
{
#ifdef ARCANE_DEBUG
  //   if (force_recompute){
  ITraceMng* msg = m_p->m_mesh->traceMng();
  msg->debug(Trace::High) << "ItemGroupImpl::invalidate(force=" << force_recompute << ")"
                          << " name=" << name();
//   }
#endif
 
  m_p->updateTimestamp();
  m_p->setNeedRecompute();
  if (force_recompute)
    checkNeedUpdate();
  m_p->notifyInvalidateObservers();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemInternalList ItemGroupImpl::
itemsInternal() const
{
  return m_p->items();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemInfoListView ItemGroupImpl::
itemInfoListView() const
{
  return m_p->itemInfoListView();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
isLocalToSubDomain() const
{
  return m_p->m_is_local_to_sub_domain;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
setLocalToSubDomain(bool v)
{
  m_p->m_is_local_to_sub_domain = v;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
addItems(Int32ConstArrayView items_local_id,bool check_if_present)
{
  ARCANE_ASSERT(( (!m_p->m_need_recompute && !isAllItems()) || (m_p->m_transaction_mode && isAllItems()) ),
    ("Operation on invalid group"));
  if (m_p->m_compute_functor && !m_p->m_transaction_mode)
    ARCANE_FATAL("Cannot add items on computed group ({0})", name());
  IMesh* amesh = mesh();
  if (!amesh)
    throw ArgumentException(A_FUNCINFO,"null group");
  ITraceMng* trace = amesh->traceMng();
 
  Integer nb_item_to_add = items_local_id.size();
  if (nb_item_to_add==0)
    return;
 
  Int32Array& items_lid = m_p->mutableItemsLocalId();
  const Integer current_size = items_lid.size();
  Integer nb_added = 0;
 
  if(isAllItems()) {
    // Ajoute les nouveaux items à la fin
    Integer nb_items_id = current_size;
    m_p->m_items_index_in_all_group.resize(m_p->maxLocalId());
    for( Integer i=0, is=nb_item_to_add; i<is; ++i ){
      Int32 local_id = items_local_id[i];
      items_lid.add(local_id);
      m_p->m_items_index_in_all_group[local_id] = nb_items_id+i;
    }
    nb_added = nb_item_to_add;
  }
  else if (check_if_present) {
    // Vérifie que les entités à ajouter ne sont pas déjà présentes
    UniqueArray<bool> presence_checks(m_p->maxLocalId());
    presence_checks.fill(false);
    for( Integer i=0, is=items_lid.size(); i<is; ++i ){
      presence_checks[items_lid[i]] = true;
    }
 
    for( Integer i=0; i<nb_item_to_add; ++i ) {
      const Integer lid = items_local_id[i];
      if (!presence_checks[lid]){
        items_lid.add(lid);
        // Met à \a true comme cela si l'entité est présente plusieurs fois
        // dans \a items_local_id cela fonctionne quand même.
        presence_checks[lid] = true;
        ++nb_added;
      }
    }
  }
  else {
    nb_added = nb_item_to_add;
    MemoryUtils::checkResizeArrayWithCapacity(items_lid,current_size+nb_added,false);
    SmallSpan<Int32> ptr = items_lid.subView(current_size,nb_added);
    MemoryUtils::copy<Int32>(ptr, items_local_id);
  }
 
  if (arcaneIsCheck()){
    trace->info(5) << "ItemGroupImpl::addItems() group <" << name() << "> "
                   << " checkpresent=" << check_if_present
                   << " nb_current=" << current_size
                   << " want_to_add=" << nb_item_to_add
                   << " effective_added=" << nb_added;
    checkValid();
  }
 
  if (nb_added!=0) {
    m_p->updateTimestamp();
    Int32ConstArrayView observation_info(nb_added, items_lid.unguardedBasePointer()+current_size);
    m_p->notifyExtendObservers(&observation_info);
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
removeItems(Int32ConstArrayView items_local_id,[[maybe_unused]] bool check_if_present)
{
  m_p->_removeItems(items_local_id);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
removeAddItems(Int32ConstArrayView removed_items_lids,
               Int32ConstArrayView added_items_lids,
               bool check_if_present)
{
  ARCANE_ASSERT(( (!m_p->m_need_recompute && !isAllItems()) || (m_p->m_transaction_mode && isAllItems()) ),
                ("Operation on invalid group"));
  if (m_p->m_compute_functor && !m_p->m_transaction_mode)
    ARCANE_FATAL("Cannot remove items on computed group ({0})", name());
  IMesh* amesh = mesh();
  if (!amesh)
    throw ArgumentException(A_FUNCINFO,"null group");
  ITraceMng* trace = amesh->traceMng();
  if (isOwn() && amesh->meshPartInfo().nbPart()!=1){
    ARCANE_THROW(NotSupportedException,"Cannot remove items if isOwn() is true");
  }
 
  Int32Array & items_lid = m_p->mutableItemsLocalId();
  
  if(isAllItems()) {
    ItemInternalList internals = itemsInternal();
    const Integer internal_size = internals.size();
    const Integer new_size = m_p->m_item_family->nbItem();
    items_lid.resize(new_size);
    m_p->m_items_index_in_all_group.resize(m_p->maxLocalId());
    if (new_size==internal_size){
      // Il n'y a pas de trous dans la numérotation
      for( Integer i=0; i<internal_size; ++i ){
        Int32 local_id = internals[i]->localId();
        items_lid[i] = local_id;
        m_p->m_items_index_in_all_group[local_id] = i;
      }
    }
    else{
      Integer index = 0;
      for( Integer i=0; i<internal_size; ++i ){
        ItemInternal* item = internals[i];
        if (!item->isSuppressed()){
          Int32 local_id = item->localId();
          items_lid[index] = local_id;
          m_p->m_items_index_in_all_group[local_id] = index;
          ++index;
        }
      }
      if (index!=new_size)
        ARCANE_FATAL("Inconsistent number of elements in the generation of the group '{0}' expected={1} present={2}",
                     name(), new_size, index);
    }
 
    // On ne peut pas savoir ce qui a changé donc dans le doute on
    // incrémente le timestamp.
    m_p->updateTimestamp();
  }
  else {
    removeItems(removed_items_lids,check_if_present);
    addItems(added_items_lids,check_if_present);
  }
  
  if(arcaneIsCheck()){
    trace->info(5) << "ItemGroupImpl::removeAddItems() group <" << name() << "> "
                   << " old_size=" << m_p->m_item_family->nbItem()
                   << " new_size=" << size()
                   << " nb_removed=" << removed_items_lids.size()
                   << " nb_added=" << added_items_lids.size();
    checkValid();
  }
}
 
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
setItems(Int32ConstArrayView items_local_id)
{
  if (m_p->m_compute_functor && !m_p->m_transaction_mode)
    ARCANE_FATAL("Cannot set items on computed group ({0})", name());
  Int32Array& buf = m_p->mutableItemsLocalId();
  buf.resize(items_local_id.size());
  Int32ArrayView buf_view(buf);
  buf_view.copy(items_local_id);
  m_p->updateTimestamp();
  m_p->m_need_recompute = false;
  if (arcaneIsCheck()){
    ITraceMng* trace = m_p->mesh()->traceMng();
    trace->debug(Trace::High) << "ItemGroupImpl::setItems() group <" << name() << "> "
                              << " size=" << size();
    checkValid();
  }
 
  // On tolére encore le setItems initial en le convertissant en un addItems
  if (size() != 0)
    m_p->notifyInvalidateObservers();
  else
    m_p->notifyExtendObservers(&items_local_id);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
class ItemGroupImpl::ItemSorter
{
 public:
  ItemSorter(ItemInfoListView items) : m_items(items){}
 public:
  void sort(Int32ArrayView local_ids) const
  {
    std::sort(std::begin(local_ids),std::end(local_ids),*this);
  }
  bool operator()(Int32 lid1,Int32 lid2) const
  {
    return m_items[lid1].uniqueId() < m_items[lid2].uniqueId();
  }
 private:
  ItemInfoListView m_items;
};
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
setItems(Int32ConstArrayView items_local_id,bool do_sort)
{
  if (!do_sort){
    setItems(items_local_id);
    return;
  }
  UniqueArray<Int32> sorted_lid(items_local_id);
  ItemSorter sorter(itemFamily()->itemInfoListView());
  sorter.sort(sorted_lid);
  setItems(sorted_lid);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
setIsAllItems()
{
  m_p->m_is_all_items = true;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
isAllItems() const
{
  return m_p->m_is_all_items;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
class ItemCheckSuppressedFunctor
{
 public:
  ItemCheckSuppressedFunctor(ItemInternalList items)
  : m_items(items)
    {
    }
 public:
  bool operator()(Integer item_lid)
    {
      return m_items[item_lid]->isSuppressed();
    }
 private:
  ItemInternalList m_items;
};
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
removeSuppressedItems()
{
  // TODO: (HP) Assertion need fix"
  // ARCANE_ASSERT((!m_p->m_need_recompute),("Operation on invalid group"));
  ITraceMng* trace = m_p->mesh()->traceMng();
  if (m_p->m_compute_functor) {
    trace->debug(Trace::High) << "ItemGroupImpl::removeSuppressedItems on " << name() << " : skip computed group";
    return;
  }
 
  ItemInternalList items = m_p->items();
  Integer nb_item = items.size();
  Int32Array& items_lid = m_p->mutableItemsLocalId();
  Integer current_size = items_lid.size();
  if (arcaneIsCheck()){
    for( Integer i=0; i<current_size; ++i ){
      if (items_lid[i]>=nb_item){
        trace->fatal() << "ItemGroupImpl::removeSuppressedItems(): bad range "
                       << " name=" << name()
                       << " i=" << i << " lid=" << items_lid[i]
                       << " max=" << nb_item;
      }
    }
  }
 
  Int32UniqueArray removed_ids, removed_lids;
  Int32ConstArrayView * observation_info  = NULL;
  Int32ConstArrayView * observation_info2 = NULL;
  Integer new_size;
  // Si le groupe posséde des observers ayant besoin d'info, il faut les calculer
  if (m_p->m_observer_need_info){
    removed_ids.reserve(current_size); // préparation à taille max
    Integer index = 0;
    for( Integer i=0; i<current_size; ++i ){
      if (!items[items_lid[i]]->isSuppressed()){
        items_lid[index] = items_lid[i];
        ++index;
      } else {
        // trace->debug(Trace::Highest) << "Remove from group " << name() << " item " << i << " " << items_lid[i] << " " << ItemPrinter(items[items_lid[i]]);
        removed_lids.add(items_lid[i]);
      }
    }
    new_size = index;
    if (new_size!=current_size){
      items_lid.resize(new_size);
      observation_info = new Int32ConstArrayView(removed_lids.size(), removed_lids.unguardedBasePointer());
    }
  }
  else{
    ItemCheckSuppressedFunctor f(m_p->items());
    auto ibegin = std::begin(items_lid);
    auto new_end = std::remove_if(ibegin,std::end(items_lid),f);
    new_size = (Integer)(new_end-ibegin);
    items_lid.resize(new_size);
  }
 
  if (arcaneIsCheck()) {
    trace->debug(Trace::High) << "ItemGroupImpl::removeSupressedItems() group <" << name()
                              << "> NEW SIZE=" << new_size << " OLD=" << current_size;
    checkValid();
  }
 
  // Ne met à jour que si le groupe a été modifié
  if (current_size != new_size) {
    m_p->updateTimestamp();
    m_p->notifyReduceObservers(observation_info);
    delete observation_info;
    delete observation_info2;
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
checkValid()
{
  m_p->checkValid();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
_checkNeedUpdate(bool do_padding)
{
  // En cas de problème sur des recalculs très imbriqués, une proposition est
  // de désactiver les lignes #A pour activer les lignes #B
  bool has_recompute = false;
  if (m_p->m_need_recompute) {
    m_p->m_need_recompute = false;
    // bool need_invalidate_on_recompute = m_p->m_need_invalidate_on_recompute; // #B
    // m_p->m_need_invalidate_on_recompute = false;                             // #B
    if (m_p->m_compute_functor) {
      m_p->m_compute_functor->executeFunctor();
    }
    //     if (need_invalidate_on_recompute) { // #B
    //       m_p->notifyInvalidateObservers(); // #B
    //     }                                   // #B
 
    if (m_p->m_need_invalidate_on_recompute) {     // #A
      m_p->m_need_invalidate_on_recompute = false; // #A
      m_p->notifyInvalidateObservers();            // #A
    }                                              // #A
    has_recompute = true;
  }
  if (do_padding)
    _checkUpdateSimdPadding();
  return has_recompute;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
_checkNeedUpdateNoPadding()
{
  return _checkNeedUpdate(false);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
_checkNeedUpdateWithPadding()
{
  return _checkNeedUpdate(true);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
checkNeedUpdate()
{
  return _checkNeedUpdate(false);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_checkUpdateSimdPadding()
{
  m_p->checkUpdateSimdPadding();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
clear()
{
  Int32Array& items_lid = m_p->mutableItemsLocalId();
  if (!items_lid.empty())
    // Incrémente uniquement si le groupe n'est pas déjà vide
    m_p->updateTimestamp();
  items_lid.clear();
  m_p->m_need_recompute = false;
  for( const auto& i : m_p->m_sub_groups )
    i.second->clear();
  m_p->notifyInvalidateObservers();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroup ItemGroupImpl::
parentGroup()
{
  if (m_p->m_parent)
    return ItemGroup(m_p->m_parent);
  return ItemGroup();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
applyOperation(IItemOperationByBasicType* operation)
{
  ARCANE_ASSERT((!m_p->m_need_recompute),("Operation on invalid group"));
  bool is_verbose = m_p->m_is_debug_apply_operation;
 
  ITraceMng* tm = m_p->mesh()->traceMng();
  if (is_verbose)
    tm->info() << "applyOperation name=" << name() << " nb_item=" << size();
  if (m_p->isUseV2ForApplyOperation()){
    if (m_p->m_children_by_type_ids.empty()){
      _initChildrenByTypeV2();
    }
    Int64 t = m_p->timestamp();
    if (is_verbose)
      tm->info() << "applyOperation timestamp=" << t << " last=" << m_p->m_children_by_type_ids_computed_timestamp;
    if (m_p->m_children_by_type_ids_computed_timestamp != t){
      _computeChildrenByTypeV2();
      m_p->m_children_by_type_ids_computed_timestamp = t;
    }
  }
  else{
    if (m_p->m_children_by_type.empty())
      _initChildrenByType();
  }
  IItemFamily* family = m_p->m_item_family;
  const bool has_only_one_type = (m_p->m_unique_children_type != IT_NullType);
  if (is_verbose)
    tm->info() << "applyOperation has_only_one_type=" << has_only_one_type << " value=" << m_p->m_unique_children_type;
  // TODO: Supprimer la macro et la remplacer par une fonction.
 
#define APPLY_OPERATION_ON_TYPE(ITEM_TYPE)      \
  if (m_p->isUseV2ForApplyOperation()){\
    Int16 type_id = IT_##ITEM_TYPE;\
    Int32ConstArrayView sub_ids = m_p->m_children_by_type_ids[type_id]; \
    if (has_only_one_type && type_id==m_p->m_unique_children_type)\
      sub_ids = itemsLocalId(); \
    if (is_verbose && sub_ids.size()>0)                                   \
      tm->info() << "Type=" << (int)IT_##ITEM_TYPE << " nb=" << sub_ids.size(); \
    if (sub_ids.size()!=0){\
      operation->apply##ITEM_TYPE(family->view(sub_ids)); \
    }\
  }\
  else {                                                     \
    ItemGroup group(m_p->m_children_by_type[IT_##ITEM_TYPE]); \
    if (!group.empty())                                       \
      operation->apply##ITEM_TYPE(group.view());              \
  }
 
  APPLY_OPERATION_ON_TYPE(Vertex);
  APPLY_OPERATION_ON_TYPE(Line2);
  APPLY_OPERATION_ON_TYPE(Triangle3);
  APPLY_OPERATION_ON_TYPE(Quad4);
  APPLY_OPERATION_ON_TYPE(Pentagon5);
  APPLY_OPERATION_ON_TYPE(Hexagon6);
  APPLY_OPERATION_ON_TYPE(Tetraedron4);
  APPLY_OPERATION_ON_TYPE(Pyramid5);
  APPLY_OPERATION_ON_TYPE(Pentaedron6);
  APPLY_OPERATION_ON_TYPE(Hexaedron8);
  APPLY_OPERATION_ON_TYPE(Heptaedron10);
  APPLY_OPERATION_ON_TYPE(Octaedron12);
  APPLY_OPERATION_ON_TYPE(HemiHexa7);
  APPLY_OPERATION_ON_TYPE(HemiHexa6);
  APPLY_OPERATION_ON_TYPE(HemiHexa5);
  APPLY_OPERATION_ON_TYPE(HemiHexa7);
  APPLY_OPERATION_ON_TYPE(AntiWedgeLeft6);
  APPLY_OPERATION_ON_TYPE(AntiWedgeRight6);
  APPLY_OPERATION_ON_TYPE(DiTetra5);
  APPLY_OPERATION_ON_TYPE(DualNode);
  APPLY_OPERATION_ON_TYPE(DualEdge);
  APPLY_OPERATION_ON_TYPE(DualFace);
  APPLY_OPERATION_ON_TYPE(DualCell);
  APPLY_OPERATION_ON_TYPE(Link);
  
#undef APPLY_OPERATION_ON_TYPE
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
needSynchronization() const
{
  return !(m_p->m_compute_functor != 0 ||
           m_p->m_is_local_to_sub_domain ||
           m_p->m_is_own);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
Int64 ItemGroupImpl::
timestamp() const
{
  return m_p->timestamp();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
attachObserver(const void * ref, IItemGroupObserver * obs)
{
  auto finder = m_p->m_observers.find(ref);
  auto end = m_p->m_observers.end();
  if (finder != end) {
    delete finder->second;
    finder->second = obs;
    return;
  }
  else {
    m_p->m_observers[ref] = obs;
  }
 
  // Mise à jour du flag de demande d'info
  _updateNeedInfoFlag(m_p->m_observer_need_info | obs->needInfo());
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
detachObserver(const void * ref)
{
  auto finder = m_p->m_observers.find(ref);
  auto end = m_p->m_observers.end();
 
  if (finder == end)
    return;
 
  IItemGroupObserver * obs = finder->second;
  delete obs;
  m_p->m_observers.erase(finder);
  // Mise à jour du flag de demande d'info
  bool new_observer_need_info = false;
  auto i = m_p->m_observers.begin();
  for( ; i != end ; ++i ) {
    IItemGroupObserver * obs = i->second;
    new_observer_need_info |= obs->needInfo();
  }
  _updateNeedInfoFlag(new_observer_need_info);
 
  // On invalide la table de hachage éventuelle
  // des variables partielles si il n'y a plus
  // de référence. 
  if(m_p->m_group_index_table.isUsed() && m_p->m_group_index_table.isUnique()) {
    m_p->m_group_index_table.reset();
    m_p->m_synchronizer.reset();
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
hasInfoObserver() const
{
  return m_p->m_observer_need_info;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
setComputeFunctor(IFunctor* functor)
{
  delete m_p->m_compute_functor;
  m_p->m_compute_functor = functor;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
hasComputeFunctor() const
{
  return (m_p->m_compute_functor);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_initChildrenByType()
{
  IItemFamily* family = m_p->m_item_family;
  ItemTypeMng* type_mng = family->mesh()->itemTypeMng();
  Integer nb_basic_item_type= ItemTypeMng::nbBasicItemType(); //NB_BASIC_ITEM_TYPE
  m_p->m_children_by_type.resize(nb_basic_item_type);
  for( Integer i=0; i<nb_basic_item_type; ++i ){
    // String child_name(name()+i);
    String child_name(type_mng->typeName(i));
    ItemGroupImpl* igi = createSubGroup(child_name,family,
                                        new ItemGroupImplItemGroupComputeFunctor(this,&ItemGroupImpl::_computeChildrenByType));
    m_p->m_children_by_type[i] = igi;
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_initChildrenByTypeV2()
{
  bool is_verbose = m_p->m_is_debug_apply_operation;
  if (is_verbose)
    m_p->mesh()->traceMng()->info() << "ItemGroupImpl::_initChildrenByTypeV2() name=" << name();
 
  Int32 nb_basic_item_type= ItemTypeMng::nbBasicItemType();
  m_p->m_children_by_type_ids.resize(nb_basic_item_type);
  for( Integer i=0; i<nb_basic_item_type; ++i ){
    m_p->m_children_by_type_ids[i] = UniqueArray<Int32>{MemoryUtils::getDefaultDataAllocator()};
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_computeChildrenByType()
{
  ItemGroup that_group(this);
  ITraceMng * trace = that_group.mesh()->traceMng();
  trace->debug(Trace::High) << "ItemGroupImpl::_computeChildrenByType for " << name();
 
  Integer nb_basic_item_type = ItemTypeMng::nbBasicItemType();
 
  UniqueArray< SharedArray<Int32> > items_by_type(nb_basic_item_type);
  for( Integer i=0; i<nb_basic_item_type; ++i ){
    ItemGroupImpl * impl = m_p->m_children_by_type[i];
    impl->beginTransaction();
  }
 
  ENUMERATE_ITEM(iitem,that_group){
    const Item& item = *iitem;
    Integer item_type = item.type();
    if (item_type<nb_basic_item_type)
      items_by_type[item_type].add(iitem.itemLocalId());
  }
 
  for( Integer i=0; i<nb_basic_item_type; ++i ){
    ItemGroupImpl * impl = m_p->m_children_by_type[i];
    impl->setItems(items_by_type[i]);
    impl->endTransaction();
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_computeChildrenByTypeV2()
{
  ItemGroup that_group(this);
  Int32 nb_item = size();
  IMesh* mesh = m_p->mesh();
  ItemTypeMng* type_mng = mesh->itemTypeMng();
  ITraceMng* trace = mesh->traceMng();
  bool is_verbose = m_p->m_is_debug_apply_operation;
  if (is_verbose)
    trace->info() << "ItemGroupImpl::_computeChildrenByTypeV2 for " << name();
 
  // Si le maillage est cartésien, on sait que toutes les entités ont le même type
  if (nb_item>0 && mesh->meshKind().meshStructure()==eMeshStructure::Cartesian){
    ItemInfoListView lv(m_p->m_item_family->itemInfoListView());
    m_p->m_unique_children_type = ItemTypeId{lv.typeId(m_p->itemsLocalId()[0])};
    return;
  }
 
  Int32 nb_basic_item_type = ItemTypeMng::nbBasicItemType();
  m_p->m_unique_children_type = ItemTypeId{IT_NullType};
 
  UniqueArray<Int32> nb_items_by_type(nb_basic_item_type);
  nb_items_by_type.fill(0);
  ENUMERATE_(Item,iitem,that_group){
    Item item = *iitem;
    Int16 item_type = item.type();
    if (item_type<nb_basic_item_type)
      ++nb_items_by_type[item_type];
  }
 
  Int32 nb_different_type = 0;
  for( Int32 i=0; i<nb_basic_item_type; ++i ){
    m_p->m_children_by_type_ids[i].clear();
    const Int32 n = nb_items_by_type[i];
    m_p->m_children_by_type_ids[i].reserve(n);
    if (n>0)
      ++nb_different_type;
    if (is_verbose)
      trace->info() << "ItemGroupImpl::_computeChildrenByTypeV2 for " << name()
                    << " type=" << type_mng->typeName(i) << " nb=" << n;
  }
  if (is_verbose)
    trace->info() << "ItemGroupImpl::_computeChildrenByTypeV2 for " << name()
                  << " nb_item=" << nb_item << " nb_different_type=" << nb_different_type;
 
  // Si nb_different_type == 1, cela signifie qu'il n'y a qu'un seul
  // type d'entité et on conserve juste ce type car dans ce cas on passera
  // directement le groupe en argument de applyOperation().
  if (nb_item>0 && nb_different_type==1){
    ItemInfoListView lv(m_p->m_item_family->itemInfoListView());
    m_p->m_unique_children_type = ItemTypeId{lv.typeId(m_p->itemsLocalId()[0])};
    if (is_verbose)
      trace->info() << "ItemGroupImpl::_computeChildrenByTypeV2 for " << name()
                    << " unique_type=" << type_mng->typeName(m_p->m_unique_children_type);
    return;
  }
 
  ENUMERATE_(Item,iitem,that_group){
    Item item = *iitem;
    Integer item_type = item.type();
    if (item_type<nb_basic_item_type)
      m_p->m_children_by_type_ids[item_type].add(iitem.itemLocalId());
  }
 
  for( Int32 i=0; i<nb_basic_item_type; ++i )
    applySimdPadding(m_p->m_children_by_type_ids[i]);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_executeExtend(const Int32ConstArrayView* info)
{
  ARCANE_UNUSED(info);
  // On ne sait encore calculer appliquer des transformations à des groupes calculés
  // On choisit l'invalidation systématique
  m_p->notifyInvalidateObservers();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_executeReduce(const Int32ConstArrayView* info)
{
  ARCANE_UNUSED(info);
  // On ne sait encore calculer appliquer des transformations à des groupes calculés
  // On choisit l'invalidation systématique
  m_p->notifyInvalidateObservers();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_executeCompact(const Int32ConstArrayView* info)
{
  ARCANE_UNUSED(info);
  // On ne sait encore calculer appliquer des transformations à des groupes calculés
  // On choisit l'invalidation systématique en différé (évaluée lors du prochain checkNeedUpdate)
  m_p->notifyInvalidateObservers();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_executeInvalidate()
{
  m_p->setNeedRecompute();
  m_p->notifyInvalidateObservers();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_updateNeedInfoFlag(const bool flag)
{
  if (m_p->m_observer_need_info == flag)
    return;
  m_p->m_observer_need_info = flag;
  // Si change, change aussi l'observer du parent pour qu'il adapte les besoins en info de transition
  ItemGroupImpl * parent = m_p->m_parent;
  if (parent) {
    parent->detachObserver(this);
    if (m_p->m_observer_need_info) {
      parent->attachObserver(this,newItemGroupObserverT(this,
                                                        &ItemGroupImpl::_executeExtend,
                                                        &ItemGroupImpl::_executeReduce,
                                                        &ItemGroupImpl::_executeCompact,
                                                        &ItemGroupImpl::_executeInvalidate));
    } else {
      parent->attachObserver(this,newItemGroupObserverT(this,
                                                        &ItemGroupImpl::_executeInvalidate));
    }
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_forceInvalidate(const bool self_invalidate)
{
  // (HP) TODO: Mettre un observer forceInvalidate pour prévenir tout le monde ?
  // avec forceInvalidate on doit invalider mais ne rien calculer
  if (self_invalidate) {
    m_p->setNeedRecompute();
    m_p->m_need_invalidate_on_recompute = true;
  }
 
  for( const auto& i :  m_p->m_sub_groups )
    i.second->_forceInvalidate(true);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
destroy()
{
  // Détache les observateurs. Cela modifie \a m_observers donc il faut
  // en faire une copie
  {
    std::vector<const void*> ptrs;
    for( const auto& i : m_p->m_observers )
      ptrs.push_back(i.first);
    for( const void* i : ptrs )
      detachObserver(i);
  }
 
  // Le groupe de toutes les entités est spécial. Il ne faut jamais le détruire
  // vraiement.
  if (m_p->m_is_all_items)
    m_p->resetSubGroups();
  else{
    delete m_p;
    m_p = new ItemGroupInternal();
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
SharedPtrT<GroupIndexTable> ItemGroupImpl::
localIdToIndex()
{ 
  if(!m_p->m_group_index_table.isUsed()) {
    m_p->m_group_index_table = SharedPtrT<GroupIndexTable>(new GroupIndexTable(this));
    ITraceMng* trace =  m_p->m_mesh->traceMng();
    trace->debug(Trace::High) << "** CREATION OF LOCAL ID TO INDEX TABLE OF GROUP : " << m_p->m_name;
    m_p->m_group_index_table->update();
  }
  return m_p->m_group_index_table;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
IVariableSynchronizer * ItemGroupImpl::
synchronizer()
{ 
  if(!m_p->m_synchronizer.get()) {
    IParallelMng* pm = m_p->m_mesh->parallelMng();
    ItemGroup this_group(this);
    m_p->m_synchronizer = ParallelMngUtils::createSynchronizerRef(pm,this_group);
    ITraceMng* trace =  m_p->m_mesh->traceMng();
    trace->debug(Trace::High) << "** CREATION OF SYNCHRONIZER OF GROUP : " << m_p->m_name;
    m_p->m_synchronizer->compute();
  }
  return m_p->m_synchronizer.get();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
hasSynchronizer()
{
  return m_p->m_synchronizer.get();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
checkIsSorted() const
{
  if (null())
    return true;
 
  // TODO: stocker l'info dans un flag et ne refaire le test que si
  // la liste des entités a changer (utiliser timestamp()).
  ItemInternalList items(m_p->items());
  Int32ConstArrayView items_lid(m_p->itemsLocalId());
  Integer nb_item = items_lid.size();
  // On est toujours trié si une seule entité ou moins.
  if (nb_item<=1)
    return true;
  // Compare chaque uniqueId() avec le précédent et vérifie qu'il est
  // supérieur.
  ItemUniqueId last_uid = items[items_lid[0]]->uniqueId();
  for( Integer i=1; i<nb_item; ++i ){
    ItemUniqueId uid = items[items_lid[i]]->uniqueId();
    if (uid<last_uid)
      return false;
    last_uid = uid;
  }
  return true;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
deleteMe()
{
  // S'il s'agit de 'shared_null'. Il faut le positionner à nullptr pour qu'il
  // soit réalloué éventuellement par _buildSharedNull().
  if (this==shared_null){
    shared_null = nullptr;
  }
  delete this;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_buildSharedNull()
{
  if (!shared_null){
    shared_null = new ItemGroupImplNull();
    shared_null->addRef();
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
_destroySharedNull()
{
  // Supprime une référence à 'shared_null'. Si le nombre de référence
  // devient égal à 0, alors l'instance 'shared_null' est détruite.
  if (shared_null)
    shared_null->removeRef();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
bool ItemGroupImpl::
isContigousLocalIds() const
{
  return m_p->isContigous();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
checkLocalIdsAreContigous() const
{
  m_p->checkIsContigous();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
Int64 ItemGroupImpl::
capacity() const
{
  Int32Array& items_lid = m_p->mutableItemsLocalId();
  return items_lid.capacity();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
void ItemGroupImpl::
shrinkMemory()
{
  if (hasComputeFunctor()){
    // Groupe calculé. On l'invalide et on supprime ses éléments
    invalidate(false);
    m_p->mutableItemsLocalId().clear();
  }
 
  if (m_p->variableItemsLocalid())
    m_p->variableItemsLocalid()->variable()->shrinkMemory();
  else
    m_p->mutableItemsLocalId().shrink();
 
  m_p->applySimdPadding();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ItemGroupImplInternal* ItemGroupImpl::
_internalApi() const
{
  return &m_p->m_internal_api;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
} // End namespace Arcane
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/