src/MEDSPLITTER/MEDSPLITTER_MESHCollectionDriver.H

   1 // Copyright (C) 2007-2012  CEA/DEN, EDF R&D
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19 #ifndef MEDSPLITTER_MESHCOLLECTIONDRIVER_H
  20 #define MEDSPLITTER_MESHCOLLECTIONDRIVER_H
  21
  22
  23 #include <libxml/tree.h>
  24 #include <libxml/parser.h>
  25 #include <libxml/xpath.h>
  26 #include <libxml/xpathInternals.h>
  27
  28 /*!reads a distributed field
  29  *
  30  * \param fields vector of fields (one field per subdomain)
  31  * \param fieldname name of the field
  32  * \param itnumber number of iteration
  33  * \param ordernumber internal number inside the iteration
  34  * */
  35 template <class T>
  36 void MESHCollectionDriver::readFields(vector<MEDMEM::FIELD<T>* >& fields,char* fieldname, int itnumber, int ordernumber)
  37 {
  38   for (int i=0; i<_collection->getMesh().size(); i++)
  39   {
  40     char filename[256];
  41     strcpy(filename,_filename[i].c_str());
  42     cout << "maillage : " << filename << " champ : " << fieldname << endl;
  43     //    MEDMEM::FIELD<T>* field = new MEDMEM::FIELD<T>(MEDMEM::MED_DRIVER,filename,fieldname,itnumber,ordernumber);
  44     fields.push_back (new MEDMEM::FIELD<T>(MEDMEM::MED_DRIVER,filename,fieldname,itnumber,ordernumber));
  45   }
  46 }
  47
  48
  49 /*!writes a distributed field
  50  *
  51  * \param fields vector of fields (one field per subdomain)
  52  * \param fieldname name of the field
  53  * */
  54 template <class T>
  55 void MESHCollectionDriver::writeFields(vector<MEDMEM::FIELD<T>* >& fields,char* fieldname)
  56 {
  57   xmlDocPtr master_doc=xmlParseFile(_master_filename.c_str());
  58
  59   if (!master_doc)
  60     throw MEDEXCEPTION("MEDSPLITTER writeFields - Master File does not exist");
  61
  62   //number of domains
  63
  64   xmlXPathContextPtr xpathCtx = xmlXPathNewContext(master_doc);
  65   xmlXPathObjectPtr xpathObj = xmlXPathEvalExpression(BAD_CAST "//mapping/mesh", xpathCtx);
  66   //assuming there is only one mesh in the XML file
  67   xmlNodePtr mesh_node= xpathObj->nodesetval->nodeTab[0];
  68   xmlNodePtr field_node = xmlNewChild(mesh_node, 0, BAD_CAST "field",0);
  69   xmlNewProp(field_node,BAD_CAST "name",BAD_CAST fieldname);
  70
  71
  72   for (int i=0; i<_collection->getMesh().size(); i++)
  73   {
  74     char filename[256];
  75     strcpy(filename,_filename[i].c_str());
  76     int driverid = fields[i]->addDriver(MEDMEM::MED_DRIVER, filename, fieldname);
  77     fields[i]->write(driverid);
  78     xmlNodePtr chunk_node= xmlNewChild(field_node,0,BAD_CAST "chunk",0);
  79     char id[8];
  80     sprintf(id,"%d",i+1);
  81     xmlNewProp(chunk_node,BAD_CAST "subdomain",BAD_CAST id);
  82     xmlNewProp(chunk_node,BAD_CAST "name", BAD_CAST fieldname);
  83
  84   }
  85   xmlSaveFormatFileEnc(_master_filename.c_str(), master_doc, "UTF-8", 1);
  86   xmlFreeDoc(master_doc);
  87   //xmlCleanupParser();
  88
  89 }
  90
  91 #endif