med.h

Aller à la documentation de ce fichier.

/* -*- mode:C; coding:utf-8 -*- */
/*  This file is part of MED.
 *
 *  COPYRIGHT (C) 1999 - 2025  EDF R&D, CEA/DEN
 *  MED is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU Lesser General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  MED is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with MED.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#ifndef MED_H
#define MED_H
 
#include <hdf5.h>
#ifdef HAVE_STDINT_H
#include <stdint.h>   // Use the C99 official header
#endif
 
#ifdef __cplusplus
extern "C" {
#endif
 
/* This is the minimum hdf version med must have  */
#define HDF_VERSION_REF         "5-1.14.0"
#define HDF_VERSION_MAJOR_REF   1
#define HDF_VERSION_MINOR_REF   14
#define HDF_VERSION_RELEASE_REF 0
#define HDF_VERSION_NUM_REF     (HDF_VERSION_MAJOR_REF * 10000 + HDF_VERSION_MINOR_REF * 100  + HDF_VERSION_RELEASE_REF )
 
#define HDF5_VERSION (H5_VERS_MAJOR * 10000 + H5_VERS_MINOR * 100  + H5_VERS_RELEASE )
 
#if HDF5_VERSION < 10607
/*TODO : A enlever après MAJ complète api hdf5-1.8 */
#define HDF5_SELECT_BUG (const hsize_t **)
#else
#define HDF5_SELECT_BUG
#endif
 
/*Le symbole H5F_LIBVER_18 n'existe pas dans hdf5-1.8.4 (uniquement à partir de la version 1.8.6)
  Pour rester dans le modèle interne hdf introduit en 1.8 quelque soit les futurs modèles internes hdf
  on impose le modèle 1.8 et non le LATEST (cf _MEDfileCreate).
  La valeur utilisée est celle définie dans les  version >= 1.8.6 (cf H5Fpublic.h) */
#ifndef H5F_LIBVER_18
 
#define H5F_LIBVER_18 1
#endif
 
/*En anglais*/
#define MED_MAJOR_NUM  6
#define MED_MINOR_NUM  0
#define MED_RELEASE_NUM 1
  
#define MED_3_LATEST_MINOR    3
#define MED_4_LATEST_MINOR    2
#define MED_40_LATEST_RELEASE 1
#define MED_41_LATEST_RELEASE 2
#define MED_42_LATEST_RELEASE @MED_42_LATEST_RELEASE@
#define MED_5_LATEST_MINOR    1
#define MED_50_LATEST_RELEASE 2
#define MED_51_LATEST_RELEASE 0
#define MED_6_LATEST_MINOR    0
#define MED_60_LATEST_RELEASE 1
 
/*En franglais*/
#define MED_NUM_MAJEUR  6
#define MED_NUM_MINEUR  0
#define MED_NUM_RELEASE 1
 
#define STR(x) #x
#define _MED_VERSION(x,y,z) STR(x) "." STR(y) "." STR(z)
#define MED_VERSION_STR _MED_VERSION(MED_NUM_MAJEUR,MED_NUM_MINEUR,MED_NUM_RELEASE)
 
 
 
#define MED_NULL       (void *) NULL
#define MED_MAX_PARA   20
 
#define MED_COMMENT_SIZE 200
#define MED_IDENT_SIZE     8
#define MED_NAME_SIZE     64
#define MED_SNAME_SIZE    16
#define MED_LNAME_SIZE    80
 
#define MED_SNAME_BLANK "                "
#define MED_NAME_BLANK  "                                                                "
/*                       0123456789012345678901234567890123456789012345678901234567890123*/
#define MED_LNAME_BLANK  MED_NAME_BLANK MED_SNAME_BLANK
 
 
#define MED_PATHNAME_SIZE 255
#define MED_MAX_CHFID_PATH MED_PATHNAME_SIZE  
/* 2.3v3.0 ok*/
typedef enum {MED_FULL_INTERLACE, 
              MED_NO_INTERLACE,    
              MED_UNDEF_INTERLACE=-1
        }  med_switch_mode;
#define med_interlacing_mode med_switch_mode
 
/* 2.3v3.0 ok*/
typedef enum { MED_UNDEF_STMODE,   
               MED_GLOBAL_STMODE,  
               MED_COMPACT_STMODE, 
               MED_GLOBAL_PFLMODE=MED_GLOBAL_STMODE,   
               MED_COMPACT_PFLMODE=MED_COMPACT_STMODE, 
               MED_UNDEF_PFLMODE=MED_UNDEF_STMODE      
             }  med_storage_mode;
 
 
