src/INTERPOLATION/test_MEDMEM_Interpolation.cxx

   1 #include "MEDMEM_Exception.hxx"
   2 #include "MEDMEM_define.hxx"
   3
   4 #include "MEDMEM_Field.hxx"
   5 #include "MEDMEM_Mesh.hxx"
   6 #include "MEDMEM_Interpolation.hxx"
   7
   8 #include "stdio.h"
   9
  10 #ifdef _DEBUG_
  11 #include "LocalTraceCollector.hxx"
  12 #endif /* ifdef _DEBUG_*/
  13
  14 main () {
  15 #ifdef _DEBUG_
  16   LocalTraceCollector::instance();
  17 #endif /* ifdef _DEBUG_*/
  18
  19   const char * fromFileName  = "fromMesh.med";
  20   const char * toFileName    = "toMesh.med";
  21   //const char * fieldName     = "fieldcelldoublevector";
  22   const char * fieldName     = "fieldnodedouble";
  23
  24   const char * fromMeshName  = "fromMesh";
  25   const char * toMeshName    = "toMesh";
  26
  27   try {
  28
  29     cout<<"Lecture du Maillage Source : "<<flush; MESH            fromMesh    (MED_DRIVER,fromFileName,fromMeshName);           cout<<"OK !"<<endl;
  30     //cout<<"Construction du support    : "<<flush; SUPPORT         fromSupport (&fromMesh,"XsupportX",MED_CELL);                 cout<<"OK !"<<endl;
  31     cout<<"Construction du support    : "<<flush; SUPPORT         fromSupport (&fromMesh,"XsupportX",MED_NODE);                 cout<<"OK !"<<endl;
  32     cout<<"Lecture du champ           : "<<flush; FIELD<double>   fromField   (&fromSupport,MED_DRIVER,fromFileName,fieldName); cout<<"OK !"<<endl;
  33     cout<<"Lecture du Mailllage Cible : "<<flush; MESH            toMesh      (MED_DRIVER,toFileName,toMeshName);               cout<<"OK !"<<endl;
  34
  35     INTERPOLATION<3> myInter (fromField,toMesh);
  36
  37     //FIELD<double>   * toField = myInter.interpolate(0,1);
  38     FIELD<double>   * toField = myInter.interpolate(1,1);
  39
  40     cout<<"Creation du driver"<<endl;
  41
  42     toField->addDriver(MED_DRIVER,toFileName,toField->getName()) ;
  43
  44     cout<<"toField->getName()                   = "<<toField->getName()                  <<endl;
  45     cout<<"toField->getDescription()            = "<<toField->getDescription()           <<endl;
  46     cout<<"toField->getNumberOfComponents()     = "<<toField->getNumberOfComponents()    <<endl;
  47     cout<<"toField->getNumberOfValues()         = "<<toField->getNumberOfValues()        <<endl;
  48     cout<<"toField->getComponentsNames()        = "<<toField->getComponentsNames()       <<endl;
  49     cout<<"toField->getComponentsDescriptions() = "<<toField->getComponentsDescriptions()<<endl;
  50     cout<<"toField->getMEDComponentsUnits()     = "<<toField->getMEDComponentsUnits()    <<endl;
  51     cout<<"toField->getIterationNumber()        = "<<toField->getIterationNumber()       <<endl;
  52     cout<<"toField->getTime()                   = "<<toField->getTime()                  <<endl;
  53     cout<<"toField->getOrderNumber()            = "<<toField->getOrderNumber()           <<endl;
  54     cout<<"toField->getValueType()              = "<<toField->getValueType()             <<endl;
  55
  56     toField->write();
  57
  58     cout<<"Fin"<<endl;
  59
  60   } catch (MEDEXCEPTION& ex){
  61     MESSAGE(ex.what()) ;
  62   }
  63 }