doc/MEDMEM/FIELDgeneral.cxx

   1 // Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   3 //
   4 using namespace std;
   5 #include "MEDMEM_Mesh.hxx"
   6 #include "MEDMEM_Field.hxx"
   7
   8 using namespace MEDMEM;
   9 using namespace MED_EN ;
  10
  11 int main (int argc, char ** argv) {
  12
  13   const string MedFile = "pointe.med" ;
  14   const string MeshName = "maa1" ;
  15   const string FieldName = "fieldcelldoublevector" ;
  16
  17   /* read MESH */
  18   MESH * myMesh = new MESH(MED_DRIVER,MedFile,MeshName) ;
  19   //  myMesh->read() ;
  20
  21   /* read FIELD */
  22   // we need a support :
  23   SUPPORT * mySupport = new SUPPORT(myMesh,"Support on all Cells",MED_CELL);
  24   FIELD<double> myField(mySupport,MED_DRIVER,MedFile,FieldName) ;
  25   //  myField.read() ;
  26
  27   /* what in Field ? */
  28   // How many components
  29   int NumberOfCompoennts = myField.getNumberOfComponents() ;
  30
  31   const string * ComponentsNames = myField.getComponentsNames();
  32   const string * ComponentsDescriptions = myField.getComponentsDescriptions();
  33   const string * ComponentsUnits = myField.getMEDComponentsUnits();
  34
  35   for(int i=0;i<NumberOfCompoennts; i++) {
  36     cout << "Component " << i << " :" <<endl ;
  37     cout << "  - name        : " << ComponentsNames[i] << endl ;
  38     cout << "  - description : " << ComponentsDescriptions[i] << endl ;
  39     cout << "  - unit        : " << ComponentsUnits[i] << endl ;
  40   }
  41
  42   // Which iteration :
  43   int IterationNumber = myField.getIterationNumber() ; // negative mean undefined
  44   int OrderNumber = myField.getOrderNumber() ;
  45   // internal iteration at this time iteration, negative mean undefined
  46   double Time = myField.getTime() ;
  47
  48   cout << "Iteration " << IterationNumber << " at time " << Time <<
  49     " (and order number " << OrderNumber << ")" << endl ;
  50
  51   // How many Value :
  52   int NumberOfValue = mySupport->getNumberOfElements(MED_ALL_ELEMENTS);
  53   // Value
  54   const double * Value = myField.getValue();
  55   for(int i=0; i<NumberOfValue; i++) {
  56     for(int j=0; j<NumberOfCompoennts; j++)
  57       cout << Value[i*NumberOfCompoennts+j] << " " ;
  58     cout << endl ;
  59   }
  60
  61   delete mySupport;
  62   delete myMesh;
  63
  64   return 0 ;
  65 }