MeshMaterialVariableScalar.cc

// -*- tab-width: 2; indent-tabs-mode: nil; coding: utf-8-with-signature -*-
//-----------------------------------------------------------------------------
// Copyright 2000-2026 CEA (www.cea.fr) IFPEN (www.ifpenergiesnouvelles.com)
// See the top-level COPYRIGHT file for details.
// SPDX-License-Identifier: Apache-2.0
//-----------------------------------------------------------------------------
/*---------------------------------------------------------------------------*/
/* MeshMaterialVariableScalar.cc                               (C) 2000-2024 */
/*                                                                           */
/* Scalar variable on a mesh material.                                       */
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#include "arcane/materials/MeshMaterialVariable.h"
 
#include "arcane/utils/NotImplementedException.h"
#include "arcane/utils/ScopedPtr.h"
#include "arcane/utils/ITraceMng.h"
#include "arcane/utils/NumericTypes.h"
#include "arcane/utils/CheckedConvert.h"
 
#include "arcane/materials/MaterialVariableBuildInfo.h"
#include "arcane/materials/IMeshMaterial.h"
#include "arcane/materials/MatItemEnumerator.h"
#include "arcane/materials/MeshMaterialVariableSynchronizerList.h"
#include "arcane/materials/ItemMaterialVariableBaseT.H"
#include "arcane/materials/IMeshMaterialVariableSynchronizer.h"
#include "arcane/materials/internal/MeshMaterialVariablePrivate.h"
 
#include "arcane/core/IItemFamily.h"
#include "arcane/core/IMesh.h"
#include "arcane/core/IParallelMng.h"
#include "arcane/core/ISubDomain.h"
#include "arcane/core/IApplication.h"
#include "arcane/core/IDataFactoryMng.h"
#include "arcane/core/IParallelExchanger.h"
#include "arcane/core/ISerializer.h"
#include "arcane/core/ISerializeMessage.h"
#include "arcane/core/internal/IVariableInternal.h"
#include "arcane/core/materials/internal/IMeshComponentInternal.h"
#include "arcane/core/VariableInfo.h"
#include "arcane/core/VariableRefArray.h"
#include "arcane/core/MeshVariable.h"
#include "arcane/core/VariableArray.h"
#include "arcane/core/ItemPrinter.h"
#include "arcane/core/IVariableSynchronizer.h"
#include "arcane/core/internal/IDataInternal.h"
#include "arcane/core/Timer.h"
#include "arcane/core/parallel/IStat.h"
#include "arcane/core/datatype/DataTypeTraits.h"
#include "arcane/core/datatype/DataStorageBuildInfo.h"
 