typedef enum {MED_ACC_RDONLY, 
              MED_ACC_RDWR,   
              MED_ACC_RDEXT,  
              MED_ACC_CREAT,  
              MED_ACC_UNDEF   
             } med_access_mode;
 
/*2.3v3.0 ok*/
typedef enum {MED_UNSTRUCTURED_MESH, MED_STRUCTURED_MESH, MED_UNDEF_MESH_TYPE=-1} med_mesh_type;
 
/* 2.3v3.0 ok*/
typedef enum {MED_CARTESIAN_GRID,
              MED_POLAR_GRID,
              MED_CURVILINEAR_GRID,
              MED_UNDEF_GRID_TYPE=-1} med_grid_type;
 
/* 2.3v3.0 ok*/
typedef enum {MED_CELL, MED_DESCENDING_FACE, MED_DESCENDING_EDGE, MED_NODE,
              MED_NODE_ELEMENT, MED_STRUCT_ELEMENT, MED_ALL_ENTITY_TYPE,
              MED_UNDEF_ENTITY_TYPE=-1 } med_entity_type;
#define MED_N_ENTITY_TYPES 6
 
/* 2.3v3.0 ok*/
typedef enum {MED_COORDINATE, MED_CONNECTIVITY, MED_NAME, MED_NUMBER,  MED_FAMILY_NUMBER, 
        MED_COORDINATE_AXIS1, MED_COORDINATE_AXIS2, MED_COORDINATE_AXIS3, 
        MED_INDEX_FACE, MED_INDEX_NODE, MED_GLOBAL_NUMBER,
        MED_VARIABLE_ATTRIBUTE, MED_COORDINATE_TRSF,MED_UNDEF_DATATYPE=-1} med_data_type;
 
/* 2.3v3.0 ok*/
typedef enum {MED_INTERNAL_FLOAT32=4, MED_INTERNAL_FLOAT64=6,
              MED_INTERNAL_INT8=20,
              MED_INTERNAL_INT32=24 , MED_INTERNAL_INT64=26, MED_INTERNAL_INT=28,
              MED_INTERNAL_NAME=30  , MED_INTERNAL_SNAME=32,
              MED_INTERNAL_LNAME=34 , MED_INTERNAL_IDENT=38,
              MED_INTERNAL_CHAR=40  , MED_INTERNAL_UNDEF =0 } med_internal_type; 
/*TODO: Réécrire les routines _MEDatt avec med_internal_type au lieu de med_field_type */
 
typedef enum {MED_DOUBLE=MED_INTERNAL_FLOAT64,  
              MED_FLOAT64=MED_INTERNAL_FLOAT64, 
              MED_FLOAT32=MED_INTERNAL_FLOAT32, 
              MED_INT32=MED_INTERNAL_INT32,     
              MED_INT64=MED_INTERNAL_INT64,     
              MED_INT=MED_INTERNAL_INT          
             }  med_field_type; 
 
typedef enum {MED_ATT_FLOAT64=MED_INTERNAL_FLOAT64,
              MED_ATT_INT=MED_INTERNAL_INT,
              MED_ATT_NAME=MED_INTERNAL_NAME,
              MED_ATT_UNDEF=MED_INTERNAL_UNDEF }  med_attribute_type; 
 
typedef med_field_type med_parameter_type; 
 
#define MED_ATT_FLOAT64_ctype   med_float       
#define MED_ATT_INT_ctype       med_int         
#define MED_ATT_NAME_ctype      char            
#define MED_ATT_UNDEF_ctype     unknown_ctype   
#define MED_CTYPE(medatttype) medatttype##_ctype 
 
/* 2.3v3.0 ok : med_type_donnee*/
typedef enum {MED_MESH, MED_FIELD, MED_LIBRARY, MED_FILE,
        MED_MESH_SUPPORT,MED_ELSTRUCT, MED_FAMILY,MED_EQUIVALENCE, MED_GROUP, MED_JOINT, 
        MED_LOCALIZATION, MED_PROFILE, MED_FILTER, MED_INTERPOLATION, MED_NUMERICAL_DATA,
        MED_LINK, MED_FAMILY_ZERO, MED_FAMILY_NODE, MED_FAMILY_ELEM,
        MED_CLASS_UNDEF=-1, MED_CLASS_ALL=-2 } med_class;
 
 
