Scatter.h

/****************************************************************************
* Copyright (c) 2026, CEA
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
* 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
* 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
* 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
* IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
* OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
*****************************************************************************/
 
#ifndef Scatter_included
#define Scatter_included
 
#include <TRUSTTabs_forward.h>
#include <Interprete.h>
#include <TRUST_Ref.h>
#include <Joint.h>
#include <Domaine_forward.h>
 
class Domaine_VF;
 
class Scatter : public Interprete
{
  Declare_instanciable(Scatter);
public:
  Entree& interpreter(Entree&) override;
  virtual void lire_domaine(Nom& fil);
  Domaine& domaine();
 
  static int Chercher_Correspondance(const DoubleTab& sommets1, const DoubleTab& sommets2, ArrOfInt& correspondance, const double epsilon);
  static void construire_correspondance_sommets_par_coordonnees(Domaine& dom, bool allow_resize=false);
  static void construire_correspondance_aretes_par_coordonnees(Domaine_VF& zvf);
  static void construire_correspondance_items_par_coordonnees(Joints& joints, const JOINT_ITEM type_item, const DoubleTab& coord_items, bool allow_resize=false);
 
  static void construire_structures_paralleles(Domaine& dom);
 
  //static void rechercher_elems_joints(Domaine & domaine);
 
  static void calculer_espace_distant(Domaine&                  domaine,
                                      const int           nb_items_reels,
                                      const ArrsOfInt& items_to_send,
                                      const JOINT_ITEM type_item);
 
  static void calculer_nb_items_virtuels(Joints& joints,
                                         const JOINT_ITEM type_item);
 
  static void calculer_renum_items_communs(Joints& joints, const JOINT_ITEM type_item);
  static void calculer_espace_distant_faces(Domaine& domaine, const int nb_faces_reelles, const IntTab& elem_faces);
  static void calculer_espace_distant_aretes(Domaine& domaine, const int nb_aretes_reelles, const IntTab& elem_aretes);
 
  static void calculer_espace_distant_elements(Domaine& dom);
  static void corriger_espace_distant_elements_perio(Domaine& dom);
 
  static void calculer_espace_distant_sommets(Domaine& dom);
  static void construire_espace_virtuel_traduction(const MD_Vector& md_indice, const MD_Vector& md_valeur, IntTab& tableau, const int error_is_fatal = 1);
 
  static void reordonner_faces_de_joint(Domaine& dom);
  static void ajouter_joints(Domaine& domaine, ArrOfInt& pe_voisins);
 
  static void trier_les_joints(Joints& joints);
  static void construire_md_vector(const Domaine&, int nb_items_reels, const JOINT_ITEM, MD_Vector&);
 
  static void check_consistancy_remote_items( Domaine& dom, const ArrOfInt& mergedDomaines );
 
  // Although Scatter is never meant to be used in 64b, those two methods are needed by both configurations:
  template <typename _SIZE_>
  static void init_sequential_domain(Domaine_32_64<_SIZE_>& dom);
  template <typename _SIZE_>
  static void uninit_sequential_domain(Domaine_32_64<_SIZE_>& dom);
 
 
protected:
  OBS_PTR(Domaine) le_domaine;
 
  void read_domain_no_comm(Entree& fic, bool& read_perio);
};
 
#endif