src/MEDWrapper/V2_1/MEDbodyFittedEcr.cxx

   1 // Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   3 //
   4 // This library is free software; you can redistribute it and/or
   5 // modify it under the terms of the GNU Lesser General Public
   6 // License as published by the Free Software Foundation; either
   7 // version 2.1 of the License.
   8 //
   9 // This library is distributed in the hope that it will be useful
  10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 // Lesser General Public License for more details.
  13 //
  14 // You should have received a copy of the GNU Lesser General Public
  15 // License along with this library; if not, write to the Free Software
  16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 //
  18 // See http://www.salome-platform.org/
  19 //
  20 /*************************************************************************
  21 * COPYRIGHT (C) 1999 - 2002  EDF R&D
  22 * THIS LIBRARY IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
  23 * IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE
  24 * AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION;
  25 * EITHER VERSION 2.1 OF THE LICENSE, OR (AT YOUR OPTION) ANY LATER VERSION.
  26 *
  27 * THIS LIBRARY IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
  28 * WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
  29 * MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
  30 * LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS.
  31 *
  32 * YOU SHOULD HAVE RECEIVED A COPY OF THE GNU LESSER GENERAL PUBLIC LICENSE
  33 * ALONG WITH THIS LIBRARY; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION,
  34 * INC., 59 TEMPLE PLACE, SUITE 330, BOSTON, MA 02111-1307 USA
  35 *
  36 *************************************************************************/
  37
  38 #include "med_outils.hxx"
  39 #include "med.hxx"
  40
  41 #if defined(IRIX64)
  42 #define MED_INTEGER MED_INT64
  43 #else
  44 #define MED_INTEGER MED_INT32
  45 #endif
  46
  47 namespace med_2_1{
  48
  49 med_err
  50 MEDbodyFittedEcr(med_idt fid, char *maa, med_int mdim, med_float *coo, med_int *nbr, med_mode_switch mode_coo,
  51                  med_repere repere, char *nomcoo, char *unicoo, med_int *fam, med_int nnoeuds, med_mode_acces mode )
  52 {
  53     /* ecriture des coordonnees */
  54     med_idt   maaid, noeid, dataset;
  55     med_size dimd[1];
  56     char     chemin[MED_TAILLE_MAA+MED_TAILLE_NOM+1];
  57     med_int  type_rep_int;
  58     int      d;
  59     char     *ds;
  60
  61     /* On inhibe le gestionnaire d'erreur HDF */
  62     _MEDmodeErreurVerrouiller();
  63
  64     /* Si le maillage n'existe pas => erreur */
  65     strcpy(chemin, MED_MAA);
  66     strcat(chemin, maa);
  67     if ((maaid = _MEDdatagroupOuvrir(fid, chemin)) < 0) {
  68         return(-1);
  69     };
  70
  71     /* Si le Data Group "NOE" n'existe pas on le cree */
  72     if ((noeid = _MEDdatagroupOuvrir(maaid, MED_NOM_NOE)) < 0) {
  73         if ((noeid = _MEDdatagroupCreer(maaid, MED_NOM_NOE)) < 0) {
  74             return(-1);
  75         };
  76     };
  77
  78     /* Creation du Data Set "BOF" */
  79     dimd[0] = nnoeuds*mdim;
  80     if (_MEDdatasetNumEcrire(noeid, MED_NOM_BOF, MED_REEL64, mode_coo, mdim, MED_ALL, MED_NOPF, 0, 0, dimd, (unsigned char*)coo, mode) < 0) {
  81         return(-1);
  82     };
  83
  84     /* On re-ouvre le Data Set "BOF" pour y placer des attributs */
  85     if ((dataset = _MEDdatasetOuvrir(noeid, MED_NOM_BOF)) < 0) {
  86         return(-1);
  87     };
  88
  89     /* Attribut NBR (nombre de noeuds) */
  90     if (_MEDattrEntierEcrire(dataset, MED_NOM_NBR, &nnoeuds, mode) < 0) {
  91         return(-1);
  92     };
  93
  94     /* L'attribut "REP" */
  95     type_rep_int = (med_int)repere;
  96     if (_MEDattrEntierEcrire(dataset, MED_NOM_REP, &type_rep_int, mode) < 0) {
  97         return(-1);
  98     };
  99
 100     /* Attribut "NOM" */
 101     if (_MEDattrStringEcrire(dataset, MED_NOM_NOM, mdim*MED_TAILLE_PNOM, nomcoo, mode) < 0) {
 102         return(-1);
 103     };
 104
 105     /* Attribut "UNI" */
 106     if (_MEDattrStringEcrire(dataset, MED_NOM_UNI, mdim*MED_TAILLE_PNOM, unicoo, mode) < 0) {
 107         return(-1);
 108     };
 109
 110     if (_MEDdatasetFermer(dataset) < 0) return(-1);
 111
 112     dimd[0] = 1;
 113     for (d=0; d<mdim; d++) {
 114         switch (d) {
 115             case 0 : {
 116                 ds = MED_NOM_IN1;
 117                 break;
 118             };
 119             case 1 : {
 120                 ds = MED_NOM_IN2;
 121                 break;
 122             };
 123             case 2 : {
 124                 ds = MED_NOM_IN3;
 125                 break;
 126             };
 127             default : {
 128                 return(-1);
 129             };
 130         };
 131
 132         /* Creation du Data Set "IN1", "IN2", "IN3" contenant la taille du bodyfitted sur cette dimension */
 133         if (_MEDdatasetNumEcrire(noeid, ds, MED_INTEGER, mode_coo, MED_DIM1, MED_ALL, MED_NOPF, 0, 0, dimd, (unsigned char*)&nbr[d], mode) < 0) {
 134             return(-1);
 135         };
 136
 137         /* On re-ouvre ce Data Set pour y placer des attributs */
 138         if ((dataset = _MEDdatasetOuvrir(noeid, ds)) < 0) {
 139             return(-1);
 140         };
 141
 142         /* Attribut NBR (nombre de noeuds dans l'une des dimensions) */
 143         if (_MEDattrEntierEcrire(dataset, MED_NOM_NBR, &nbr[d], mode) < 0) {
 144             return(-1);
 145         };
 146
 147         if (_MEDdatasetFermer(dataset) < 0) return(-1);
 148     };
 149
 150     /* Ecriture des numeros de familles */
 151     if (MEDfamEcr(fid, maa, fam, nnoeuds, mode, MED_NOEUD, MED_POINT1) < 0) {
 152         return(-1);
 153     };
 154
 155     /* On ferme tout */
 156     if (_MEDdatagroupFermer(noeid) < 0) return(-1);
 157     if (_MEDdatagroupFermer(maaid) < 0) return(-1);
 158     return(0);
 159 }
 160
 161 }