typedef int med_geometry_type;
 
/* Mailles dont la géométrie à une connectivité fixe*/
 
#define MED_POINT1 001
 
#define MED_SEG2 102
#define MED_SEG3 103
#define MED_SEG4 104
#define MED_TRIA3 203
#define MED_QUAD4 204
#define MED_TRIA6 206
#define MED_TRIA7 207
#define MED_QUAD8 208
#define MED_QUAD9 209
 
#define MED_TETRA4  304
#define MED_PYRA5   305
#define MED_PENTA6  306
#define MED_HEXA8   308
#define MED_TETRA10 310
#define MED_OCTA12  312
#define MED_PYRA13  313
#define MED_PENTA15 315
#define MED_PENTA18 318
#define MED_HEXA20  320
#define MED_HEXA27  327
 
/* Mailles dont la géométrie à une connectivité variable*/
#define MED_POLYGON    400
#define MED_POLYGON2   420
#define MED_POLYHEDRON 500
 
/* Mailles dont la géométrie est dynamique (découverte du modèle dans le fichier) */
#define MED_STRUCT_GEO_INTERNAL     600
#define MED_STRUCT_GEO_SUP_INTERNAL 700
 
#define MED_NONE 0
#define MED_NO_GEOTYPE MED_NONE
#define MED_UNDEF_GEOTYPE MED_NONE
#define MED_UNDEF_GEOMETRY_TYPE MED_NONE
#define MED_ALL_GEOTYPE -1
#define MED_GEO_ALL MED_ALL_GEOTYPE
 
#define MED_N_CELL_GEO            25 
#define MED_N_CELL_FIXED_GEO      24 
#define MED_N_CELL_GEO_FIXED_CON  21 
#define MED_N_FACE_GEO               8
#define MED_N_FACE_FIXED_GEO         8
#define MED_N_FACE_GEO_FIXED_CON     6
 
#define MED_N_EDGE_TYPES             3
#define MED_N_EDGE_FIXED_GEO         3
#define MED_N_EDGE_GEO_FIXED_CON     3
 
#define MED_N_NODE_GEO               1
#define MED_N_NODE_FIXED_GEO         1
#define MED_N_NODE_GEO_FIXED_CON     1
 
/*2.3v3.0 ok*/
typedef enum {MED_NODAL, MED_DESCENDING, MED_UNDEF_CONNECTIVITY_MODE=-1, MED_NO_CMODE=-1} med_connectivity_mode; 
 
/*2.3v3.0 ok*/
typedef enum {MED_CARTESIAN, MED_CYLINDRICAL, MED_SPHERICAL, MED_UNDEF_AXIS_TYPE=-1} med_axis_type;
 
typedef enum {MED_FALSE, MED_TRUE} med_bool ;
 
#define MED_GAUSS_ELNO "MED_GAUSS_ELNO" 
#define MED_IPOINT_ELNO MED_GAUSS_ELNO
 
 
#define MED_NO_NAME ""                  
#define MED_NO_MESHNAME ""              
#define MED_NO_MESH     ""              
#define MED_NO_MESH_SUPPORT ""          
#define MED_NO_LOCALIZATION ""          
#define MED_NO_INTERPOLATION ""         
#define MED_NO_IPOINT_INTERNAL 1       
#define MED_NO_PROFILE ""                     
#define MED_NO_GROUP ""                 
#define MED_ALLENTITIES_PROFILE MED_NO_PROFILE 
#define MED_NO_PROFILE_INTERNAL   "MED_NO_PROFILE_INTERNAL"   
#define MED_SAME_PROFILE_INTERNAL "MED_SAME_PROFILE_INTERNAL" 
#define MED_ALL_CONSTITUENT 0                 
#define MED_UNDEF_CONSTITUENT -1              
#define MED_UNDEF_SIZE 0
#define MED_NO_PROFILE_SIZE 0  
typedef enum {MED_SORT_DTIT,MED_SORT_ITDT,MED_SORT_UNDEF=-1 } med_sorting_type; 
#define MED_ALL_DT   -2 
#define MED_ALL_IT   -2 
#define MED_NO_DT    -1 
#define MED_NO_IT    -1 
#define MED_UNDEF_DT  0 
#define MED_ATT_NOT_FILLED -1
 
 
/* correspondance des types avec HDF 5 */
 
