1 // Copyright (C) 2007-2013 CEA/DEN, EDF R&D
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.
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.
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
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
19 #ifndef MEDSPLITTER_MESHCOLLECTIONDRIVER_H
20 #define MEDSPLITTER_MESHCOLLECTIONDRIVER_H
23 #include <libxml/tree.h>
24 #include <libxml/parser.h>
25 #include <libxml/xpath.h>
26 #include <libxml/xpathInternals.h>
28 /*!reads a distributed field
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
36 void MESHCollectionDriver::readFields(vector<MEDMEM::FIELD<T>* >& fields,char* fieldname, int itnumber, int ordernumber)
38 for (int i=0; i<_collection->getMesh().size(); i++)
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));
49 /*!writes a distributed field
51 * \param fields vector of fields (one field per subdomain)
52 * \param fieldname name of the field
55 void MESHCollectionDriver::writeFields(vector<MEDMEM::FIELD<T>* >& fields,char* fieldname)
57 xmlDocPtr master_doc=xmlParseFile(_master_filename.c_str());
60 throw MEDEXCEPTION("MEDSPLITTER writeFields - Master File does not exist");
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);
72 for (int i=0; i<_collection->getMesh().size(); i++)
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);
81 xmlNewProp(chunk_node,BAD_CAST "subdomain",BAD_CAST id);
82 xmlNewProp(chunk_node,BAD_CAST "name", BAD_CAST fieldname);
85 xmlSaveFormatFileEnc(_master_filename.c_str(), master_doc, "UTF-8", 1);
86 xmlFreeDoc(master_doc);