src/ParaMEDLoader/ParaMEDFileMesh.cxx

   1 // Copyright (C) 2007-2014  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, or (at your option) any later version.
   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 // Author : Anthony Geay (EDF R&D)
  20
  21 #include "ParaMEDFileMesh.hxx"
  22 #include "MEDCouplingAutoRefCountObjectPtr.hxx"
  23 #include "MEDFileMesh.hxx"
  24 #include "MEDFileMeshLL.hxx"
  25 #include "MEDLoader.hxx"
  26
  27 using namespace ParaMEDMEM;
  28
  29 MEDFileMesh *ParaMEDFileMesh::New(int iPart, int nbOfParts, const std::string& fileName, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs)
  30 {
  31   MEDFileUtilities::CheckFileForRead(fileName);
  32   ParaMEDMEM::MEDCouplingMeshType meshType;
  33   MEDFileUtilities::AutoFid fid(MEDfileOpen(fileName.c_str(),MED_ACC_RDONLY));
  34   int dummy0,dummy1;
  35   std::string dummy2;
  36   MEDFileMeshL2::GetMeshIdFromName(fid,mName,meshType,dummy0,dummy1,dummy2);
  37   switch(meshType)
  38   {
  39     case UNSTRUCTURED:
  40       {
  41         return ParaMEDFileUMesh::New(iPart,nbOfParts,fileName,mName,dt,it,mrs);
  42       }
  43     default:
  44       throw INTERP_KERNEL::Exception("ParaMEDFileMesh::New : only unstructured mesh supported for the moment !");
  45   }
  46 }
  47
  48 MEDFileMesh *ParaMEDFileMesh::ParaNew(int iPart, int nbOfParts, const MPI_Comm com, const MPI_Info nfo, const std::string& fileName, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs)
  49 {
  50   MEDFileUtilities::CheckFileForRead(fileName);
  51   ParaMEDMEM::MEDCouplingMeshType meshType;
  52   MEDFileUtilities::AutoFid fid(MEDfileOpen(fileName.c_str(),MED_ACC_RDONLY));
  53   int dummy0,dummy1;
  54   std::string dummy2;
  55   MEDFileMeshL2::GetMeshIdFromName(fid,mName,meshType,dummy0,dummy1,dummy2);
  56   switch(meshType)
  57   {
  58     case UNSTRUCTURED:
  59       {
  60         return ParaMEDFileUMesh::ParaNew(iPart,nbOfParts,com,nfo,fileName,mName,dt,it,mrs);
  61       }
  62     default:
  63       throw INTERP_KERNEL::Exception("ParaMEDFileMesh::ParaNew : only unstructured mesh supported for the moment !");
  64   }
  65 }
  66
  67 MEDFileUMesh *ParaMEDFileUMesh::New(int iPart, int nbOfParts, const std::string& fileName, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs)
  68 {
  69   MEDCouplingAutoRefCountObjectPtr<MEDFileUMesh> ret;
  70   MEDFileUtilities::AutoFid fid(MEDfileOpen(fileName.c_str(),MED_ACC_RDONLY));
  71   //
  72   int meshDim, spaceDim, numberOfNodes;
  73   std::vector< std::vector< std::pair<INTERP_KERNEL::NormalizedCellType,int> > > typesDistrib(MEDLoader::GetUMeshGlobalInfo(fileName,mName,meshDim,spaceDim,numberOfNodes));
  74   std::vector<INTERP_KERNEL::NormalizedCellType> types;
  75   std::vector<int> distrib;
  76   for(std::vector< std::vector< std::pair<INTERP_KERNEL::NormalizedCellType,int> > >::const_iterator it0=typesDistrib.begin();it0!=typesDistrib.end();it0++)
  77     for(std::vector< std::pair<INTERP_KERNEL::NormalizedCellType,int> >::const_iterator it1=(*it0).begin();it1!=(*it0).end();it1++)
  78       {
  79         types.push_back((*it1).first);
  80         int tmp[3];
  81         DataArray::GetSlice(0,(*it1).second,1,iPart,nbOfParts,tmp[0],tmp[1]);
  82         tmp[2]=1;
  83         distrib.insert(distrib.end(),tmp,tmp+3);
  84       }
  85   ret=MEDFileUMesh::LoadPartOf(fileName,mName,types,distrib,dt,it,mrs);
  86   return ret.retn();
  87 }
  88
  89 MEDFileUMesh *ParaMEDFileUMesh::ParaNew(int iPart, int nbOfParts, const MPI_Comm com, const MPI_Info nfo, const std::string& fileName, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs)
  90 {
  91   MEDCouplingAutoRefCountObjectPtr<MEDFileUMesh> ret(MEDFileUMesh::New());
  92   MEDFileUtilities::AutoFid fid(MEDparFileOpen(fileName.c_str(),MED_ACC_RDONLY,com,nfo));
  93   //
  94   int meshDim, spaceDim, numberOfNodes;
  95   std::vector< std::vector< std::pair<INTERP_KERNEL::NormalizedCellType,int> > > typesDistrib(MEDLoader::GetUMeshGlobalInfo(fileName,mName,meshDim,spaceDim,numberOfNodes));
  96   return ret.retn();
  97 }
  98
  99 MEDFileMeshes *ParaMEDFileMeshes::New(int iPart, int nbOfParts, const std::string& fileName)
 100 {
 101   std::vector<std::string> ms(MEDLoader::GetMeshNames(fileName));
 102   MEDCouplingAutoRefCountObjectPtr<MEDFileMeshes> ret(MEDFileMeshes::New());
 103   for(std::vector<std::string>::const_iterator it=ms.begin();it!=ms.end();it++)
 104     {
 105       MEDCouplingAutoRefCountObjectPtr<MEDFileMesh> mesh(ParaMEDFileMesh::New(iPart,nbOfParts,fileName,(*it)));
 106       ret->pushMesh(mesh);
 107     }
 108   return ret.retn();
 109 }
 110
 111 MEDFileMeshes *ParaMEDFileMeshes::ParaNew(int iPart, int nbOfParts, const MPI_Comm com, const MPI_Info nfo, const std::string& fileName)
 112 {
 113   std::vector<std::string> ms(MEDLoader::GetMeshNames(fileName));
 114   MEDCouplingAutoRefCountObjectPtr<MEDFileMeshes> ret(MEDFileMeshes::New());
 115   for(std::vector<std::string>::const_iterator it=ms.begin();it!=ms.end();it++)
 116     {
 117       MEDCouplingAutoRefCountObjectPtr<MEDFileMesh> mesh(ParaMEDFileMesh::ParaNew(iPart,nbOfParts,com,nfo,fileName,(*it)));
 118       ret->pushMesh(mesh);
 119     }
 120   return ret.retn();
 121 }