typedef hsize_t        med_size;  
typedef hssize_t       med_ssize; 
typedef hid_t          med_idt;   
typedef herr_t         med_err;   
/* types elementaires */
 
typedef double         med_float; 
typedef double         med_float64; 
typedef double         med_double;
typedef float          med_float32;
 
/*La ligne suivante est déterminée à l'étape de configuration*/
typedef int  med_int; 
typedef int32_t  med_int32; 
typedef int64_t  med_int64; 
#define MED_MAX_NAME_IN_PATH 5
#define MED_MAX_DATAGROUPNAME_IN_PATH 2*MED_MAX_NAME_IN_PATH*MED_NAME_SIZE 
typedef struct {
  /*PUBLIC IN attributes*/
  med_class           medclass;
  med_int             nname;                      
  const char *        name[MED_MAX_NAME_IN_PATH]; 
  /*PRIVATE attributes*/
  char                _datagroupname[MED_MAX_DATAGROUPNAME_IN_PATH+1]; 
} med_path;
 
#define INIT2X( x ) x, x,
#define INIT4X( x ) x, x, x, x
#define INIT10X( x ) x, x, x, x, x, x, x, x, x, x
#define MED_MAX_FILTER_SPACES 4100
#define MED_MAX_FILTER_SPACES_INIT {INIT4X(INIT10X(INIT10X(INIT10X(0)))), INIT10X(INIT10X(0)) }
 
typedef enum {MED_FILTER_LIST       = 1,
              MED_FILTER_BLOCK      = 2,
              MED_FILTER_TYPE_UNDEF =-1 } med_filter_type; 
 
typedef struct {
  /*PRIVATE attributes*/
  med_int          nspaces;
  med_idt          memspace[MED_MAX_FILTER_SPACES];
  med_idt          diskspace[MED_MAX_FILTER_SPACES];
  med_filter_type  filtertype;
  /*PUBLIC IN attributes*/
  med_int          nentity;
  med_int          nvaluesperentity;
  med_int          nconstituentpervalue;
  med_int          constituentselect;
  med_switch_mode  switchmode;
  med_switch_mode  diskswitchmode;
  /* Sélection par liste d'entités */
  med_int          filterarraysize;
  const med_int *  filterarray;
  med_size*        filterarray23v30;
  med_int          profilearraysize;/*PRIVATE attribute en IN et OUT positionnée dans MEDfilterCr*/
  /* profil utilisé et mode de stockage en mémoire */
  med_storage_mode storagemode;
  char             profilename[MED_NAME_SIZE+1];
  /* Sélection par block d'entités */
  med_size         start        ;
  med_size         stride       ;
  med_size         count        ;
  med_size         blocksize    ;
  med_size         lastblocksize;  
} med_filter; 
 
 
#define MED_FILTER_INIT { 0, MED_MAX_FILTER_SPACES_INIT, MED_MAX_FILTER_SPACES_INIT,         \
                          MED_FILTER_TYPE_UNDEF, 0, 0, 0, 0,                                 \
                          MED_UNDEF_INTERLACE,MED_UNDEF_INTERLACE,MED_NO_FILTER_SIZE,        \
      NULL, NULL, MED_NO_PROFILE_SIZE,MED_UNDEF_PFLMODE,MED_NO_PROFILE,\
                          0, 0, 0, 0, 0 }
 
#define MED_NO_FILTER NULL
#define MED_NO_FILTER_SIZE 0
#define MED_NO_PROFILE_F "                                "
 
typedef struct {
  med_int majeur; med_int mineur; med_int release;
} med_file_version;
#define MED_FILE_VERSION_INIT {0, 0, 0};
 
typedef struct {
  void   *app_image_ptr;  
  size_t  app_image_size; 
  int     ref_count; 
  void           *fapl_image_ptr;  
  size_t          fapl_image_size; 
  int             fapl_ref_count;  
  void           *vfd_image_ptr;   
  size_t          vfd_image_size;  
  int             vfd_ref_count;   
  med_access_mode flags; 
} med_memfile ;
 
#define MED_MEMFILE_INIT {NULL, 0, 1, NULL, 0, 0, NULL, 0, 0, MED_ACC_RDONLY}
 
