LireMED.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 LireMED_included
#define LireMED_included
 
/*! @brief Classe LireMED Lire un fichier MED
 *
 *     Structure du jeu de donnee (en dimension 2) :
 *     LireMED dom medfile
 */
#include <Interprete_geometrique_base.h>
#include <TRUSTTabs_forward.h>
#include <Domaine_forward.h>
#include <med++.h>
#include <medcoupling++.h>
 
#ifdef MEDCOUPLING_
#include <MEDLoader.hxx>
#include <MEDFileMesh.hxx>
#endif
 
class Nom;
 
template <typename _SIZE_>
class LireMED_32_64 : public Interprete_geometrique_base_32_64<_SIZE_>
{
  Declare_instanciable_32_64(LireMED_32_64);
public :
  using int_t = _SIZE_;
  using ArrOfInt_t = ArrOfInt_T<_SIZE_>;
  using IntTab_t = IntTab_T<_SIZE_>;
  using BigArrOfInt_ = TRUSTArray<int, _SIZE_>; // Always int stored inside (family numbers typically)
  using DoubleTab_t = DoubleTab_T<_SIZE_>;
 
 
  using Domaine_t = Domaine_32_64<_SIZE_>;
  using Elem_geom_t = OWN_PTR(Elem_geom_base_32_64<_SIZE_>);
 
  LireMED_32_64(const Nom& file_name, const Nom& mesh_name);
  Entree& interpreter_(Entree&) override;
  void lire_geom(bool subDom=true);
  void retrieve_MC_objects();
 
#ifdef MEDCOUPLING_
  inline const MEDCoupling::MEDCouplingUMesh* get_mc_mesh() const
  {
    // really using the cast operator of MEDCoupling, not a wild cast:
    return (const MEDCoupling::MEDCouplingUMesh*)mcumesh_;
  }
#endif
 
protected:
  Nom nom_fichier_;              ///< Name of the MED file to read
  Nom nom_mesh_ = "--any--";     ///< Name of the mesh in the MED file to read
  bool axi1d_ = false;           ///< Are we in Axi1D
  bool convertAllToPoly_ = false;///< Should the mesh be converted to polygons/polyedrons
  bool isVEFForce_ = false;      ///< Should we force element type to be VEF compatible
  int space_dim_ = -1;           ///< Space dimension read in the MED file
  Elem_geom_t type_elem_;          ///< Highest dimension element type
  Nom type_face_;                ///< Boundary element type
  Noms noms_bords_;              ///< Names of the boundaries
  Noms exclude_grps_;            ///< Names of the (face) groups to skip when reading the file
  Noms restrict_ssz_;           ///< Names of the subzones to keep only in the .geo file
  Noms internal_face_grps_;      ///< Names of the internals face groups to read in the file
#ifdef MEDCOUPLING_
  MEDCoupling::MEDCouplingAxisType axis_type_  = MEDCoupling::MEDCouplingAxisType::AX_CART ;
  MEDCoupling::MCAuto<MEDCoupling::MEDCouplingUMesh> mcumesh_;
  MEDCoupling::MCAuto<MEDCoupling::MEDFileUMesh> mfumesh_;
#endif
 
  Nom type_medcoupling_to_type_geo_trio(int type_cell, bool cell_from_boundary) const;
  void prepare_som_and_elem(DoubleTab_t& sommets, IntTab_t& les_elems);
  void finalize_sommets(const DoubleTab_t& sommets2, DoubleTab_t& sommets) const;
  void write_sub_dom_datasets() const;
  void read_boundaries(BigArrOfInt_& familles, IntTab_t& all_faces_bords);
  void fill_frontieres(const BigArrOfInt_& familles, const IntTab_t& all_faces_bords);
};
 
using LireMED = LireMED_32_64<int>;
using LireMED_64 = LireMED_32_64<trustIdType>;
 
#endif