+// HOMARD HOMARD : implementaion of HOMARD idl descriptions
+//
+// Copyright (C) 2011 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.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
+
#include <string.h>
#include <unistd.h>
#include <sys/stat.h>
std::string dir;
std::string executable;
if ( !(dirchar = getenv("HOMARD_REP_EXE")) )
- dir = "/logiciels/HOMARD/HOMARD_PUBLIC/Linux64" ;
+ dir = "/local00/HOMARD_SVN/trunk/bin" ;
else
dir = std::string(dirchar);
if ( !(execchar = getenv("HOMARD_EXE")) )
- executable = "HOMARD";
+ executable = "HOMARD.out";
else
executable = std::string(execchar);
_HOMARD_Exec = dir + "/" + executable ;
- MESSAGE("Dans HomardDriver::HomardDriver, _HOMARD_Exec ="<<_HOMARD_Exec);
+ MESSAGE("Dans HomardDriver, _HOMARD_Exec ="<<_HOMARD_Exec);
//
_siter = siter ;
_siterp1 = siterp1 ;
////=============================================================================
void HomardDriver::TexteInit( const std::string DirCompute, const std::string DirComputePa, const std::string MessFile )
{
- MESSAGE("Dans HomardDriver::TexteInit, MessFile ="<<MessFile);
+ MESSAGE("Dans TexteInit, MessFile ="<<MessFile);
//
_Texte = "Action homa\n" ;
_Texte += "CCAssoci med\n" ;
////=============================================================================
void HomardDriver::TexteMaillage( const std::string NomMesh, const std::string MeshFile, int apres )
{
- MESSAGE("Dans HomardDriver::TexteMaillage, NomMesh ="<<NomMesh<<", MeshFile ="<<MeshFile);
+ MESSAGE("Dans TexteMaillage, NomMesh ="<<NomMesh<<", MeshFile ="<<MeshFile);
std::string saux ;
saux = "P1" ;
if ( apres < 1 ) { saux = "__" ; }
////=============================================================================
void HomardDriver::TexteConfRaffDera( int ConfType, int TypeAdap, int TypeRaff, int TypeDera )
{
- MESSAGE("Dans HomardDriver::TexteConfRaffDera, ConfType ="<<ConfType);
- MESSAGE("Dans HomardDriver::TexteConfRaffDera, TypeAdap ="<<TypeAdap<<", TypeRaff ="<<TypeRaff<<", TypeDera ="<<TypeDera);
+ MESSAGE("Dans TexteConfRaffDera, ConfType ="<<ConfType);
+ MESSAGE("Dans TexteConfRaffDera, TypeAdap ="<<TypeAdap<<", TypeRaff ="<<TypeRaff<<", TypeDera ="<<TypeDera);
//
// Type de conformite
//
////=============================================================================
void HomardDriver::TexteCompo( int NumeComp, const std::string NompCompo)
{
- MESSAGE("Dans HomardDriver::TexteCompo, NumeComp = "<<NumeComp<<", NompCompo = "<<NompCompo);
+ MESSAGE("Dans TexteCompo, NumeComp = "<<NumeComp<<", NompCompo = "<<NompCompo);
_Texte +="CCCoChaI \"" + NompCompo + "\"\n" ;
}
////=============================================================================
-void HomardDriver::TexteZone( int NumeZone, int ZoneType, double x0, double x1, double x2, double x3, double x4, double x5 )
+void HomardDriver::TexteZone( int NumeZone, int ZoneType, double x0, double x1, double x2, double x3, double x4, double x5, double x6, double x7, double x8 )
{
- MESSAGE("Dans HomardDriver::TexteZone, NumeZone = "<<NumeZone<<", ZoneType = "<<ZoneType);
- MESSAGE("Dans HomardDriver::TexteZone, coor = "<< x0<<","<<x1<< ","<< x2<< ","<< x3<<","<<x4<<","<<x5);
+ MESSAGE("TexteZone, NumeZone = "<<NumeZone<<", ZoneType = "<<ZoneType);
+ MESSAGE("TexteZone, coor = "<< x0<<","<<x1<< ","<< x2<< ","<< x3<<","<<x4<<","<<x5<<","<<x6<<","<<x7<<","<<x8);
//
std::string saux, saux2 ;
//
// Type de zones
+//
+ int ZoneTypeHOMARD ;
+ if ( ZoneType >= 11 and ZoneType <= 13 ) { ZoneTypeHOMARD = 1 ; }
+ else if ( ZoneType >= 31 and ZoneType <= 33 ) { ZoneTypeHOMARD = 3 ; }
+ else if ( ZoneType >= 61 and ZoneType <= 63 ) { ZoneTypeHOMARD = 6 ; }
+ else { ZoneTypeHOMARD = ZoneType ; }
//
std::stringstream saux1 ;
saux1 << NumeZone ;
saux = "#\n# Zone numero " + saux2 + "\n" ;
//
{ std::stringstream saux1 ;
- saux1 << NumeZone << " " << ZoneType ;
+ saux1 << NumeZone << " " << ZoneTypeHOMARD ;
saux2 = saux1.str() ;
saux += "ZoRaType " + saux2 + "\n" ;
}
//
-// Cas du rectangle/parallelepipede
+// Cas du rectangle
+//
+ if ( ZoneType == 11 ) // Z est constant X Homard <=> X Salome
+// Y Homard <=> Y Salome
+ {
+ saux += "#Rectangle\n" ;
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x0 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXmin " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x1 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXmax " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x2 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYmin " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x3 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYmax " + saux2 + "\n" ;
+ }
+ }
+//
+ else if ( ZoneType == 12 ) // X est constant X Homard <=> Y Salome
+// Y Homard <=> Z Salome
+ {
+ saux += "#Rectangle\n" ;
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x2 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXmin " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x3 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXmax " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x4 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYmin " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x5 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYmax " + saux2 + "\n" ;
+ }
+ }
+//
+ else if ( ZoneType == 13 ) // Y est constant X Homard <=> X Salome
+// Y Homard <=> Z Salome
+ {
+ saux += "#Rectangle\n" ;
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x0 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXmin " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x1 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXmax " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x4 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYmin " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x5 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYmax " + saux2 + "\n" ;
+ }
+ }
+//
+// Cas du parallelepipede
//
- if ( ZoneType == 2 )
+ else if ( ZoneType == 2 )
{
saux += "# Boite\n" ;
{ std::stringstream saux1 ;
}
}
//
-// Cas du cercle/sphere
+// Cas du disque
+//
+ else if ( ZoneType == 31 or ZoneType == 61 )
+ {
+ saux += "# Sphere\n" ;
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x0 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXCen " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x1 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYCen " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x6 ;
+ saux2 = saux1.str() ;
+ if ( ZoneType == 61 ) { saux += "ZoRaRayE " + saux2 + "\n" ; }
+ else { saux += "ZoRaRayo " + saux2 + "\n" ; }
+ }
+ if ( ZoneType == 61 )
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x8 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaRayI " + saux2 + "\n" ;
+ }
+ }
+ else if ( ZoneType == 32 or ZoneType == 62 )
+ {
+ saux += "# Sphere\n" ;
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x1 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXCen " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x2 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYCen " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x6 ;
+ saux2 = saux1.str() ;
+ if ( ZoneType == 62 ) { saux += "ZoRaRayE " + saux2 + "\n" ; }
+ else { saux += "ZoRaRayo " + saux2 + "\n" ; }
+ }
+ if ( ZoneType == 62 )
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x8 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaRayI " + saux2 + "\n" ;
+ }
+ }
+ else if ( ZoneType == 33 or ZoneType == 63 )
+ {
+ saux += "# Sphere\n" ;
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x0 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXCen " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x2 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYCen " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x6 ;
+ saux2 = saux1.str() ;
+ if ( ZoneType == 63 ) { saux += "ZoRaRayE " + saux2 + "\n" ; }
+ else { saux += "ZoRaRayo " + saux2 + "\n" ; }
+ }
+ if ( ZoneType == 63 )
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x8 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaRayI " + saux2 + "\n" ;
+ }
+ }
+//
+// Cas de la sphere
//
else if ( ZoneType == 4 )
{
saux += "ZoRaRayo " + saux2 + "\n" ;
}
}
+//
+// Cas du cylindre ou du tuyau
+//
+ else if ( ZoneType == 5 or ZoneType == 7 )
+ {
+ if ( ZoneType == 5 ) { saux += "# Cylindre\n" ; }
+ else { saux += "# Tuyau\n" ; }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x0 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXBas " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x1 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYBas " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x2 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaZBas " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x3 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaXAxe " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x4 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaYAxe " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x5 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaZAxe " + saux2 + "\n" ;
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x6 ;
+ saux2 = saux1.str() ;
+ if ( ZoneType == 5 ) { saux += "ZoRaRayo " + saux2 + "\n" ; }
+ else { saux += "ZoRaRayE " + saux2 + "\n" ; }
+ }
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x7 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaHaut " + saux2 + "\n" ;
+ }
+ if ( ZoneType == 7 )
+ { std::stringstream saux1 ;
+ saux1 << NumeZone << " " << x8 ;
+ saux2 = saux1.str() ;
+ saux += "ZoRaRayI " + saux2 + "\n" ;
+ }
+ }
//
_Texte += saux + "#\n" ;
//
////=============================================================================
void HomardDriver::TexteField( const std::string FieldName, const std::string FieldFile, int TimeStep, int Rank,
- int TypeThR, double ThreshR, int TypeThC, double ThreshC, int UsCmpI )
+ int TypeThR, double ThreshR, int TypeThC, double ThreshC,
+ int UsField, int UsCmpI )
{
- MESSAGE("Dans HomardDriver::TexteField, FieldName = "<<FieldName<<", FieldFile = "<<FieldFile);
- MESSAGE("Dans HomardDriver::TexteField, TimeStep = "<<TimeStep<<", Rank = "<<Rank);
+ MESSAGE("Dans TexteField, FieldName = "<<FieldName<<", FieldFile = "<<FieldFile);
+ MESSAGE("Dans TexteField, TimeStep = "<<TimeStep<<", Rank = "<<Rank);
std::string saux, saux2 ;
//
saux2 = saux1.str() ;
_Texte += "Seuil" + saux + " " + saux2 + "\n" ;
}
+//
+ saux = " " ;
+ if ( UsField == 0 )
+ { saux = "MAILLE" ; }
+ if ( UsField == 1 )
+ { saux = "SAUT" ; }
+ if ( saux != " " )
+ {
+ _Texte += "CCModeFI " + saux + "\n" ;
+ }
//
saux = " " ;
if ( UsCmpI == 0 )
////=============================================================================
void HomardDriver::TexteGroup( const std::string GroupName )
{
- MESSAGE("Dans HomardDriver::TexteGroup, GroupName = "<<GroupName);
+ MESSAGE("Dans TexteGroup, GroupName = "<<GroupName);
//
_Texte += "CCGroAda " + GroupName + "\n" ;
//
////=============================================================================
void HomardDriver::TexteBoundaryOption( int BoundaryOption )
{
- MESSAGE("Dans HomardDriver::TexteBoundaryOption, BoundaryOption = "<<BoundaryOption);
+ MESSAGE("Dans TexteBoundaryOption, BoundaryOption = "<<BoundaryOption);
//
// Type de suivi de frontiere
//
//
}
////=============================================================================
-void HomardDriver::TexteBoundaryDi( const std::string MeshName, const std::string MeshFile, const std::string GroupName )
+void HomardDriver::TexteBoundaryDi( const std::string MeshName, const std::string MeshFile )
{
- MESSAGE("Dans HomardDriver::TexteBoundaryDi, MeshName = "<<MeshName);
- MESSAGE("Dans HomardDriver::TexteBoundaryDi, MeshFile = "<<MeshFile);
- MESSAGE("Dans HomardDriver::TexteBoundaryDi, GroupName = "<<GroupName);
+ MESSAGE("Dans TexteBoundaryDi, MeshName = "<<MeshName);
+ MESSAGE("Dans TexteBoundaryDi, MeshFile = "<<MeshFile);
//
_Texte += "CCNoMFro " + MeshName + "\n" ;
_Texte += "CCFronti " + MeshFile + "\n" ;
- if ( GroupName.size() > 0 ) _Texte += "CCGroFro " + GroupName + "\n" ;
//
}
////=============================================================================
-void HomardDriver::TexteBoundaryAn( int NumeBoundary, int BoundaryType, const std::string Group, double x0, double x1, double x2, double x3, double x4, double x5, double x6 )
+void HomardDriver::TexteBoundaryDiGr( const std::string GroupName )
{
- MESSAGE("Dans HomardDriver::TexteBoundaryAn, NumeBoundary = "<<NumeBoundary);
- MESSAGE("Dans HomardDriver::TexteBoundaryAn, BoundaryType = "<<BoundaryType);
- MESSAGE("Dans HomardDriver::TexteBoundaryAn, coor = "<< x0<<","<<x1<< ","<< x2<< ","<< x3<<","<<x4<<","<<x5<<","<<x6);
+ MESSAGE("Dans TexteBoundaryDiGr, GroupName = "<<GroupName);
+//
+ _Texte += "CCGroFro " + GroupName + "\n" ;
+//
+}
+////=============================================================================
+void HomardDriver::TexteBoundaryAn( const std::string NameBoundary, int NumeBoundary, int BoundaryType, double x0, double x1, double x2, double x3, double x4, double x5, double x6 )
+{
+ MESSAGE("Dans TexteBoundaryAn, NameBoundary = "<<NameBoundary);
+// MESSAGE("Dans TexteBoundaryAn, NumeBoundary = "<<NumeBoundary);
+ MESSAGE("Dans TexteBoundaryAn, BoundaryType = "<<BoundaryType);
+// MESSAGE("Dans TexteBoundaryAn, coor = "<< x0<<","<<x1<< ","<< x2<< ","<< x3<<","<<x4<<","<<x5<<","<<x6);
//
std::string saux, saux2 ;
//
-// Type de zones
+// Commentaires
//
std::stringstream saux1 ;
saux1 << NumeBoundary ;
{ saux += "# Cylindre\n" ; }
if ( BoundaryType == 2 )
{ saux += "# Sphere\n" ; }
+//
+// Le nom de la frontiere
//
{ std::stringstream saux1 ;
- saux1 << NumeBoundary << " " << BoundaryType ;
+ saux1 << NumeBoundary ;
saux2 = saux1.str() ;
- saux += "FAType " + saux2 + "\n" ;
+ saux += "FANom " + saux2 + " \"" + NameBoundary + "\"\n" ;
}
//
-// Le groupe
+// Type de frontiere
//
{ std::stringstream saux1 ;
- saux1 << NumeBoundary ;
+ saux1 << NumeBoundary << " " << BoundaryType ;
saux2 = saux1.str() ;
- saux += "FAGroupe " + saux2 + " '" + Group + "'\n" ;
+ saux += "FAType " + saux2 + "\n" ;
}
//
// Cas du cylindre
//
}
////=============================================================================
+void HomardDriver::TexteBoundaryAnGr( const std::string NameBoundary, int NumeBoundary, const std::string GroupName )
+{
+ MESSAGE("Dans TexteBoundaryAnGr, NameBoundary = "<<NameBoundary);
+// MESSAGE("Dans TexteBoundaryAnGr, NumeBoundary = "<<NumeBoundary);
+// MESSAGE("Dans TexteBoundaryAnGr, GroupName = "<<GroupName);
+//
+// Commentaires
+//
+ std::string saux, saux2 ;
+ std::stringstream saux1 ;
+ saux1 << NumeBoundary ;
+ saux2 = saux1.str() ;
+ saux = "#\n# Lien Frontiere/Groupe numero " + saux2 + "\n" ;
+//
+ saux += "FGNomFro " + saux2 + " \"" + NameBoundary + "\"\n" ;
+ saux += "FGNomGro " + saux2 + " \"" + GroupName + "\"\n" ;
+//
+ _Texte += saux + "#\n" ;
+//
+}
+////=============================================================================
void HomardDriver::TexteFieldInterp( int TypeFieldInterp, const std::string FieldFile, const std::string MeshFile, int TimeStep, int Rank )
{
- MESSAGE("Dans HomardDriver::TexteFieldInterp, TypeFieldInterp = "<<TypeFieldInterp);
- MESSAGE("Dans HomardDriver::TexteFieldInterp, FieldFile = "<<FieldFile<<", MeshFile = "<<MeshFile);
- MESSAGE("Dans HomardDriver::TexteFieldInterp, TimeStep = "<<TimeStep<<", Rank = "<<Rank);
+ MESSAGE("Dans TexteFieldInterp, TypeFieldInterp = "<<TypeFieldInterp);
+ MESSAGE("Dans TexteFieldInterp, FieldFile = "<<FieldFile<<", MeshFile = "<<MeshFile);
+ MESSAGE("Dans TexteFieldInterp, TimeStep = "<<TimeStep<<", Rank = "<<Rank);
//
// Type d'interpolation
//
////=============================================================================
void HomardDriver::TexteFieldInterpName( int NumeChamp, const std::string FieldName)
{
- MESSAGE("Dans HomardDriver::TexteFieldInterpName, NumeChamp = "<<NumeChamp<<", FieldName = "<<FieldName);
+ MESSAGE("Dans TexteFieldInterpName, NumeChamp = "<<NumeChamp<<", FieldName = "<<FieldName);
std::stringstream saux1 ;
saux1 << NumeChamp+1 ;
std::string saux = saux1.str() ;
}
}
}
+////=============================================================================
+void HomardDriver::TexteAdvanced( int NivMax, double DiamMin)
+{
+ MESSAGE("Dans TexteAdvanced, NivMax ="<<NivMax<<", DiamMin ="<<DiamMin);
+ std::string saux ;
+ _Texte += "# Niveaux extremes\n" ;
+ { std::stringstream saux1 ;
+ saux1 << NivMax ;
+ std::string saux2 = saux1.str() ;
+ _Texte += "NiveauMa " + saux2 + "\n" ;
+ }
+ {
+ std::stringstream saux1 ;
+ saux1 << DiamMin ;
+ std::string saux2 = saux1.str() ;
+ _Texte += "DiametMi " + saux2 + "\n" ;
+ }
+}
////=============================================================================
void HomardDriver::CreeFichier( )