#define AFF_MEMFILE \
XSCRUTE(memfile->app_image_ptr);\
ISCRUTE_long(memfile->app_image_size);\
ISCRUTE_int(memfile->ref_count);\
XSCRUTE(memfile->fapl_image_ptr);\
ISCRUTE_long(memfile->fapl_image_size);\
ISCRUTE_int(memfile->fapl_ref_count);\
XSCRUTE(memfile->vfd_image_ptr);\
ISCRUTE_long(memfile->vfd_image_size);\
ISCRUTE_int(memfile->vfd_ref_count);\
ISCRUTE_int(memfile->flags);
 
 
#include "med_proto.h"
 
/*En cas de modifications/ajouts de type d'entité/géométrique mettre à jour les itérateurs suivants : */
 
extern MEDC_EXPORT med_entity_type MED_GET_ENTITY_TYPE[MED_N_ENTITY_TYPES+2];
extern MEDC_EXPORT const char * const MED_GET_ENTITY_TYPENAME[MED_N_ENTITY_TYPES+2];
extern MEDC_EXPORT med_geometry_type MED_GET_CELL_GEOMETRY_TYPE[MED_N_CELL_FIXED_GEO+2];
extern MEDC_EXPORT const char * const MED_GET_CELL_GEOMETRY_TYPENAME[MED_N_CELL_FIXED_GEO+2];
extern MEDC_EXPORT med_geometry_type MED_GET_FACE_GEOMETRY_TYPE[MED_N_FACE_FIXED_GEO+2];
extern MEDC_EXPORT const char * const MED_GET_FACE_GEOMETRY_TYPENAME[MED_N_FACE_FIXED_GEO+2];
extern MEDC_EXPORT med_geometry_type MED_GET_EDGE_GEOMETRY_TYPE[MED_N_EDGE_FIXED_GEO+2];
extern MEDC_EXPORT const char * MED_GET_EDGE_GEOMETRY_TYPENAME[MED_N_EDGE_FIXED_GEO+2];
extern MEDC_EXPORT med_geometry_type MED_GET_NODE_GEOMETRY_TYPE[MED_N_NODE_FIXED_GEO+2];
extern MEDC_EXPORT const char * MED_GET_NODE_GEOMETRY_TYPENAME[MED_N_NODE_FIXED_GEO+2];
extern MEDC_EXPORT const med_geometry_type * const MED_GET_CELL_GEOMETRY_TYPE_FROM_ENTITY_TYPE[MED_N_ENTITY_TYPES+2];
#define MED_GET_GEOMETRY_TYPE_IT_FROM_ENTITY_TYPE_IT MED_GET_CELL_GEOMETRY_TYPE_FROM_ENTITY_TYPE
extern MEDC_EXPORT const int MED_GET_N_FIXED_GEO_FOR_ENTITY_TYPE_IT[MED_N_ENTITY_TYPES+2];
 
  
/* Itérateurs sur les types d'entité et types géométriques      */
/* TODO : Ajouter un module de documentation pour ces fonctions */
extern MEDC_EXPORT int                 MEDgetEntityTypeIt         (med_entity_type   entitytype);
extern MEDC_EXPORT int                 MEDgetCellGeometryTypeIt   (med_geometry_type geotype   );
extern MEDC_EXPORT int                 MEDgetFaceGeometryTypeIt   (med_geometry_type geotype   );
extern MEDC_EXPORT int                 MEDgetEdgeGeometryTypeIt   (med_geometry_type geotype   );
extern MEDC_EXPORT int                 MEDgetNodeGeometryTypeIt   (med_geometry_type geotype   );
extern MEDC_EXPORT int                 MEDgetEntityGeometryTypeIt (med_entity_type   entitytype, med_geometry_type geotype);
 
extern MEDC_EXPORT const char * const  MEDgetGeometryTypeName     (med_geometry_type geotype   );
extern MEDC_EXPORT const char * const  MEDgetEntityTypeName       (med_entity_type   entitytype);
extern MEDC_EXPORT const char * const  MEDgetInterlacingModeName(med_interlacing_mode interlacing);
 
/*      Predefined names for struct elements            */
#define MED_PARTICLE_NAME   "MED_PARTICLE"
#define MED_BALL_NAME       "MED_BALL"
#define MED_BEAM_NAME       "MED_BEAM"
 
/*      Predefined attributes names for struct elements */
#define MED_PARTICLE_LABEL  "MED_PARTICLE_LABEL"
#define MED_BALL_DIAMETER   "MED_BALL_DIAMETER"
#define MED_BEAM_THICKNESS  "MED_BEAM_THICKNESS"
 
#ifdef __cplusplus
}
#endif
 
#endif  /* MED_H */