-// Copyright (C) 2011-2012 CEA/DEN, EDF R&D
+// Copyright (C) 2011-2014 CEA/DEN, EDF R&D
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
-// version 2.1 of the License.
+// version 2.1 of the License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
#include "HOMARD_Iteration_i.hxx"
#include "HOMARD_Boundary_i.hxx"
#include "HOMARD_Zone_i.hxx"
+#include "HOMARD_YACS_i.hxx"
#include "HomardDriver.hxx"
#include "HOMARD_DriverTools.hxx"
#include "HomardMedCommun.h"
+#include "YACSDriver.hxx"
+#include "HOMARD.hxx"
#include "HOMARD_version.h"
#include CORBA_CLIENT_HEADER(SALOME_ModuleCatalog)
#include CORBA_CLIENT_HEADER(SMESH_Gen)
+#include <cmath>
#include <stdlib.h>
#include <sys/stat.h>
+#ifndef WIN32
#include <dirent.h>
+#endif
#include <string>
+#include <cstring>
#include <iostream>
+#include <fstream>
#include <iomanip>
-#include <sys/stat.h>
#include <set>
#include <vector>
#include <stdio.h>
+#ifdef WIN32
+#include <direct.h>
+#endif
+
using namespace std;
//=============================================================================
_NS = SINGLETON_<SALOME_NamingService>::Instance();
ASSERT(SINGLETON_<SALOME_NamingService>::IsAlreadyExisting());
_NS->init_orb(_orb);
+
+ _tag_gene = 0 ;
+ _tag_boun = 0 ;
+ _tag_hypo = 0 ;
+ _tag_yacs = 0 ;
+ _tag_zone = 0 ;
+
+ SetPreferences( ) ;
}
//=================================
/*!
// Utilitaires pour l'iteration
//=============================================================================
//=============================================================================
-void HOMARD_Gen_i::SetEtatIter(const char* nomIter, const CORBA::Boolean EtatCalcul)
+void HOMARD_Gen_i::SetEtatIter(const char* nomIter, const CORBA::Long Etat)
//=====================================================================================
{
- MESSAGE( "SetEtatIter : affectation de l'etat '" << EtatCalcul << "' a l'iteration " << nomIter );
+ MESSAGE( "SetEtatIter : affectation de l'etat " << Etat << " a l'iteration " << nomIter );
HOMARD::HOMARD_Iteration_var myIteration = myContextMap[GetCurrentStudyID()]._mesIterations[nomIter];
if (CORBA::is_nil(myIteration))
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Iteration ";
+ es.text = "Invalid iteration";
throw SALOME::SALOME_Exception(es);
return ;
};
- myIteration->SetEtat(EtatCalcul);
+ myIteration->SetState(Etat);
SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
SALOMEDS::SObject_var aIterSO = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myIteration)));
- int number = myIteration->GetNumber() ;
- const char* icone ;
- if ( number == 0 )
+ std::string icone ;
+ if ( Etat <= 0 )
icone = "iter0.png" ;
- else if (EtatCalcul)
+ else if ( Etat == 2 )
icone = "iter_calculee.png" ;
else
icone = "iter_non_calculee.png" ;
- PublishInStudyAttr(aStudyBuilder, aIterSO, NULL , NULL, icone, NULL) ;
+ PublishInStudyAttr(aStudyBuilder, aIterSO, NULL , NULL, icone.c_str(), NULL) ;
aStudyBuilder->CommitCommand();
HOMARD::HOMARD_Boundary_var myBoundary = myContextMap[GetCurrentStudyID()]._mesBoundarys[BoundaryName];
if (CORBA::is_nil(myBoundary))
{
-// const char* message = "Invalid Boundary " ;
-// SALOMEException(message) ;
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Boundary ";
+ es.text = "Invalid boundary";
throw SALOME::SALOME_Exception(es);
return 1 ;
};
return 0 ;
}
//=============================================================================
-CORBA::Long HOMARD_Gen_i::DeleteCase(const char* nomCas)
+CORBA::Long HOMARD_Gen_i::DeleteCase(const char* nomCas, CORBA::Long Option)
{
// Pour detruire un cas
- MESSAGE ( "DeleteCase : nomCas = " << nomCas );
+ MESSAGE ( "DeleteCase : nomCas = " << nomCas << ", avec option = " << Option );
HOMARD::HOMARD_Cas_var myCase = myContextMap[GetCurrentStudyID()]._mesCas[nomCas];
if (CORBA::is_nil(myCase))
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Case Context ";
+ es.text = "Invalid case context";
throw SALOME::SALOME_Exception(es);
return 1;
};
// On commence par detruire toutes les iterations en partant de l'initiale et y compris elle
CORBA::String_var nomIter = myCase->GetIter0Name();
CORBA::Long Option1 = 0 ;
- CORBA::Long Option2 = 1 ;
- if ( DeleteIterationOption(nomIter, Option1, Option2) != 0 )
+ if ( DeleteIterationOption(nomIter, Option1, Option) != 0 )
{
return 2;
};
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Hypothesis ";
+ es.text = "Invalid hypothesis";
throw SALOME::SALOME_Exception(es);
return 1 ;
};
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "This hypothesis is used and cannot be deleted.";
+ es.text = "This hypothesis is used in an iteration and cannot be deleted.";
throw SALOME::SALOME_Exception(es);
return 2 ;
};
{
// Option = 0 : On ne supprime pas le fichier du maillage associe
// Option = 1 : On supprime le fichier du maillage associe
- // Option = 2 : On supprime le fichier du maillage associe, sauf si c'est la derniere iteration
// Pour detruire une iteration courante
MESSAGE ( "DeleteIteration : nomIter = " << nomIter << ", avec option = " << Option );
CORBA::Long Option1 = 1 ;
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Iteration ";
+ es.text = "Invalid iteration";
throw SALOME::SALOME_Exception(es);
return 1 ;
};
int numero = myIteration->GetNumber();
MESSAGE ( "DeleteIterationOption : numero = " << numero );
- if ( numero == 0 and Option1 == 1 )
+ if ( numero == 0 && Option1 == 1 )
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
// On arrive ici pour une iteration sans fille
MESSAGE ( "Destruction effective de " << nomIter );
// On commence par invalider l'iteration pour faire le menage des dependances
- // et eventeullement du maillage associe
- InvalideIterOption(nomIter, Option2) ;
+ // et eventuellement du maillage associe
+ int option ;
+ if ( numero == 0 ) { option = 0 ; }
+ else { option = Option2 ; }
+ InvalideIterOption(nomIter, option) ;
// Retrait dans la descendance de l'iteration parent
if ( numero > 0 )
{
- const char* nomIterationParent = myIteration->GetIterParentName();
+ std::string nomIterationParent = myIteration->GetIterParentName();
MESSAGE ( "Retrait dans la descendance de nomIterationParent " << nomIterationParent );
HOMARD::HOMARD_Iteration_var myIterationParent = myContextMap[GetCurrentStudyID()]._mesIterations[nomIterationParent];
if (CORBA::is_nil(myIterationParent))
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Iteration ";
+ es.text = "Invalid iteration";
throw SALOME::SALOME_Exception(es);
return 3 ;
};
// suppression du lien avec l'hypothese
if ( numero > 0 )
{
- const char* nomHypo = myIteration->GetHypoName();
+ std::string nomHypo = myIteration->GetHypoName();
HOMARD::HOMARD_Hypothesis_var myHypo = myContextMap[GetCurrentStudyID()]._mesHypotheses[nomHypo];
ASSERT(!CORBA::is_nil(myHypo));
myHypo->UnLinkIteration(nomIter);
return 0 ;
}
//=============================================================================
+CORBA::Long HOMARD_Gen_i::DeleteYACS(const char* nomYACS, CORBA::Long Option)
+{
+ // Option = 0 : On ne supprime pas le fichier du schema associe
+ // Option = 1 : On supprime le fichier du schema associe
+ MESSAGE ( "DeleteYACS : nomYACS = " << nomYACS << ", avec option = " << Option );
+ HOMARD::HOMARD_YACS_var myYACS = myContextMap[GetCurrentStudyID()]._mesYACSs[nomYACS];
+ if (CORBA::is_nil(myYACS))
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid schema YACS";
+ throw SALOME::SALOME_Exception(es);
+ return 1 ;
+ };
+ // Suppression eventuelle du fichier XML
+ if ( Option == 1 )
+ {
+ std::string nomFichier = myYACS->GetXMLFile();
+ std::string commande = "rm -rf " + nomFichier ;
+ MESSAGE ( "commande = " << commande );
+ if ((system(commande.c_str())) != 0)
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The xml file for the schema YACS cannot be removed." ;
+ throw SALOME::SALOME_Exception(es);
+ return 2 ;
+ }
+ }
+ // comme on a un _var comme pointeur CORBA, on ne se preoccupe pas du delete
+ myContextMap[GetCurrentStudyID()]._mesYACSs.erase(nomYACS);
+ SALOMEDS::Study::ListOfSObject_var listSO = myCurrentStudy->FindObjectByName(nomYACS, ComponentDataType());
+ SALOMEDS::SObject_var aSO =listSO[0];
+ SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
+ myCurrentStudy->NewBuilder()->RemoveObjectWithChildren(aSO);
+
+ return 0 ;
+}
+//=============================================================================
CORBA::Long HOMARD_Gen_i::DeleteZone(const char* nomZone)
{
MESSAGE ( "DeleteZone : nomZone = " << nomZone );
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Zone ";
+ es.text = "Invalid zone";
throw SALOME::SALOME_Exception(es);
return 1 ;
};
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "This zone is used and cannot be deleted.";
+ es.text = "This zone is used in an hypothesis and cannot be deleted.";
throw SALOME::SALOME_Exception(es);
return 2 ;
};
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Boundary ";
+ es.text = "Invalid boundary";
throw SALOME::SALOME_Exception(es);
return ;
}
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Hypothesis ";
+ es.text = "Invalid hypothesis";
throw SALOME::SALOME_Exception(es);
return ;
};
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Iteration ";
+ es.text = "Invalid iteration";
throw SALOME::SALOME_Exception(es);
return ;
};
if (!so->FindAttribute(anAttr, "AttributeComment")) continue;
SALOMEDS::AttributeComment_var aCommentAttr = SALOMEDS::AttributeComment::_narrow(anAttr);
std::string value (aCommentAttr->Value());
+ if(value == std::string("IterationHomard")) continue;
if(value == std::string("HypoHomard")) continue;
SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
aStudyBuilder->RemoveObject(so);
}
- int number = myIteration->GetNumber();
- if ( number > 0 )
+ int etat = myIteration->GetState();
+ if ( etat > 0 )
{
- SetEtatIter(nomIter,false);
+ SetEtatIter(nomIter,1);
const char * nomCas = myIteration->GetCaseName();
HOMARD::HOMARD_Cas_var myCase = myContextMap[GetCurrentStudyID()]._mesCas[nomCas];
if (CORBA::is_nil(myCase))
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Case Context ";
+ es.text = "Invalid case context";
throw SALOME::SALOME_Exception(es);
return ;
};
- const char* dirCase = myCase->GetDirName();
- const char* nomDir = myIteration->GetDirName();
- const char* nomFichier = myIteration->GetMeshFile();
- std::string commande= "rm -rf " + std::string(dirCase) + "/" + std::string(nomDir);
+ std::string nomDir = myIteration->GetDirName();
+ std::string nomFichier = myIteration->GetMeshFile();
+ std::string commande = "rm -rf " + std::string(nomDir);
if ( Option == 1 ) { commande = commande + ";rm -rf " + std::string(nomFichier) ; }
MESSAGE ( "commande = " << commande );
if ((system(commande.c_str())) != 0)
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Menage du repertoire de calcul impossible" ;
- throw SALOME::SALOME_Exception(es);
- return ;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The directory for the calculation cannot be cleared." ;
+ throw SALOME::SALOME_Exception(es);
+ return ;
}
// Suppression du maillage publie dans SMESH
- const char* MeshName = myIteration->GetMeshName() ;
+ std::string MeshName = myIteration->GetMeshName() ;
DeleteResultInSmesh(nomFichier, MeshName) ;
};
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Iteration ";
+ es.text = "Invalid iteration";
throw SALOME::SALOME_Exception(es);
return ;
};
SALOMEDS::AttributeComment_var aCommentAttr = SALOMEDS::AttributeComment::_narrow(anAttr);
std::string value (aCommentAttr->Value());
/* MESSAGE("... value = " << value);*/
- if( (value == std::string("logInfo")) or ( value == std::string("SummaryInfo")) )
+ if( (value == std::string("logInfo")) || ( value == std::string("SummaryInfo")) )
{
SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
aStudyBuilder->RemoveObject(so);
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Case Context ";
+ es.text = "Invalid case context";
throw SALOME::SALOME_Exception(es);
return ;
};
- const char* dirCase = myCase->GetDirName();
const char* nomDir = myIteration->GetDirName();
- std::string commande = "rm -f " + std::string(dirCase) + "/" + std::string(nomDir) + "/info* " ;
- commande += std::string(dirCase) + "/" + std::string(nomDir) + "/Liste.*info" ;
+ std::string commande = "rm -f " + std::string(nomDir) + "/info* " ;
+ commande += std::string(nomDir) + "/Liste.*info" ;
/* MESSAGE ( "commande = " << commande );*/
if ((system(commande.c_str())) != 0)
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Menage du repertoire de calcul impossible" ;
- throw SALOME::SALOME_Exception(es);
- return ;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The directory for the calculation cannot be cleared." ;
+ throw SALOME::SALOME_Exception(es);
+ return ;
+ }
+}
+//=============================================================================
+void HOMARD_Gen_i::InvalideYACS(const char* YACSName)
+{
+ MESSAGE( "InvalideYACS : YACSName = " << YACSName );
+ HOMARD::HOMARD_YACS_var myYACS = myContextMap[GetCurrentStudyID()]._mesYACSs[YACSName];
+ if (CORBA::is_nil(myYACS))
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid schema YACS";
+ throw SALOME::SALOME_Exception(es);
+ return ;
+ };
+ //
+ SALOMEDS::SObject_var aYACSSO = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myYACS)));
+ SALOMEDS::ChildIterator_var aYACS = myCurrentStudy->NewChildIterator(aYACSSO);
+ for (; aYACS->More(); aYACS->Next())
+ {
+ SALOMEDS::SObject_var so = aYACS->Value();
+ SALOMEDS::GenericAttribute_var anAttr;
+ if (!so->FindAttribute(anAttr, "AttributeComment")) continue;
+ SALOMEDS::AttributeComment_var aCommentAttr = SALOMEDS::AttributeComment::_narrow(anAttr);
+ std::string value (aCommentAttr->Value());
+ if( value == std::string("xml") )
+ {
+ SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
+ aStudyBuilder->RemoveObject(so);
+ }
+ }
+ std::string nomFichier = myYACS->GetXMLFile();
+ std::string commande = "rm -rf " + std::string(nomFichier) ;
+ MESSAGE ( "commande = " << commande );
+ if ((system(commande.c_str())) != 0)
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The XML file for the schema YACS cannot be removed." ;
+ throw SALOME::SALOME_Exception(es);
+ return ;
}
-
}
//=============================================================================
void HOMARD_Gen_i::InvalideZone(const char* ZoneName)
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Zone ";
+ es.text = "Invalid zone";
throw SALOME::SALOME_Exception(es);
return ;
};
//=============================================================================
void HOMARD_Gen_i::AssociateCaseIter(const char* nomCas, const char* nomIter, const char* labelIter)
{
- MESSAGE( "AssociateCaseIter : " << nomCas << " ," << nomIter << "," << labelIter );
+ MESSAGE( "AssociateCaseIter : " << nomCas << ", " << nomIter << ", " << labelIter );
HOMARD::HOMARD_Cas_var myCase = myContextMap[GetCurrentStudyID()]._mesCas[nomCas];
if (CORBA::is_nil(myCase))
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Case ";
- throw SALOME::SALOME_Exception(es);
- return ;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid case";
+ throw SALOME::SALOME_Exception(es);
+ return ;
};
HOMARD::HOMARD_Iteration_var myIteration = myContextMap[GetCurrentStudyID()]._mesIterations[nomIter];
if (CORBA::is_nil(myIteration))
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid iteration ";
- throw SALOME::SALOME_Exception(es);
- return ;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid iteration";
+ throw SALOME::SALOME_Exception(es);
+ return ;
};
SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
SALOMEDS::SObject_var aCasSO = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myCase)));
if (CORBA::is_nil(aCasSO))
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Case ";
- throw SALOME::SALOME_Exception(es);
- return ;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid case";
+ throw SALOME::SALOME_Exception(es);
+ return ;
};
aStudyBuilder->NewCommand();
{
MESSAGE("AssociateIterHypo : nomHypo = " << nomHypo << " nomIter = " << nomIter);
+ // Verification de l'existence de l'hypothese
HOMARD::HOMARD_Hypothesis_var myHypo = myContextMap[GetCurrentStudyID()]._mesHypotheses[nomHypo];
ASSERT(!CORBA::is_nil(myHypo));
SALOMEDS::SObject_var aHypoSO = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myHypo)));
ASSERT(!CORBA::is_nil(aHypoSO));
+ // Verification de l'existence de l'iteration
HOMARD::HOMARD_Iteration_var myIteration = myContextMap[GetCurrentStudyID()]._mesIterations[nomIter];
ASSERT(!CORBA::is_nil(myIteration));
SALOMEDS::SObject_var aIterSO = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myIteration)));
ASSERT(!CORBA::is_nil(aIterSO));
+ // Gestion de l'arbre d'etudes
SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
-
aStudyBuilder->NewCommand();
-
SALOMEDS::SObject_var aSubSO = aStudyBuilder->NewObject(aIterSO);
aStudyBuilder->Addreference(aSubSO, aHypoSO);
-
aStudyBuilder->CommitCommand();
+ // Liens reciproques
myIteration->SetHypoName(nomHypo);
myHypo->LinkIteration(nomIter);
+
+ // On stocke les noms des champ a interpoler pour le futur controle de la donnee des pas de temps
+ myIteration->SupprFieldInterps() ;
+ HOMARD::listeFieldInterpsHypo* ListField = myHypo->GetFieldInterps();
+ int numberOfFieldsx2 = ListField->length();
+ for (int iaux = 0; iaux< numberOfFieldsx2; iaux++)
+ {
+ std::string FieldName = std::string((*ListField)[iaux]) ;
+ myIteration->SetFieldInterp(FieldName.c_str()) ;
+ iaux++ ;
+ }
};
//=============================================================================
//=============================================================================
// Recuperation des listes
//=============================================================================
//=============================================================================
+HOMARD::listeBoundarys* HOMARD_Gen_i::GetAllBoundarysName()
+{
+ MESSAGE("GetAllBoundarysName");
+ IsValidStudy () ;
+
+ HOMARD::listeBoundarys_var ret = new HOMARD::listeBoundarys;
+ ret->length(myContextMap[GetCurrentStudyID()]._mesBoundarys.size());
+ std::map<std::string, HOMARD::HOMARD_Boundary_var>::const_iterator it;
+ int i = 0;
+ for (it = myContextMap[GetCurrentStudyID()]._mesBoundarys.begin();
+ it != myContextMap[GetCurrentStudyID()]._mesBoundarys.end(); it++)
+ {
+ ret[i++] = CORBA::string_dup((*it).first.c_str());
+ }
+
+ return ret._retn();
+}
+//=============================================================================
HOMARD::listeCases* HOMARD_Gen_i::GetAllCasesName()
{
MESSAGE("GetAllCasesName");
return ret._retn();
}
//=============================================================================
-HOMARD::listeZones* HOMARD_Gen_i::GetAllZonesName()
+HOMARD::listeIterations* HOMARD_Gen_i::GetAllIterationsName()
{
- MESSAGE("GetAllZonesName");
+ MESSAGE("GetAllIterationsName");
IsValidStudy () ;
- HOMARD::listeZones_var ret = new HOMARD::listeZones;
- ret->length(myContextMap[GetCurrentStudyID()]._mesZones.size());
- std::map<std::string, HOMARD::HOMARD_Zone_var>::const_iterator it;
+ HOMARD::listeIterations_var ret = new HOMARD::listeIterations;
+ ret->length(myContextMap[GetCurrentStudyID()]._mesIterations.size());
+ std::map<std::string, HOMARD::HOMARD_Iteration_var>::const_iterator it;
int i = 0;
- for (it = myContextMap[GetCurrentStudyID()]._mesZones.begin();
- it != myContextMap[GetCurrentStudyID()]._mesZones.end(); it++)
+ for (it = myContextMap[GetCurrentStudyID()]._mesIterations.begin();
+ it != myContextMap[GetCurrentStudyID()]._mesIterations.end(); it++)
{
ret[i++] = CORBA::string_dup((*it).first.c_str());
}
return ret._retn();
}
//=============================================================================
-HOMARD::listeIterations* HOMARD_Gen_i::GetAllIterationsName()
+HOMARD::listeYACSs* HOMARD_Gen_i::GetAllYACSsName()
{
- MESSAGE("GetAllIterationsName");
+ MESSAGE("GetAllYACSsName");
IsValidStudy () ;
- HOMARD::listeIterations_var ret = new HOMARD::listeIterations;
- ret->length(myContextMap[GetCurrentStudyID()]._mesIterations.size());
- std::map<std::string, HOMARD::HOMARD_Iteration_var>::const_iterator it;
+ HOMARD::listeYACSs_var ret = new HOMARD::listeYACSs;
+ ret->length(myContextMap[GetCurrentStudyID()]._mesYACSs.size());
+ std::map<std::string, HOMARD::HOMARD_YACS_var>::const_iterator it;
int i = 0;
- for (it = myContextMap[GetCurrentStudyID()]._mesIterations.begin();
- it != myContextMap[GetCurrentStudyID()]._mesIterations.end(); it++)
+ for (it = myContextMap[GetCurrentStudyID()]._mesYACSs.begin();
+ it != myContextMap[GetCurrentStudyID()]._mesYACSs.end(); it++)
{
ret[i++] = CORBA::string_dup((*it).first.c_str());
}
return ret._retn();
}
//=============================================================================
-HOMARD::listeBoundarys* HOMARD_Gen_i::GetAllBoundarysName()
+HOMARD::listeZones* HOMARD_Gen_i::GetAllZonesName()
{
- MESSAGE("GetAllBoundarysName");
+ MESSAGE("GetAllZonesName");
IsValidStudy () ;
- HOMARD::listeBoundarys_var ret = new HOMARD::listeBoundarys;
- ret->length(myContextMap[GetCurrentStudyID()]._mesBoundarys.size());
- std::map<std::string, HOMARD::HOMARD_Boundary_var>::const_iterator it;
+ HOMARD::listeZones_var ret = new HOMARD::listeZones;
+ ret->length(myContextMap[GetCurrentStudyID()]._mesZones.size());
+ std::map<std::string, HOMARD::HOMARD_Zone_var>::const_iterator it;
int i = 0;
- for (it = myContextMap[GetCurrentStudyID()]._mesBoundarys.begin();
- it != myContextMap[GetCurrentStudyID()]._mesBoundarys.end(); it++)
+ for (it = myContextMap[GetCurrentStudyID()]._mesZones.begin();
+ it != myContextMap[GetCurrentStudyID()]._mesZones.end(); it++)
{
ret[i++] = CORBA::string_dup((*it).first.c_str());
}
// Recuperation des structures identifiees par leurs noms
//=============================================================================
//=============================================================================
+HOMARD::HOMARD_Boundary_ptr HOMARD_Gen_i::GetBoundary(const char* nomBoundary)
+{
+ HOMARD::HOMARD_Boundary_var myBoundary = myContextMap[GetCurrentStudyID()]._mesBoundarys[nomBoundary];
+ ASSERT(!CORBA::is_nil(myBoundary));
+ return HOMARD::HOMARD_Boundary::_duplicate(myBoundary);
+}
+//=============================================================================
HOMARD::HOMARD_Cas_ptr HOMARD_Gen_i::GetCase(const char* nomCas)
{
HOMARD::HOMARD_Cas_var myCase = myContextMap[GetCurrentStudyID()]._mesCas[nomCas];
return HOMARD::HOMARD_Cas::_duplicate(myCase);
}
//=============================================================================
-HOMARD::HOMARD_Zone_ptr HOMARD_Gen_i::GetZone(const char* ZoneName)
-{
- HOMARD::HOMARD_Zone_var myZone = myContextMap[GetCurrentStudyID()]._mesZones[ZoneName];
- ASSERT(!CORBA::is_nil(myZone));
- return HOMARD::HOMARD_Zone::_duplicate(myZone);
-}
-//=============================================================================
HOMARD::HOMARD_Hypothesis_ptr HOMARD_Gen_i::GetHypothesis(const char* nomHypothesis)
{
HOMARD::HOMARD_Hypothesis_var myHypothesis = myContextMap[GetCurrentStudyID()]._mesHypotheses[nomHypothesis];
return HOMARD::HOMARD_Iteration::_duplicate(myIteration);
}
//=============================================================================
-HOMARD::HOMARD_Boundary_ptr HOMARD_Gen_i::GetBoundary(const char* nomBoundary)
+HOMARD::HOMARD_YACS_ptr HOMARD_Gen_i::GetYACS(const char* nomYACS)
{
- HOMARD::HOMARD_Boundary_var myBoundary = myContextMap[GetCurrentStudyID()]._mesBoundarys[nomBoundary];
- ASSERT(!CORBA::is_nil(myBoundary));
- return HOMARD::HOMARD_Boundary::_duplicate(myBoundary);
+ HOMARD::HOMARD_YACS_var myYACS = myContextMap[GetCurrentStudyID()]._mesYACSs[nomYACS];
+ ASSERT(!CORBA::is_nil(myYACS));
+ return HOMARD::HOMARD_YACS::_duplicate(myYACS);
+}
+//=============================================================================
+HOMARD::HOMARD_Zone_ptr HOMARD_Gen_i::GetZone(const char* ZoneName)
+{
+ HOMARD::HOMARD_Zone_var myZone = myContextMap[GetCurrentStudyID()]._mesZones[ZoneName];
+ ASSERT(!CORBA::is_nil(myZone));
+ return HOMARD::HOMARD_Zone::_duplicate(myZone);
}
//=============================================================================
//=============================================================================
//=============================================================================
void HOMARD_Gen_i::MeshInfo(const char* nomCas, const char* MeshName, const char* MeshFile, const char* DirName, CORBA::Long Qual, CORBA::Long Diam, CORBA::Long Conn, CORBA::Long Tail, CORBA::Long Inte)
{
- MESSAGE ( "MeshInfo : nomCas = " << nomCas << ", MeshName = " << MeshName << ", MeshFile = " << MeshFile );
- MESSAGE ( "Qual = " << Qual << ", Diam = " << Diam << ", Conn = " << Conn << ", Tail = " << Tail << ", Inte = " << Inte );
+ INFOS ( "MeshInfo : nomCas = " << nomCas << ", MeshName = " << MeshName << ", MeshFile = " << MeshFile );
+ INFOS ( "Qual = " << Qual << ", Diam = " << Diam << ", Conn = " << Conn << ", Tail = " << Tail << ", Inte = " << Inte );
IsValidStudy () ;
// Creation du cas
- HOMARD::HOMARD_Cas_ptr myCase = CreateCase(nomCas, MeshName, MeshFile) ;
+ int option = 1 ;
+ if ( _PublisMeshIN != 0 ) option = 2 ;
+ HOMARD::HOMARD_Cas_ptr myCase = CreateCase0(nomCas, MeshName, MeshFile, 1, 0, option) ;
myCase->SetDirName(DirName) ;
// Analyse
myCase->MeshInfo(Qual, Diam, Conn, Tail, Inte) ;
return aZone._retn();
}
//=============================================================================
+HOMARD::HOMARD_YACS_ptr HOMARD_Gen_i::newYACS()
+{
+ HOMARD::HOMARD_Gen_var engine = POA_HOMARD::HOMARD_Gen::_this();
+ HOMARD_YACS_i* aServant = new HOMARD_YACS_i(_orb, engine);
+ HOMARD::HOMARD_YACS_var aYACS = HOMARD::HOMARD_YACS::_narrow(aServant->_this());
+ return aYACS._retn();
+}
+//=============================================================================
//=============================================================================
//=============================================================================
//=============================================================================
//=============================================================================
HOMARD::HOMARD_Cas_ptr HOMARD_Gen_i::CreateCase(const char* nomCas, const char* MeshName, const char* MeshFile)
+//
+// Creation d'un cas initial
+// nomCas : nom du cas a creer
+// MeshName, MeshFile : nom et fichier du maillage correspondant
+//
{
- MESSAGE ( "CreateCase : nomCas = " << nomCas << ", MeshName = " << MeshName << ", MeshFile = " << MeshFile );
- IsValidStudy () ;
+ INFOS ( "CreateCase : nomCas = " << nomCas << ", MeshName = " << MeshName << ", MeshFile = " << MeshFile );
- if ((myContextMap[GetCurrentStudyID()]._mesCas).find(nomCas)!=(myContextMap[GetCurrentStudyID()]._mesCas).end())
+ int option = 1 ;
+ if ( _PublisMeshIN != 0 ) option = 2 ;
+ HOMARD::HOMARD_Cas_ptr myCase = CreateCase0(nomCas, MeshName, MeshFile, 0, 0, option) ;
+
+// Valeurs par defaut des filtrages
+ myCase->SetPyram(0);
+
+ return HOMARD::HOMARD_Cas::_duplicate(myCase);
+}
+//=============================================================================
+HOMARD::HOMARD_Cas_ptr HOMARD_Gen_i::CreateCaseFromIteration(const char* nomCas, const char* DirNameStart)
+//
+// nomCas : nom du cas a creer
+// DirNameStart : nom du repertoire contenant l'iteration de reprise
+//
+{
+ INFOS ( "CreateCaseFromIteration : nomCas = " << nomCas << ", DirNameStart = " << DirNameStart );
+ std::string nomDirWork = getenv("PWD") ;
+ int codret ;
+
+ // A. Decodage du point de reprise
+ // A.1. Controle du repertoire de depart de l'iteration
+ codret = CHDIR(DirNameStart) ;
+ if ( codret != 0 )
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "This case has already been defined";
+ es.text = "The directory of the iteration does not exist.";
throw SALOME::SALOME_Exception(es);
return 0;
};
-
- int existe = MEDFileExist ( MeshFile ) ;
- if ( existe == 0 )
+ // A.2. Reperage des fichiers du repertoire de reprise
+ std::string file_configuration = "" ;
+ std::string file_maillage_homard = "" ;
+ int bilan ;
+#ifndef WIN32
+ DIR *dp;
+ struct dirent *dirp;
+ dp = opendir(DirNameStart);
+ while ( (dirp = readdir(dp)) != NULL )
+ {
+ std::string file_name(dirp->d_name);
+// MESSAGE ( file_name );
+ bilan = file_name.find("HOMARD.Configuration.") ;
+ if ( bilan != string::npos ) { file_configuration = file_name ; }
+ bilan = file_name.find("maill.") ;
+ if ( bilan != string::npos )
+ {
+ bilan = file_name.find(".hom.med") ;
+ if ( bilan != string::npos ) { file_maillage_homard = file_name ; }
+ }
+ }
+ closedir(dp);
+#else
+ HANDLE hFind = INVALID_HANDLE_VALUE;
+ WIN32_FIND_DATA ffd;
+ hFind = FindFirstFile(DirNameStart, &ffd);
+ if (INVALID_HANDLE_VALUE != hFind) {
+ while (FindNextFile(hFind, &ffd) != 0) {
+ if (ffd.dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) continue; //skip directories
+ std::string file_name(ffd.cFileName);
+ bilan = file_name.find("HOMARD.Configuration.") ;
+ if ( bilan != string::npos ) { file_configuration = file_name ; }
+ bilan = file_name.find("maill.") ;
+ if ( bilan != string::npos )
+ {
+ bilan = file_name.find(".hom.med") ;
+ if ( bilan != string::npos ) { file_maillage_homard = file_name ; }
+ }
+ }
+ FindClose(hFind);
+ }
+#endif
+ MESSAGE ( "==> file_configuration : " << file_configuration ) ;
+ MESSAGE ( "==> file_maillage_homard : " << file_maillage_homard ) ;
+ // A.3. Controle
+ if ( ( file_configuration == "" ) || ( file_maillage_homard == "" ) )
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "The mesh file does not exist.";
+ std::string text ;
+ if ( file_configuration == "" ) { text = "The configuration file cannot be found." ; }
+ else { text = "The HOMARD mesh file cannot be found." ; }
+ es.text = CORBA::string_dup(text.c_str());
throw SALOME::SALOME_Exception(es);
- return 0;
}
- HOMARD::HOMARD_Cas_var myCase = newCase();
- myCase->SetName(nomCas);
- SALOMEDS::SObject_var aSO;
- PublishInStudy(myCurrentStudy, aSO, myCase, nomCas);
- myContextMap[GetCurrentStudyID()]._mesCas[nomCas] = myCase;
+ // B. Lecture du fichier de configuration
+ // ATTENTION : on doit veiller a la coherence entre HomardDriver et CreateCaseFromIteration
+ int NumeIter ;
+ int TypeConf = 0 ;
+ int Pyram = 0 ;
+ char* MeshName ;
+ char* MeshFile ;
+ // le constructeur de ifstream permet d'ouvrir un fichier en lecture
+ std::ifstream fichier( file_configuration.c_str() );
+ if ( fichier ) // ce test échoue si le fichier n'est pas ouvert
+ {
+ std::string ligne; // variable contenant chaque ligne lue
+ std::string mot_cle;
+ std::string argument;
+ int decalage;
+ // cette boucle sur les lignes s'arrête dès qu'une erreur de lecture survient
+ while ( std::getline( fichier, ligne ) )
+ {
+ // B.1. Pour la ligne courante, on identifie le premier mot : le mot-cle
+ std::istringstream ligne_bis(ligne); // variable contenant chaque ligne sous forme de flux
+ ligne_bis >> mot_cle ;
+ // B.2. Des valeurs entieres : le second bloc de la ligne
+ if ( mot_cle == "NumeIter" )
+ {
+ ligne_bis >> NumeIter ;
+ NumeIter += 1 ;
+ }
+ // B.3. Des valeurs caracteres brutes : le second bloc de la ligne est la valeur
+ else if ( ( mot_cle == "TypeConf" ) || ( mot_cle == "TypeElem" ) )
+ {
+ ligne_bis >> argument ;
- std::vector<double> LesExtremes =GetBoundingBoxInMedFile(MeshFile);
- HOMARD::extrema_var aSeq = new HOMARD::extrema();
- if (LesExtremes.size()!=10) { return false; }
- aSeq->length(10);
- for (int i =0; i< LesExtremes.size(); i++)
- aSeq[i]=LesExtremes[i];
- myCase->SetBoundingBox(aSeq);
-
- std::set<std::string> LesGroupes =GetListeGroupesInMedFile(MeshFile);
- HOMARD::ListGroupType_var aSeqGroupe = new HOMARD::ListGroupType;
- aSeqGroupe->length(LesGroupes.size());
- std::set<std::string>::const_iterator it;
- int i = 0;
- for (it=LesGroupes.begin() ; it != LesGroupes.end(); it++)
- aSeqGroupe[i++]=(*it).c_str();
- myCase->SetGroups(aSeqGroupe);
+ if ( mot_cle == "TypeConf" )
+ {
+ if ( argument == "conforme" ) { TypeConf = 1 ; }
+ else if ( argument == "non_conforme_1_noeud" ) { TypeConf = 2 ; }
+ else if ( argument == "non_conforme_1_arete" ) { TypeConf = 3 ; }
+ else if ( argument == "non_conforme_indicateur" ) { TypeConf = 4 ; }
+ }
+ else if ( mot_cle == "TypeElem" )
+ {
+ if ( argument == "ignore_pyra" ) { Pyram = 1 ; }
+ else if ( argument == "HOMARD" ) { Pyram = 0 ; }
+ }
+ }
+ // B.4. Des valeurs caracteres : le deuxieme bloc de la ligne peut etre encadre par des quotes :
+ // il faut les supprimer
+ else if ( ( mot_cle == "CCNoMNP1" ) || ( mot_cle == "CCMaiNP1" ) )
+ {
+ ligne_bis >> argument ;
+ if ( argument[0] == '"' ) { decalage = 1 ; }
+ else { decalage = 0 ; }
+ size_t size = argument.size() + 1 - 2*decalage ;
-// Recherche d'un nom pour l'iteration 0. Par defaut, on prend le nom
-// du maillage du cas. Si ce nom existe deja, on incremente avec 0, 1, 2, etc.
- int monNum=0;
- std::string NomIteration = std::string(MeshName) ;
- while ((myContextMap[GetCurrentStudyID()]._mesIterations).find(NomIteration) != (myContextMap[GetCurrentStudyID()]._mesIterations.end()))
+ if ( mot_cle == "CCNoMNP1" )
+ {
+ MeshName = new char[ size ];
+ strncpy( MeshName, argument.c_str()+decalage, size );
+ MeshName[size-1] = '\0' ;
+ }
+ else if ( mot_cle == "CCMaiNP1" )
+ {
+ MeshFile = new char[ size ];
+ strncpy( MeshFile, argument.c_str()+decalage, size );
+ MeshFile[size-1] = '\0' ;
+ }
+ }
+ }
+ MESSAGE ( "==> TypeConf : " << TypeConf ) ;
+ MESSAGE ( "==> MeshName : " << MeshName ) ;
+ MESSAGE ( "==> MeshFile : " << MeshFile ) ;
+ MESSAGE ( "==> NumeIter : " << NumeIter ) ;
+ MESSAGE ( "==> Pyram : " << Pyram ) ;
+ }
+ else
{
- std::ostringstream nom;
- nom << MeshName << monNum;
- NomIteration=nom.str();
- monNum = monNum+1;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ std::string text = "The configuration file cannot be read." ;
+ es.text = CORBA::string_dup(text.c_str());
+ throw SALOME::SALOME_Exception(es);
}
- MESSAGE ( "CreateCase : NomIteration = " << NomIteration );
-
- HOMARD::HOMARD_Iteration_var anIter = newIteration();
- myContextMap[GetCurrentStudyID()]._mesIterations[NomIteration] = anIter;
- std::ostringstream DirName;
- DirName << "I00";
- anIter->SetDirName(DirName.str().c_str());
- anIter->SetName(NomIteration.c_str());
- anIter->SetMeshFile(MeshFile);
- anIter->SetMeshName(MeshName);
- anIter->SetNumber(0);
+ // C. Creation effective du cas
+
+ int option = 1 ;
+ if ( _PublisMeshIN != 0 ) option = 2 ;
+ HOMARD::HOMARD_Cas_ptr myCase = CreateCase0(nomCas, MeshName, MeshFile, 1, NumeIter, option) ;
+
+ // D. Parametrages lus dans le fichier de configuration
+
+ myCase->SetConfType (TypeConf) ;
+ myCase->SetPyram (Pyram) ;
+
+ // E. Copie du fichier de maillage homard
+ // E.1. Repertoire associe au cas
+ char* nomDirCase = myCase->GetDirName() ;
+ // E.2. Repertoire associe a l'iteration de ce cas
+ char* IterName ;
+ IterName = myCase->GetIter0Name() ;
+ HOMARD::HOMARD_Iteration_var Iter = GetIteration(IterName) ;
+ char* nomDirIter = CreateDirNameIter(nomDirCase, 0 );
+ Iter->SetDirNameLoc(nomDirIter);
+ std::string nomDirIterTotal ;
+ nomDirIterTotal = std::string(nomDirCase) + "/" + std::string(nomDirIter) ;
+#ifndef WIN32
+ if (mkdir(nomDirIterTotal.c_str(), S_IRWXU|S_IRGRP|S_IXGRP) != 0)
+#else
+ if (_mkdir(nomDirIterTotal.c_str()) != 0)
+#endif
+ {
+ MESSAGE ( "nomDirIterTotal : " << nomDirIterTotal ) ;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ std::string text = "The directory for the computation cannot be created." ;
+ es.text = CORBA::string_dup(text.c_str());
+ throw SALOME::SALOME_Exception(es);
+ }
+ // E.3. Copie du maillage HOMARD au format MED
+ codret = CHDIR(DirNameStart) ;
+ std::string commande = "cp " + file_maillage_homard + " " + nomDirIterTotal ;
+ MESSAGE ( "commande : " << commande ) ;
+ codret = system(commande.c_str()) ;
+ MESSAGE ( "codret : " << codret ) ;
+ if ( codret != 0 )
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The starting point for the case cannot be copied into the working directory.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
- AssociateCaseIter (nomCas,NomIteration.c_str(),"IterationHomard");
- SetEtatIter(NomIteration.c_str(),true);
-//
- PublishResultInSmesh(MeshFile, 0);
+ // F. Menage
-// Valeurs par defaut des filtrages
- myCase->SetPyram(0);
+ delete[] MeshName ;
+ delete[] MeshFile ;
+ CHDIR(nomDirWork.c_str());
return HOMARD::HOMARD_Cas::_duplicate(myCase);
}
//=============================================================================
-HOMARD::HOMARD_Hypothesis_ptr HOMARD_Gen_i::CreateHypothesis(const char* nomHypothesis)
+HOMARD::HOMARD_Cas_ptr HOMARD_Gen_i::CreateCaseFromCaseLastIteration(const char* nomCas, const char* DirNameStart)
+//
+// nomCas : nom du cas a creer
+// DirNameStart : nom du repertoire du cas contenant l'iteration de reprise
+//
+{
+ INFOS ( "CreateCaseFromCaseLastIteration : nomCas = " << nomCas << ", DirNameStart = " << DirNameStart );
+
+ std::string DirNameStartIter = CreateCase1(DirNameStart, -1) ;
+
+ DirNameStartIter = string(DirNameStart) + "/" + DirNameStartIter ;
+ HOMARD::HOMARD_Cas_ptr myCase = CreateCaseFromIteration(nomCas, DirNameStartIter.c_str()) ;
+
+ return HOMARD::HOMARD_Cas::_duplicate(myCase);
+}
+//=============================================================================
+HOMARD::HOMARD_Cas_ptr HOMARD_Gen_i::CreateCaseFromCaseIteration(const char* nomCas, const char* DirNameStart, CORBA::Long Number)
+//
+// nomCas : nom du cas a creer
+// DirNameStart : nom du repertoire du cas contenant l'iteration de reprise
+// Number : numero de l'iteration de depart
+//
+{
+ INFOS ( "CreateCaseFromCaseIteration : nomCas = " << nomCas << ", DirNameStart = " << DirNameStart << ", Number = " << Number );
+ if ( Number < 0 )
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The number of iteration must be positive.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+
+ std::string DirNameStartIter = CreateCase1(DirNameStart, Number) ;
+
+ DirNameStartIter = string(DirNameStart) + "/" + DirNameStartIter ;
+ HOMARD::HOMARD_Cas_ptr myCase = CreateCaseFromIteration(nomCas, DirNameStartIter.c_str()) ;
+
+ return HOMARD::HOMARD_Cas::_duplicate(myCase);
+}
+//=============================================================================
+std::string HOMARD_Gen_i::CreateCase1(const char* DirNameStart, CORBA::Long Number)
+//
+// Retourne le nom du repertoire ou se trouve l'iteration voulue.
+// DirNameStart : nom du repertoire du cas contenant l'iteration de reprise
+// Number : numero de l'iteration de depart ou -1 si on cherche la derniere
+//
+{
+ MESSAGE ( "CreateCase1 : DirNameStart = " << DirNameStart << ", Number = " << Number );
+ std::string nomDirWork = getenv("PWD") ;
+ std::string DirNameStartIter ;
+ int codret ;
+ int NumeIterMax = -1 ;
+
+ // A.1. Controle du repertoire de depart du cas
+ codret = CHDIR(DirNameStart) ;
+ if ( codret != 0 )
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The directory of the case for the pursuit does not exist.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+ // A.2. Reperage des sous-repertoire du repertoire de reprise
+ bool existe = false ;
+#ifndef WIN32
+ DIR *dp;
+ struct dirent *dirp;
+ dp = opendir(DirNameStart);
+ while ( (dirp = readdir(dp)) != NULL ) {
+ std::string DirName_1(dirp->d_name);
+#else
+ HANDLE hFind = INVALID_HANDLE_VALUE;
+ WIN32_FIND_DATA ffd;
+ hFind = FindFirstFile(DirNameStart, &ffd);
+ if (INVALID_HANDLE_VALUE != hFind) {
+ while (FindNextFile(hFind, &ffd) != 0) {
+ std::string DirName_1 = "";
+ if (ffd.dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) {
+ DirName_1 = std::string(ffd.cFileName);
+ }
+#endif
+ if ( ( DirName_1 != "." ) && ( DirName_1 != ".." ) )
+ {
+ if ( CHDIR(DirName_1.c_str()) == 0 )
+ {
+// On cherche le fichier de configuration dans ce sous-repertoire
+ codret = CHDIR(DirNameStart);
+#ifndef WIN32
+ DIR *dp_1;
+ struct dirent *dirp_1;
+ dp_1 = opendir(DirName_1.c_str()) ;
+ while ( (dirp_1 = readdir(dp_1)) != NULL )
+ {
+ std::string file_name_1(dirp_1->d_name);
+#else
+ HANDLE hFind1 = INVALID_HANDLE_VALUE;
+ WIN32_FIND_DATA ffd1;
+ hFind1 = FindFirstFile(DirName_1.c_str(), &ffd1);
+ while (FindNextFile(hFind1, &ffd1) != 0)
+ {
+ if (ffd1.dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) continue; //skip directories
+ std::string file_name_1(ffd1.cFileName);
+#endif
+ int bilan = file_name_1.find("HOMARD.Configuration.") ;
+ if ( bilan != string::npos )
+ {
+ // Decodage du fichier pour trouver le numero d'iteration
+ CHDIR(DirName_1.c_str()) ;
+
+ std::ifstream fichier( file_name_1.c_str() );
+ if ( fichier ) // ce test échoue si le fichier n'est pas ouvert
+ {
+ int NumeIter ;
+ std::string ligne; // variable contenant chaque ligne lue
+ std::string mot_cle;
+ // cette boucle sur les lignes s'arrête dès qu'une erreur de lecture survient
+ while ( std::getline( fichier, ligne ) )
+ {
+ // B.1. Pour la ligne courante, on identifie le premier mot : le mot-cle
+ std::istringstream ligne_bis(ligne); // variable contenant chaque ligne sous forme de flux
+ ligne_bis >> mot_cle ;
+ if ( mot_cle == "NumeIter" )
+ {
+ ligne_bis >> NumeIter ;
+ NumeIter += 1 ;
+// MESSAGE ( "==> NumeIter : " << NumeIter ) ;
+ if ( Number == - 1 )
+ {
+ if ( NumeIter >= NumeIterMax )
+ {
+ NumeIterMax = NumeIter ;
+ DirNameStartIter = DirName_1 ;
+ }
+ }
+ else
+ {
+ if ( NumeIter == Number )
+ {
+ DirNameStartIter = DirName_1 ;
+ existe = true ;
+ break ;
+ }
+ }
+ }
+ }
+ }
+ else
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ std::string text = "The configuration file cannot be read." ;
+ es.text = CORBA::string_dup(text.c_str());
+ throw SALOME::SALOME_Exception(es);
+ }
+ CHDIR(DirNameStart) ;
+ }
+ if ( existe ) { break ; }
+ }
+#ifndef WIN32
+ closedir(dp_1);
+#else
+ FindClose(hFind1);
+#endif
+ if ( existe ) { break ; }
+ }
+ }
+ }
+#ifndef WIN32
+ closedir(dp);
+#else
+ FindClose(hFind);
+#endif
+ CHDIR(nomDirWork.c_str());
+
+ if ( ( Number >= 0 && ( !existe ) ) || ( Number < 0 && ( NumeIterMax == -1 ) ) )
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The directory of the iteration does not exist.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+
+ return DirNameStartIter ;
+}
+//=============================================================================
+HOMARD::HOMARD_Cas_ptr HOMARD_Gen_i::CreateCase0(const char* nomCas, const char* MeshName, const char* MeshFile, CORBA::Long MeshOption, CORBA::Long NumeIter, CORBA::Long Option)
+//
+// nomCas : nom du cas a creer
+// MeshName, MeshFile : nom et fichier du maillage correspondant
+// MeshOption : 0 : le maillage fourni est obligatoirement present ==> erreur si absent
+// 1 : le maillage fourni peut ne pas exister ==> on continue si absent
+// -1 : le maillage n'est pas fourni
+// NumeIter : numero de l'iteration correspondante : 0, pour un depart, n>0 pour une poursuite
+// Option : multiple de nombres premiers
+// 1 : aucune option
+// x2 : publication du maillage dans SMESH
{
- MESSAGE ( "CreateHypothesis : nomHypothesis = " << nomHypothesis );
+ MESSAGE ( "CreateCase0 : nomCas = " << nomCas );
+ MESSAGE ( "CreateCase0 : MeshName = " << MeshName << ", MeshFile = " << MeshFile << ", MeshOption = " << MeshOption );
+ MESSAGE ( "CreateCase0 : NumeIter = " << NumeIter << ", Option = " << Option );
+//
+ // A. Controles
+ // A.1. L'etude
IsValidStudy () ;
- if ((myContextMap[GetCurrentStudyID()]._mesHypotheses).find(nomHypothesis) != (myContextMap[GetCurrentStudyID()]._mesHypotheses).end())
+ // A.2. Controle du nom :
+ if ((myContextMap[GetCurrentStudyID()]._mesCas).find(nomCas)!=(myContextMap[GetCurrentStudyID()]._mesCas).end())
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This case has already been defined.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+
+ // A.3. Controle du fichier du maillage
+ int existeMeshFile ;
+ if ( MeshOption >= 0 )
{
+ existeMeshFile = MEDFileExist ( MeshFile ) ;
+ MESSAGE ( "CreateCase0 : existeMeshFile = " << existeMeshFile );
+ if ( ( existeMeshFile == 0 ) && ( MeshOption == 0 ) )
+ {
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "This hypothesis is already defined.";
+ es.text = "The mesh file does not exist.";
throw SALOME::SALOME_Exception(es);
return 0;
}
+ }
+ else { existeMeshFile = 0 ; }
+
+ // B. Creation de l'objet cas et publication
+// MESSAGE ( "CreateCase0 : Creation de l'objet" );
+ HOMARD::HOMARD_Cas_var myCase = newCase();
+ myCase->SetName(nomCas);
+ SALOMEDS::SObject_var aSO;
+ SALOMEDS::SObject_var aResultSO=PublishInStudy(myCurrentStudy, aSO, myCase, nomCas);
+ myContextMap[GetCurrentStudyID()]._mesCas[nomCas] = myCase;
+
+ // C. Caracteristiques du maillage
+ if ( existeMeshFile != 0 )
+ {
+ // Les valeurs extremes des coordonnees
+// MESSAGE ( "CreateCase0 : Les valeurs extremes des coordonnees" );
+ std::vector<double> LesExtremes =GetBoundingBoxInMedFile(MeshFile) ;
+ HOMARD::extrema_var aSeq = new HOMARD::extrema() ;
+ if (LesExtremes.size()!=10) { return false; }
+ aSeq->length(10) ;
+ for (int i =0 ; i< LesExtremes.size() ; i++)
+ aSeq[i]=LesExtremes[i] ;
+ myCase->SetBoundingBox(aSeq) ;
+ // Les groupes
+// MESSAGE ( "CreateCase0 : Les groupes" );
+ std::set<std::string> LesGroupes =GetListeGroupesInMedFile(MeshFile) ;
+ HOMARD::ListGroupType_var aSeqGroupe = new HOMARD::ListGroupType ;
+ aSeqGroupe->length(LesGroupes.size());
+ std::set<std::string>::const_iterator it ;
+ int i = 0 ;
+ for (it=LesGroupes.begin() ; it != LesGroupes.end() ; it++)
+ aSeqGroupe[i++]=(*it).c_str() ;
+ myCase->SetGroups(aSeqGroupe) ;
+ }
+
+ // D. L'iteration initiale du cas
+ MESSAGE ( "CreateCase0 : iteration initiale du cas" );
+ // D.1. Recherche d'un nom : par defaut, on prend le nom du maillage correspondant.
+ // Si ce nom d'iteration existe deja, on incremente avec 0, 1, 2, etc.
+ int monNum = 0;
+ std::string NomIteration = std::string(MeshName) ;
+ while ( (myContextMap[GetCurrentStudyID()]._mesIterations).find(NomIteration) != (myContextMap[GetCurrentStudyID()]._mesIterations.end()) )
+ {
+ std::ostringstream nom;
+ nom << MeshName << monNum;
+ NomIteration = nom.str();
+ monNum += 1;
+ }
+ MESSAGE ( "CreateCas0 : ==> NomIteration = " << NomIteration );
+
+ // D.2. Creation de l'iteration
+ HOMARD::HOMARD_Iteration_var anIter = newIteration();
+ myContextMap[GetCurrentStudyID()]._mesIterations[NomIteration] = anIter;
+ anIter->SetName(NomIteration.c_str());
+ AssociateCaseIter (nomCas, NomIteration.c_str(), "IterationHomard");
+
+ // D.4. Maillage correspondant
+ if ( existeMeshFile != 0 )
+ {
+ anIter->SetMeshFile(MeshFile);
+ if ( Option % 2 == 0 ) { PublishResultInSmesh(MeshFile, 0); }
+ }
+ anIter->SetMeshName(MeshName);
+
+ // D.5. Numero d'iteration
+ anIter->SetNumber(NumeIter);
+
+ // D.6. Etat
+ SetEtatIter(NomIteration.c_str(), -NumeIter);
+//
+
+ return HOMARD::HOMARD_Cas::_duplicate(myCase);
+}
+//=============================================================================
+HOMARD::HOMARD_Hypothesis_ptr HOMARD_Gen_i::CreateHypothesis(const char* nomHypothesis)
+{
+ INFOS ( "CreateHypothesis : nomHypothesis = " << nomHypothesis );
+ IsValidStudy () ;
+
+ // A. Controle du nom :
+ if ((myContextMap[GetCurrentStudyID()]._mesHypotheses).find(nomHypothesis) != (myContextMap[GetCurrentStudyID()]._mesHypotheses).end())
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This hypothesis has already been defined.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ }
+ // B. Creation de l'objet
HOMARD::HOMARD_Hypothesis_var myHypothesis = newHypothesis();
+ if (CORBA::is_nil(myHypothesis))
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Unable to create the hypothesis";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
myHypothesis->SetName(nomHypothesis);
+
+ // C. Enregistrement
myContextMap[GetCurrentStudyID()]._mesHypotheses[nomHypothesis] = myHypothesis;
+
SALOMEDS::SObject_var aSO;
- PublishInStudy(myCurrentStudy, aSO, myHypothesis, nomHypothesis);
+ SALOMEDS::SObject_var aResultSO=PublishInStudy(myCurrentStudy, aSO, myHypothesis, nomHypothesis);
-// Valeurs par defaut des options avancees
+ // D. Valeurs par defaut des options avancees
myHypothesis->SetNivMax(-1);
myHypothesis->SetDiamMin(-1.0);
myHypothesis->SetAdapInit(0);
- myHypothesis->SetLevelOutput(0);
+ myHypothesis->SetExtraOutput(1);
return HOMARD::HOMARD_Hypothesis::_duplicate(myHypothesis);
}
HOMARD::HOMARD_Iteration_ptr HOMARD_Gen_i::CreateIteration(const char* NomIteration, const char* nomIterParent)
//=============================================================================
{
- MESSAGE ("CreateIteration : NomIteration = " << NomIteration << ", nomIterParent = " << nomIterParent);
+ INFOS ("CreateIteration : NomIteration = " << NomIteration << ", nomIterParent = " << nomIterParent);
IsValidStudy () ;
HOMARD::HOMARD_Iteration_var myIterationParent = myContextMap[GetCurrentStudyID()]._mesIterations[nomIterParent];
if (CORBA::is_nil(myIterationParent))
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "The parent iteration is not defined.";
- throw SALOME::SALOME_Exception(es);
- return 0;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The parent iteration is not defined.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
};
const char* nomCas = myIterationParent->GetCaseName();
HOMARD::HOMARD_Cas_var myCase = myContextMap[GetCurrentStudyID()]._mesCas[nomCas];
if (CORBA::is_nil(myCase))
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Case Context ";
- throw SALOME::SALOME_Exception(es);
- return 0;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid case context";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
};
+ const char* nomDirCase = myCase->GetDirName();
+ // Controle du nom :
if ((myContextMap[GetCurrentStudyID()]._mesIterations).find(NomIteration)!=(myContextMap[GetCurrentStudyID()]._mesIterations).end())
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "This iteration is already defined. ";
- throw SALOME::SALOME_Exception(es);
- return 0;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This iteration has already been defined.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
};
HOMARD::HOMARD_Iteration_var myIteration = newIteration();
if (CORBA::is_nil(myIteration))
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Unable to create Iteration ";
- throw SALOME::SALOME_Exception(es);
- return 0;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Unable to create the iteration";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
};
myContextMap[GetCurrentStudyID()]._mesIterations[std::string(NomIteration)] = myIteration;
// Nom de l'iteration et du maillage
myIteration->SetName(NomIteration);
myIteration->SetMeshName(NomIteration);
+ myIteration->SetState(1);
int numero = myIterationParent->GetNumber() + 1;
myIteration->SetNumber(numero);
// Nombre d'iterations deja connues pour le cas, permettant
// la creation d'un sous-repertoire unique
- int nbitercase = myCase->GetNumber();
- std::ostringstream iaux ;
- iaux << std::setw(2) << std::setfill('0') << nbitercase ;
- std::stringstream DirName;
- DirName << "I" << iaux.str();
- myIteration->SetDirName(DirName.str().c_str());
+ int nbitercase = myCase->GetNumberofIter();
+ char* nomDirIter = CreateDirNameIter(nomDirCase, nbitercase );
+ myIteration->SetDirNameLoc(nomDirIter);
// Le nom du fichier du maillage MED est indice par le nombre d'iterations du cas.
// Si on a une chaine unique depuis le depart, ce nombre est le meme que le
// situation la plus frequente.
// Si on a plusieurs branches, donc des iterations du meme niveau d'adaptation, utiliser
// le nombre d'iterations du cas permet d'eviter les collisions.
- std::stringstream MeshFile;
- const char* nomDir = myCase->GetDirName();
- MeshFile << nomDir << "/maill." << iaux.str() << ".med";
- myIteration->SetMeshFile(MeshFile.str().c_str());
+ int jaux ;
+ if ( nbitercase < 100 ) { jaux = 2 ; }
+ else if ( nbitercase < 1000 ) { jaux = 3 ; }
+ else if ( nbitercase < 10000 ) { jaux = 4 ; }
+ else if ( nbitercase < 100000 ) { jaux = 5 ; }
+ else { jaux = 9 ; }
+ std::ostringstream iaux ;
+ iaux << std::setw(jaux) << std::setfill('0') << nbitercase ;
+ std::stringstream MeshFile;
+ MeshFile << nomDirCase << "/maill." << iaux.str() << ".med";
+ myIteration->SetMeshFile(MeshFile.str().c_str());
// Association avec le cas
std::string label = "IterationHomard_" + std::string(nomIterParent);
// Lien avec l'iteration precedente
myIterationParent->LinkNextIteration(NomIteration);
myIteration->SetIterParentName(nomIterParent);
+ // Gestion de l'arbre d'etudes
+ SALOMEDS::SObject_var aIterSOParent = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myIterationParent)));
+ SALOMEDS::SObject_var aIterSO = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myIteration)));
+ SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
+ aStudyBuilder->NewCommand();
+ SALOMEDS::SObject_var aSubSO = aStudyBuilder->NewObject(aIterSO);
+ aStudyBuilder->Addreference(aSubSO, aIterSOParent);
+ aStudyBuilder->CommitCommand();
return HOMARD::HOMARD_Iteration::_duplicate(myIteration);
}
MESSAGE ("CreateBoundary : BoundaryName = " << BoundaryName << ", BoundaryType = " << BoundaryType);
IsValidStudy () ;
+ // Controle du nom :
if ((myContextMap[GetCurrentStudyID()]._mesBoundarys).find(BoundaryName)!=(myContextMap[GetCurrentStudyID()]._mesBoundarys).end())
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "This boundary has already been defined";
- throw SALOME::SALOME_Exception(es);
- return 0;
+ MESSAGE ("CreateBoundary : la frontiere " << BoundaryName << " existe deja");
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This boundary has already been defined";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
};
HOMARD::HOMARD_Boundary_var myBoundary = newBoundary();
//=============================================================================
HOMARD::HOMARD_Boundary_ptr HOMARD_Gen_i::CreateBoundaryDi(const char* BoundaryName, const char* MeshName, const char* MeshFile)
{
- MESSAGE ("CreateBoundaryDi : BoundaryName = " << BoundaryName << "MeshName = " << MeshName );
+ INFOS ("CreateBoundaryDi : BoundaryName = " << BoundaryName << ", MeshName = " << MeshName );
HOMARD::HOMARD_Boundary_var myBoundary = CreateBoundary(BoundaryName, 0);
myBoundary->SetMeshFile( MeshFile ) ;
myBoundary->SetMeshName( MeshName ) ;
//=============================================================================
HOMARD::HOMARD_Boundary_ptr HOMARD_Gen_i::CreateBoundaryCylinder(const char* BoundaryName,
CORBA::Double Xcentre, CORBA::Double Ycentre, CORBA::Double Zcentre,
- CORBA::Double Xaxis, CORBA::Double Yaxis, CORBA::Double Zaxis,
+ CORBA::Double Xaxe, CORBA::Double Yaxe, CORBA::Double Zaxe,
CORBA::Double Rayon)
{
- MESSAGE ("CreateBoundaryCylinder : BoundaryName = " << BoundaryName ) ;
+ INFOS ("CreateBoundaryCylinder : BoundaryName = " << BoundaryName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Rayon <= 0.0 )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ double daux = fabs(Xaxe) + fabs(Yaxe) + fabs(Zaxe) ;
+ if ( daux < 0.0000001 )
+ { es.text = "The axis must be a non 0 vector." ;
+ error = 2 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
HOMARD::HOMARD_Boundary_var myBoundary = CreateBoundary(BoundaryName, 1) ;
- myBoundary->SetCylinder( Xcentre, Ycentre, Zcentre, Xaxis, Yaxis, Zaxis, Rayon ) ;
+ myBoundary->SetCylinder( Xcentre, Ycentre, Zcentre, Xaxe, Yaxe, Zaxe, Rayon ) ;
return HOMARD::HOMARD_Boundary::_duplicate(myBoundary) ;
}
CORBA::Double Xcentre, CORBA::Double Ycentre, CORBA::Double Zcentre,
CORBA::Double Rayon)
{
- MESSAGE ("CreateBoundarySphere : BoundaryName = " << BoundaryName ) ;
+ INFOS ("CreateBoundarySphere : BoundaryName = " << BoundaryName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Rayon <= 0.0 )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
HOMARD::HOMARD_Boundary_var myBoundary = CreateBoundary(BoundaryName, 2) ;
myBoundary->SetSphere( Xcentre, Ycentre, Zcentre, Rayon ) ;
return HOMARD::HOMARD_Boundary::_duplicate(myBoundary) ;
}
//=============================================================================
+HOMARD::HOMARD_Boundary_ptr HOMARD_Gen_i::CreateBoundaryConeA(const char* BoundaryName,
+ CORBA::Double Xaxe, CORBA::Double Yaxe, CORBA::Double Zaxe, CORBA::Double Angle,
+ CORBA::Double Xcentre, CORBA::Double Ycentre, CORBA::Double Zcentre)
+{
+ INFOS ("CreateBoundaryConeA : BoundaryName = " << BoundaryName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Angle <= 0.0 || Angle >= 90.0 )
+ { es.text = "The angle must be included higher than 0 degree and lower than 90 degrees." ;
+ error = 1 ; }
+ double daux = fabs(Xaxe) + fabs(Yaxe) + fabs(Zaxe) ;
+ if ( daux < 0.0000001 )
+ { es.text = "The axis must be a non 0 vector." ;
+ error = 2 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
+ HOMARD::HOMARD_Boundary_var myBoundary = CreateBoundary(BoundaryName, 3) ;
+ myBoundary->SetConeA( Xaxe, Yaxe, Zaxe, Angle, Xcentre, Ycentre, Zcentre ) ;
+
+ return HOMARD::HOMARD_Boundary::_duplicate(myBoundary) ;
+}
+//=============================================================================
+HOMARD::HOMARD_Boundary_ptr HOMARD_Gen_i::CreateBoundaryConeR(const char* BoundaryName,
+ CORBA::Double Xcentre1, CORBA::Double Ycentre1, CORBA::Double Zcentre1, CORBA::Double Rayon1,
+ CORBA::Double Xcentre2, CORBA::Double Ycentre2, CORBA::Double Zcentre2, CORBA::Double Rayon2)
+{
+ INFOS ("CreateBoundaryConeR : BoundaryName = " << BoundaryName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Rayon1 < 0.0 || Rayon2 < 0.0 )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ double daux = fabs(Rayon2-Rayon1) ;
+ if ( daux < 0.0000001 )
+ { es.text = "The radius must be different." ;
+ error = 2 ; }
+ daux = fabs(Xcentre2-Xcentre1) + fabs(Ycentre2-Ycentre1) + fabs(Zcentre2-Zcentre1) ;
+ if ( daux < 0.0000001 )
+ { es.text = "The centers must be different." ;
+ error = 3 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
+ HOMARD::HOMARD_Boundary_var myBoundary = CreateBoundary(BoundaryName, 4) ;
+ myBoundary->SetConeR( Xcentre1, Ycentre1, Zcentre1, Rayon1, Xcentre2, Ycentre2, Zcentre2, Rayon2 ) ;
+
+ return HOMARD::HOMARD_Boundary::_duplicate(myBoundary) ;
+}
+//=============================================================================
+HOMARD::HOMARD_Boundary_ptr HOMARD_Gen_i::CreateBoundaryTorus(const char* BoundaryName,
+ CORBA::Double Xcentre, CORBA::Double Ycentre, CORBA::Double Zcentre,
+ CORBA::Double Xaxe, CORBA::Double Yaxe, CORBA::Double Zaxe,
+ CORBA::Double RayonRev, CORBA::Double RayonPri)
+{
+ INFOS ("CreateBoundaryTorus : BoundaryName = " << BoundaryName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( ( RayonRev <= 0.0 ) || ( RayonPri <= 0.0 ) )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ double daux = fabs(Xaxe) + fabs(Yaxe) + fabs(Zaxe) ;
+ if ( daux < 0.0000001 )
+ { es.text = "The axis must be a non 0 vector." ;
+ error = 2 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
+ HOMARD::HOMARD_Boundary_var myBoundary = CreateBoundary(BoundaryName, 5) ;
+ myBoundary->SetTorus( Xcentre, Ycentre, Zcentre, Xaxe, Yaxe, Zaxe, RayonRev, RayonPri ) ;
+
+ return HOMARD::HOMARD_Boundary::_duplicate(myBoundary) ;
+}
+//=============================================================================
HOMARD::HOMARD_Zone_ptr HOMARD_Gen_i::CreateZone(const char* ZoneName, CORBA::Long ZoneType)
{
MESSAGE ("CreateZone : ZoneName = " << ZoneName << ", ZoneType = " << ZoneType);
IsValidStudy () ;
+ // Controle du nom :
if ((myContextMap[GetCurrentStudyID()]._mesZones).find(ZoneName)!=(myContextMap[GetCurrentStudyID()]._mesZones).end())
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "This zone has already been defined";
- throw SALOME::SALOME_Exception(es);
- return 0;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This zone has already been defined";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
};
HOMARD::HOMARD_Zone_var myZone = newZone();
CORBA::Double Ymini, CORBA::Double Ymaxi,
CORBA::Double Zmini, CORBA::Double Zmaxi)
{
- MESSAGE ("CreateZoneBox : ZoneName = " << ZoneName ) ;
+ INFOS ("CreateZoneBox : ZoneName = " << ZoneName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Xmini > Xmaxi )
+ { es.text = "The X coordinates are not coherent." ;
+ error = 1 ; }
+ if ( Ymini > Ymaxi )
+ { es.text = "The Y coordinates are not coherent." ;
+ error = 2 ; }
+ if ( Zmini > Zmaxi )
+ { es.text = "The Z coordinates are not coherent." ;
+ error = 3 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
HOMARD::HOMARD_Zone_var myZone = CreateZone(ZoneName, 2) ;
myZone->SetBox ( Xmini, Xmaxi, Ymini, Ymaxi, Zmini, Zmaxi) ;
HOMARD::HOMARD_Zone_ptr HOMARD_Gen_i::CreateZoneSphere(const char* ZoneName,
CORBA::Double Xcentre, CORBA::Double Ycentre, CORBA::Double Zcentre, CORBA::Double Rayon)
{
- MESSAGE ("CreateZoneSphere : ZoneName = " << ZoneName ) ;
+ INFOS ("CreateZoneSphere : ZoneName = " << ZoneName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Rayon <= 0.0 )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
HOMARD::HOMARD_Zone_var myZone = CreateZone(ZoneName, 4) ;
myZone->SetSphere( Xcentre, Ycentre, Zcentre, Rayon ) ;
CORBA::Double Xaxe, CORBA::Double Yaxe, CORBA::Double Zaxe,
CORBA::Double Rayon, CORBA::Double Haut)
{
- MESSAGE ("CreateZoneCylinder : ZoneName = " << ZoneName ) ;
+ INFOS ("CreateZoneCylinder : ZoneName = " << ZoneName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Rayon <= 0.0 )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ double daux = fabs(Xaxe) + fabs(Yaxe) + fabs(Zaxe) ;
+ if ( daux < 0.0000001 )
+ { es.text = "The axis must be a non 0 vector." ;
+ error = 2 ; }
+ if ( Haut <= 0.0 )
+ { es.text = "The height must be positive." ;
+ error = 3 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
HOMARD::HOMARD_Zone_var myZone = CreateZone(ZoneName, 5) ;
myZone->SetCylinder( Xcentre, Ycentre, Zcentre, Xaxe, Yaxe, Zaxe, Rayon, Haut ) ;
CORBA::Double Xaxe, CORBA::Double Yaxe, CORBA::Double Zaxe,
CORBA::Double Rayon, CORBA::Double Haut, CORBA::Double Rayonint)
{
- MESSAGE ("CreateZonePipe : ZoneName = " << ZoneName ) ;
+ INFOS ("CreateZonePipe : ZoneName = " << ZoneName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Rayon <= 0.0 || Rayonint <= 0.0 )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ double daux = fabs(Xaxe) + fabs(Yaxe) + fabs(Zaxe) ;
+ if ( daux < 0.0000001 )
+ { es.text = "The axis must be a non 0 vector." ;
+ error = 2 ; }
+ if ( Haut <= 0.0 )
+ { es.text = "The height must be positive." ;
+ error = 3 ; }
+ if ( Rayon <= Rayonint )
+ { es.text = "The external radius must be higher than the internal radius." ;
+ error = 4 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
HOMARD::HOMARD_Zone_var myZone = CreateZone(ZoneName, 7) ;
myZone->SetPipe( Xcentre, Ycentre, Zcentre, Xaxe, Yaxe, Zaxe, Rayon, Haut, Rayonint ) ;
CORBA::Double Vmini, CORBA::Double Vmaxi,
CORBA::Long Orient)
{
- MESSAGE ("CreateZoneBox2D : ZoneName = " << ZoneName ) ;
+ INFOS ("CreateZoneBox2D : ZoneName = " << ZoneName ) ;
// MESSAGE ("Umini = " << Umini << ", Umaxi =" << Umaxi ) ;
// MESSAGE ("Vmini = " << Vmini << ", Vmaxi =" << Vmaxi ) ;
// MESSAGE ("Orient = " << Orient ) ;
-
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Umini > Umaxi )
+ { es.text = "The first coordinates are not coherent." ;
+ error = 1 ; }
+ if ( Vmini > Vmaxi )
+ { es.text = "The second coordinates are not coherent." ;
+ error = 2 ; }
+ if ( Orient < 1 || Orient > 3 )
+ { es.text = "The orientation must be 1, 2 or 3." ;
+ error = 3 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
double Xmini, Xmaxi ;
double Ymini, Ymaxi ;
double Zmini, Zmaxi ;
Ymaxi = 0. ;
Zmini = Umini ;
Zmaxi = Umaxi ; }
- else { ASSERT( Orient >= 1 and Orient <= 3 ) ; }
+ else { VERIFICATION( (Orient>=1) && (Orient<=3) ) ; }
HOMARD::HOMARD_Zone_var myZone = CreateZone(ZoneName, 10+Orient) ;
myZone->SetBox ( Xmini, Xmaxi, Ymini, Ymaxi, Zmini, Zmaxi) ;
CORBA::Double Rayon,
CORBA::Long Orient)
{
- MESSAGE ("CreateZoneDisk : ZoneName = " << ZoneName ) ;
+ INFOS ("CreateZoneDisk : ZoneName = " << ZoneName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Rayon <= 0.0 )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ if ( Orient < 1 || Orient > 3 )
+ { es.text = "The orientation must be 1, 2 or 3." ;
+ error = 3 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
double Xcentre ;
double Ycentre ;
double Zcentre ;
{ Xcentre = Vcentre ;
Ycentre = 0. ;
Zcentre = Ucentre ; }
- else { ASSERT( Orient >= 1 and Orient <= 3 ) ; }
+ else { VERIFICATION( (Orient>=1) && (Orient<=3) ) ; }
HOMARD::HOMARD_Zone_var myZone = CreateZone(ZoneName, 30+Orient) ;
myZone->SetCylinder( Xcentre, Ycentre, Zcentre, 0., 0., 1., Rayon, 1. ) ;
CORBA::Double Rayon, CORBA::Double Rayonint,
CORBA::Long Orient)
{
- MESSAGE ("CreateZoneDiskWithHole : ZoneName = " << ZoneName ) ;
+ INFOS ("CreateZoneDiskWithHole : ZoneName = " << ZoneName ) ;
+//
+ SALOME::ExceptionStruct es;
+ int error = 0 ;
+ if ( Rayon <= 0.0 || Rayonint <= 0.0 )
+ { es.text = "The radius must be positive." ;
+ error = 1 ; }
+ if ( Orient < 1 || Orient > 3 )
+ { es.text = "The orientation must be 1, 2 or 3." ;
+ error = 3 ; }
+ if ( Rayon <= Rayonint )
+ { es.text = "The external radius must be higher than the internal radius." ;
+ error = 4 ; }
+ if ( error != 0 )
+ {
+ es.type = SALOME::BAD_PARAM;
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+//
double Xcentre ;
double Ycentre ;
double Zcentre ;
{ Xcentre = Vcentre ;
Ycentre = 0. ;
Zcentre = Ucentre ; }
- else { ASSERT( Orient >= 1 and Orient <= 3 ) ; }
+ else { VERIFICATION( (Orient>=1) && (Orient<=3) ) ; }
HOMARD::HOMARD_Zone_var myZone = CreateZone(ZoneName, 60+Orient) ;
myZone->SetPipe( Xcentre, Ycentre, Zcentre, 0., 0., 1., Rayon, 1., Rayonint ) ;
// etatMenage = 1 : destruction du repertoire d'execution
// modeHOMARD = 1 : adaptation
// != 1 : information avec les options modeHOMARD
-// Option >0 : appel depuis python
-// <0 : appel depuis GUI
-//=============================================================================
+// Option1 >0 : appel depuis python
+// <0 : appel depuis GUI
+// Option2 : multiple de nombres premiers
+// 1 : aucune option
+// x2 : publication du maillage dans SMESH
//=============================================================================
-CORBA::Long HOMARD_Gen_i::Compute(const char* NomIteration, CORBA::Long etatMenage, CORBA::Long modeHOMARD, CORBA::Long Option)
+CORBA::Long HOMARD_Gen_i::Compute(const char* NomIteration, CORBA::Long etatMenage, CORBA::Long modeHOMARD, CORBA::Long Option1, CORBA::Long Option2)
{
- MESSAGE ( "Compute : traitement de " << NomIteration << ", avec modeHOMARD = " << modeHOMARD << ", avec Option = " << Option );
+ INFOS ( "Compute : traitement de " << NomIteration << ", avec modeHOMARD = " << modeHOMARD << ", Option1 = " << Option1 << ", Option2 = " << Option2 );
// A. Prealable
int codret = 0;
HOMARD::HOMARD_Iteration_var myIteration = myContextMap[GetCurrentStudyID()]._mesIterations[NomIteration];
ASSERT(!CORBA::is_nil(myIteration));
- // A.2. Numero de l'iteration
+ // A.2. Controle de la possibilite d'agir
+ // A.2.1. Etat de l'iteration
+ int etat = myIteration->GetState();
+ MESSAGE ( "etat = "<<etat );
+ // A.2.2. On ne calcule pas l'iteration initiale, ni une iteration deja calculee
+ if ( modeHOMARD == 1 )
+ {
+ if ( etat <= 0 )
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This iteration is the first of the case and cannot be computed.";
+ throw SALOME::SALOME_Exception(es);
+ return 1 ;
+ }
+ else if ( ( etat == 2 ) & ( modeHOMARD == 1 ) )
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This iteration is already computed.";
+ throw SALOME::SALOME_Exception(es);
+ return 1 ;
+ }
+ }
+ // A.2.3. On n'analyse pas une iteration non calculee
+ else
+ {
+ if ( etat == 1 )
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This iteration is not computed.";
+ throw SALOME::SALOME_Exception(es);
+ return 1 ;
+ }
+ }
+
+ // A.3. Numero de l'iteration
// siterp1 : numero de l'iteration a traiter
// Si adaptation :
// siter : numero de l'iteration parent, ou 0 si deja au debut mais cela ne servira pas !
else
{ siter = siterp1 ; }
- // A.3. Le cas
+ // A.4. Le cas
const char* nomCas = myIteration->GetCaseName();
HOMARD::HOMARD_Cas_var myCase = myContextMap[GetCurrentStudyID()]._mesCas[nomCas];
ASSERT(!CORBA::is_nil(myCase));
// B. Les repertoires
// B.1. Le repertoire courant
- char* nomDirWork = getenv("PWD") ;
+ std::string nomDirWork = getenv("PWD") ;
// B.2. Le sous-repertoire de l'iteration a traiter
- char* DirCompute = ComputeDir(myCase, myIteration, etatMenage);
+ char* DirCompute = ComputeDirManagement(myCase, myIteration, etatMenage);
MESSAGE( ". DirCompute = " << DirCompute );
// C. Le fichier des messages
// D. On passe dans le repertoire de l'iteration a calculer
MESSAGE ( ". On passe dans DirCompute = " << DirCompute );
- chdir(DirCompute) ;
+ CHDIR(DirCompute);
// E. Les donnees de l'execution HOMARD
// E.1. L'objet du texte du fichier de configuration
HomardDriver* myDriver = new HomardDriver(siter, siterp1);
- myDriver->TexteInit(DirCompute, LogFile);
+ myDriver->TexteInit(DirCompute, LogFile, _Langue);
// E.2. Le maillage associe a l'iteration
const char* NomMesh = myIteration->GetMeshName();
iaux = 1 ;
myDriver->TexteMaillageHOMARD( DirCompute, siterp1, iaux ) ;
myDriver->TexteMaillage(NomMesh, MeshFile, 1);
- codret = ComputeAdap(myCase, myIteration, etatMenage, myDriver, Option) ;
+ codret = ComputeAdap(myCase, myIteration, etatMenage, myDriver, Option1, Option2) ;
}
else
{
DriverTexteBoundary(myCase, myDriver) ;
// E.5. Ecriture du texte dans le fichier
+ MESSAGE ( ". Ecriture du texte dans le fichier de configuration ; codret = "<<codret );
if (codret == 0)
{ myDriver->CreeFichier(); }
int codretexec = 12 ;
if (codret == 0)
{
- codretexec = myDriver->ExecuteHomard(Option);
+ codretexec = myDriver->ExecuteHomard(Option1);
//
MESSAGE ( "Erreur en executant HOMARD : " << codretexec );
+ // En mode adaptation, on ajuste l'etat de l'iteration
if ( modeHOMARD == 1 )
{
- if (codretexec == 0) { SetEtatIter(NomIteration,true); }
- else { SetEtatIter(NomIteration,false); }
+ if (codretexec == 0) { SetEtatIter(NomIteration,2); }
+ else { SetEtatIter(NomIteration,1); }
// GERALD -- QMESSAGE BOX
}
}
if ( modeHOMARD == 1 )
{
std::stringstream saux0 ;
- Commentaire = "Iteration" ;
- Commentaire += " " + siter ;
+ Commentaire = "Mesh" ;
+ Commentaire += " " + siterp1 ;
PublishFileUnderIteration(NomIteration, MeshFile, Commentaire.c_str());
- PublishResultInSmesh(MeshFile, 1);
+ if ( Option2 % 2 == 0 ) { PublishResultInSmesh(MeshFile, 1); }
}
}
- // H.3 Message d'erreur en cas de probleme
- else
+ // H.3 Message d'erreur
+ if (codretexec != 0)
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- std::string text = "Error during the adaptation.\n" ;
- try
+ std::string text = "" ;
+ // Message d'erreur en cas de probleme en adaptation
+ if ( modeHOMARD == 1 )
{
- ifstream fichier(LogFile.c_str(), ios::in);
- string ligne;
- while(getline(fichier, ligne) and (ligne != "===== HOMARD ===== STOP ====="));
- while (getline(fichier, ligne)) { text += ligne+ "\n";};
- }
- catch (...) {
- text += "no log file ....";
+ text = "Error during the adaptation.\n" ;
+ bool stopvu = false ;
+ std::ifstream fichier( LogFile.c_str() );
+ if ( fichier ) // ce test échoue si le fichier n'est pas ouvert
+ {
+ std::string ligne; // variable contenant chaque ligne lue
+ while ( std::getline( fichier, ligne ) )
+ {
+// INFOS(ligne);
+ if ( stopvu )
+ { text += ligne+ "\n"; }
+ else
+ {
+ int position = ligne.find( "===== HOMARD ===== STOP =====" ) ;
+ if ( position > 0 ) { stopvu = true ; }
+ }
+ }
+ }
}
+ text += "\n\nSee the file " + LogFile + "\n" ;
+ INFOS ( text ) ;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
es.text = CORBA::string_dup(text.c_str());
throw SALOME::SALOME_Exception(es);
+//
+ // En mode information, on force le succes pour pouvoir consulter le fichier log
+ if ( modeHOMARD != 1 ) { codretexec = 0 ; }
}
}
{
delete myDriver;
MESSAGE ( ". On retourne dans nomDirWork = " << nomDirWork );
- chdir(nomDirWork);
+
+ CHDIR(nomDirWork.c_str());
}
return codretexec ;
//=============================================================================
// Calcul d'une iteration : partie specifique a l'adaptation
//=============================================================================
-CORBA::Long HOMARD_Gen_i::ComputeAdap(HOMARD::HOMARD_Cas_var myCase, HOMARD::HOMARD_Iteration_var myIteration, CORBA::Long etatMenage, HomardDriver* myDriver, CORBA::Long Option)
+CORBA::Long HOMARD_Gen_i::ComputeAdap(HOMARD::HOMARD_Cas_var myCase, HOMARD::HOMARD_Iteration_var myIteration, CORBA::Long etatMenage, HomardDriver* myDriver, CORBA::Long Option1, CORBA::Long Option2)
{
MESSAGE ( "ComputeAdap" );
std::string siter = saux0.str() ;
if (NumeIter < 11) { siter = "0" + siter ; }
- // A.2. On ne calcule pas l iteration 0
- if ( NumeIter == 0 )
- {
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "This iteration is the first of the case and cannot be computed.";
- throw SALOME::SALOME_Exception(es);
- return 1;
- };
-
- // A.3. On verifie qu il y a une hypothese (erreur improbable);
+ // A.2. On verifie qu il y a une hypothese (erreur improbable);
const char* nomHypo = myIteration->GetHypoName();
if (std::string(nomHypo) == std::string(""))
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text= "This iteration does not have any associated hypothesis.";
+ es.text = "This iteration does not have any associated hypothesis.";
throw SALOME::SALOME_Exception(es);
return 2;
};
HOMARD::HOMARD_Iteration_var myIterationParent = myContextMap[GetCurrentStudyID()]._mesIterations[nomIterationParent];
ASSERT(!CORBA::is_nil(myIterationParent));
// Si l'iteration parent n'est pas calculee, on le fait (recursivite amont)
- if ( ! myIterationParent->GetEtat() )
+ if ( myIterationParent->GetState() == 1 )
{
- int iaux = 1 ;
- int codret = Compute(nomIterationParent, etatMenage, iaux, Option);
- if (codret != 0)
- {
- // GERALD -- QMESSAGE BOX
- ASSERT("Pb au calcul de l'iteration precedente" == 0);
- }
+ int iaux = 1 ;
+ int codret = Compute(nomIterationParent, etatMenage, iaux, Option1, Option2);
+ if (codret != 0)
+ {
+ // GERALD -- QMESSAGE BOX
+ VERIFICATION("Pb au calcul de l'iteration precedente" == 0);
+ }
};
// C. Le sous-repertoire de l'iteration precedente
- char* DirComputePa = ComputeDirPa(myCase, myIteration);
+ char* DirComputePa = ComputeDirPaManagement(myCase, myIteration);
MESSAGE( ". DirComputePa = " << DirComputePa );
// D. Les donnees de l'adaptation HOMARD
fclose(file);
if (etatMenage == 0)
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- std::string text = "MeshFile : " + std::string(MeshFile) + " already exists ";
- es.text = CORBA::string_dup(text.c_str());
- throw SALOME::SALOME_Exception(es);
- return 4;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ std::string text = "MeshFile : " + std::string(MeshFile) + " already exists ";
+ es.text = CORBA::string_dup(text.c_str());
+ throw SALOME::SALOME_Exception(es);
+ return 4;
}
else
{
- std::string commande = "rm -f " + std::string(MeshFile);
- codret = system(commande.c_str());
- if (codret != 0)
- {
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- std::string text = "PB with meshfile destruction ";
- es.text = CORBA::string_dup(text.c_str());
- throw SALOME::SALOME_Exception(es);
- return 5;
- }
+ std::string commande = "rm -f " + std::string(MeshFile);
+ codret = system(commande.c_str());
+ if (codret != 0)
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The mesh file cannot be deleted.";
+ throw SALOME::SALOME_Exception(es);
+ return 5;
}
+ }
}
// D.4. Les types de raffinement et de deraffinement
int TypeAdap = (*ListTypes)[0];
int TypeRaff = (*ListTypes)[1];
int TypeDera = (*ListTypes)[2];
+// MESSAGE ( ". TypeAdap = " << TypeAdap << ", TypeRaff = " << TypeRaff << ", TypeDera = " << TypeDera );
// E. Texte du fichier de configuration
// E.1. Incontournables du texte
MESSAGE ( ". DiamMin = " << DiamMin );
int AdapInit = myHypo->GetAdapInit();
MESSAGE ( ". AdapInit = " << AdapInit );
- int LevelOutput = myHypo->GetLevelOutput();
- MESSAGE ( ". LevelOutput = " << LevelOutput );
- myDriver->TexteAdvanced(Pyram, NivMax, DiamMin, AdapInit, LevelOutput);
+ int ExtraOutput = myHypo->GetExtraOutput();
+ MESSAGE ( ". ExtraOutput = " << ExtraOutput );
+ myDriver->TexteAdvanced(Pyram, NivMax, DiamMin, AdapInit, ExtraOutput);
+
+ // E.7. Ajout des informations sur le deroulement de l'execution
+ int MessInfo = myIteration->GetInfoCompute();
+ MESSAGE ( ". MessInfo = " << MessInfo );
+ myDriver->TexteInfoCompute(MessInfo);
return codret ;
}
//=============================================================================
+// Creation d'un nom de sous-repertoire pour l'iteration au sein d'un repertoire parent
+// nomrep : nom du repertoire parent
+// num : le nom du sous-repertoire est sous la forme 'In', n est >= num
+//=============================================================================
+char* HOMARD_Gen_i::CreateDirNameIter(const char* nomrep, CORBA::Long num )
+{
+ MESSAGE ( "CreateDirNameIter : nomrep ="<< nomrep << ", num = "<<num);
+ // On verifie que le repertoire parent existe
+ int codret = CHDIR(nomrep) ;
+ if ( codret != 0 )
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "The directory of the case does not exist.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+ std::string nomDirActuel = getenv("PWD") ;
+ std::string DirName ;
+ // On boucle sur tous les noms possibles jusqu'a trouver un nom correspondant a un repertoire inconnu
+ bool a_chercher = true ;
+ while ( a_chercher )
+ {
+ // On passe dans le repertoire parent
+
+ CHDIR(nomrep);
+ // On recherche un nom sous la forme Iabc, avec abc representant le numero
+ int jaux ;
+ if ( num < 100 ) { jaux = 2 ; }
+ else if ( num < 1000 ) { jaux = 3 ; }
+ else if ( num < 10000 ) { jaux = 4 ; }
+ else if ( num < 100000 ) { jaux = 5 ; }
+ else { jaux = 9 ; }
+ std::ostringstream iaux ;
+ iaux << std::setw(jaux) << std::setfill('0') << num ;
+ std::ostringstream DirNameA ;
+ DirNameA << "I" << iaux.str();
+ // Si on ne pas peut entrer dans le repertoire, on doit verifier
+ // que c'est bien un probleme d'absence
+ if ( CHDIR(DirNameA.str().c_str()) != 0 )
+ {
+ bool existe = false ;
+#ifndef WIN32
+ DIR *dp;
+ struct dirent *dirp;
+ dp = opendir(nomrep);
+ while ( (dirp = readdir(dp)) != NULL )
+ {
+ std::string file_name(dirp->d_name);
+#else
+ HANDLE hFind = INVALID_HANDLE_VALUE;
+ WIN32_FIND_DATA ffd;
+ hFind = FindFirstFile(nomrep, &ffd);
+ if (INVALID_HANDLE_VALUE != hFind) {
+ while (FindNextFile(hFind, &ffd) != 0) {
+ if (ffd.dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) continue; //skip directories
+ std::string file_name(ffd.cFileName);
+#endif
+ if ( file_name == DirNameA.str() ) { existe = true ; }
+ }
+#ifndef WIN32
+ closedir(dp);
+#else
+ FindClose(hFind);
+#endif
+ if ( !existe )
+ {
+ DirName = DirNameA.str() ;
+ a_chercher = false ;
+ break ;
+ }
+ }
+ num += 1;
+ }
+
+ MESSAGE ( "==> DirName = " << DirName);
+ MESSAGE ( ". On retourne dans nomDirActuel = " << nomDirActuel );
+ CHDIR(nomDirActuel.c_str());
+ return CORBA::string_dup( DirName.c_str() );
+}
+//=============================================================================
// Calcul d'une iteration : gestion du repertoire de calcul
// Si le sous-repertoire existe :
// etatMenage = 0 : on sort en erreur si le repertoire n'est pas vide
// etatMenage = 1 : on fait le menage du repertoire
// etatMenage = -1 : on ne fait rien
//=============================================================================
-char* HOMARD_Gen_i::ComputeDir(HOMARD::HOMARD_Cas_var myCase, HOMARD::HOMARD_Iteration_var myIteration, CORBA::Long etatMenage)
+char* HOMARD_Gen_i::ComputeDirManagement(HOMARD::HOMARD_Cas_var myCase, HOMARD::HOMARD_Iteration_var myIteration, CORBA::Long etatMenage)
{
- MESSAGE ( "ComputeDir : repertoires pour le calcul" );
+ MESSAGE ( "ComputeDirManagement : repertoires pour le calcul" );
// B.2. Le repertoire du cas
const char* nomDirCase = myCase->GetDirName();
MESSAGE ( ". nomDirCase = " << nomDirCase );
// B.3. Le sous-repertoire de l'iteration a calculer, puis le repertoire complet a creer
// B.3.1. Le nom du sous-repertoire
- const char* nomDirIt = myIteration->GetDirName();
+ const char* nomDirIt = myIteration->GetDirNameLoc();
// B.3.2. Le nom complet du sous-repertoire
std::stringstream DirCompute ;
MESSAGE (". DirCompute = " << DirCompute.str() );
// B.3.3. Si le sous-repertoire n'existe pas, on le cree
- if (chdir(DirCompute.str().c_str()) != 0)
+ if (CHDIR(DirCompute.str().c_str()) != 0)
{
-// Creation du repertoire car il n'existe pas :
+#ifndef WIN32
if (mkdir(DirCompute.str().c_str(), S_IRWXU|S_IRGRP|S_IXGRP) != 0)
+#else
+ if (_mkdir(DirCompute.str().c_str()) != 0)
+#endif
{
// GERALD -- QMESSAGE BOX
std::cerr << "Pb Creation du repertoire DirCompute = " << DirCompute.str() << std::endl;
- ASSERT("Pb a la creation du repertoire" == 0);
+ VERIFICATION("Pb a la creation du repertoire" == 0);
}
}
else
// On demande de faire le menage de son contenu :
if (etatMenage == 1)
{
- MESSAGE (". Menage du repertoire DirCompute = " << DirCompute.str());
- std::string commande= "rm -rf " + DirCompute.str()+"/*" ;
- int codret = system(commande.c_str());
- if (codret != 0)
- {
- // GERALD -- QMESSAGE BOX
- std::cerr << ". Menage du repertoire de calcul" << DirCompute.str() << std::endl;
- ASSERT("Pb au menage du repertoire de calcul" == 0);
- }
+ MESSAGE (". Menage du repertoire DirCompute = " << DirCompute.str());
+ std::string commande = "rm -rf " + DirCompute.str()+"/*" ;
+ int codret = system(commande.c_str());
+ if (codret != 0)
+ {
+ // GERALD -- QMESSAGE BOX
+ std::cerr << ". Menage du repertoire de calcul" << DirCompute.str() << std::endl;
+ VERIFICATION("Pb au menage du repertoire de calcul" == 0);
+ }
}
// On n'a pas demande de faire le menage de son contenu : on sort en erreur :
else
{
if (etatMenage == 0)
{
+#ifndef WIN32
DIR *dp;
struct dirent *dirp;
dp = opendir(DirCompute.str().c_str());
bool result = true;
while ((dirp = readdir(dp)) != NULL && result )
{
- std::string file_name(dirp->d_name);
- result = file_name.empty() || file_name == "." || file_name == ".."; //if any file - break and return false
+ std::string file_name(dirp->d_name);
+ result = file_name.empty() || file_name == "." || file_name == ".."; //if any file - break and return false
}
closedir(dp);
+#else
+ HANDLE hFind = INVALID_HANDLE_VALUE;
+ WIN32_FIND_DATA ffd;
+ hFind = FindFirstFile(DirCompute.str().c_str(), &ffd);
+ bool result = true;
+ if (INVALID_HANDLE_VALUE != hFind) {
+ while (FindNextFile(hFind, &ffd) != 0) {
+ if (ffd.dwFileAttributes & FILE_ATTRIBUTE_DIRECTORY) continue; //skip directories
+ std::string file_name(ffd.cFileName);
+ result = file_name.empty() || file_name == "." || file_name == ".."; //if any file - break and return false
+ }
+ }
+ FindClose(hFind);
+#endif
if ( result == false)
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- std::string text = "Directory : " + DirCompute.str() + "is not empty";
- es.text = CORBA::string_dup(text.c_str());
- throw SALOME::SALOME_Exception(es);
- ASSERT("Directory is not empty" == 0);
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ std::string text = "Directory : " + DirCompute.str() + " is not empty";
+ es.text = CORBA::string_dup(text.c_str());
+ throw SALOME::SALOME_Exception(es);
+ VERIFICATION("Directory is not empty" == 0);
}
}
}
//=============================================================================
// Calcul d'une iteration : gestion du repertoire de calcul de l'iteration parent
//=============================================================================
-char* HOMARD_Gen_i::ComputeDirPa(HOMARD::HOMARD_Cas_var myCase, HOMARD::HOMARD_Iteration_var myIteration)
+char* HOMARD_Gen_i::ComputeDirPaManagement(HOMARD::HOMARD_Cas_var myCase, HOMARD::HOMARD_Iteration_var myIteration)
{
- MESSAGE ( "ComputeDirPa : repertoires pour le calcul" );
+ MESSAGE ( "ComputeDirPaManagement : repertoires pour le calcul" );
// Le repertoire du cas
const char* nomDirCase = myCase->GetDirName();
MESSAGE ( ". nomDirCase = " << nomDirCase );
const char* nomIterationParent = myIteration->GetIterParentName();
HOMARD::HOMARD_Iteration_var myIterationParent = myContextMap[GetCurrentStudyID()]._mesIterations[nomIterationParent];
- const char* nomDirItPa = myIterationParent->GetDirName();
+ const char* nomDirItPa = myIterationParent->GetDirNameLoc();
std::stringstream DirComputePa ;
DirComputePa << nomDirCase << "/" << nomDirItPa;
MESSAGE( ". nomDirItPa = " << nomDirItPa);
MESSAGE ( "... ZoneType = " << ZoneType << ", TypeUse = "<<TypeUse);
NumZone = iaux/2 + 1 ;
HOMARD::double_array* zone = myZone->GetCoords();
- if ( ZoneType == 2 or ( ZoneType>=11 and ZoneType <=13 ) ) // Cas d un parallelepipede ou d'un rectangle
- {
- myDriver->TexteZone(NumZone, ZoneType, TypeUse, (*zone)[0], (*zone)[1], (*zone)[2], (*zone)[3], (*zone)[4], (*zone)[5], 0., 0., 0.);
- }
+ if ( ZoneType == 2 || ( ZoneType>=11 && ZoneType <=13 ) ) // Cas d un parallelepipede ou d'un rectangle
+ { myDriver->TexteZone(NumZone, ZoneType, TypeUse, (*zone)[0], (*zone)[1], (*zone)[2], (*zone)[3], (*zone)[4], (*zone)[5], 0., 0., 0.); }
else if ( ZoneType == 4 ) // Cas d une sphere
- {
- myDriver->TexteZone(NumZone, ZoneType, TypeUse, (*zone)[0], (*zone)[1], (*zone)[2], (*zone)[3], 0., 0., 0., 0., 0.);
- }
- else if ( ZoneType == 5 or ( ZoneType>=31 and ZoneType <=33 ) ) // Cas d un cylindre ou d'un disque
- {
- myDriver->TexteZone(NumZone, ZoneType, TypeUse, (*zone)[0], (*zone)[1], (*zone)[2], (*zone)[3], (*zone)[4], (*zone)[5], (*zone)[6], (*zone)[7], 0.);
- }
- else if ( ZoneType == 7 or ( ZoneType>=61 and ZoneType <=63 ) ) // Cas d un tuyau ou disque perce
- {
- myDriver->TexteZone(NumZone, ZoneType, TypeUse, (*zone)[0], (*zone)[1], (*zone)[2], (*zone)[3], (*zone)[4], (*zone)[5], (*zone)[6], (*zone)[7], (*zone)[8]);
- }
- else { ASSERT("ZoneType est incorrect." == 0) ; }
+ { myDriver->TexteZone(NumZone, ZoneType, TypeUse, (*zone)[0], (*zone)[1], (*zone)[2], (*zone)[3], 0., 0., 0., 0., 0.); }
+ else if ( ZoneType == 5 || ( ZoneType>=31 && ZoneType <=33 ) ) // Cas d un cylindre ou d'un disque
+ { myDriver->TexteZone(NumZone, ZoneType, TypeUse, (*zone)[0], (*zone)[1], (*zone)[2], (*zone)[3], (*zone)[4], (*zone)[5], (*zone)[6], (*zone)[7], 0.); }
+ else if ( ZoneType == 7 || ( ZoneType>=61 && ZoneType <=63 ) ) // Cas d un tuyau ou disque perce
+ { myDriver->TexteZone(NumZone, ZoneType, TypeUse, (*zone)[0], (*zone)[1], (*zone)[2], (*zone)[3], (*zone)[4], (*zone)[5], (*zone)[6], (*zone)[7], (*zone)[8]); }
+ else { VERIFICATION("ZoneType est incorrect." == 0) ; }
iaux += 1 ;
}
return ;
{
// GERALD -- QMESSAGE BOX
std::cerr << "Le fichier du champ n'a pas ete fourni." << std::endl;
- ASSERT("The file for the field is not given." == 0);
+ VERIFICATION("The file for the field is not given." == 0);
}
-// Les caracteristiques d'instants
+// Les caracteristiques d'instants du champ de pilotage
int TimeStep = myIteration->GetTimeStep();
MESSAGE( ". TimeStep = " << TimeStep );
int Rank = myIteration->GetRank();
MESSAGE( ". Rank = " << Rank );
-// Les informations sur les champ
+// Les informations sur les champs
HOMARD::InfosHypo* aInfosHypo = myHypo->GetField();
// Le nom
const char* FieldName = aInfosHypo->FieldName;
myDriver->TexteField(FieldName, FieldFile, TimeStep, Rank, TypeThR, ThreshR, TypeThC, ThreshC, UsField, UsCmpI);
//
// Les composantes
- HOMARD::listeComposantsHypo* mescompo = myHypo->GetListComp();
+ HOMARD::listeComposantsHypo* mescompo = myHypo->GetComps();
int numberOfCompos = mescompo->length();
MESSAGE( ". numberOfCompos = " << numberOfCompos );
for (int NumeComp = 0; NumeComp< numberOfCompos; NumeComp++)
}
//=============================================================================
// Calcul d'une iteration : ecriture des frontieres dans le fichier de configuration
+// On ecrit dans l'ordre :
+// 1. la definition des frontieres
+// 2. les liens avec les groupes
+// 3. un entier resumant le type de comportement pour les frontieres
//=============================================================================
void HOMARD_Gen_i::DriverTexteBoundary(HOMARD::HOMARD_Cas_var myCase, HomardDriver* myDriver)
{
MESSAGE ( "... DriverTexteBoundary" );
- // On ecrit d'abord la definition des frontieres, puis les liens avec les groupes
+ // 1. Recuperation des frontieres
std::list<std::string> ListeBoundaryTraitees ;
HOMARD::ListBoundaryGroupType* ListBoundaryGroupType = myCase->GetBoundaryGroup();
int numberOfitems = ListBoundaryGroupType->length();
MESSAGE ( "... number of string for Boundary+Group = " << numberOfitems);
int BoundaryOption = 1 ;
+ // 2. Parcours des frontieres pour ecrire leur description
int NumBoundaryAnalytical = 0 ;
for (int NumBoundary = 0; NumBoundary< numberOfitems; NumBoundary=NumBoundary+2)
{
std::string BoundaryName = std::string((*ListBoundaryGroupType)[NumBoundary]);
MESSAGE ( "... BoundaryName = " << BoundaryName);
+ // 2.1. La frontiere a-t-elle deja ete ecrite ?
+ // Cela arrive quand elle est liee a plusieurs groupes. Il ne faut l'ecrire que la premiere fois
int A_faire = 1 ;
std::list<std::string>::const_iterator it = ListeBoundaryTraitees.begin();
while (it != ListeBoundaryTraitees.end())
{
- MESSAGE ( "... BoundaryNameTraitee = " << *it);
+ MESSAGE ( "..... BoundaryNameTraitee = " << *it);
if ( BoundaryName == *it ) { A_faire = 0 ; }
it++;
}
+ // 2.2. Ecriture de la frontiere
if ( A_faire == 1 )
{
-// Caracteristiques de la frontiere
+ // 2.2.1. Caracteristiques de la frontiere
HOMARD::HOMARD_Boundary_var myBoundary = myContextMap[GetCurrentStudyID()]._mesBoundarys[BoundaryName];
ASSERT(!CORBA::is_nil(myBoundary));
int BoundaryType = myBoundary->GetType();
MESSAGE ( "... BoundaryType = " << BoundaryType );
-// Ecriture selon le type
- if (BoundaryType == 0) // Cas d une frontiere discrete
+ // 2.2.2. Ecriture selon le type
+ // 2.2.2.1. Cas d une frontiere discrete
+ if (BoundaryType == 0)
{
const char* MeshName = myBoundary->GetMeshName() ;
+ MESSAGE ( ". MeshName = " << MeshName );
const char* MeshFile = myBoundary->GetMeshFile() ;
+ MESSAGE ( ". MeshFile = " << MeshFile );
myDriver->TexteBoundaryDi( MeshName, MeshFile);
- BoundaryOption = BoundaryOption*2 ;
+ if ( BoundaryOption % 2 != 0 ) { BoundaryOption = BoundaryOption*2 ; }
}
- else // Cas d une frontiere analytique
+ // 2.2.2.1. Cas d une frontiere analytique
+ else
{
NumBoundaryAnalytical++ ;
HOMARD::double_array* coor = myBoundary->GetCoords();
if (BoundaryType == 1) // Cas d un cylindre
{
- myDriver->TexteBoundaryAn(BoundaryName, NumBoundaryAnalytical, BoundaryType, (*coor)[0], (*coor)[1], (*coor)[2], (*coor)[3], (*coor)[4], (*coor)[5], (*coor)[6]);
- BoundaryOption = BoundaryOption*3 ;
+ myDriver->TexteBoundaryAn(BoundaryName, NumBoundaryAnalytical, BoundaryType, (*coor)[0], (*coor)[1], (*coor)[2], (*coor)[3], (*coor)[4], (*coor)[5], (*coor)[6], 0.);
+ if ( BoundaryOption % 3 != 0 ) { BoundaryOption = BoundaryOption*3 ; }
}
else if (BoundaryType == 2) // Cas d une sphere
{
- myDriver->TexteBoundaryAn(BoundaryName, NumBoundaryAnalytical, BoundaryType, (*coor)[0], (*coor)[1], (*coor)[2], (*coor)[3], 0., 0., 0.);
- BoundaryOption = BoundaryOption*3 ;
+ myDriver->TexteBoundaryAn(BoundaryName, NumBoundaryAnalytical, BoundaryType, (*coor)[0], (*coor)[1], (*coor)[2], (*coor)[3], 0., 0., 0., 0.);
+ if ( BoundaryOption % 3 != 0 ) { BoundaryOption = BoundaryOption*3 ; }
+ }
+ else if (BoundaryType == 3) // Cas d un cone defini par un axe et un angle
+ {
+ myDriver->TexteBoundaryAn(BoundaryName, NumBoundaryAnalytical, BoundaryType, (*coor)[0], (*coor)[1], (*coor)[2], (*coor)[3], (*coor)[4], (*coor)[5], (*coor)[6], 0.);
+ if ( BoundaryOption % 3 != 0 ) { BoundaryOption = BoundaryOption*3 ; }
+ }
+ else if (BoundaryType == 4) // Cas d un cone defini par les 2 rayons
+ {
+ myDriver->TexteBoundaryAn(BoundaryName, NumBoundaryAnalytical, BoundaryType, (*coor)[0], (*coor)[1], (*coor)[2], (*coor)[3], (*coor)[4], (*coor)[5], (*coor)[6], (*coor)[7]);
+ if ( BoundaryOption % 3 != 0 ) { BoundaryOption = BoundaryOption*3 ; }
+ }
+ else if (BoundaryType == 5) // Cas d un tore
+ {
+ myDriver->TexteBoundaryAn(BoundaryName, NumBoundaryAnalytical, BoundaryType, (*coor)[0], (*coor)[1], (*coor)[2], (*coor)[3], (*coor)[4], (*coor)[5], (*coor)[6], (*coor)[7]);
+ if ( BoundaryOption % 3 != 0 ) { BoundaryOption = BoundaryOption*3 ; }
}
}
-// Memorisation du traitement
+ // 2.2.3. Memorisation du traitement
ListeBoundaryTraitees.push_back( BoundaryName );
}
}
+ // 3. Parcours des frontieres pour ecrire les liens avec les groupes
NumBoundaryAnalytical = 0 ;
for (int NumBoundary = 0; NumBoundary< numberOfitems; NumBoundary=NumBoundary+2)
{
ASSERT(!CORBA::is_nil(myBoundary));
int BoundaryType = myBoundary->GetType();
MESSAGE ( "... BoundaryType = " << BoundaryType );
-// Recuperation du nom du groupe
+ // 3.1. Recuperation du nom du groupe
std::string GroupName = std::string((*ListBoundaryGroupType)[NumBoundary+1]);
MESSAGE ( "... GroupName = " << GroupName);
- if (BoundaryType == 0) // Cas d une frontiere discrete
+ // 3.2. Cas d une frontiere discrete
+ if ( BoundaryType == 0 )
{
if ( GroupName.size() > 0 ) { myDriver->TexteBoundaryDiGr ( GroupName ) ; }
}
- else // Cas d une frontiere analytique
+ // 3.3. Cas d une frontiere analytique
+ else
{
NumBoundaryAnalytical++ ;
myDriver->TexteBoundaryAnGr ( BoundaryName, NumBoundaryAnalytical, GroupName ) ;
}
}
+ // 4. Ecriture de l'option finale
myDriver->TexteBoundaryOption(BoundaryOption);
+//
return ;
}
//=============================================================================
{
MESSAGE ( "... DriverTexteFieldInterp" );
int TypeFieldInterp = myHypo->GetTypeFieldInterp();
+ MESSAGE ( "... TypeFieldInterp = " << TypeFieldInterp);
if (TypeFieldInterp != 0)
{
// Le fichier des champs
if (strlen(FieldFile) == 0)
{
// GERALD -- QMESSAGE BOX
- std::cerr << "Le fichier du champ n'a pas ete fourni." << std::endl;
- ASSERT("The file for the field is not given." == 0);
- }
- // Les caracteristiques d'instants
- int TimeStep = myIteration->GetTimeStep();
- MESSAGE( ". TimeStep = " << TimeStep );
- int Rank = myIteration->GetRank();
- MESSAGE( ". Rank = " << Rank );
+ VERIFICATION("The file for the field is not given." == 0);
+ }
//
const char* MeshFile = myIteration->GetMeshFile();
- myDriver->TexteFieldInterp(TypeFieldInterp, FieldFile, MeshFile, TimeStep, Rank);
- // Les champs
- if (TypeFieldInterp == 2)
+ myDriver->TexteFieldInterp(FieldFile, MeshFile);
+
+ // Les champs
+ // Interpolation de tous les champs
+ if ( TypeFieldInterp == 1 )
{
- HOMARD::listFieldInterpHypo* meschamps = myHypo->GetListFieldInterp();
- int numberOfFields = meschamps->length();
- MESSAGE( ". numberOfFields = " << numberOfFields );
- for (int NumeChamp = 0; NumeChamp< numberOfFields; NumeChamp++)
+ myDriver->TexteFieldInterpAll();
+ }
+ // Interpolation de certains champs
+ else if (TypeFieldInterp == 2)
+ {
+ // Les champs et leurs instants pour l'iteration
+ HOMARD::listeFieldInterpTSRsIter* ListFieldTSR = myIteration->GetFieldInterpsTimeStepRank();
+ int numberOfFieldsx3 = ListFieldTSR->length();
+ MESSAGE( ". pour iteration, numberOfFields = " << numberOfFieldsx3/3 );
+ // Les champs pour l'hypothese
+ HOMARD::listeFieldInterpsHypo* ListField = myHypo->GetFieldInterps();
+ int numberOfFieldsx2 = ListField->length();
+ MESSAGE( ". pour hypothese, numberOfFields = " << numberOfFieldsx2/2 );
+ // On parcourt tous les champs de l'hypothese
+ int NumField = 0 ;
+ for (int iaux = 0; iaux< numberOfFieldsx2; iaux++)
{
- std::string nomChamp = std::string((*meschamps)[NumeChamp]);
- MESSAGE( "... nomChamp = " << nomChamp );
- myDriver->TexteFieldInterpName(NumeChamp, nomChamp);
+ // Le nom du champ
+ std::string FieldName = std::string((*ListField)[iaux]) ;
+ // Le type d'interpolation
+ std::string TypeInterpstr = std::string((*ListField)[iaux+1]) ;
+ MESSAGE( "... FieldName = " << FieldName << ", TypeInterp = " << TypeInterpstr );
+ // On cherche à savoir si des instants ont été précisés pour cette itération
+ int tsrvu = 0;
+ for (int jaux = 0; jaux< numberOfFieldsx3; jaux++)
+ {
+ // Le nom du champ
+ std::string FieldName2 = std::string((*ListFieldTSR)[jaux]) ;
+ MESSAGE( "..... FieldName2 = " << FieldName2 );
+ // Quand c'est le bon champ, on ecrit le pas de temps
+ if ( FieldName == FieldName2 )
+ {
+ tsrvu = 1 ;
+ // Le pas de temps
+ std::string TimeStepstr = std::string((*ListFieldTSR)[jaux+1]) ;
+ // Le numero d'ordre
+ std::string Rankstr = std::string((*ListFieldTSR)[jaux+2]) ;
+ MESSAGE( "..... TimeStepstr = " << TimeStepstr <<", Rankstr = "<<Rankstr );
+ NumField += 1 ;
+ int TimeStep = atoi( TimeStepstr.c_str() );
+ int Rank = atoi( Rankstr.c_str() );
+ myDriver->TexteFieldInterpNameType(NumField, FieldName, TypeInterpstr, TimeStep, Rank);
+ }
+ jaux += 2 ;
+ }
+ // Si aucun instant n'a été défini
+ if ( tsrvu == 0 )
+ {
+ NumField += 1 ;
+ myDriver->TexteFieldInterpNameType(NumField, FieldName, TypeInterpstr, -1, -1);
+ }
+ iaux++ ;
}
}
}
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Study Context ";
+ es.text = "Invalid study context";
throw SALOME::SALOME_Exception(es);
return 0;
};
// Recuperation de l'objet correspondant, en essayant chacun des types possibles
// Rq : Iteration est publiee ailleurs
+ HOMARD::HOMARD_Boundary_var aBoundary = HOMARD::HOMARD_Boundary::_narrow(theObject);
HOMARD::HOMARD_Cas_var aCase = HOMARD::HOMARD_Cas::_narrow(theObject);
HOMARD::HOMARD_Hypothesis_var aHypo = HOMARD::HOMARD_Hypothesis::_narrow(theObject);
+ HOMARD::HOMARD_YACS_var aYACS = HOMARD::HOMARD_YACS::_narrow(theObject);
HOMARD::HOMARD_Zone_var aZone = HOMARD::HOMARD_Zone::_narrow(theObject);
- HOMARD::HOMARD_Boundary_var aBoundary = HOMARD::HOMARD_Boundary::_narrow(theObject);
addInStudy(theStudy);
// Controle de la non publication d'un objet de meme nom
- if ((!aHypo->_is_nil()) or (!aZone->_is_nil()) or (!aBoundary->_is_nil()))
+ if ( (!aBoundary->_is_nil()) || (!aHypo->_is_nil()) || (!aYACS->_is_nil()) || (!aZone->_is_nil()) )
+ {
+ SALOMEDS::Study::ListOfSObject_var listSO = theStudy->FindObjectByName(theName, ComponentDataType());
+ if (listSO->length() >= 1)
{
- SALOMEDS::Study::ListOfSObject_var listSO = theStudy->FindObjectByName(theName, ComponentDataType());
- if (listSO->length() >= 1)
- {
- MESSAGE("This name "<<theName<<" is already used "<<listSO->length()<<" time(s)");
- std::cerr <<"This name "<<theName<<" is already used "<<listSO->length()<<" time(s)" << std::endl;
- aResultSO = listSO[0];
- return aResultSO._retn();
- }
+ MESSAGE("This name "<<theName<<" is already used "<<listSO->length()<<" time(s)");
+ aResultSO = listSO[0];
+ return aResultSO._retn();
}
+ }
// Caracteristiques de l'etude
SALOMEDS::StudyBuilder_var aStudyBuilder = theStudy->NewBuilder();
aStudyBuilder->NewCommand();
- if(!aCase->_is_nil())
+ if(!aBoundary->_is_nil())
+ aResultSO = PublishBoundaryInStudy(theStudy, aStudyBuilder, aBoundary, theName);
+ else if(!aCase->_is_nil())
aResultSO = PublishCaseInStudy(theStudy, aStudyBuilder, aCase, theName);
else if(!aHypo->_is_nil())
aResultSO = PublishHypotheseInStudy(theStudy, aStudyBuilder, aHypo, theName);
+ else if(!aYACS->_is_nil())
+ aResultSO = PublishYACSInStudy(theStudy, aStudyBuilder, aYACS, theName);
else if(!aZone->_is_nil())
aResultSO = PublishZoneInStudy(theStudy, aStudyBuilder, aZone, theName);
- else if(!aBoundary->_is_nil())
- aResultSO = PublishBoundaryInStudy(theStudy, aStudyBuilder, aBoundary, theName);
- aStudyBuilder->CommitCommand();
+ aStudyBuilder->CommitCommand();
+
return aResultSO._retn();
};
//=============================================================================
+SALOMEDS::SObject_ptr HOMARD_Gen_i::PublishBoundaryInStudy(SALOMEDS::Study_ptr theStudy,
+ SALOMEDS::StudyBuilder_var aStudyBuilder,
+ HOMARD::HOMARD_Boundary_ptr theObject, const char* theName)
+{
+ MESSAGE("PublishBoundaryStudy pour "<<theName);
+ SALOMEDS::SObject_var aResultSO;
+
+ // Caracteristique de la Boundary
+ HOMARD::HOMARD_Boundary_var myBoundary = myContextMap[GetCurrentStudyID()]._mesBoundarys[theName];
+
+ // On recupere le module pere dans l etude
+ SALOMEDS::SComponent_var theFatherHomard = theStudy->FindComponent(ComponentDataType());
+ if (theFatherHomard->_is_nil())
+ {
+ MESSAGE("theFatherHomard->_is_nil()");
+ return aResultSO._retn();
+ }
+
+ // On ajoute la categorie des boundarys dans l etude si necessaire
+ if ( _tag_boun == 0 )
+ {
+ _tag_gene += 1 ;
+ _tag_boun = _tag_gene ;
+ }
+ MESSAGE("PublishBoundaryInStudy _tag_gene = "<<_tag_gene << ", _tag_boun = "<<_tag_boun );
+ SALOMEDS::SObject_var aSObject;
+ if (!theFatherHomard->FindSubObject(_tag_boun, aSObject))
+ {
+ MESSAGE("Ajout de la categorie des boundarys");
+ aSObject = aStudyBuilder->NewObjectToTag(theFatherHomard, _tag_boun);
+ PublishInStudyAttr(aStudyBuilder, aSObject, "Boundaries", "BoundList", "zone_icone_2.png", NULL ) ;
+ }
+ else { MESSAGE("La categorie des boundarys existe deja."); }
+
+ CORBA::Long BoundaryType = myBoundary->GetType();
+// MESSAGE("BoundaryType : "<<BoundaryType);
+ std::string icone ;
+ std::string value ;
+ switch (BoundaryType)
+ {
+ case 0 :
+ { value = "BoundaryDiHomard" ;
+ icone = "mesh_tree_mesh.png" ;
+ break;
+ }
+ case 1 :
+ { value = "BoundaryAnHomard" ;
+ icone = "cylinderpointvector_2.png" ;
+ break;
+ }
+ case 2 :
+ { value = "BoundaryAnHomard" ;
+ icone = "spherepoint_2.png" ;
+ break;
+ }
+ case 3 :
+ { value = "BoundaryAnHomard" ;
+ icone = "conepointvector.png" ;
+ break;
+ }
+ case 4 :
+ { value = "BoundaryAnHomard" ;
+ icone = "conedxyz.png" ;
+ break;
+ }
+ case 5 :
+ { value = "BoundaryAnHomard" ;
+ icone = "toruspointvector.png" ;
+ break;
+ }
+ }
+ aResultSO = aStudyBuilder->NewObject(aSObject);
+ PublishInStudyAttr(aStudyBuilder, aResultSO, theName, value.c_str(), icone.c_str(), _orb->object_to_string(theObject));
+ return aResultSO._retn();
+}
+//=============================================================================
SALOMEDS::SObject_ptr HOMARD_Gen_i::PublishCaseInStudy(SALOMEDS::Study_ptr theStudy,
SALOMEDS::StudyBuilder_var aStudyBuilder,
HOMARD::HOMARD_Cas_ptr theObject, const char* theName)
{
MESSAGE("PublishCaseInStudy pour "<<theName);
SALOMEDS::SObject_var aResultSO;
- SALOMEDS::GenericAttribute_var anAttr;
if (CORBA::is_nil(theObject)) {
MESSAGE("HOMARD_Gen_i::theObject->_is_nil()");
return aResultSO._retn();
}
+ // On recupere le module pere dans l etude
SALOMEDS::SComponent_var theFatherHomard = theStudy->FindComponent(ComponentDataType());
if (theFatherHomard->_is_nil())
{
return aResultSO._retn();
}
+ _tag_gene += 1 ;
+ MESSAGE("PublishCaseInStudy _tag_gene = "<<_tag_gene );
aResultSO = aStudyBuilder->NewObject(theFatherHomard);
PublishInStudyAttr(aStudyBuilder, aResultSO, theName, "CasHomard", "cas_calcule.png",
_orb->object_to_string(theObject) ) ;
+
return aResultSO._retn();
}
-
//=============================================================================
-SALOMEDS::SObject_ptr HOMARD_Gen_i::PublishZoneInStudy(SALOMEDS::Study_ptr theStudy,
- SALOMEDS::StudyBuilder_var aStudyBuilder,
- HOMARD::HOMARD_Zone_ptr theObject, const char* theName)
+SALOMEDS::SObject_ptr HOMARD_Gen_i::PublishHypotheseInStudy(SALOMEDS::Study_ptr theStudy,
+ SALOMEDS::StudyBuilder_var aStudyBuilder,
+ HOMARD::HOMARD_Hypothesis_ptr theObject, const char* theName)
{
- MESSAGE("PublishZoneStudy pour "<<theName);
+ MESSAGE("PublishHypotheseInStudy pour "<<theName);
SALOMEDS::SObject_var aResultSO;
- SALOMEDS::GenericAttribute_var anAttr;
- if (CORBA::is_nil(theObject))
+ // On recupere le module pere dans l etude
+ // On ajoute la categorie des hypotheses dans l etude si necessaire
+ SALOMEDS::SComponent_var theFatherHomard = theStudy->FindComponent(ComponentDataType());
+ if (theFatherHomard->_is_nil())
{
- MESSAGE("PublishZoneInStudy : theObject->_is_nil()");
+ MESSAGE("theFatherHomard->_is_nil()");
return aResultSO._retn();
}
- if (theStudy->_is_nil())
- {
+
+ // On ajoute la categorie des hypotheses dans l etude si necessaire
+ SALOMEDS::SObject_var aSObject;
+ if ( _tag_hypo == 0 )
+ {
+ _tag_gene += 1 ;
+ _tag_hypo = _tag_gene ;
+ }
+ MESSAGE("PublishHypotheseInStudy _tag_gene = "<<_tag_gene << ", _tag_hypo = "<<_tag_hypo );
+ if (!theFatherHomard->FindSubObject(_tag_hypo, aSObject))
+ {
+ MESSAGE("Ajout de la categorie des hypotheses");
+ aSObject = aStudyBuilder->NewObjectToTag(theFatherHomard, _tag_hypo);
+ PublishInStudyAttr(aStudyBuilder, aSObject, "Hypothesis", "HypoList", "hypotheses.png", NULL);
+ }
+ else { MESSAGE("La categorie des hypotheses existe deja."); }
+
+// Creation du resultat dans l'etude
+ aResultSO = aStudyBuilder->NewObject(aSObject);
+ PublishInStudyAttr(aStudyBuilder, aResultSO, theName, "HypoHomard", NULL, _orb->object_to_string(theObject) ) ;
+
+ return aResultSO._retn();
+}
+//=============================================================================
+SALOMEDS::SObject_ptr HOMARD_Gen_i::PublishYACSInStudy(SALOMEDS::Study_ptr theStudy,
+ SALOMEDS::StudyBuilder_var aStudyBuilder,
+ HOMARD::HOMARD_YACS_ptr theObject, const char* theName)
+{
+ MESSAGE("PublishYACSInStudy pour "<<theName);
+ SALOMEDS::SObject_var aResultSO;
+
+ // On recupere le module pere dans l etude
+ // On ajoute la categorie des schemas YACS dans l etude si necessaire
+ SALOMEDS::SComponent_var theFatherHomard = theStudy->FindComponent(ComponentDataType());
+ if (theFatherHomard->_is_nil())
+ {
+ MESSAGE("theFatherHomard->_is_nil()");
+ return aResultSO._retn();
+ }
+ // On ajoute la categorie des schemas YACS dans l etude si necessaire
+ if ( _tag_yacs == 0 )
+ {
+ _tag_gene += 1 ;
+ _tag_yacs = _tag_gene ;
+ }
+ MESSAGE("PublishZoneStudy _tag_gene = "<<_tag_gene << ", _tag_yacs = "<<_tag_yacs );
+ SALOMEDS::SObject_var aSObject;
+ if (!theFatherHomard->FindSubObject(_tag_yacs, aSObject))
+ {
+ MESSAGE("Ajout de la categorie des schemas YACS");
+ aSObject = aStudyBuilder->NewObjectToTag(theFatherHomard, _tag_yacs);
+ PublishInStudyAttr(aStudyBuilder, aSObject, "YACS", "YACSList", "full_view.png", NULL);
+ }
+ else { MESSAGE("La categorie des schemas YACS existe deja."); }
+
+// Creation du resultat dans l'etude
+ aResultSO = aStudyBuilder->NewObject(aSObject);
+ PublishInStudyAttr(aStudyBuilder, aResultSO, theName, "YACSHomard", "schema.png", _orb->object_to_string(theObject) ) ;
+
+ return aResultSO._retn();
+}
+
+//=============================================================================
+SALOMEDS::SObject_ptr HOMARD_Gen_i::PublishZoneInStudy(SALOMEDS::Study_ptr theStudy,
+ SALOMEDS::StudyBuilder_var aStudyBuilder,
+ HOMARD::HOMARD_Zone_ptr theObject, const char* theName)
+{
+ MESSAGE("PublishZoneStudy pour "<<theName);
+ SALOMEDS::SObject_var aResultSO;
+ if (CORBA::is_nil(theObject))
+ {
+ MESSAGE("PublishZoneInStudy : theObject->_is_nil()");
+ return aResultSO._retn();
+ }
+ if (theStudy->_is_nil())
+ {
MESSAGE("PublishZoneInStudy : theStudy->_is_nil()");
return aResultSO._retn();
}
CORBA::Long ZoneType = myZone->GetType();
// On ajoute la categorie des zones dans l etude si necessaire
- SALOMEDS::SObject_var aSOZone;
- if (!theFatherHomard->FindSubObject(100, aSOZone))
+ if ( _tag_zone == 0 )
+ {
+ _tag_gene += 1 ;
+ _tag_zone = _tag_gene ;
+ }
+ MESSAGE("PublishZoneStudy _tag_gene = "<<_tag_gene << ", _tag_zone = "<<_tag_zone );
+ SALOMEDS::SObject_var aSObject;
+ if (!theFatherHomard->FindSubObject(_tag_zone, aSObject))
{
MESSAGE("Ajout de la categorie des zones");
- aSOZone = aStudyBuilder->NewObjectToTag(theFatherHomard, 100);
- PublishInStudyAttr(aStudyBuilder, aSOZone, "Zones", "ZoneList", "zone_icone_2.png", NULL ) ;
+ aSObject = aStudyBuilder->NewObjectToTag(theFatherHomard, _tag_zone);
+ PublishInStudyAttr(aStudyBuilder, aSObject, "Zones", "ZoneList", "zone_icone_2.png", NULL ) ;
}
else { MESSAGE("La categorie des zones existe deja."); }
- aResultSO = aStudyBuilder->NewObject(aSOZone);
- const char* icone ;
+ aResultSO = aStudyBuilder->NewObject(aSObject);
+ std::string icone ;
switch (ZoneType)
{
case 11 :
break ;
}
}
- PublishInStudyAttr(aStudyBuilder, aResultSO, theName, "ZoneHomard", icone, _orb->object_to_string(theObject) ) ;
-
- return aResultSO._retn();
-}
-//=============================================================================
-SALOMEDS::SObject_ptr HOMARD_Gen_i::PublishBoundaryInStudy(SALOMEDS::Study_ptr theStudy,
- SALOMEDS::StudyBuilder_var aStudyBuilder,
- HOMARD::HOMARD_Boundary_ptr theObject, const char* theName)
-{
- MESSAGE("PublishBoundaryStudy pour "<<theName);
- SALOMEDS::SObject_var aResultSO;
- SALOMEDS::GenericAttribute_var anAttr;
-
- // Caracteristique de la Boundary
- HOMARD::HOMARD_Boundary_var myBoundary = myContextMap[GetCurrentStudyID()]._mesBoundarys[theName];
-
- // On recupere le module pere dans l etude
- SALOMEDS::SComponent_var theFatherHomard = theStudy->FindComponent(ComponentDataType());
- if (theFatherHomard->_is_nil())
- {
- MESSAGE("theFatherHomard->_is_nil()");
- return aResultSO._retn();
- }
-
- // On ajoute la categorie des boundarys dans l etude si necessaire
- SALOMEDS::SObject_var aSOBoundary;
- if (!theFatherHomard->FindSubObject(101, aSOBoundary))
- {
- MESSAGE("Ajout de la categorie des boundarys");
- aSOBoundary = aStudyBuilder->NewObjectToTag(theFatherHomard, 101);
- PublishInStudyAttr(aStudyBuilder, aSOBoundary, "Boundaries", "BoundList", "zone_icone_2.png", NULL ) ;
- }
- else { MESSAGE("La categorie des boundarys existe deja."); }
-
- aResultSO = aStudyBuilder->NewObject(aSOBoundary);
- CORBA::Long BoundaryType = myBoundary->GetType();
-// MESSAGE("BoundaryType : "<<BoundaryType);
- const char* icone ;
- const char* value ;
- switch (BoundaryType)
- {
- case 0 :
- { value = "BoundaryDiHomard" ;
- icone = "mesh_tree_mesh.png" ;
- break;
- }
- case 1 :
- { value = "BoundaryAnHomard" ;
- icone = "cylinderpointvector_2.png" ;
- break;
- }
- case 2 :
- { value = "BoundaryAnHomard" ;
- icone = "spherepoint_2.png" ;
- break;
- }
- }
- PublishInStudyAttr(aStudyBuilder, aResultSO, theName, value, icone, _orb->object_to_string(theObject));
- return aResultSO._retn();
-}
-
-//=============================================================================
-SALOMEDS::SObject_ptr HOMARD_Gen_i::PublishHypotheseInStudy(SALOMEDS::Study_ptr theStudy,
- SALOMEDS::StudyBuilder_var aStudyBuilder,
- HOMARD::HOMARD_Hypothesis_ptr theObject, const char* theName)
-{
- MESSAGE("PublishHypotheseInStudy pour "<<theName);
- SALOMEDS::SObject_var aResultSO;
- SALOMEDS::GenericAttribute_var anAttr;
-
- // On recupere le module pere dans l etude
- // On ajoute la categorie des hypotheses dans l etude si necessaire
- SALOMEDS::SComponent_var theFatherHomard = theStudy->FindComponent(ComponentDataType());
- if (theFatherHomard->_is_nil())
- {
- MESSAGE("theFatherHomard->_is_nil()");
- return aResultSO._retn();
- }
- SALOMEDS::SObject_var aSOHypothese;
- if (!theFatherHomard->FindSubObject(0, aSOHypothese))
- {
- MESSAGE("Ajout de la categorie des hypotheses");
- aSOHypothese = aStudyBuilder->NewObjectToTag(theFatherHomard, 0);
- PublishInStudyAttr(aStudyBuilder, aSOHypothese, "Hypothesis", "HypoList","hypotheses.png", NULL);
- }
- else { MESSAGE("La categorie des hypotheses existe deja."); }
-
-// Creation du resultat dans l'etude
- aResultSO = aStudyBuilder->NewObject(aSOHypothese);
- PublishInStudyAttr(aStudyBuilder, aResultSO, theName, "HypoHomard", NULL, _orb->object_to_string(theObject) ) ;
+ PublishInStudyAttr(aStudyBuilder, aResultSO, theName, "ZoneHomard", icone.c_str(), _orb->object_to_string(theObject) ) ;
return aResultSO._retn();
}
//===========================================================================
void HOMARD_Gen_i::PublishInStudyAttr(SALOMEDS::StudyBuilder_var aStudyBuilder,
SALOMEDS::SObject_var aResultSO,
- const char* name, const char* value, const char* icone, const char* ior)
+ const char* name, const char* comment, const char* icone, const char* ior)
{
+ MESSAGE("PublishInStudyAttr pour name = "<<name<<", comment = "<<comment);
+// MESSAGE("icone = "<<icone);
+// MESSAGE("ior = "<<ior);
SALOMEDS::GenericAttribute_var anAttr ;
// Ajout du nom
if ( name != NULL )
}
// Ajout du commentaire
- if ( value != NULL )
+ if ( comment != NULL )
{
anAttr = aStudyBuilder->FindOrCreateAttribute(aResultSO, "AttributeComment");
SALOMEDS::AttributeComment_var aCommentAttrib = SALOMEDS::AttributeComment::_narrow(anAttr);
- aCommentAttrib->SetValue(value);
+ aCommentAttrib->SetValue(comment);
}
// Ajout de l'icone
if ( icone != NULL )
{
- anAttr = aStudyBuilder->FindOrCreateAttribute(aResultSO,"AttributePixMap");
+ anAttr = aStudyBuilder->FindOrCreateAttribute(aResultSO, "AttributePixMap");
SALOMEDS::AttributePixMap_var aPixmap = SALOMEDS::AttributePixMap::_narrow(anAttr);
aPixmap->SetPixMap(icone);
}
SALOMEDS::SObject_var aSubSO = aStudyBuilder->NewObject(aCaseSO);
aStudyBuilder->Addreference(aSubSO, aBoundarySO);
+// aStudyBuilder->RemoveReference(aSubSO);
+
+ aStudyBuilder->CommitCommand();
+
+};
+//=====================================================================================
+void HOMARD_Gen_i::PublishCaseUnderYACS(const char* nomYACS, const char* CaseName)
+{
+ MESSAGE ( "PublishCaseUnderYACS : nomYACS = " << nomYACS << ", CaseName= " << CaseName );
+
+ HOMARD::HOMARD_YACS_var myYACS = myContextMap[GetCurrentStudyID()]._mesYACSs[nomYACS];
+ ASSERT(!CORBA::is_nil(myYACS));
+ SALOMEDS::SObject_var aYACSSO = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myYACS)));
+ ASSERT(!CORBA::is_nil(aYACSSO));
+
+ HOMARD::HOMARD_Cas_var myCase = myContextMap[GetCurrentStudyID()]._mesCas[CaseName];
+ ASSERT(!CORBA::is_nil(myCase));
+ SALOMEDS::SObject_var aCaseSO = SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myCase)));
+ ASSERT(!CORBA::is_nil(aCaseSO));
+
+ SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
+
+ aStudyBuilder->NewCommand();
+
+ SALOMEDS::SObject_var aSubSO = aStudyBuilder->NewObject(aYACSSO);
+ aStudyBuilder->Addreference(aSubSO, aCaseSO);
aStudyBuilder->CommitCommand();
MESSAGE( "PublishResultInSmesh " << NomFich << ", avec Option = " << Option);
if (CORBA::is_nil(myCurrentStudy))
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Study Context ";
- throw SALOME::SALOME_Exception(es);
- return ;
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid study context";
+ throw SALOME::SALOME_Exception(es);
+ return ;
};
// Le module SMESH est-il actif ?
SALOMEDS::ChildIterator_var aIter = myCurrentStudy->NewChildIterator(aSmeshSO);
for (; aIter->More(); aIter->Next())
{
- SALOMEDS::SObject_var aSO = aIter->Value();
- SALOMEDS::GenericAttribute_var aGAttr;
- if (aSO->FindAttribute(aGAttr,"AttributeExternalFileDef"))
- {
- SALOMEDS::AttributeExternalFileDef_var anAttr = SALOMEDS::AttributeExternalFileDef::_narrow(aGAttr);
- CORBA::String_var value=anAttr->Value();
- if (strcmp((const char*)value,NomFich) == 0)
- {
- MESSAGE ( "PublishResultInSmesh : le fichier " << NomFich << " est deja publie." );
- // Pour un fichier importe, on ne republie pas
- if ( Option == 0 )
- {
- return;
- }
- // Pour un fichier calcule, on commence par faire la depublication
- else
- {
- MESSAGE ( "PublishResultInSmesh : depublication" );
- SALOMEDS::AttributeName_var anAttr2 = SALOMEDS::AttributeName::_narrow(aGAttr);
- CORBA::String_var value2=anAttr2->Value();
- const char* MeshName = value2 ;
- MESSAGE ( "PublishResultInSmesh : depublication de " << MeshName );
- DeleteResultInSmesh(NomFich, MeshName) ;
- }
- }
- }
- }
+ SALOMEDS::SObject_var aSO = aIter->Value();
+ SALOMEDS::GenericAttribute_var aGAttr;
+ if (aSO->FindAttribute(aGAttr,"AttributeExternalFileDef"))
+ {
+ SALOMEDS::AttributeExternalFileDef_var anAttr = SALOMEDS::AttributeExternalFileDef::_narrow(aGAttr);
+ CORBA::String_var value=anAttr->Value();
+ if (strcmp((const char*)value,NomFich) == 0)
+ {
+ MESSAGE ( "PublishResultInSmesh : le fichier " << NomFich << " est deja publie." );
+ // Pour un fichier importe, on ne republie pas
+ if ( Option == 0 ) { return; }
+ // Pour un fichier calcule, on commence par faire la depublication
+ else
+ {
+ MESSAGE ( "PublishResultInSmesh : depublication" );
+ SALOMEDS::AttributeName_var anAttr2 = SALOMEDS::AttributeName::_narrow(aGAttr);
+ CORBA::String_var value2=anAttr2->Value();
+ std::string MeshName = string(value2) ;
+ MESSAGE ( "PublishResultInSmesh : depublication de " << MeshName );
+ DeleteResultInSmesh(NomFich, MeshName) ;
+ }
+ }
+ }
+ }
}
// On enregistre le fichier
SMESH::mesh_array* mesMaillages=aSmeshEngine->CreateMeshesFromMED(NomFich, theStatus);
for (int i = 0; i < mesMaillages->length(); i++)
{
- MESSAGE( ". Mise a jour des attributs du maillage");
- SMESH::SMESH_Mesh_var monMaillage= (*mesMaillages)[i];
- SALOMEDS::SObject_var aSO=SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(monMaillage)));
- SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
- SALOMEDS::GenericAttribute_var aGAttr = aStudyBuilder->FindOrCreateAttribute(aSO, "AttributeExternalFileDef");
- SALOMEDS::AttributeExternalFileDef_var anAttr = SALOMEDS::AttributeExternalFileDef::_narrow(aGAttr);
- anAttr->SetValue(NomFich);
- SALOMEDS::GenericAttribute_var aPixMap = aStudyBuilder->FindOrCreateAttribute(aSO, "AttributePixMap" );
- SALOMEDS::AttributePixMap_var anAttr2 = SALOMEDS::AttributePixMap::_narrow(aPixMap);
- const char* icone ;
- if ( Option == 0 ) { icone = "mesh_tree_importedmesh.png" ; }
- else { icone = "mesh_tree_mesh.png" ; }
- anAttr2->SetPixMap( icone );
+ MESSAGE( ". Mise a jour des attributs du maillage");
+ SMESH::SMESH_Mesh_var monMaillage= (*mesMaillages)[i];
+ SALOMEDS::SObject_var aSO=SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(monMaillage)));
+ SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
+ SALOMEDS::GenericAttribute_var aGAttr = aStudyBuilder->FindOrCreateAttribute(aSO, "AttributeExternalFileDef");
+ SALOMEDS::AttributeExternalFileDef_var anAttr = SALOMEDS::AttributeExternalFileDef::_narrow(aGAttr);
+ anAttr->SetValue(NomFich);
+ SALOMEDS::GenericAttribute_var aPixMap = aStudyBuilder->FindOrCreateAttribute(aSO, "AttributePixMap" );
+ SALOMEDS::AttributePixMap_var anAttr2 = SALOMEDS::AttributePixMap::_narrow(aPixMap);
+ const char* icone ;
+ if ( Option == 0 ) { icone = "mesh_tree_importedmesh.png" ; }
+ else { icone = "mesh_tree_mesh.png" ; }
+ anAttr2->SetPixMap( icone );
+ }
+
+}
+//=============================================================================
+void HOMARD_Gen_i::DeleteResultInSmesh(std::string NomFich, std::string MeshName)
+{
+ MESSAGE ("DeleteResultInSmesh pour le maillage " << MeshName << " dans le fichier " << NomFich );
+ if (CORBA::is_nil(myCurrentStudy))
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid study context";
+ throw SALOME::SALOME_Exception(es);
+ return ;
+ };
+
+// Le module SMESH est-il actif ?
+ SALOMEDS::SObject_var aSmeshSO = myCurrentStudy->FindComponent("SMESH");
+//
+ if (CORBA::is_nil(aSmeshSO))
+ {
+ return ;
+ };
+// On verifie que le fichier est deja publie
+ SALOMEDS::StudyBuilder_var myBuilder = myCurrentStudy->NewBuilder();
+ SALOMEDS::ChildIterator_var aIter = myCurrentStudy->NewChildIterator(aSmeshSO);
+ for (; aIter->More(); aIter->Next())
+ {
+ SALOMEDS::SObject_var aSO = aIter->Value();
+ SALOMEDS::GenericAttribute_var aGAttr;
+ if (aSO->FindAttribute(aGAttr,"AttributeExternalFileDef"))
+ {
+ SALOMEDS::AttributeExternalFileDef_var anAttr = SALOMEDS::AttributeExternalFileDef::_narrow(aGAttr);
+ CORBA::String_var value=anAttr->Value();
+ if (strcmp((const char*)value,NomFich.c_str()) == 0)
+ {
+ if (aSO->FindAttribute(aGAttr,"AttributeName"))
+ {
+ SALOMEDS::AttributeName_var anAttr2 = SALOMEDS::AttributeName::_narrow(aGAttr);
+ CORBA::String_var value2=anAttr2->Value();
+ if (strcmp((const char*)value2,MeshName.c_str()) == 0)
+ {
+ myBuilder->RemoveObjectWithChildren( aSO ) ;
+ }
+ }
+ }
+ }
+ }
+ return ;
+}
+//=============================================================================
+void HOMARD_Gen_i::PublishMeshIterInSmesh(const char* NomIter)
+{
+ MESSAGE( "PublishMeshIterInSmesh " << NomIter);
+ HOMARD::HOMARD_Iteration_var myIteration = myContextMap[GetCurrentStudyID()]._mesIterations[NomIter];
+
+ SALOMEDS::SObject_var aIterSO=SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myIteration)));
+ if (CORBA::is_nil(myIteration))
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid iterationStudy Object";
+ throw SALOME::SALOME_Exception(es);
+ return ;
+ };
+ const char* MeshFile = myIteration->GetMeshFile() ;
+ const char* MeshName = myIteration->GetMeshName() ;
+ CORBA::Long Option = -1 ;
+ int etat = myIteration->GetState();
+// Iteration initiale
+ if ( etat <= 0 ) { Option = 0 ; }
+// ou iteration calculee
+ else if ( etat == 2 ) { Option = 1 ; }
+// Publication effective apres menage eventuel
+ if ( Option >= 0 )
+ {
+ DeleteResultInSmesh(MeshFile, MeshName) ;
+ PublishResultInSmesh(MeshFile, Option) ;
+ }
+
+}
+//=============================================================================
+void HOMARD_Gen_i::PublishFileUnderIteration(const char* NomIter, const char* NomFich, const char* Commentaire)
+{
+// MESSAGE ("PublishFileUnderIteration pour l'iteration " << NomIter << " du fichier " << NomFich << " avec le commentaire " << Commentaire );
+ HOMARD::HOMARD_Iteration_var myIteration = myContextMap[GetCurrentStudyID()]._mesIterations[NomIter];
+
+ SALOMEDS::SObject_var aIterSO=SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myIteration)));
+ if (CORBA::is_nil(myIteration))
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid iterationStudy Object";
+ throw SALOME::SALOME_Exception(es);
+ return ;
+ };
+
+ SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
+
+ aStudyBuilder->NewCommand();
+
+ SALOMEDS::SObject_var aSubSO = aStudyBuilder->NewObject(aIterSO);
+// Pour les fichiers med, on affiche une icone de maillage
+// Pour les fichiers qui sont du texte, on affiche une icone de fichier texte 'texte'
+// Le reperage se fait par la 1ere lettre du commentaire : I pour Iteration n
+ std::string icone ;
+ std::string ior = " " ;
+ if ( Commentaire[0] == 'I' )
+ { icone = "med.png" ; }
+ else
+ { icone = "texte_2.png" ; }
+ PublishInStudyAttr(aStudyBuilder, aSubSO, NomFich, Commentaire, icone.c_str(), ior.c_str() ) ;
+
+ aStudyBuilder->CommitCommand();
+}
+//
+//=============================================================================
+void HOMARD_Gen_i::PublishFileUnderYACS(const char* nomYACS, const char* NomFich, const char* Commentaire)
+{
+// MESSAGE ("PublishFileUnderYACS pour le schema " << nomYACS << " du fichier " << NomFich << " avec le commentaire " << Commentaire );
+ HOMARD::HOMARD_YACS_var myYACS = myContextMap[GetCurrentStudyID()]._mesYACSs[nomYACS];
+
+ SALOMEDS::SObject_var aYACSSO=SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myYACS)));
+ if (CORBA::is_nil(myYACS))
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Invalid YACSStudy Object";
+ throw SALOME::SALOME_Exception(es);
+ return ;
+ };
+
+ SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
+
+ aStudyBuilder->NewCommand();
+
+ SALOMEDS::SObject_var aSubSO = aStudyBuilder->NewObject(aYACSSO);
+ std::string icone = "texte_2.png" ;
+ std::string ior = " " ;
+ PublishInStudyAttr(aStudyBuilder, aSubSO, NomFich, Commentaire, icone.c_str(), ior.c_str() ) ;
+
+ aStudyBuilder->CommitCommand();
+}
+//
+//=============================================================================
+//=============================================================================
+// YACS
+//=============================================================================
+//=============================================================================
+//=============================================================================
+// Creation d'un schema YACS
+// nomCas : nom du cas a traiter
+// FileName : nom du fichier contenant le script de lancement du calcul
+// DirName : le repertoire de lancement des calculs du schéma
+//=============================================================================
+HOMARD::HOMARD_YACS_ptr HOMARD_Gen_i::CreateYACSSchema (const char* nomYACS, const char* nomCas, const char* ScriptFile, const char* DirName, const char* MeshFile)
+{
+ INFOS ( "CreateYACSSchema : Schema YACS " << nomYACS );
+ INFOS ( ". nomCas : " << nomCas);
+ INFOS ( ". ScriptFile : " << ScriptFile);
+ INFOS ( ". DirName : " << DirName);
+ INFOS ( ". MeshFile : " << MeshFile);
+
+ // A. Controle du nom :
+ if ((myContextMap[GetCurrentStudyID()]._mesYACSs).find(nomYACS) != (myContextMap[GetCurrentStudyID()]._mesYACSs).end())
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "This schema YACS has already been defined.";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ }
+
+ // B. Creation de l'objet
+ HOMARD::HOMARD_YACS_var myYACS = newYACS();
+ if (CORBA::is_nil(myYACS))
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ es.text = "Unable to create the schema YACS";
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ };
+ myYACS->SetName( nomYACS ) ;
+
+ // C. Enregistrement
+ myContextMap[GetCurrentStudyID()]._mesYACSs[nomYACS] = myYACS;
+
+ SALOMEDS::SObject_var aSO;
+ SALOMEDS::SObject_var aResultSO=PublishInStudy(myCurrentStudy, aSO, myYACS, nomYACS);
+
+ PublishCaseUnderYACS(nomYACS, nomCas);
+
+ // D. Caracterisation
+ // D.1. Options
+ myYACS->SetDirName( DirName ) ;
+ myYACS->SetMeshFile( MeshFile ) ;
+ myYACS->SetScriptFile( ScriptFile ) ;
+ myYACS->SetCaseName( nomCas ) ;
+ // D.2. Defaut
+ int defaut_i ;
+ // D.2.1. Type constant
+ myYACS->SetType( 1 ) ;
+ // D.2.2. Convergence
+ defaut_i = GetYACSMaxIter() ;
+ myYACS->SetMaxIter( defaut_i ) ;
+ defaut_i = GetYACSMaxNode() ;
+ myYACS->SetMaxNode( defaut_i ) ;
+ defaut_i = GetYACSMaxElem() ;
+ myYACS->SetMaxElem( defaut_i ) ;
+ // D.3. Fichier de sauvegarde dans le repertoire du cas
+ HOMARD::HOMARD_Cas_ptr caseyacs = GetCase(nomCas) ;
+ std::string dirnamecase = caseyacs->GetDirName() ;
+ std::string XMLFile ;
+ XMLFile = dirnamecase + "/schema.xml" ;
+ myYACS->SetXMLFile( XMLFile.c_str() ) ;
+
+ return HOMARD::HOMARD_YACS::_duplicate(myYACS);
+}
+//=============================================================================
+// Ecriture d'un schema YACS
+//=============================================================================
+CORBA::Long HOMARD_Gen_i::YACSWrite(const char* nomYACS)
+{
+ INFOS ( "YACSWrite : Ecriture de " << nomYACS );
+// Le repertoire du cas
+ HOMARD::HOMARD_YACS_var myYACS = myContextMap[GetCurrentStudyID()]._mesYACSs[nomYACS];
+ ASSERT(!CORBA::is_nil(myYACS));
+// Le nom du fichier du schema
+ std::string XMLFile ;
+ XMLFile = myYACS->GetXMLFile() ;
+
+ int codret = YACSWriteOnFile(nomYACS, XMLFile.c_str()) ;
+
+ return codret ;
+}
+//=============================================================================
+// Ecriture d'un schema YACS sur un fichier donne
+//=============================================================================
+CORBA::Long HOMARD_Gen_i::YACSWriteOnFile(const char* nomYACS, const char* XMLFile)
+{
+ INFOS ( "YACSWriteOnFile : Ecriture de " << nomYACS << " sur " << XMLFile );
+
+ // A. Prealable
+ int codret = 0;
+
+ // B. L'objet YACS
+ // B.1. L'objet
+ HOMARD::HOMARD_YACS_var myYACS = myContextMap[GetCurrentStudyID()]._mesYACSs[nomYACS];
+ ASSERT(!CORBA::is_nil(myYACS));
+ // B.2. Les caracteristiques
+ std::string DirName = myYACS->GetDirName() ;
+ std::string MeshFile = myYACS->GetMeshFile() ;
+ std::string ScriptFile = myYACS->GetScriptFile() ;
+ // B.3. Les caracteristiques de convergence
+ int MaxIter = myYACS->GetMaxIter() ;
+ int MaxNode = myYACS->GetMaxNode() ;
+ int MaxElem = myYACS->GetMaxElem() ;
+
+ // C. Le cas
+ // C.1. L'objet cas
+ const char* nomCas = myYACS->GetCaseName();
+ HOMARD::HOMARD_Cas_var myCase = myContextMap[GetCurrentStudyID()]._mesCas[nomCas];
+ ASSERT(!CORBA::is_nil(myCase));
+ // C.2. Les instructions python associees au cas
+ CORBA::String_var dumpCorbaCase = myCase->GetDumpPython();
+ std::string pythonCas = dumpCorbaCase.in();
+ MESSAGE ("pythonCas :\n"<<pythonCas<<"\n");
+
+ // D. Les iterations
+ // D.1. L'iteration initiale
+ HOMARD::HOMARD_Iteration_var Iter0 = myCase->GetIter0() ;
+ std::string Iter0Name = myCase->GetIter0Name() ;
+ MESSAGE (". Iter0Name = " << Iter0Name);
+ std::string MeshName = Iter0->GetMeshName();
+ MESSAGE (". MeshName = " << MeshName);
+ // D.2. L'iteration numero 1
+ HOMARD::listeIterFilles* maListe = Iter0->GetIterations();
+ int numberOfIter = maListe->length();
+ if ( numberOfIter > 1 )
+ {
+ MESSAGE (". numberOfIter = " << numberOfIter);
+ SALOME::ExceptionStruct es ;
+ es.type = SALOME::BAD_PARAM;
+ std::string text = "Une seule iteration est permise." ;
+ es.text = CORBA::string_dup(text.c_str());
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ }
+ std::string Iter1Name = std::string((*maListe)[0]);
+ MESSAGE ("... Iter1Name = " << Iter1Name);
+ HOMARD::HOMARD_Iteration_var Iter1 = GetIteration(Iter1Name.c_str()) ;
+ // D.3. Les instructions python associees a l'iteration
+ CORBA::String_var dumpCorbaIter = Iter1->GetDumpPython();
+ std::string pythonIter = dumpCorbaIter.in();
+ MESSAGE ("pythonIter :\n"<<pythonIter<<"\n");
+
+ // E. L'hypothese pour passer de l'iteration initiale a la suivante
+ // E.1. La structure
+ std::string nomHypo = Iter1->GetHypoName();
+ MESSAGE (". nomHypo = " << nomHypo);
+ HOMARD::HOMARD_Hypothesis_var myHypo = myContextMap[GetCurrentStudyID()]._mesHypotheses[nomHypo];
+ ASSERT(!CORBA::is_nil(myHypo));
+ // E.2. Les caracteristiques de l'adaptation
+ HOMARD::listeTypes* ListTypes = myHypo->GetAdapRefinUnRef();
+ ASSERT(ListTypes->length() == 3);
+ int TypeAdap = (*ListTypes)[0];
+// int TypeRaff = (*ListTypes)[1];
+// int TypeDera = (*ListTypes)[2];
+ // E.3. Les instructions python associees a l'hypothese
+ CORBA::String_var dumpCorbaHypo = myHypo->GetDumpPython();
+ std::string pythonHypo = dumpCorbaHypo.in();
+ MESSAGE ("pythonHypo :\n"<<pythonHypo<<"\n");
+
+ // F. Le fichier du schema de reference
+ // HOMARD_ROOT_DIR : repertoire ou se trouve le module HOMARD
+ std::string XMLFile_base ;
+ if ( getenv("HOMARD_ROOT_DIR") != NULL ) { XMLFile_base = getenv("HOMARD_ROOT_DIR") ; }
+ else
+ {
+ SALOME::ExceptionStruct es ;
+ es.type = SALOME::BAD_PARAM;
+ std::string text = "HOMARD_ROOT_DIR est inconnu." ;
+ es.text = CORBA::string_dup(text.c_str());
+ throw SALOME::SALOME_Exception(es);
+ return 0;
+ }
+ XMLFile_base += "/share/salome/resources/homard/yacs_01." + _LangueShort + ".xml" ;
+// if ( _Langue ==
+ MESSAGE("XMLFile_base ="<<XMLFile_base);
+
+ // G. Lecture du schema de reference et insertion des donnees propres au fil de la rencontre des mots-cles
+ YACSDriver* myDriver = new YACSDriver(XMLFile, DirName);
+ std::ifstream fichier( XMLFile_base.c_str() );
+ if ( fichier ) // ce test échoue si le fichier n'est pas ouvert
+ {
+ // G.1. Lecture du schema de reference et insertion des donnees propres au fil de la rencontre des mots-cles
+ std::string ligne; // variable contenant chaque ligne lue
+ std::string mot_cle;
+ while ( std::getline( fichier, ligne ) )
+ {
+ // G.1.1. Pour la ligne courante, on identifie le premier mot : le mot-cle eventuel
+ std::istringstream ligne_bis(ligne); // variable contenant chaque ligne sous forme de flux
+ ligne_bis >> mot_cle ;
+ // G.1.2. Le maillage initial
+ if ( mot_cle == "DataInit_MeshFile" )
+ { myDriver->Texte_DataInit_MeshFile(MeshFile); }
+ // G.1.3. Le script de lancement
+ else if ( mot_cle == "Alternance_Calcul_HOMARD_Calcul" )
+ { myDriver->Texte_Alternance_Calcul_HOMARD_Calcul(ScriptFile); }
+ // G.1.4. Les options du cas
+ else if ( mot_cle == "Iter_1_Case_Options" )
+ { myDriver->Texte_Iter_1_Case_Options(pythonCas); }
+ // G.1.5. Execution de HOMARD : le repertoire du cas
+ else if ( mot_cle == "HOMARD_Exec_DirName" )
+ { myDriver->Texte_HOMARD_Exec_DirName(); }
+ // G.1.6. Execution de HOMARD : le nom du maillage
+ else if ( mot_cle == "HOMARD_Exec_MeshName" )
+ {
+ myDriver->Texte_HOMARD_Exec_MeshName(MeshName);
+ std::string node = "Boucle_de_convergence.Alternance_Calcul_HOMARD.Adaptation.p0_Adaptation_HOMARD.HOMARD_Initialisation.p1_Iter_1.CreateCase" ;
+ myDriver->TexteParametre( node, "MeshName", "string", MeshName ) ;
+ }
+ // G.1.7. Execution de HOMARD : les options de l'hypothese
+ else if ( mot_cle == "HOMARD_Exec_Hypo_Options" )
+ { myDriver->Texte_python_1( pythonHypo, 3, "Hypo" ) ; }
+ // G.1.8. Execution de HOMARD : les options de l'iteration
+ else if ( mot_cle == "HOMARD_Exec_Iter_Options" )
+ { myDriver->Texte_python_2( pythonIter, "TimeStep", "Iter" ) ; }
+ // G.1.9. a. Creation eventuelles des zones et frontieres
+ // b. Enchainement
+ else if ( mot_cle == "Iter_1" )
+ {
+ std::string texte_control = "" ;
+ if ( TypeAdap == 0 ) { texte_control += YACSDriverTexteZone( myHypo, myDriver ) ; }
+ texte_control += YACSDriverTexteBoundary( myCase, myDriver ) ;
+ texte_control += myDriver->Texte_Iter_1_control() ;
+ myDriver->TexteAdd(texte_control);
+ }
+ // G.1.10. Les tests de convergence
+ else if ( mot_cle == "Analyse_Test_Convergence" )
+ { myDriver->TexteAnalyse_Test_Convergence(MaxIter, MaxNode, MaxElem); }
+ // G.1.11. Les parametres
+ else if ( mot_cle == "PARAMETRES" )
+ { myDriver->TexteAddParametres(); }
+ // G.1.n. La ligne est recopiee telle quelle
+ else { myDriver->TexteAdd(ligne); }
+ }
+ // G.2. Ecriture du texte dans le fichier
+ if ( codret == 0 )
+ { myDriver->CreeFichier(); }
+ }
+ else
+ {
+ SALOME::ExceptionStruct es;
+ es.type = SALOME::BAD_PARAM;
+ std::string text = "The reference file for the YACS schema cannot be read." ;
+ es.text = CORBA::string_dup(text.c_str());
+ throw SALOME::SALOME_Exception(es);
}
+ delete myDriver;
+
+ // H. Publication du fichier dans l'arbre
+
+ std::string Commentaire = "xml" ;
+ PublishFileUnderYACS(nomYACS, XMLFile, Commentaire.c_str());
+
+ return codret ;
}
//=============================================================================
-void HOMARD_Gen_i::DeleteResultInSmesh(const char* NomFich, const char* MeshName)
+// Ecriture d'un schema YACS : ecriture des zones associees a une hypothese
+//=============================================================================
+std::string HOMARD_Gen_i::YACSDriverTexteZone(HOMARD::HOMARD_Hypothesis_var myHypo, YACSDriver* myDriver)
{
- MESSAGE (" DeleteResultInSmesh pour le maillage " << MeshName << " dans le fichier " << NomFich );
- if (CORBA::is_nil(myCurrentStudy))
- {
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Study Context ";
- throw SALOME::SALOME_Exception(es);
- return ;
- };
+ MESSAGE ( "YACSDriverTexteZone" );
+ // A. Les zones associees a cette hypothese
+ HOMARD::listeZonesHypo* ListZone = myHypo->GetZones();
+ int numberOfZonesx2 = ListZone->length();
-// Le module SMESH est-il actif ?
- SALOMEDS::SObject_var aSmeshSO = myCurrentStudy->FindComponent("SMESH");
-//
- if (CORBA::is_nil(aSmeshSO))
- {
- return ;
- };
-// On verifie que le fichier est deja publie
- SALOMEDS::StudyBuilder_var myBuilder = myCurrentStudy->NewBuilder();
- SALOMEDS::ChildIterator_var aIter = myCurrentStudy->NewChildIterator(aSmeshSO);
- for (; aIter->More(); aIter->Next())
+ // B. Parcours des zones
+ std::string texte_control ;
+ for (int iaux = 0; iaux< numberOfZonesx2; iaux++)
{
- SALOMEDS::SObject_var aSO = aIter->Value();
- SALOMEDS::GenericAttribute_var aGAttr;
- if (aSO->FindAttribute(aGAttr,"AttributeExternalFileDef"))
- {
- SALOMEDS::AttributeExternalFileDef_var anAttr = SALOMEDS::AttributeExternalFileDef::_narrow(aGAttr);
- CORBA::String_var value=anAttr->Value();
- if (strcmp((const char*)value,NomFich) == 0)
- {
- if (aSO->FindAttribute(aGAttr,"AttributeName"))
- {
- SALOMEDS::AttributeName_var anAttr2 = SALOMEDS::AttributeName::_narrow(aGAttr);
- CORBA::String_var value2=anAttr2->Value();
- if (strcmp((const char*)value2,MeshName) == 0)
- {
- myBuilder->RemoveObjectWithChildren( aSO ) ;
- }
- }
- }
- }
+ // 1. Reperage de la zone
+ std::string ZoneName = std::string((*ListZone)[iaux]);
+ MESSAGE ( "\n. ZoneName = " << ZoneName << " - " <<iaux);
+ HOMARD::HOMARD_Zone_var myZone = myContextMap[GetCurrentStudyID()]._mesZones[ZoneName];
+ ASSERT(!CORBA::is_nil(myZone));
+ // 2. Les instructions python associees a la zone
+ // La premiere ligne est un commentaire a eliminer
+ // La seconde ligne est l'instruction a proprement parler ; on ne garde que ce qui suit le "."
+ CORBA::String_var dumpCorbaZone = myZone->GetDumpPython();
+ std::string pythonStructure_0 = dumpCorbaZone.in();
+ MESSAGE ("pythonStructure_0 :"<<pythonStructure_0);
+ std::istringstream tout (pythonStructure_0) ;
+ std::string ligne ;
+ std::string pythonStructure ;
+ while ( std::getline( tout, ligne ) )
+ { pythonStructure = GetStringInTexte ( ligne, ".", 1 ) ; }
+ MESSAGE ("pythonStructure :\n"<<pythonStructure);
+ // 3. Decodage du nom du service
+ std::string methode = GetStringInTexte ( pythonStructure, "(", 0 ) ;
+ MESSAGE ( "... methode = " << methode);
+ // 4. Mise en place des instructions
+ int ZoneType = myZone->GetType();
+ MESSAGE ( "... ZoneType = " << ZoneType);
+ std::string texte_control_0 ;
+ texte_control_0 = myDriver->Texte_Iter_1_Zone(ZoneType, pythonStructure, methode, ZoneName );
+ texte_control += texte_control_0 ;
+ // 5. Decalage
+ iaux ++ ;
}
+ return texte_control ;
}
//=============================================================================
-void HOMARD_Gen_i::PublishFileUnderIteration(const char* NomIter, const char* NomFich, const char* Commentaire)
+// Ecriture d'un schema YACS : ecriture des frontieres associees au cas
+//=============================================================================
+std::string HOMARD_Gen_i::YACSDriverTexteBoundary(HOMARD::HOMARD_Cas_var myCase, YACSDriver* myDriver)
{
-// MESSAGE (" PublishFileUnderIteration pour l'iteration " << NomIter << " du fichier " << NomFich << " avec le commentaire " << Commentaire );
- if (CORBA::is_nil(myCurrentStudy))
- {
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Study Context ";
- throw SALOME::SALOME_Exception(es);
- return ;
- };
+ MESSAGE ( "YACSDriverTexteBoundary" );
+ // A. Les frontieres associees au cas
+ HOMARD::ListBoundaryGroupType* ListBoundaryGroupType = myCase->GetBoundaryGroup();
+ int numberOfitems = ListBoundaryGroupType->length();
- HOMARD::HOMARD_Iteration_var myIteration = myContextMap[GetCurrentStudyID()]._mesIterations[NomIter];
- if (CORBA::is_nil(myIteration))
- {
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Iteration ";
- throw SALOME::SALOME_Exception(es);
- return ;
- };
- SALOMEDS::SObject_var aIterSO=SALOMEDS::SObject::_narrow(myCurrentStudy->FindObjectIOR(_orb->object_to_string(myIteration)));
- if (CORBA::is_nil(myIteration))
+ // B. Parcours des frontieres
+ std::string texte_control ;
+ std::list<std::string> ListeBoundaryTraitees ;
+ for (int NumBoundary = 0; NumBoundary< numberOfitems; NumBoundary=NumBoundary+2)
{
- SALOME::ExceptionStruct es;
- es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Iteration Study Object";
- throw SALOME::SALOME_Exception(es);
- return ;
- };
-
- SALOMEDS::StudyBuilder_var aStudyBuilder = myCurrentStudy->NewBuilder();
-
- aStudyBuilder->NewCommand();
-
- SALOMEDS::SObject_var aSubSO = aStudyBuilder->NewObject(aIterSO);
-// Pour les fichiers med, on affiche une icone de maillage
-// Pour les fichiers qui sont texte, on affiche une icone de fichier texte 'texte'
-// Le reperage se fait par la 1ere lettre du commentaire : I pour Iteration n
- const char* icone ;
- const char* ior = " " ;
- if ( Commentaire[0] == 'I' )
- { icone = "med.png" ; }
- else
- { icone = "texte_2.png" ; }
- PublishInStudyAttr(aStudyBuilder, aSubSO, NomFich, Commentaire, icone, ior ) ;
+ std::string BoundaryName = std::string((*ListBoundaryGroupType)[NumBoundary]);
+ MESSAGE ( "... BoundaryName = " << BoundaryName);
+ // Attention a n'ecrire la definition qu'une seule fois car elle peut se trouver
+ // plusieurs fois dans la definition du cas, d'ou la liste ListeBoundaryTraitees
+ int A_faire = 1 ;
+ std::list<std::string>::const_iterator it = ListeBoundaryTraitees.begin();
+ while (it != ListeBoundaryTraitees.end())
+ {
+ MESSAGE ( "..... BoundaryNameTraitee = " << *it);
+ if ( BoundaryName == *it ) { A_faire = 0 ; }
+ it++;
+ }
+ if ( A_faire == 1 )
+ {
+ // 1. Caracteristiques de la frontiere
+ HOMARD::HOMARD_Boundary_var myBoundary = myContextMap[GetCurrentStudyID()]._mesBoundarys[BoundaryName];
+ ASSERT(!CORBA::is_nil(myBoundary));
+ // 2. Les instructions python associees a la frontiere
+ // La premiere ligne est un commentaire a eliminer
+ // La seconde ligne est l'instruction a proprement parler ; on ne garde que ce qui suit le "."
+ CORBA::String_var dumpCorbaBoundary = myBoundary->GetDumpPython();
+ std::string pythonStructure_0 = dumpCorbaBoundary.in();
+ MESSAGE ("pythonStructure_0 :"<<pythonStructure_0);
+ std::istringstream tout (pythonStructure_0) ;
+ std::string ligne ;
+ std::string pythonStructure ;
+ while ( std::getline( tout, ligne ) )
+ { pythonStructure = GetStringInTexte ( ligne, ".", 1 ) ; }
+ MESSAGE ("pythonStructure :\n"<<pythonStructure);
+ // 3. Decodage du nom du service
+ std::string methode = GetStringInTexte ( pythonStructure, "(", 0 ) ;
+ MESSAGE ( "... methode = " << methode);
+ // 4. Mise en place des instructions
+ int BoundaryType = myBoundary->GetType();
+ MESSAGE ( "... BoundaryType = " << BoundaryType);
+ const char* MeshName ;
+ const char* MeshFile ;
+ if (BoundaryType == 0)
+ {
+ MeshName = myBoundary->GetMeshName() ;
+ MESSAGE ( ". MeshName = " << MeshName );
+ MeshFile = myBoundary->GetMeshFile() ;
+ MESSAGE ( ". MeshFile = " << MeshFile );
+ }
+ std::string texte_control_0 ;
+ texte_control_0 = myDriver->Texte_Iter_1_Boundary(BoundaryType, pythonStructure, methode, BoundaryName, MeshName, MeshFile );
+ texte_control += texte_control_0 ;
+ // 5. Memorisation du traitement
+ ListeBoundaryTraitees.push_back( BoundaryName );
+ }
+ }
- aStudyBuilder->CommitCommand();
+ return texte_control ;
}
-//=============================================================================
-//=============================================================================
//
//=============================================================================
//=============================================================================
const char* theURL,
CORBA::Boolean isMultiFile)
{
- MESSAGE (" Save for theURL = "<< theURL);
+ MESSAGE ("Save for theURL = "<< theURL);
SALOMEDS::TMPFile_var aStreamFile;
// get temporary directory name
const char* theURL,
CORBA::Boolean isMultiFile)
{
- MESSAGE (" Load pour theURL = "<< theURL);
+ MESSAGE ("Load pour theURL = "<< theURL);
SALOMEDS::Study_var aStudy = theComponent->GetStudy();
// set current study
while (f) {
std::getline(f, line);
+ std::string bounSignature = HOMARD::GetSignature(HOMARD::Boundary);
std::string caseSignature = HOMARD::GetSignature(HOMARD::Case);
- std::string zoneSignature = HOMARD::GetSignature(HOMARD::Zone);
- std::string iterSignature = HOMARD::GetSignature(HOMARD::Iteration);
std::string hypoSignature = HOMARD::GetSignature(HOMARD::Hypothesis);
- std::string bounSignature = HOMARD::GetSignature(HOMARD::Boundary);
- if (line.substr(0, caseSignature.size()) == caseSignature) {
+ std::string iterSignature = HOMARD::GetSignature(HOMARD::Iteration);
+ std::string zoneSignature = HOMARD::GetSignature(HOMARD::Zone);
+ std::string yacsSignature = HOMARD::GetSignature(HOMARD::YACS);
+ if (line.substr(0, bounSignature.size()) == bounSignature) {
+ // re-create boundary
+ MESSAGE ("Recreation de la frontiere" );
+ HOMARD::HOMARD_Boundary_var aBoundary = newBoundary();
+ PortableServer::ServantBase_var aServant = GetServant(aBoundary);
+ HOMARD_Boundary_i* aBoundaryServant = dynamic_cast<HOMARD_Boundary_i*>(aServant.in());
+ if (aBoundaryServant && aBoundaryServant->Restore(line.substr(bounSignature.size()))) {
+ context._mesBoundarys[aBoundary->GetName()] = aBoundary;
+ context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aBoundaryServant);
+ }
+ }
+ else if (line.substr(0, caseSignature.size()) == caseSignature) {
// re-create case
- MESSAGE (" Recreation du cas" );
+ MESSAGE ("Recreation du cas" );
HOMARD::HOMARD_Cas_var aCase = newCase();
PortableServer::ServantBase_var aServant = GetServant(aCase);
HOMARD_Cas_i* aCaseServant = dynamic_cast<HOMARD_Cas_i*>(aServant.in());
context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aCaseServant);
}
}
- else if (line.substr(0, zoneSignature.size()) == zoneSignature) {
- MESSAGE (" Recreation de la zone" );
- // re-create zone
- HOMARD::HOMARD_Zone_var aZone = newZone();
- PortableServer::ServantBase_var aServant = GetServant(aZone);
- HOMARD_Zone_i* aZoneServant = dynamic_cast<HOMARD_Zone_i*>(aServant.in());
- if (aZoneServant && aZoneServant->Restore(line.substr(zoneSignature.size()))) {
- context._mesZones[aZone->GetName()] = aZone;
- context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aZoneServant);
+ else if (line.substr(0, hypoSignature.size()) == hypoSignature) {
+ // re-create hypothesis
+ MESSAGE ("Recreation de l hypothese" );
+ HOMARD::HOMARD_Hypothesis_var aHypo = newHypothesis();
+ PortableServer::ServantBase_var aServant = GetServant(aHypo);
+ HOMARD_Hypothesis_i* aHypoServant = dynamic_cast<HOMARD_Hypothesis_i*>(aServant.in());
+ if (aHypoServant && aHypoServant->Restore(line.substr(hypoSignature.size()))) {
+ context._mesHypotheses[aHypo->GetName()] = aHypo;
+ context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aHypoServant);
}
}
else if (line.substr(0, iterSignature.size()) == iterSignature) {
// re-create iteration
- MESSAGE (" Recreation de l iteration" );
+ MESSAGE ("Recreation de l'iteration" );
HOMARD::HOMARD_Iteration_var aIter = newIteration();
PortableServer::ServantBase_var aServant = GetServant(aIter);
HOMARD_Iteration_i* aIterServant = dynamic_cast<HOMARD_Iteration_i*>(aServant.in());
context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aIterServant);
}
}
- else if (line.substr(0, hypoSignature.size()) == hypoSignature) {
- // re-create hypothesis
- MESSAGE (" Recreation de l hypothese" );
- HOMARD::HOMARD_Hypothesis_var aHypo = newHypothesis();
- PortableServer::ServantBase_var aServant = GetServant(aHypo);
- HOMARD_Hypothesis_i* aHypoServant = dynamic_cast<HOMARD_Hypothesis_i*>(aServant.in());
- if (aHypoServant && aHypoServant->Restore(line.substr(hypoSignature.size()))) {
- context._mesHypotheses[aHypo->GetName()] = aHypo;
- context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aHypoServant);
+ else if (line.substr(0, zoneSignature.size()) == zoneSignature) {
+ MESSAGE ("Recreation de la zone" );
+ // re-create zone
+ HOMARD::HOMARD_Zone_var aZone = newZone();
+ PortableServer::ServantBase_var aServant = GetServant(aZone);
+ HOMARD_Zone_i* aZoneServant = dynamic_cast<HOMARD_Zone_i*>(aServant.in());
+ if (aZoneServant && aZoneServant->Restore(line.substr(zoneSignature.size()))) {
+ context._mesZones[aZone->GetName()] = aZone;
+ context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aZoneServant);
}
}
- else if (line.substr(0, bounSignature.size()) == bounSignature) {
- // re-create boundary
- MESSAGE (" Recreation de la frontiere" );
- HOMARD::HOMARD_Boundary_var aBoundary = newBoundary();
- PortableServer::ServantBase_var aServant = GetServant(aBoundary);
- HOMARD_Boundary_i* aBoundaryServant = dynamic_cast<HOMARD_Boundary_i*>(aServant.in());
- if (aBoundaryServant && aBoundaryServant->Restore(line.substr(bounSignature.size()))) {
- context._mesBoundarys[aBoundary->GetName()] = aBoundary;
- context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aBoundaryServant);
+ else if (line.substr(0, zoneSignature.size()) == yacsSignature) {
+ MESSAGE ("Recreation du schema YACS" );
+ // re-create YACS
+ HOMARD::HOMARD_YACS_var aYACS = newYACS();
+ PortableServer::ServantBase_var aServant = GetServant(aYACS);
+ HOMARD_YACS_i* aYACSServant = dynamic_cast<HOMARD_YACS_i*>(aServant.in());
+ if (aYACSServant && aYACSServant->Restore(line.substr(yacsSignature.size()))) {
+ context._mesYACSs[aYACS->GetName()] = aYACS;
+ context._idmap[id] = dynamic_cast<PortableServer::ServantBase*>(aYACSServant);
}
}
id++;
std::string aScript = "\"\"\"\n";
aScript += "Python script for HOMARD\n";
- aScript += "Copyright EDF-R&D 2013\n";
+ aScript += "Copyright EDF-R&D 1996, 2011, 2014\n";
aScript += "\"\"\"\n";
aScript += "__revision__ = \"V1.2\"\n";
aScript += "import HOMARD\n";
for ( it_zone = myContextMap[GetCurrentStudyID()]._mesZones.begin();
it_zone != myContextMap[GetCurrentStudyID()]._mesZones.end(); ++it_zone)
{
- HOMARD::HOMARD_Zone_var maZone = (*it_zone).second;
- CORBA::String_var dumpCorbaZone = maZone->GetDumpPython();
+ HOMARD::HOMARD_Zone_var myZone = (*it_zone).second;
+ CORBA::String_var dumpCorbaZone = myZone->GetDumpPython();
std::string dumpZone = dumpCorbaZone.in();
MESSAGE (dumpZone<<"\n");
aScript += dumpZone;
std::string dumpCas = std::string("\n# Creation of the case ") ;
dumpCas += nomCas + std::string("\n");
dumpCas += std::string("\t") + nomCas;
- dumpCas += std::string(" = homard.CreateCase('") + nomCas + std::string("', '");
+ dumpCas += std::string(" = homard.CreateCase(\"") + nomCas + std::string("\", \"");
HOMARD::HOMARD_Cas_var myCase = (*it_cas).second;
CORBA::String_var cIter0= myCase->GetIter0Name();
std::string mesh0 = cMesh0.in();
CORBA::String_var cMeshName0= myIteration->GetMeshName();
std::string meshName0 = cMeshName0.in();
- dumpCas += meshName0 + std::string("', '")+ mesh0 + std::string("')\n");
+ dumpCas += meshName0 + std::string("\", \"")+ mesh0 + std::string("\")\n");
CORBA::String_var dumpCorbaCase = myCase->GetDumpPython();
std::string dumpCas2= dumpCorbaCase.in();
aScript += dumpIter;
}
+
+ if (myContextMap[GetCurrentStudyID()]._mesYACSs.size() > 0)
+ {
+ MESSAGE (". Ecritures des schemas YACS");
+ aScript += "#\n# Creation of the schemas YACS";
+ aScript += "\n# ============================";
+ }
+ std::map<std::string, HOMARD::HOMARD_YACS_var>::const_iterator it_yacs;
+ for ( it_yacs = myContextMap[GetCurrentStudyID()]._mesYACSs.begin();
+ it_yacs != myContextMap[GetCurrentStudyID()]._mesYACSs.end(); ++it_yacs)
+ {
+ HOMARD::HOMARD_YACS_var myYACS = (*it_yacs).second;
+ CORBA::String_var dumpCorbaYACS = myYACS->GetDumpPython();
+ std::string dumpYACS = dumpCorbaYACS.in();
+ MESSAGE (dumpYACS<<"\n");
+ aScript += dumpYACS;
+ }
+
MESSAGE (". Ecritures finales");
if( isMultiFile )
aScript += "\n\tpass";
{
SALOME::ExceptionStruct es;
es.type = SALOME::BAD_PARAM;
- es.text = "Invalid Study Context";
+ es.text = "Invalid study context";
throw SALOME::SALOME_Exception(es);
};
return ;
}
//=============================================================================
-CORBA::Boolean HOMARD_Gen_i::VerifieDir(const char* nomDir)
+char* HOMARD_Gen_i::VerifieDir(const char* nomDir)
{
+ std::string casename = std::string("") ;
std::map<std::string, HOMARD::HOMARD_Cas_var>::const_iterator it;
for (it = myContextMap[GetCurrentStudyID()]._mesCas.begin();
it != myContextMap[GetCurrentStudyID()]._mesCas.end(); it++)
{
- if (std::string(nomDir) == std::string(it->second->GetDirName())) return false;
+ if (std::string(nomDir) == std::string(it->second->GetDirName()))
+ {
+ casename = std::string(it->second->GetName()) ;
+ break ;
+ }
}
- return true;
+ return CORBA::string_dup( casename.c_str() );
}
/*//=============================================================================
void SALOMEException( std::string message )
return CORBA::string_dup(HOMARD_VERSION_STR);
#endif
}
+//===============================================================================
+// Recuperation de la chaine de caracteres par rapport a l'apparition d'un caractere
+// ligne : la ligne a manipuler
+// caractere : le caractere a reperer
+// option : 0 : la chaine avant la premiere apparition du caractere
+// 1 : la chaine apres la premiere apparition du caractere
+// 2 : la chaine avant la derniere apparition du caractere
+// 3 : la chaine apres la derniere apparition du caractere
+// Si le caractere est absent, on retourne la chaine totale
+//===============================================================================
+std::string HOMARD_Gen_i::GetStringInTexte( const std::string ligne, const std::string caractere, int option )
+{
+// MESSAGE("GetStringInTexte, recherche de '"<<caractere<<"' dans '"<<ligne<<"'"<<", option = "<<option);
+//
+ std::string chaine = ligne ;
+ int position ;
+ if ( option < 2 ) { position = ligne.find_first_of( caractere ) ; }
+ else { position = ligne.find_last_of( caractere ) ; }
+// MESSAGE("position = "<<position);
+// MESSAGE("a = "<<ligne.substr( 0, position ).c_str());
+// MESSAGE("b = "<<ligne.substr( position+1 ).c_str());
+//
+ if ( position != std::string::npos )
+ {
+ if ( ( option == 0 ) || ( option == 2 ) ) { chaine = ligne.substr( 0, position ) ; }
+ else { chaine = ligne.substr( position+1 ) ; }
+ }
+ return chaine ;
+//
+}
+//=============================================================================
+//=============================================================================
+// Gestion des preferences
+//=============================================================================
+//=============================================================================
+// Decodage du fichier d'arcihvage des preferences
+//
+void HOMARD_Gen_i::SetPreferences( )
+{
+ MESSAGE ( "SetPreferences" );
+
+ std::string ligne, mot_cle, salome_version ;
+ bool ok = true ;
+
+ // A. Les valeurs par defaut ; elles doivent etre coherentes
+ std::string LanguageShort = "en" ;
+ int PublisMeshIN = 0 ;
+ int PublisMeshOUT = 0 ;
+ int YACSMaxIter = 0 ;
+ int YACSMaxNode = 0 ;
+ int YACSMaxElem = 0 ;
+ std::string YACSTypeTestchaine = "None" ;
+
+ // B. La version de salome
+ // Cela se presente sous la forme :
+ // [SALOME KERNEL] : 7.3.0
+ std::string File ;
+ File = getenv("KERNEL_ROOT_DIR") ;
+ File += "/bin/salome/VERSION" ;
+ MESSAGE ( "File = "<<File ) ;
+ std::ifstream fichier0( File.c_str() ) ;
+ if ( fichier0 ) // ce test échoue si le fichier n'est pas ouvert
+ {
+ std::string ligne; // variable contenant chaque ligne lue
+ while ( std::getline( fichier0, ligne ) )
+ {
+ std::istringstream ligne_bis(ligne); // variable contenant chaque ligne sous forme de flux
+ ligne_bis >> mot_cle ;
+ if ( mot_cle == "[SALOME" )
+ {
+ salome_version = GetStringInTexte ( ligne, " ", 3 ) ;
+// MESSAGE ( "salome_version = "<<salome_version<<"|||");
+ break ;
+ }
+ }
+ }
+ else { ok = false ; }
+
+ // B. Decodage du fichier de preferences
+ if ( ok )
+ {
+ std::string PrefFile ;
+ PrefFile = getenv("HOME") ;
+ PrefFile += "/.config/salome/SalomeApprc." + salome_version ;
+ MESSAGE ( "PrefFile = "<<PrefFile ) ;
+
+ std::ifstream fichier( PrefFile.c_str() );
+ if ( fichier ) // ce test échoue si le fichier n'est pas ouvert
+ {
+ bool section_langue = false ;
+ bool section_homard = false ;
+ while ( std::getline( fichier, ligne ) )
+ {
+ std::string chaine ;
+ // 1. Pour la ligne courante, on identifie le premier mot : le mot-cle eventuel
+ std::istringstream ligne_bis(ligne); // variable contenant chaque ligne sous forme de flux
+ ligne_bis >> mot_cle ;
+
+ // 2. Les sections
+ // 2.1. Debut d'une section
+ // MESSAGE(mot_cle);
+ if ( mot_cle == "<section" )
+ { /*MESSAGE ( "Debut de la section : "<< ligne);*/
+ ligne_bis >> mot_cle ;
+ chaine = GetStringInTexte ( mot_cle, "\"", 1 ) ;
+ chaine = GetStringInTexte ( chaine, "\"", 0 ) ;
+ if ( chaine == "language" ) { section_langue = true ; }
+ if ( chaine == "HOMARD" ) { section_homard = true ; }
+ // MESSAGE ( "section_langue = "<<section_langue<<", section_homard = "<<section_homard);
+ }
+ // 2.2. Fin d'une section
+ else if ( mot_cle == "</section>" )
+ { /*MESSAGE ( "Fin de la section : "<< ligne<<", section_langue = "<<section_langue<<", section_homard = "<<section_homard);*/
+ section_langue = false ;
+ section_homard = false ; }
+
+ // 3. Parametres
+ // 3.1. La langue
+ else if ( section_langue || section_homard )
+ { MESSAGE ( "a decoder : "<< ligne);
+// La valeur : entre les deux premieres quotes
+ chaine = GetStringInTexte ( ligne, "\"", 1 ) ;
+// MESSAGE("chaine 1 = |"<<chaine<<"|");
+ chaine = GetStringInTexte ( chaine, "\"", 0 ) ;
+// MESSAGE("chaine = |"<<chaine<<"|");
+// Le mot_cle : entre les deux dernieres quotes
+ std::string chaine2 = GetStringInTexte ( ligne, "\"", 2 ) ;
+// MESSAGE("chaine2 1 = |"<<chaine2<<"|");
+ chaine2 = GetStringInTexte ( chaine2, "\"", 3 ) ;
+// MESSAGE("chaine2 = |"<<chaine2<<"|");
+ // 3.1. La langue
+ if ( section_langue )
+ { if ( chaine2 == "language" ) { LanguageShort = chaine ; } }
+ // 3.2. HOMARD
+ if ( section_homard )
+ {
+ std::istringstream chainebis( chaine ) ;
+ // 3.2.1. Les publications
+ if ( chaine2 == "publish_mesh_in" ) { chainebis >> PublisMeshIN ; }
+ if ( chaine2 == "publish_mesh_out" ) { chainebis >> PublisMeshOUT ; }
+ // 3.2.2. Les maximum pour YACS
+ if ( chaine2 == "yacs_max_iter" ) { chainebis >> YACSMaxIter ; }
+ if ( chaine2 == "yacs_max_node" ) { chainebis >> YACSMaxNode ; }
+ if ( chaine2 == "yacs_max_elem" ) { chainebis >> YACSMaxElem ; }
+ if ( chaine2 == "yacs_type_test" ) { YACSTypeTestchaine = chaine ; }
+ }
+ }
+ }
+ }
+ }
+
+ // C. Enregistrements
+ MESSAGE ("Enregistrement de LanguageShort = " << LanguageShort );
+ SetLanguageShort( LanguageShort.c_str() ) ;
+
+ MESSAGE ("Enregistrement de PublisMeshIN = " << PublisMeshIN<<", PublisMeshOUT = "<< PublisMeshOUT);
+ SetPublisMesh(PublisMeshIN, PublisMeshOUT) ;
+
+ MESSAGE ("Enregistrement de YACSMaxIter = " << YACSMaxIter<<", YACSMaxNode = "<< YACSMaxNode<<", YACSMaxElem = "<< YACSMaxElem);
+ SetYACSMaximum(YACSMaxIter, YACSMaxNode, YACSMaxElem) ;
+
+ MESSAGE ("Enregistrement de TypeTest = " << YACSTypeTestchaine.c_str() );
+ int YACSTypeTest ;
+ if ( ( YACSTypeTestchaine == "VTest > VRef" ) || ( YACSTypeTestchaine == "VTest > VRef" ) ) { YACSTypeTest = 1 ; }
+ else if ( ( YACSTypeTestchaine == "VTest < VRef" ) || ( YACSTypeTestchaine == "VTest < VRef" ) ) { YACSTypeTest = 2 ; }
+ else { YACSTypeTest = 0 ; }
+ MESSAGE ("==> TypeTest = " << YACSTypeTest );
+ SetYACSConvergenceType( YACSTypeTest ) ;
+
+ return ;
+}
+//===============================================================================
+// Langue de SALOME
+//===============================================================================
+void HOMARD_Gen_i::SetLanguageShort(const char* LanguageShort)
+{
+// MESSAGE ("SetLanguageShort pour LanguageShort = " << LanguageShort );
+ _LangueShort = LanguageShort ;
+ if ( _LangueShort == "fr" ) { _Langue = "Francais" ; }
+ else { _Langue = "English" ; }
+ return ;
+}
+char* HOMARD_Gen_i::GetLanguageShort()
+{
+// MESSAGE ("GetLanguageShort");
+ return CORBA::string_dup( _LangueShort.c_str() );
+}
+//===============================================================================
+// Options de publications
+//===============================================================================
+void HOMARD_Gen_i::SetPublisMesh(CORBA::Long PublisMeshIN, CORBA::Long PublisMeshOUT)
+{
+ _PublisMeshIN = PublisMeshIN ;
+ _PublisMeshOUT = PublisMeshOUT ;
+ return ;
+}
+CORBA::Long HOMARD_Gen_i::GetPublisMeshIN()
+{
+ return _PublisMeshIN ;
+}
+CORBA::Long HOMARD_Gen_i::GetPublisMeshOUT()
+{
+ return _PublisMeshOUT ;
+}
+//===============================================================================
+// YACS - test de convergence
+//===============================================================================
+void HOMARD_Gen_i::SetYACSMaximum(CORBA::Long YACSMaxIter, CORBA::Long YACSMaxNode, CORBA::Long YACSMaxElem)
+{
+ _YACSMaxIter = YACSMaxIter ;
+ _YACSMaxNode = YACSMaxNode ;
+ _YACSMaxElem = YACSMaxElem ;
+ return ;
+}
+CORBA::Long HOMARD_Gen_i::GetYACSMaxIter()
+{
+ return _YACSMaxIter ;
+}
+CORBA::Long HOMARD_Gen_i::GetYACSMaxNode()
+{
+ return _YACSMaxNode ;
+}
+CORBA::Long HOMARD_Gen_i::GetYACSMaxElem()
+{
+ return _YACSMaxElem ;
+}
+void HOMARD_Gen_i::SetYACSConvergenceType(CORBA::Long YACSTypeTest)
+{
+ _YACSTypeTest = YACSTypeTest ;
+ return ;
+}
+CORBA::Long HOMARD_Gen_i::GetYACSConvergenceType()
+{
+ return _YACSTypeTest ;
+}
//=============================================================================
extern "C"
{
+ HOMARDENGINE_EXPORT
PortableServer::ObjectId* HOMARDEngine_factory(CORBA::ORB_ptr orb,
PortableServer::POA_ptr poa,
PortableServer::ObjectId* contId,
