src/INTERPOLATION/UseCaseInterpolationwots.cxx

   1 //  Copyright (C) 2007-2008  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 //  Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 //  This library is free software; you can redistribute it and/or
   7 //  modify it under the terms of the GNU Lesser General Public
   8 //  License as published by the Free Software Foundation; either
   9 //  version 2.1 of the License.
  10 //
  11 //  This library is distributed in the hope that it will be useful,
  12 //  but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 //  Lesser General Public License for more details.
  15 //
  16 //  You should have received a copy of the GNU Lesser General Public
  17 //  License along with this library; if not, write to the Free Software
  18 //  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 //  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22 #include "MEDMEM_Exception.hxx"
  23 #include "MEDMEM_define.hxx"
  24
  25 #include "MEDMEM_Field.hxx"
  26 #include "MEDMEM_Mesh.hxx"
  27 #include "MEDMEM_Interpolation.hxx"
  28
  29 #include "stdio.h"
  30
  31 using namespace MEDMEM;
  32 using namespace MED_EN;
  33
  34 int main () {
  35   const char * fromFileName  = "fromMesh.med";
  36   const char * toFileName    = "toMesh.med";
  37   const char * fieldName     = "fieldnodedouble";
  38
  39   const char * fromMeshName  = "fromMesh";
  40   const char * toMeshName    = "toMesh";
  41
  42   const int flagConvex = 1;
  43   const int interpolationType = 1;
  44
  45   try {
  46
  47     MESH            fromMesh    (MED_DRIVER,fromFileName,fromMeshName);
  48     SUPPORT         fromSupport (&fromMesh,"XsupportX",MED_NODE);
  49     FIELD<double>   fromField   (&fromSupport,MED_DRIVER,fromFileName,fieldName);
  50     MESH            toMesh      (MED_DRIVER,toFileName,toMeshName);
  51
  52     INTERPOLATION<3> myInter (fromField,toMesh);
  53
  54     FIELD<double>   * toField = myInter.interpolate(interpolationType,flagConvex);
  55
  56     toField->addDriver(MED_DRIVER,toFileName,toField->getName()) ;
  57
  58     toField->write();
  59
  60
  61   } catch (MEDEXCEPTION& ex){
  62     MESSAGE_MED(ex.what()) ;
  63   }
  64 }