src/HOMARD_I/HOMARD_Cas_i.cxx

   1 // Copyright (C) 2011-2016  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
  20 // Remarques :
  21 // L'ordre de description des fonctions est le meme dans tous les fichiers
  22 // HOMARD_aaaa.idl, HOMARD_aaaa.hxx, HOMARD_aaaa.cxx, HOMARD_aaaa_i.hxx, HOMARD_aaaa_i.cxx :
  23 // 1. Les generalites : Name, Delete, DumpPython, Dump, Restore
  24 // 2. Les caracteristiques
  25 // 3. Le lien avec les autres structures
  26 //
  27 // Quand les 2 fonctions Setxxx et Getxxx sont presentes, Setxxx est decrit en premier
  28 //
  29
  30 #include "HOMARD_Cas_i.hxx"
  31 #include "HOMARD_Gen_i.hxx"
  32 #include "HOMARD_Cas.hxx"
  33 #include "HOMARD_DriverTools.hxx"
  34 #include "HOMARD.hxx"
  35
  36 #include "utilities.h"
  37 #include <vector>
  38 #include <sys/stat.h>
  39
  40 #ifdef WIN32
  41 #include <direct.h>
  42 #endif
  43
  44 //=============================================================================
  45 /*!
  46  *  standard constructor
  47  */
  48 //=============================================================================
  49 HOMARD_Cas_i::HOMARD_Cas_i()
  50 {
  51   MESSAGE( "Default constructor, not for use" );
  52   ASSERT( 0 );
  53 }
  54
  55 //=============================================================================
  56 /*!
  57  *  standard constructor
  58  */
  59 //=============================================================================
  60 HOMARD_Cas_i::HOMARD_Cas_i( CORBA::ORB_ptr orb,
  61                             HOMARD::HOMARD_Gen_var engine )
  62 {
  63   MESSAGE( "HOMARD_Cas_i" );
  64   _gen_i = engine;
  65   _orb = orb;
  66   myHomardCas = new ::HOMARD_Cas();
  67   ASSERT( myHomardCas );
  68 }
  69
  70 //=============================================================================
  71 /*!
  72  *  standard destructor
  73  */
  74 //=============================================================================
  75 HOMARD_Cas_i::~HOMARD_Cas_i()
  76 {
  77 }
  78 //=============================================================================
  79 //=============================================================================
  80 // Generalites
  81 //=============================================================================
  82 //=============================================================================
  83 void HOMARD_Cas_i::SetName( const char* Name )
  84 {
  85   ASSERT( myHomardCas );
  86   myHomardCas->SetName( Name );
  87 }
  88 //=============================================================================
  89 char* HOMARD_Cas_i::GetName()
  90 {
  91   ASSERT( myHomardCas );
  92   return CORBA::string_dup( myHomardCas->GetName().c_str() );
  93 }
  94 //=============================================================================
  95 CORBA::Long  HOMARD_Cas_i::Delete( CORBA::Long Option )
  96 {
  97   ASSERT( myHomardCas );
  98   char* CaseName = GetName() ;
  99   MESSAGE ( "Delete : destruction du cas " << CaseName << ", Option = " << Option );
 100   return _gen_i->DeleteCase(CaseName, Option) ;
 101 }
 102 //=============================================================================
 103 char* HOMARD_Cas_i::GetDumpPython()
 104 {
 105   ASSERT( myHomardCas );
 106   return CORBA::string_dup( myHomardCas->GetDumpPython().c_str() );
 107 }
 108 //=============================================================================
 109 std::string HOMARD_Cas_i::Dump() const
 110 {
 111   return HOMARD::Dump( *myHomardCas );
 112 }
 113 //=============================================================================
 114 bool HOMARD_Cas_i::Restore( const std::string& stream )
 115 {
 116   return HOMARD::Restore( *myHomardCas, stream );
 117 }
 118 //=============================================================================
 119 //=============================================================================
 120 // Caracteristiques
 121 //=============================================================================
 122 //=============================================================================
 123 void HOMARD_Cas_i::SetDirName( const char* NomDir )
 124 {
 125   ASSERT( myHomardCas );
 126   int codret ;
 127   // A. recuperation du nom ; on ne fait rien si c'est le meme
 128   char* oldrep = GetDirName() ;
 129   if ( strcmp(oldrep,NomDir) == 0 )
 130   {
 131    return ;
 132   }
 133   MESSAGE ( "SetDirName : passage de oldrep = "<< oldrep << " a NomDir = "<<NomDir);
 134   // B. controle de l'usage du repertoire
 135   char* CaseName = GetName() ;
 136   char* casenamedir = _gen_i->VerifieDir(NomDir) ;
 137   if ( ( std::string(casenamedir).size() > 0 ) & ( strcmp(CaseName,casenamedir)!=0 ) )
 138   {
 139     INFOS ( "Le repertoire " << NomDir << " est deja utilise pour le cas "<< casenamedir );
 140     SALOME::ExceptionStruct es;
 141     es.type = SALOME::BAD_PARAM;
 142     std::string text ;
 143     text = "The directory " + std::string(NomDir) + " is already used for the case " + std::string(casenamedir) ;
 144     es.text = CORBA::string_dup(text.c_str());
 145     throw SALOME::SALOME_Exception(es);
 146   }
 147   // C. Changement/creation du repertoire
 148   codret = myHomardCas->SetDirName( NomDir );
 149   if ( codret != 0 )
 150   {
 151     SALOME::ExceptionStruct es;
 152     es.type = SALOME::BAD_PARAM;
 153     std::string text ;
 154     if ( codret == 1 ) { text = "The directory for the case cannot be modified because some iterations are already defined." ; }
 155     else               { text = "The directory for the case cannot be reached." ; }
 156     es.text = CORBA::string_dup(text.c_str());
 157     throw SALOME::SALOME_Exception(es);
 158   }
 159   // D. En cas de reprise, deplacement du point de depart
 160   if ( GetState() != 0 )
 161   {
 162     MESSAGE ( "etat : " << GetState() ) ;
 163     // D.1. Nom local du repertoire de l'iteration de depart dans le repertoire actuel du cas
 164     HOMARD::HOMARD_Iteration_ptr Iter = GetIter0() ;
 165     char* DirNameIter = Iter->GetDirNameLoc() ;
 166     MESSAGE ( "SetDirName : nom actuel pour le repertoire de l iteration, DirNameIter = "<< DirNameIter);
 167     // D.2. Recherche d'un nom local pour l'iteration de depart dans le futur repertoire du cas
 168     char* nomDirIter = _gen_i->CreateDirNameIter(NomDir, 0 );
 169     MESSAGE ( "SetDirName : nom futur pour le repertoire de l iteration, nomDirIter = "<< nomDirIter);
 170     // D.3. Creation du futur repertoire local pour l'iteration de depart
 171     std::string nomDirIterTotal ;
 172     nomDirIterTotal = std::string(NomDir) + "/" + std::string(nomDirIter) ;
 173 #ifndef WIN32
 174     if (mkdir(nomDirIterTotal.c_str(), S_IRWXU|S_IRGRP|S_IXGRP) != 0)
 175 #else
 176     if (_mkdir(nomDirIterTotal.c_str()) != 0)
 177 #endif
 178     {
 179       MESSAGE ( "nomDirIterTotal : " << nomDirIterTotal ) ;
 180       SALOME::ExceptionStruct es;
 181       es.type = SALOME::BAD_PARAM;
 182       std::string text = "The directory for the starting iteration cannot be created." ;
 183       es.text = CORBA::string_dup(text.c_str());
 184       throw SALOME::SALOME_Exception(es);
 185     }
 186     // D.4. Deplacement du contenu du repertoire
 187     std::string oldnomDirIterTotal ;
 188     oldnomDirIterTotal = std::string(oldrep) + "/" + std::string(DirNameIter) ;
 189     std::string commande = "mv " + std::string(oldnomDirIterTotal) + "/*" + " " + std::string(nomDirIterTotal) ;
 190     codret = system(commande.c_str()) ;
 191     if ( codret != 0 )
 192     {
 193       SALOME::ExceptionStruct es;
 194       es.type = SALOME::BAD_PARAM;
 195       std::string text = "The starting point for the case cannot be moved into the new directory." ;
 196       es.text = CORBA::string_dup(text.c_str());
 197       throw SALOME::SALOME_Exception(es);
 198     }
 199     commande = "rm -rf " + std::string(oldnomDirIterTotal) ;
 200     codret = system(commande.c_str()) ;
 201     if ( codret != 0 )
 202     {
 203       SALOME::ExceptionStruct es;
 204       es.type = SALOME::BAD_PARAM;
 205       std::string text = "The starting point for the case cannot be deleted." ;
 206       es.text = CORBA::string_dup(text.c_str());
 207       throw SALOME::SALOME_Exception(es);
 208     }
 209     // D.5. Memorisation du nom du repertoire de l'iteration
 210     Iter->SetDirNameLoc(nomDirIter) ;
 211   }
 212   return ;
 213 }
 214 //=============================================================================
 215 char* HOMARD_Cas_i::GetDirName()
 216 {
 217   ASSERT( myHomardCas );
 218   return CORBA::string_dup( myHomardCas->GetDirName().c_str() );
 219 }
 220 //=============================================================================
 221 CORBA::Long HOMARD_Cas_i::GetState()
 222 {
 223   ASSERT( myHomardCas );
 224 // Nom de l'iteration initiale
 225   char* Iter0Name = GetIter0Name() ;
 226   HOMARD::HOMARD_Iteration_ptr Iter = _gen_i->GetIteration(Iter0Name) ;
 227   int state = Iter->GetNumber() ;
 228   return state ;
 229 }
 230 //=============================================================================
 231 CORBA::Long HOMARD_Cas_i::GetNumberofIter()
 232 {
 233   ASSERT( myHomardCas );
 234   return myHomardCas->GetNumberofIter();
 235 }
 236 //=============================================================================
 237 void HOMARD_Cas_i::SetConfType( CORBA::Long ConfType )
 238 {
 239   ASSERT( myHomardCas );
 240 //   VERIFICATION( (ConfType>=-2) && (ConfType<=3) );
 241   myHomardCas->SetConfType( ConfType );
 242 }
 243 //=============================================================================
 244 CORBA::Long HOMARD_Cas_i::GetConfType()
 245 {
 246   ASSERT( myHomardCas );
 247   return myHomardCas->GetConfType();
 248 }
 249 //=============================================================================
 250 void HOMARD_Cas_i::SetExtType( CORBA::Long ExtType )
 251 {
 252   ASSERT( myHomardCas );
 253 //   VERIFICATION( (ExtType>=0) && (ExtType<=2) );
 254   myHomardCas->SetExtType( ExtType );
 255 }
 256 //=============================================================================
 257 CORBA::Long HOMARD_Cas_i::GetExtType()
 258 {
 259   ASSERT( myHomardCas );
 260   return myHomardCas->GetExtType();
 261 }
 262 //=============================================================================
 263 void HOMARD_Cas_i::SetBoundingBox( const HOMARD::extrema& LesExtrema )
 264 {
 265   ASSERT( myHomardCas );
 266   std::vector<double> VExtrema;
 267   ASSERT( LesExtrema.length() == 10 );
 268   VExtrema.resize( LesExtrema.length() );
 269   for ( int i = 0; i < LesExtrema.length(); i++ )
 270   {
 271     VExtrema[i] = LesExtrema[i];
 272   }
 273   myHomardCas->SetBoundingBox( VExtrema );
 274 }
 275 //=============================================================================
 276 HOMARD::extrema* HOMARD_Cas_i::GetBoundingBox()
 277 {
 278   ASSERT(myHomardCas );
 279   HOMARD::extrema_var aResult = new HOMARD::extrema();
 280   std::vector<double> LesExtremes = myHomardCas->GetBoundingBox();
 281   ASSERT( LesExtremes.size() == 10 );
 282   aResult->length( 10 );
 283   for ( int i = 0; i < LesExtremes.size(); i++ )
 284   {
 285     aResult[i] = LesExtremes[i];
 286   }
 287   return aResult._retn();
 288 }
 289 //=============================================================================
 290 void HOMARD_Cas_i::AddGroup( const char* Group)
 291 {
 292   ASSERT( myHomardCas );
 293   myHomardCas->AddGroup( Group );
 294 }
 295 //=============================================================================
 296 void HOMARD_Cas_i::SetGroups( const HOMARD::ListGroupType& ListGroup )
 297 {
 298   ASSERT( myHomardCas );
 299   std::list<std::string> ListString ;
 300   for ( int i = 0; i < ListGroup.length(); i++ )
 301   {
 302     ListString.push_back(std::string(ListGroup[i]));
 303   }
 304   myHomardCas->SetGroups( ListString );
 305 }
 306 //=============================================================================
 307 HOMARD::ListGroupType* HOMARD_Cas_i::GetGroups()
 308 {
 309   ASSERT(myHomardCas );
 310   const std::list<std::string>& ListString = myHomardCas->GetGroups();
 311   HOMARD::ListGroupType_var aResult = new HOMARD::ListGroupType();
 312   aResult->length( ListString.size() );
 313   std::list<std::string>::const_iterator it;
 314   int i = 0;
 315   for ( it = ListString.begin(); it != ListString.end(); it++ )
 316   {
 317     aResult[i++] = CORBA::string_dup( (*it).c_str() );
 318   }
 319   return aResult._retn();
 320 }
 321 //=============================================================================
 322 void HOMARD_Cas_i::AddBoundaryGroup( const char* BoundaryName, const char* Group)
 323 {
 324   MESSAGE ("AddBoundaryGroup : BoundaryName = "<< BoundaryName << ", Group = " << Group );
 325   ASSERT( myHomardCas );
 326   // A. La liste des frontiere+groupes
 327   const std::list<std::string>& ListBoundaryGroup = myHomardCas->GetBoundaryGroup();
 328   std::list<std::string>::const_iterator it;
 329   // B. La frontiere
 330   // B.1. La frontiere est-elle deja enregistree pour ce cas ?
 331   bool existe = false ;
 332   for ( it = ListBoundaryGroup.begin(); it != ListBoundaryGroup.end(); it++ )
 333   {
 334 //     MESSAGE ("..  Frontiere : "<< *it );
 335     if ( *it == BoundaryName ) { existe = true ; }
 336     it++ ;
 337   }
 338   // B.2. Pour une nouvelle frontiere, publication dans l'arbre d'etudes sous le cas
 339   if ( !existe )
 340   {
 341     char* CaseName = GetName() ;
 342     MESSAGE ( "AddBoundaryGroup : insertion de la frontiere dans l'arbre de " << CaseName );
 343     _gen_i->PublishBoundaryUnderCase(CaseName, BoundaryName) ;
 344   }
 345   // C. Le groupe est-il deja enregistre pour une frontiere de ce cas ?
 346   for ( it = ListBoundaryGroup.begin(); it != ListBoundaryGroup.end(); it++ )
 347   {
 348     std::string boun = *it ;
 349     it++ ;
 350 //     MESSAGE ("..  Group : "<< *it );
 351     if ( *it == Group )
 352     { INFOS ("Frontiere " << boun << " Un groupe est deja associe " << Group ) ;
 353       SALOME::ExceptionStruct es;
 354       es.type = SALOME::BAD_PARAM;
 355       es.text = "Invalid AddBoundaryGroup";
 356       throw SALOME::SALOME_Exception(es);
 357       return ;
 358     }
 359   }
 360   // D. Enregistrement du couple (frontiere,groupe) dans la reference du cas
 361   myHomardCas->AddBoundaryGroup( BoundaryName, Group );
 362 }
 363 //=============================================================================
 364 HOMARD::ListBoundaryGroupType* HOMARD_Cas_i::GetBoundaryGroup()
 365 {
 366   MESSAGE ("GetBoundaryGroup");
 367   ASSERT(myHomardCas );
 368   const std::list<std::string>& ListBoundaryGroup = myHomardCas->GetBoundaryGroup();
 369   HOMARD::ListBoundaryGroupType_var aResult = new HOMARD::ListBoundaryGroupType();
 370   aResult->length( ListBoundaryGroup.size() );
 371   std::list<std::string>::const_iterator it;
 372   int i = 0;
 373   for ( it = ListBoundaryGroup.begin(); it != ListBoundaryGroup.end(); it++ )
 374   {
 375     aResult[i++] = CORBA::string_dup( (*it).c_str() );
 376   }
 377   return aResult._retn();
 378 }
 379 //=============================================================================
 380 void HOMARD_Cas_i::SupprBoundaryGroup()
 381 {
 382   MESSAGE ("SupprBoundaryGroup");
 383   ASSERT(myHomardCas );
 384   myHomardCas->SupprBoundaryGroup();
 385 }
 386 //=============================================================================
 387 void HOMARD_Cas_i::SetPyram( CORBA::Long Pyram )
 388 {
 389   MESSAGE ("SetPyram, Pyram = " << Pyram );
 390   ASSERT( myHomardCas );
 391   myHomardCas->SetPyram( Pyram );
 392 }
 393 //=============================================================================
 394 CORBA::Long HOMARD_Cas_i::GetPyram()
 395 {
 396   MESSAGE ("GetPyram");
 397   ASSERT( myHomardCas );
 398   return myHomardCas->GetPyram();
 399 }
 400 //=============================================================================
 401 void HOMARD_Cas_i::MeshInfo(CORBA::Long Qual, CORBA::Long Diam, CORBA::Long Conn, CORBA::Long Tail, CORBA::Long Inte)
 402 {
 403   MESSAGE ( "MeshInfo : information sur le maillage initial du cas" );
 404   ASSERT( myHomardCas );
 405 //
 406 // Nom de l'iteration
 407   char* IterName = GetIter0Name() ;
 408   CORBA::Long etatMenage = -1 ;
 409   CORBA::Long modeHOMARD = 7 ;
 410   CORBA::Long Option1 = 1 ;
 411   CORBA::Long Option2 = 1 ;
 412   if ( Qual != 0 ) { modeHOMARD = modeHOMARD*5 ; }
 413   if ( Diam != 0 ) { modeHOMARD = modeHOMARD*19 ; }
 414   if ( Conn != 0 ) { modeHOMARD = modeHOMARD*11 ; }
 415   if ( Tail != 0 ) { modeHOMARD = modeHOMARD*13 ; }
 416   if ( Inte != 0 ) { modeHOMARD = modeHOMARD*3 ; }
 417   CORBA::Long codret = _gen_i->Compute(IterName, etatMenage, modeHOMARD, Option1, Option2) ;
 418   MESSAGE ( "MeshInfo : codret = " << codret );
 419   return ;
 420 }
 421 //=============================================================================
 422 //=============================================================================
 423 // Liens avec les autres structures
 424 //=============================================================================
 425 //=============================================================================
 426 char* HOMARD_Cas_i::GetIter0Name()
 427 {
 428   ASSERT( myHomardCas );
 429   return CORBA::string_dup( myHomardCas->GetIter0Name().c_str() );
 430 }
 431 //=============================================================================
 432 HOMARD::HOMARD_Iteration_ptr HOMARD_Cas_i::GetIter0()
 433 {
 434 // Nom de l'iteration initiale
 435   char* Iter0Name = GetIter0Name() ;
 436   MESSAGE ( "GetIter0 : Iter0Name      = " << Iter0Name );
 437   return _gen_i->GetIteration(Iter0Name) ;
 438 }
 439 //=============================================================================
 440 HOMARD::HOMARD_Iteration_ptr HOMARD_Cas_i::NextIteration( const char* IterName )
 441 {
 442 // Nom de l'iteration parent
 443   char* NomIterParent = GetIter0Name() ;
 444   MESSAGE ( "NextIteration : IterName      = " << IterName );
 445   MESSAGE ( "NextIteration : NomIterParent = " << NomIterParent );
 446   return _gen_i->CreateIteration(IterName, NomIterParent) ;
 447 }
 448 //=============================================================================
 449 HOMARD::HOMARD_Iteration_ptr HOMARD_Cas_i::LastIteration( )
 450 {
 451   HOMARD::HOMARD_Iteration_ptr Iter ;
 452   HOMARD::listeIterFilles_var ListeIterFilles ;
 453   char* IterName ;
 454 // Iteration initiale du cas
 455   IterName = GetIter0Name() ;
 456 // On va explorer la descendance de cette iteration initiale
 457 // jusqu'a trouver celle qui n'a pas de filles
 458   int nbiterfilles = 1 ;
 459   while ( nbiterfilles == 1 )
 460   {
 461 // L'iteration associee
 462 //     MESSAGE ( ".. IterName = " << IterName );
 463     Iter = _gen_i->GetIteration(IterName) ;
 464 // Les filles de cette iteration
 465     ListeIterFilles = Iter->GetIterations() ;
 466     nbiterfilles = ListeIterFilles->length() ;
 467 //     MESSAGE ( ".. nbiterfilles = " << nbiterfilles );
 468 // S'il y a au moins 2 filles, arret : on ne sait pas faire
 469     VERIFICATION( nbiterfilles <= 1 ) ;
 470 // S'il y a une fille unique, on recupere le nom de la fille et on recommence
 471     if ( nbiterfilles == 1 )
 472     { IterName = ListeIterFilles[0] ; }
 473   }
 474 //
 475   return Iter ;
 476 }
 477 //=============================================================================
 478 void HOMARD_Cas_i::AddIteration( const char* NomIteration )
 479 {
 480   ASSERT( myHomardCas );
 481   myHomardCas->AddIteration( NomIteration );
 482 }
 483 //=============================================================================
 484 //=============================================================================
 485 // YACS
 486 //=============================================================================
 487 //=============================================================================
 488 //=============================================================================
 489 // Creation d'un schema YACS
 490 // YACSName : nom du schema
 491 // ScriptFile : nom du fichier contenant le script de lancement du calcul
 492 // DirName : le repertoire de lancement des calculs du sch?ma
 493 // MeshFile : nom du fichier contenant le maillage pour le premier calcul
 494 //=============================================================================
 495 HOMARD::HOMARD_YACS_ptr HOMARD_Cas_i::CreateYACSSchema( const char* YACSName, const char* ScriptFile, const char* DirName, const char* MeshFile )
 496 {
 497 // Nom du cas
 498   const char* CaseName = GetName() ;
 499   MESSAGE ( "CreateYACSSchema : Schema YACS pour le cas " << YACSName);
 500   MESSAGE ( "nomCas     : " << CaseName);
 501   MESSAGE ( "ScriptFile : " << ScriptFile);
 502   MESSAGE ( "DirName    : " << DirName);
 503   MESSAGE ( "MeshFile   : " << MeshFile);
 504   return _gen_i->CreateYACSSchema(YACSName, CaseName, ScriptFile, DirName, MeshFile) ;
 505 }