1 // Copyright (C) 2007-2012 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 // Author : Anthony Geay (CEA/DEN)
21 #include "MEDFileData.hxx"
23 using namespace ParaMEDMEM;
25 MEDFileData *MEDFileData::New(const char *fileName) throw(INTERP_KERNEL::Exception)
27 return new MEDFileData(fileName);
30 MEDFileData *MEDFileData::New()
32 return new MEDFileData;
35 MEDFileData *MEDFileData::deepCpy() const throw(INTERP_KERNEL::Exception)
37 MEDCouplingAutoRefCountObjectPtr<MEDFileFields> fields;
38 if((const MEDFileFields *)_fields)
39 fields=_fields->deepCpy();
40 MEDCouplingAutoRefCountObjectPtr<MEDFileMeshes> meshes;
41 if((const MEDFileMeshes *)_meshes)
42 meshes=_meshes->deepCpy();
43 MEDCouplingAutoRefCountObjectPtr<MEDFileData> ret=MEDFileData::New();
44 ret->_fields=fields; ret->_meshes=meshes;
48 std::size_t MEDFileData::getHeapMemorySize() const
51 if((const MEDFileFields *)_fields)
52 ret+=_fields->getHeapMemorySize();
53 if((const MEDFileMeshes *)_meshes)
54 ret+=_meshes->getHeapMemorySize();
58 MEDFileFields *MEDFileData::getFields() const
60 return const_cast<MEDFileFields *>(static_cast<const MEDFileFields *>(_fields));
63 MEDFileMeshes *MEDFileData::getMeshes() const
65 return const_cast<MEDFileMeshes *>(static_cast<const MEDFileMeshes *>(_meshes));
68 void MEDFileData::setFields(MEDFileFields *fields) throw(INTERP_KERNEL::Exception)
71 throw INTERP_KERNEL::Exception("MEDFileData::setFields : input pointer is null !");
76 void MEDFileData::setMeshes(MEDFileMeshes *meshes) throw(INTERP_KERNEL::Exception)
79 throw INTERP_KERNEL::Exception("MEDFileData::setMeshes : input pointer is null !");
84 int MEDFileData::getNumberOfFields() const throw(INTERP_KERNEL::Exception)
86 const MEDFileFields *f=_fields;
88 throw INTERP_KERNEL::Exception("MEDFileData::getNumberOfFields : no fields set !");
89 return f->getNumberOfFields();
92 int MEDFileData::getNumberOfMeshes() const throw(INTERP_KERNEL::Exception)
94 const MEDFileMeshes *m=_meshes;
96 throw INTERP_KERNEL::Exception("MEDFileData::getNumberOfMeshes : no meshes set !");
97 return m->getNumberOfMeshes();
100 std::string MEDFileData::simpleRepr() const
102 std::ostringstream oss;
103 oss << "(***************)\n(* MEDFileData *)\n(***************)\n\nFields part :\n*************\n\n";
104 const MEDFileFields *tmp=_fields;
107 tmp->simpleRepr(0,oss);
111 oss << "No fields set !!!\n\n";
112 oss << "Meshes part :\n*************\n\n";
113 const MEDFileMeshes *tmp2=_meshes;
116 tmp2->simpleReprWithoutHeader(oss);
119 oss << "No meshes set !!!\n";
123 bool MEDFileData::changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab) throw(INTERP_KERNEL::Exception)
126 MEDFileFields *fields=_fields;
128 ret=fields->changeMeshNames(modifTab) || ret;
129 MEDFileMeshes *meshes=_meshes;
131 ret=meshes->changeNames(modifTab) || ret;
135 bool MEDFileData::changeMeshName(const char *oldMeshName, const char *newMeshName) throw(INTERP_KERNEL::Exception)
137 std::string oldName(oldMeshName);
138 std::vector< std::pair<std::string,std::string> > v(1);
139 v[0].first=oldName; v[0].second=newMeshName;
140 return changeMeshNames(v);
144 * This method performs unpolyzation in meshes in \a this and if it leads to a modification to one or more than one meshes in \a this
145 * the fields are automatically renumbered (for those impacted, that is to say here fields on cells and fields on gauss points on impacted fields)
147 * \return If true is returned it means that some meshes in \a this has been modified and eventually fields have been renumbered.
148 * \n If false \a this remains unchanged.
150 bool MEDFileData::unPolyzeMeshes() throw(INTERP_KERNEL::Exception)
152 MEDFileMeshes *ms=_meshes;
155 std::vector< MEDFileMesh * > meshesImpacted;
156 std::vector< DataArrayInt * > renumParamsOfMeshImpacted;//same size as meshesImpacted
157 std::vector< std::vector<int> > oldCodeOfMeshImpacted,newCodeOfMeshImpacted;//same size as meshesImpacted
158 std::vector<MEDCouplingAutoRefCountObjectPtr<DataArrayInt> > memSaverIfThrow;//same size as meshesImpacted
159 for(int i=0;i<ms->getNumberOfMeshes();i++)
161 MEDFileMesh *m=ms->getMeshAtPos(i);
164 std::vector<int> oldCode,newCode;
165 DataArrayInt *o2nRenumCell=0;
166 bool modif=m->unPolyze(oldCode,newCode,o2nRenumCell);
169 renumParamsOfMeshImpacted.push_back(o2nRenumCell); memSaverIfThrow.push_back(o2nRenumCell);
170 oldCodeOfMeshImpacted.push_back(oldCode);
171 newCodeOfMeshImpacted.push_back(newCode);
172 meshesImpacted.push_back(m);
175 if(!meshesImpacted.empty())
177 MEDFileFields *fs=_fields;
179 for(std::size_t i=0;i<meshesImpacted.size();i++)
180 fs->renumberEntitiesLyingOnMesh(meshesImpacted[i]->getName(),oldCodeOfMeshImpacted[i],newCodeOfMeshImpacted[i],renumParamsOfMeshImpacted[i]);
182 return !meshesImpacted.empty();
185 MEDFileData::MEDFileData()
189 MEDFileData::MEDFileData(const char *fileName) throw(INTERP_KERNEL::Exception)
192 _fields=MEDFileFields::New(fileName);
193 _meshes=MEDFileMeshes::New(fileName);
195 catch(INTERP_KERNEL::Exception& e)
200 void MEDFileData::write(const char *fileName, int mode) const throw(INTERP_KERNEL::Exception)
202 med_access_mode medmod=MEDFileUtilities::TraduceWriteMode(mode);
203 MEDFileUtilities::AutoFid fid=MEDfileOpen(fileName,medmod);
204 const MEDFileMeshes *ms=_meshes;
207 const MEDFileFields *fs=_fields;