Format_Post_CGNS.cpp

/****************************************************************************
* 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.
*
*****************************************************************************/
 
#include <Format_Post_CGNS.h>
#include <My_Comm_Group.h>
#include <Option_CGNS.h>
#include <PE_Groups.h>
#include <EChaine.h>
#include <Param.h>
 
Implemente_instanciable_sans_constructeur(Format_Post_CGNS, "Format_Post_CGNS", Format_Post_base);
 
static void init_my_com_group()
{
  if (!PE_Groups::has_user_defined_group())
    {
      const int nb_nodes = PE_Groups::get_node_group().get_number_of_nodes();
      if (nb_nodes > 1)
        {
          Option_CGNS::LINKED_FILES_PER_COMM_GROUP = true;
          Option_CGNS::USE_LINKS = true;
          Option_CGNS::PARALLEL_OVER_ZONE = false;
          Nom ech = "{ Group_nb ";
          ech += Nom(nb_nodes);
          ech += " }";
 
          My_Comm_Group comm;
          EChaine ech2(ech);
 
          ech2 >> comm;
        }
    }
}
 
Format_Post_CGNS::Format_Post_CGNS()
{
  reset();
}
 
void Format_Post_CGNS::reset()
{
  cgns_basename_ = "??";
}
 
Sortie& Format_Post_CGNS::printOn(Sortie& os) const
{
  Process::exit("Format_Post_CGNS::printOn - Should not be called ! ");
  return os;
}
 
Entree& Format_Post_CGNS::readOn(Entree& is)
{
  Process::exit("Format_Post_CGNS::readOn - Should not be called ! ");
  verify_if_cgns(__func__);
  init_my_com_group();
  return Format_Post_base::readOn(is);
}
 
void Format_Post_CGNS::set_param(Param& param) const
{
  Process::exit("Format_Post_CGNS::readOn - Should not be called ! ");
  verify_if_cgns(__func__);
  param.ajouter("nom_fichier", &cgns_basename_, Param::REQUIRED);
}
 
void Format_Post_CGNS::resetTime(double t, const std::string dirname)
{
  if (dirname.empty())
    Process::exit("\nError in Format_Post_CGNS::resetTime. The name of the directory is empty. This is forbidden !!! \n");
 
#ifdef HAS_CGNS
  cgns_writer_.cgns_resetTime(t, dirname, cgns_basename_);
#endif
}
 
int Format_Post_CGNS::initialize_by_default(const Nom& file_basename)
{
  verify_if_cgns(__func__);
  init_my_com_group();
  cgns_basename_ = file_basename;
#ifdef HAS_CGNS
  cgns_writer_.cgns_set_base_name(cgns_basename_);
#endif
  return 1;
}
 
int Format_Post_CGNS::initialize(const Nom& file_basename, const int format, const Nom& option_para)
{
  verify_if_cgns(__func__);
  init_my_com_group();
  cgns_basename_ = file_basename;
#ifdef HAS_CGNS
  cgns_writer_.cgns_set_base_name(cgns_basename_);
#endif
  return 1;
}
 
int Format_Post_CGNS::ecrire_entete(const double temps_courant,const int reprise,const int est_le_premier_post)
{
  verify_if_cgns(__func__);
 
  if (est_le_premier_post)
    {
#ifdef HAS_CGNS
      cgns_writer_.cgns_open_file();
#endif
    }
  return 1;
}
 
void Format_Post_CGNS::set_discr_type(const Nom& type)
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_set_discr_type(type);
#endif
}
 
void Format_Post_CGNS::set_postraiter_domain()
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_set_postraiter_domain();
#endif
}
 
void Format_Post_CGNS::set_deformable_domain()
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_set_is_deformable_domain();
#endif
}
 
void Format_Post_CGNS::set_lagrangian_domain()
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_set_lagrangian_domain();
#endif
}
 
void Format_Post_CGNS::set_loc_vector(const std::vector<std::string>& vec)
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_set_loc_vector(vec);
#endif
}
 
int Format_Post_CGNS::ecrire_temps(const double t)
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_add_time(t);
#endif
  return 1;
}
 
int Format_Post_CGNS::finir_ecriture(double temps_courant)
{
#ifdef HAS_CGNS
  cgns_writer_.finir_ecriture(temps_courant);
#endif
  return 1;
}
 
int Format_Post_CGNS::finir(const int est_le_dernier_post)
{
  if (est_le_dernier_post)
    {
#ifdef HAS_CGNS
      cgns_writer_.cgns_finir();
#endif
    }
  return 1;
}
 
void Format_Post_CGNS::ecrire_domaine_dual(const Domaine& domaine, const int est_le_premier_post)
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_write_domaine_dual(domaine.le_nom(), est_le_premier_post);
#endif
}
 
int Format_Post_CGNS::ecrire_domaine_low_level(const Domaine * dom,const Nom& nom_dom, const DoubleTab& som, const IntTab& elem, const Motcle& type_e)
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_write_domaine(dom, nom_dom, som, elem, type_e);
#endif
  return 1;
}
 
int Format_Post_CGNS::ecrire_domaine_dis(const Domaine& domaine,const OBS_PTR(Domaine_dis_base)& domaine_dis_base,const int est_le_premier_post)
{
  if (domaine_dis_base)
    {
      domaine_dis_ = domaine_dis_base;
#ifdef HAS_CGNS
      cgns_writer_.cgns_associer_domaine_dis(domaine_dis_base.valeur());
#endif
    }
  return ecrire_domaine(domaine, est_le_premier_post);
}
 
int Format_Post_CGNS::ecrire_domaine(const Domaine& domaine, const int est_le_premier_post)
{
  Motcle type_e = domaine.type_elem()->que_suis_je();
  ecrire_domaine_low_level(&domaine, domaine.le_nom(), domaine.les_sommets(), domaine.les_elems(), type_e);
 
  // Si on a des frontieres domaine, on les ecrit egalement
  const LIST(OBS_PTR(Domaine)) &bords = domaine.domaines_frontieres();
  for (auto &itr : bords) ecrire_domaine(itr.valeur(), est_le_premier_post);
  return 1;
}
 
int Format_Post_CGNS::ecrire_champ(const Domaine& domaine, const Noms& unite_, const Noms& noms_compo, int ncomp, double temps,
                                   const Nom& id_du_champ, const Nom& id_du_domaine, const Nom& localisation,
                                   const Nom& nature, const DoubleTab& valeurs)
{
#ifdef HAS_CGNS
  cgns_writer_.cgns_write_field(domaine, noms_compo, temps, id_du_champ, id_du_domaine, localisation, valeurs );
#endif
  return 1;
}