#include <vector>
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
namespace Arcane::Materials
{
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
MaterialVariableScalarTraits<DataType>::
saveData(IMeshComponent* component, IData* data,
         Array<ContainerViewType>& cviews)
{
  ConstArrayView<ArrayView<DataType>> views = cviews;
  auto* true_data = dynamic_cast<ValueDataType*>(data);
  ARCANE_CHECK_POINTER(true_data);
  ContainerType& values = true_data->_internal()->_internalDeprecatedValue();
  ENUMERATE_COMPONENTCELL (icell, component) {
    MatVarIndex mvi = icell._varIndex();
    values.add(views[mvi.arrayIndex()][mvi.valueIndex()]);
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
MaterialVariableScalarTraits<DataType>::
copyTo(SmallSpan<const DataType> input, SmallSpan<const Int32> input_indexes,
       SmallSpan<DataType> output, SmallSpan<const Int32> output_indexes,
       const RunQueue& queue)
{
  // TODO: check identical index sizes
  Integer nb_value = input_indexes.size();
  auto command = makeCommand(queue);
  ARCANE_CHECK_ACCESSIBLE_POINTER(queue, output.data());
  ARCANE_CHECK_ACCESSIBLE_POINTER(queue, input.data());
  ARCANE_CHECK_ACCESSIBLE_POINTER(queue, input_indexes.data());
  ARCANE_CHECK_ACCESSIBLE_POINTER(queue, output_indexes.data());
  command << RUNCOMMAND_LOOP1(iter, nb_value)
  {
    auto [i] = iter();
    output[output_indexes[i]] = input[input_indexes[i]];
  };
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
MaterialVariableScalarTraits<DataType>::
resizeAndFillWithDefault(ValueDataType* data, ContainerType& container, Integer dim1_size)
{
  ARCANE_UNUSED(data);
  // TODO: make an Array2 version that specifies a value to give
  // for initialization during a resize() (like for Array::resize()).
  container.resize(dim1_size, DataType());
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
MaterialVariableScalarTraits<DataType>::
resizeWithReserve(PrivatePartType* var, Integer dim1_size, Real reserve_ratio)
{
  // To avoid reallocating every time there is an increase in
  // the number of material cells, allocate a little more than necessary.
  // By default, we allocate 5% more.
  Int32 nb_add = static_cast<Int32>(dim1_size * reserve_ratio);
  var->_internalApi()->resize(VariableResizeArgs(dim1_size, nb_add, true));
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> ItemMaterialVariableScalar<DataType>::
ItemMaterialVariableScalar(const MaterialVariableBuildInfo& v,
                           PrivatePartType* global_var,
                           VariableRef* global_var_ref, MatVarSpace mvs)
: BaseClass(v, global_var, global_var_ref, mvs)
{
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*!
 * \brief Fills the variable values for a material from an array.
 *
 * This method performs the following operation:
 \code
 * Integer index=0;
 * ENUMERATE_MATCELL(imatcell,mat){
 *  matvar[imatcell] = values[index];
 *  ++index;
 * }
 \endcode
*/
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
fillFromArray(IMeshMaterial* mat, ConstArrayView<DataType> values)
{
  // TODO: make a version with IMeshComponent
  Integer index = 0;
  ENUMERATE_COMPONENTITEM (MatCell, imatcell, mat) {
    MatCell mc = *imatcell;
    MatVarIndex mvi = mc._varIndex();
    setValue(mvi, values[index]);
    ++index;
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*!
 * \brief Fills the variable values for a material from an array.
 *
 * This method performs the following operation:
 \code
 * Integer index=0;
 * ENUMERATE_MATCELL(imatcell,mat){
 *  matvar[imatcell] = values[index];
 *  ++index;
 * }
 \endcode
*/
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
fillFromArray(IMeshMaterial* mat, ConstArrayView<DataType> values,
              Int32ConstArrayView indexes)
{
  ConstArrayView<MatVarIndex> mat_indexes = mat->_internalApi()->variableIndexer()->matvarIndexes();
  Integer nb_index = indexes.size();
  for (Integer i = 0; i < nb_index; ++i) {
    MatVarIndex mvi = mat_indexes[indexes[i]];
    setValue(mvi, values[i]);
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*!
 * \brief Fills an array from the variable values for a material.
 *
 * This method performs the following operation:
 \code
 * Integer index=0;
 * ENUMERATE_MATCELL(imatcell,mat){
 *  values[index] = matvar[imatcell];
 *  ++index;
 * }
 \endcode
*/
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
fillToArray(IMeshMaterial* mat, ArrayView<DataType> values)
{
  Integer index = 0;
  ENUMERATE_COMPONENTITEM (MatCell, imatcell, mat) {
    MatCell mc = *imatcell;
    MatVarIndex mvi = mc._varIndex();
    values[index] = this->operator[](mvi);
    ++index;
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*!
 * \brief Fills an array from the variable values for a material.
 *
 * This method performs the following operation:
 \code
 * Integer index=0;
 * ENUMERATE_MATCELL(imatcell,mat){
 *  values[index] = matvar[imatcell];
 *  ++index;
 * }
 \endcode
*/
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
fillToArray(IMeshMaterial* mat, ArrayView<DataType> values, Int32ConstArrayView indexes)
{
  ConstArrayView<MatVarIndex> mat_indexes = mat->_internalApi()->variableIndexer()->matvarIndexes();
  Integer nb_index = indexes.size();
  for (Integer i = 0; i < nb_index; ++i) {
    MatVarIndex mvi = mat_indexes[indexes[i]];
    values[i] = this->operator[](mvi);
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*!
 * \brief Fills the partial values with the value \a value.
 */
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
fillPartialValues(const DataType& value)
{
  // The index 0 variable corresponding to the global variable.
  // Therefore, it should not be taken into account.
  Integer nb_var = m_vars.size();
  for (Integer i = 1; i < nb_var; ++i) {
    if (m_vars[i])
      m_vars[i]->fill(value);
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> Int32
ItemMaterialVariableScalar<DataType>::
dataTypeSize() const
{
  return (Int32)sizeof(DataType);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
_copyToBufferLegacy(SmallSpan<const MatVarIndex> matvar_indexes, Span<std::byte> bytes) const
{
  // TODO: Check that the size is a multiple of sizeof(DataType) and that
  // the alignment is correct.
  const Integer value_size = arcaneCheckArraySize(bytes.size() / sizeof(DataType));
  ArrayView<DataType> values(value_size, reinterpret_cast<DataType*>(bytes.data()));
  for (Integer z = 0; z < value_size; ++z) {
    values[z] = this->operator[](matvar_indexes[z]);
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
copyToBuffer(ConstArrayView<MatVarIndex> matvar_indexes, ByteArrayView bytes) const
{
  auto* ptr = reinterpret_cast<std::byte*>(bytes.data());
  return _copyToBufferLegacy(matvar_indexes, { ptr, bytes.size() });
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
_copyFromBufferLegacy(SmallSpan<const MatVarIndex> matvar_indexes, Span<const std::byte> bytes)
{
  // TODO: Check that the size is a multiple of sizeof(DataType) and that
  // the alignment is correct.
  const Int32 value_size = CheckedConvert::toInt32(bytes.size() / sizeof(DataType));
  ConstArrayView<DataType> values(value_size, reinterpret_cast<const DataType*>(bytes.data()));
  for (Integer z = 0; z < value_size; ++z) {
    setValue(matvar_indexes[z], values[z]);
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
copyFromBuffer(ConstArrayView<MatVarIndex> matvar_indexes, ByteConstArrayView bytes)
{
  auto* ptr = reinterpret_cast<const std::byte*>(bytes.data());
  return _copyFromBufferLegacy(matvar_indexes, { ptr, bytes.size() });
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
synchronize()
{
  IParallelMng* pm = m_p->materialMng()->mesh()->parallelMng();
  Timer timer(pm->timerMng(), "MatTimer", Timer::TimerReal);
  Int64 message_size = 0;
  {
    Timer::Sentry ts(&timer);
    message_size = _synchronize2();
  }
  pm->stat()->add("MaterialSync", timer.lastActivationTime(), message_size);
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> Int64
ItemMaterialVariableScalar<DataType>::
_synchronize2()
{
  Int64 message_size = 0;
  Integer sync_version = m_p->materialMng()->synchronizeVariableVersion();
  // Only versions 6 and later are available for medium variables.
  if (m_p->space() == MatVarSpace::Environment) {
    if (sync_version < 6)
      sync_version = 6;
  }
  if (sync_version >= 6) {
    MeshMaterialVariableSynchronizerList mmvsl(m_p->materialMng());
    mmvsl.add(this);
    mmvsl.apply();
    message_size = mmvsl.totalMessageSize();
  }
  else if (sync_version == 5 || sync_version == 4 || sync_version == 3) {
    _synchronizeV5();
  }
  else if (sync_version == 2) {
    _synchronizeV2();
  }
  else
    _synchronizeV1();
  return message_size;
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
synchronize(MeshMaterialVariableSynchronizerList& sync_list)
{
  Integer sync_version = m_p->materialMng()->synchronizeVariableVersion();
  if (sync_version >= 6) {
    sync_list.add(this);
  }
  else
    this->synchronize();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
_synchronizeV1()
{
  // Synchronization:
  // For now, an unoptimized algorithm that performs several
  // syncs. We use the global variable for this and must
  // therefore save it at the beginning of the function and restore it at the end.
  // The principle is simple: for each material and medium, we
  // copy the corresponding partial value into the global variable
  // and perform a sync of the global variable. We then
  // copy from the global value into the partial value.
  // This means we have as many syncs as there are materials and media
 
  // TODO: check that the list of materials is consistent between
  // the subdomains. This is not necessary with this simplified algorithm
  // but it will be with the final algorithm. For this test, it is enough to put
  // a special value in the global variable for each synchronization
  // (for example MIN_FLOAT) and after synchronization, look at each
  // cell whose value is different from MIN_FLOAT
  // if it is in our material
  IMeshMaterialMng* material_mng = m_p->materialMng();
  IMesh* mesh = material_mng->mesh();
  //TODO: Use a type other than cellFamily() to allow for other kinds
  // of elements.
  IItemFamily* family = mesh->cellFamily();
  IParallelMng* pm = mesh->parallelMng();
  if (!pm->isParallel())
    return;
  ItemGroup all_items = family->allItems();
 
  UniqueArray<DataType> saved_values;
  // Save the values of the global variable because we are going to change them
  {
    ConstArrayView<DataType> var_values = m_global_variable->valueView();
    saved_values.resize(var_values.size());
    saved_values.copy(var_values);
  }
 
  ENUMERATE_ENV (ienv, material_mng) {
    IMeshEnvironment* env = *ienv;
    ENUMERATE_MAT (imat, env) {
      IMeshMaterial* mat = *imat;
 
      {
        ArrayView<DataType> var_values = m_global_variable->valueView();
        // Copy material values into the global variable then synchronize
        ENUMERATE_MATCELL (imatcell, mat) {
          MatCell mc = *imatcell;
          Cell c = mc.globalCell();
          MatVarIndex mvi = mc._varIndex();
          var_values[c.localId()] = this->operator[](mvi);
        }
      }
      m_global_variable->synchronize();
      {
        ConstArrayView<DataType> var_values = m_global_variable->valueView();
        // Copy values from the global variable to the material part
        ENUMERATE_MATCELL (imatcell, mat) {
          MatCell mc = *imatcell;
          Cell c = mc.globalCell();
          MatVarIndex mvi = mc._varIndex();
          setValue(mvi, var_values[c.localId()]);
        }
      }
    }
 
    // Do the same for the environment
    {
      ArrayView<DataType> var_values = m_global_variable->valueView();
      // Copy environment values into the global variable then synchronize
      ENUMERATE_ENVCELL (ienvcell, env) {
        EnvCell ec = *ienvcell;
        Cell c = ec.globalCell();
        MatVarIndex mvi = ec._varIndex();
        var_values[c.localId()] = this->operator[](mvi);
      }
    }
    m_global_variable->synchronize();
    {
      ConstArrayView<DataType> var_values = m_global_variable->valueView();
      // Copy values from the global variable to the environment part
      ENUMERATE_ENVCELL (ienvcell, env) {
        EnvCell ec = *ienvcell;
        Cell c = ec.globalCell();
        MatVarIndex mvi = ec._varIndex();
        setValue(mvi, var_values[c.localId()]);
      }
    }
  }
 
  // Restore the values of the global variable.
  {
    ArrayView<DataType> var_values = m_global_variable->valueView();
    var_values.copy(saved_values);
  }
  m_global_variable->synchronize();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
_synchronizeV2()
{
  // Synchronization:
  // For now, non-optimized algorithm that performs several
  // synchronizations. We use the global variable for this and
  // therefore must save it at the beginning of the function and restore it at the end.
  // The principle is simple: for each material and environment, we
  // copy the corresponding partial value into the global variable
  // and perform a synchronization of the global variable. We then
  // copy from the global value back into the partial value.
  // This means we have as many synchronizations as there are materials and environments
 
  // TODO: check that the list of materials is consistent between
  // the subdomains. This is not necessary with this simplified algorithm
  // but it will be with the final algorithm. For this test, it is enough to put
  // a special value in the global variable for each synchronization
  // (for example MIN_FLOAT) and after synchronization, look at each
  // cell whose value is different from MIN_FLOAT
  // if it is in our material
  IMeshMaterialMng* material_mng = m_p->materialMng();
  IMesh* mesh = material_mng->mesh();
  //TODO: Use a type other than cellFamily() to allow for other kinds
  // of elements.
  IItemFamily* family = mesh->cellFamily();
  IParallelMng* pm = mesh->parallelMng();
  if (!pm->isParallel())
    return;
  ItemGroup all_items = family->allItems();
  ITraceMng* tm = pm->traceMng();
  tm->info(4) << "MAT_SYNCHRONIZE_V2 name=" << this->name();
  IDataFactoryMng* df = m_global_variable->dataFactoryMng();
 
  DataStorageTypeInfo storage_type_info(VariableInfo::_internalGetStorageTypeInfo(m_global_variable->dataType(), 2, 0));
  DataStorageBuildInfo storage_build_info(tm);
  String storage_full_type = storage_type_info.fullName();
 
  Ref<IData> xdata(df->createSimpleDataRef(storage_full_type, storage_build_info));
  auto* data = dynamic_cast<IArray2DataT<DataType>*>(xdata.get());
  if (!data)
    ARCANE_FATAL("Bad type");
 
  ConstArrayView<MeshMaterialVariableIndexer*> indexers = material_mng->_internalApi()->variablesIndexer();
  Integer nb_indexer = indexers.size();
  data->_internal()->_internalDeprecatedValue().resize(family->maxLocalId(), nb_indexer + 1);
  Array2View<DataType> values(data->view());
 
  // Copy partial values into the array.
  for (MeshMaterialVariableIndexer* indexer : indexers) {
    ConstArrayView<MatVarIndex> matvar_indexes = indexer->matvarIndexes();
    ConstArrayView<Int32> local_ids = indexer->localIds();
    for (Integer j = 0, n = matvar_indexes.size(); j < n; ++j) {
      MatVarIndex mvi = matvar_indexes[j];
      values[local_ids[j]][mvi.arrayIndex()] = this->operator[](mvi);
    }
  }
  family->allItemsSynchronizer()->synchronizeData(data);
 
  // Copy the synchronized array back into the partial values.
  for (MeshMaterialVariableIndexer* indexer : indexers) {
    ConstArrayView<MatVarIndex> matvar_indexes = indexer->matvarIndexes();
    ConstArrayView<Int32> local_ids = indexer->localIds();
    for (Integer j = 0, n = matvar_indexes.size(); j < n; ++j) {
      MatVarIndex mvi = matvar_indexes[j];
      setValue(mvi, values[local_ids[j]][mvi.arrayIndex()]);
    }
  }
 
  m_global_variable->synchronize();
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
_synchronizeV5()
{
  // Synchronization version that sends only
  // the values of materials and environments for shared cells.
  // NOTE: This version requires that the materials are correctly
  // synchronized between the subdomains.
 
  // This version is similar to V4 in principle but
  // performs send/receive directly without using serialization.
  IMeshMaterialMng* material_mng = m_p->materialMng();
 
  IMeshMaterialVariableSynchronizer* mmvs = material_mng->_internalApi()->allCellsMatEnvSynchronizer();
  IVariableSynchronizer* var_syncer = mmvs->variableSynchronizer();
  IParallelMng* pm = var_syncer->parallelMng();
 
  if (!pm->isParallel())
    return;
 
  mmvs->checkRecompute();
 
  //ItemGroup all_items = family->allItems();
  ITraceMng* tm = pm->traceMng();
  tm->info(4) << "MAT_SYNCHRONIZE_V5 name=" << this->name();
  //tm->info() << "STACK="<< platform::getStackTrace();
  const Integer data_type_size = (Integer)sizeof(DataType);
  {
    Int32ConstArrayView ranks = var_syncer->communicatingRanks();
    Integer nb_rank = ranks.size();
    std::vector<UniqueArray<DataType>> shared_values(nb_rank);
    std::vector<UniqueArray<DataType>> ghost_values(nb_rank);
 
    UniqueArray<Parallel::Request> requests;
 
    // Post the receives.
    Int32UniqueArray recv_ranks(nb_rank);
    for (Integer i = 0; i < nb_rank; ++i) {
      Int32 rank = ranks[i];
      ConstArrayView<MatVarIndex> ghost_matcells(mmvs->ghostItems(i));
      Integer total = ghost_matcells.size();
      ghost_values[i].resize(total);
      Integer total_byte = CheckedConvert::multiply(total, data_type_size);
      ByteArrayView bytes(total_byte, (Byte*)(ghost_values[i].unguardedBasePointer()));
      requests.add(pm->recv(bytes, rank, false));
    }
 
    // Post the sends
    for (Integer i = 0; i < nb_rank; ++i) {
      Int32 rank = ranks[i];
      ConstArrayView<MatVarIndex> shared_matcells(mmvs->sharedItems(i));
      Integer total = shared_matcells.size();
      shared_values[i].resize(total);
      ArrayView<DataType> values(shared_values[i]);
      for (Integer z = 0; z < total; ++z) {
        values[z] = this->operator[](shared_matcells[z]);
      }
      Integer total_byte = CheckedConvert::multiply(total, data_type_size);
      ByteArrayView bytes(total_byte, (Byte*)(values.unguardedBasePointer()));
      requests.add(pm->send(bytes, rank, false));
    }
 
    pm->waitAllRequests(requests);
 
    // Copy the received data into the ghost cells.
    for (Integer i = 0; i < nb_rank; ++i) {
      ConstArrayView<MatVarIndex> ghost_matcells(mmvs->ghostItems(i));
      Integer total = ghost_matcells.size();
      ConstArrayView<DataType> values(ghost_values[i].constView());
      for (Integer z = 0; z < total; ++z) {
        setValue(ghost_matcells[z], values[z]);
      }
    }
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
dumpValues(std::ostream& ostr, AllEnvCellVectorView view)
{
  ostr << "Dumping values for material variable name=" << this->name() << '\n';
  ENUMERATE_ALLENVCELL (iallenvcell, view) {
    AllEnvCell all_env_cell = *iallenvcell;
    ostr << "Cell uid=" << ItemPrinter(all_env_cell.globalCell()) << " v=" << value(all_env_cell._varIndex()) << '\n';
    for (CellComponentCellEnumerator ienvcell(all_env_cell); ienvcell.hasNext(); ++ienvcell) {
      MatVarIndex evi = ienvcell._varIndex();
      ostr << "env_value=" << value(evi) << ", mvi=" << evi << '\n';
      for (CellComponentCellEnumerator imatcell(*ienvcell); imatcell.hasNext(); ++imatcell) {
        MatVarIndex mvi = imatcell._varIndex();
        ostr << "mat_value=" << value(mvi) << ", mvi=" << mvi << '\n';
      }
    }
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
dumpValues(std::ostream& ostr)
{
  ostr << "Dumping values for material variable name=" << this->name() << '\n';
  IItemFamily* family = m_global_variable->itemFamily();
  if (!family)
    return;
  IMeshMaterialMng* material_mng = m_p->materialMng();
  dumpValues(ostr, material_mng->view(family->allItems()));
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename DataType> void
ItemMaterialVariableScalar<DataType>::
serialize(ISerializer* sbuf, Int32ConstArrayView ids)
{
  IItemFamily* family = m_global_variable->itemFamily();
  if (!family)
    return;
  ITraceMng* tm = family->traceMng();
  IMeshMaterialMng* mat_mng = m_p->materialMng();
  const Integer nb_count = DataTypeTraitsT<DataType>::nbBasicType();
  typedef typename DataTypeTraitsT<DataType>::BasicType BasicType;
  const eDataType data_type = DataTypeTraitsT<BasicType>::type();
  ItemVectorView ids_view(family->view(ids));
  bool has_mat = this->space() != MatVarSpace::Environment;
  switch (sbuf->mode()) {
  case ISerializer::ModeReserve: {
    Integer nb_val = 0;
    ENUMERATE_ALLENVCELL (iallenvcell, mat_mng, ids_view) {
      ENUMERATE_CELL_ENVCELL (ienvcell, (*iallenvcell)) {
        EnvCell envcell = *ienvcell;
        ++nb_val; // 1 value for the environment
        if (has_mat)
          nb_val += envcell.nbMaterial(); // 1 value per environment material.
      }
    }
    tm->info() << "RESERVE: nb_value=" << 1 << " size=" << (nb_val * nb_count);
    sbuf->reserve(DT_Int64, 1); // For the number of values.
    sbuf->reserveSpan(data_type, nb_val * nb_count); // For the number of values.
  } break;
  case ISerializer::ModePut: {
    UniqueArray<DataType> values;
    ENUMERATE_ALLENVCELL (iallenvcell, mat_mng, ids_view) {
      ENUMERATE_CELL_ENVCELL (ienvcell, (*iallenvcell)) {
        values.add(value(ienvcell._varIndex()));
        if (has_mat) {
          ENUMERATE_CELL_MATCELL (imatcell, (*ienvcell)) {
            values.add(value(imatcell._varIndex()));
          }
        }
      }
    }
    Integer nb_value = values.size();
    ConstArrayView<BasicType> basic_values(nb_value * nb_count, reinterpret_cast<BasicType*>(values.data()));
    tm->info() << "PUT: nb_value=" << nb_value << " size=" << basic_values.size();
    sbuf->putInt64(nb_value);
    sbuf->putSpan(basic_values);
  } break;
  case ISerializer::ModeGet: {
    UniqueArray<BasicType> basic_values;
    Int64 nb_value = sbuf->getInt64();
    basic_values.resize(nb_value * nb_count);
    sbuf->getSpan(basic_values);
    Span<const DataType> data_values(reinterpret_cast<DataType*>(basic_values.data()), nb_value);
    Integer index = 0;
    ENUMERATE_ALLENVCELL (iallenvcell, mat_mng, ids_view) {
      ENUMERATE_CELL_ENVCELL (ienvcell, (*iallenvcell)) {
        EnvCell envcell = *ienvcell;
        setValue(ienvcell._varIndex(), data_values[index]);
        ++index;
        if (has_mat) {
          ENUMERATE_CELL_MATCELL (imatcell, envcell) {
            setValue(imatcell._varIndex(), data_values[index]);
            ++index;
          }
        }
      }
    }
  } break;
  default:
    throw NotSupportedException(A_FUNCINFO, "Invalid serialize");
  }
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename ItemType, typename DataType>
MeshMaterialVariableScalar<ItemType, DataType>::
MeshMaterialVariableScalar(const MaterialVariableBuildInfo& v, PrivatePartType* global_var,
                           VariableRefType* global_var_ref, MatVarSpace mvs)
: ItemMaterialVariableScalar<DataType>(v, global_var, global_var_ref, mvs)
, m_true_global_variable_ref(global_var_ref) // Will be destroyed by the base class
{
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
template <typename ItemType, typename DataType>
MeshMaterialVariableScalar<ItemType, DataType>::
~MeshMaterialVariableScalar()
{
}
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#define ARCANE_INSTANTIATE_MAT(type) \
  template class ItemMaterialVariableBase<MaterialVariableScalarTraits<type>>; \
  template class ItemMaterialVariableScalar<type>; \
  template class MeshMaterialVariableScalar<Cell, type>; \
  template class MeshMaterialVariableCommonStaticImpl<MeshMaterialVariableScalar<Cell, type>>
 
ARCANE_INSTANTIATE_MAT(Byte);
ARCANE_INSTANTIATE_MAT(Int8);
ARCANE_INSTANTIATE_MAT(Int16);
ARCANE_INSTANTIATE_MAT(Int32);
ARCANE_INSTANTIATE_MAT(Int64);
ARCANE_INSTANTIATE_MAT(BFloat16);
ARCANE_INSTANTIATE_MAT(Float16);
ARCANE_INSTANTIATE_MAT(Float32);
ARCANE_INSTANTIATE_MAT(Real);
ARCANE_INSTANTIATE_MAT(Real2);
ARCANE_INSTANTIATE_MAT(Real3);
ARCANE_INSTANTIATE_MAT(Real2x2);
ARCANE_INSTANTIATE_MAT(Real3x3);
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
} // namespace Arcane::Materials
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/