1 # -*- coding: utf-8 -*-
2 # Copyright (C) 2011-2016 CEA/DEN, EDF R&D
4 # This library is free software; you can redistribute it and/or
5 # modify it under the terms of the GNU Lesser General Public
6 # License as published by the Free Software Foundation; either
7 # version 2.1 of the License, or (at your option) any later version.
9 # This library is distributed in the hope that it will be useful,
10 # but WITHOUT ANY WARRANTY; without even the implied warranty of
11 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 # Lesser General Public License for more details.
14 # You should have received a copy of the GNU Lesser General Public
15 # License along with this library; if not, write to the Free Software
16 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 Python script for HOMARD
22 Test tutorial_6 associe au tutorial 6
24 __revision__ = "V4.04"
26 #========================================================================
27 TEST_NAME = "tutorial_6"
30 #========================================================================
36 # ==================================
37 PATH_HOMARD = os.getenv('HOMARD_ROOT_DIR')
38 # Repertoire des scripts utilitaires
39 REP_PYTHON = os.path.join(PATH_HOMARD, "bin", "salome", "test", "HOMARD")
40 REP_PYTHON = os.path.normpath(REP_PYTHON)
41 sys.path.append(REP_PYTHON)
42 from test_util import get_dir
43 from test_util import get_dir_tutorial
44 from test_util import test_results
45 # ==================================
46 # Répertoires pour ce test
47 REP_DATA, DIRCASE = get_dir(PATH_HOMARD, TEST_NAME, DEBUG)
48 DATA_TUTORIAL = get_dir_tutorial(PATH_HOMARD)
49 # ==================================
50 sys.path.append(DATA_TUTORIAL)
51 from tutorial_util import gzip_gunzip
52 # ==================================
53 gzip_gunzip(DATA_TUTORIAL, 4, -1)
54 gzip_gunzip(DATA_TUTORIAL, 6, -1)
55 # ==================================
59 IPAR = iparameters.IParameters(salome.myStudy.GetCommonParameters("Interface Applicative", 1))
60 IPAR.append("AP_MODULES_LIST", "Homard")
63 #========================= Debut de la fonction ==================================
65 def homard_exec(nom, ficmed, nomfr, ficfrmed, verbose=False):
67 Python script for HOMARD
74 HOMARD.SetCurrentStudy(salome.myStudy)
80 boun_6_1 = HOMARD.CreateBoundaryDi('intersection', nomfr, ficfrmed)
82 boun_6_2 = HOMARD.CreateBoundaryCylinder('cyl_1_ext', 0.0, 25., -25., 25., 50., 75., 100.)
84 boun_6_3 = HOMARD.CreateBoundaryCylinder('cyl_2_ext', 17.5, -2.5, -12.5, -100., -75., -25., 50.)
86 boun_6_6 = HOMARD.CreateBoundaryCylinder('cyl_1_int', 0.0, 25., -25., 25., 50., 75., 75.)
88 boun_6_5 = HOMARD.CreateBoundaryCylinder('cyl_2_int', 17.5, -2.5, -12.5, -100., -75., -25., 25.)
94 # Creation of the hypothesis hypo_0_1
95 l_hypothese_0_1 = HOMARD.CreateHypothesis('hypo_6_0_1')
96 l_hypothese_0_1.SetUnifRefinUnRef(1)
97 l_hypothese_0_1.AddGroup('IN1')
98 l_hypothese_0_1.AddGroup('IN2')
99 l_hypothese_0_1.AddGroup('T1_INT_I')
100 l_hypothese_0_1.AddGroup('T1_INT_O')
101 l_hypothese_0_1.AddGroup('T2_INT')
102 # Creation of the hypothesis hypo_1_2
103 l_hypothese_1_2 = HOMARD.CreateHypothesis('hypo_6_1_2')
104 l_hypothese_1_2.SetUnifRefinUnRef(1)
105 l_hypothese_1_2.AddGroup('T1_EXT_I')
106 l_hypothese_1_2.AddGroup('T1_EXT_O')
107 l_hypothese_1_2.AddGroup('T2_EXT')
108 # Creation of the hypothesis hypo_2_3
109 l_hypothese_2_3 = HOMARD.CreateHypothesis('hypo_6_2_3')
110 l_hypothese_2_3.SetUnifRefinUnRef(1)
111 l_hypothese_2_3.AddGroup('INT_I')
112 l_hypothese_2_3.AddGroup('INT_E')
113 l_hypothese_2_3.AddGroup('IN1')
114 l_hypothese_2_3.AddGroup('IN2')
120 le_cas = HOMARD.CreateCase('case_'+nom, nom, ficmed)
121 le_cas.SetDirName(DIRCASE)
122 le_cas.AddBoundary( 'intersection' )
123 le_cas.AddBoundaryGroup( 'cyl_1_int', 'T1_INT_I' )
124 le_cas.AddBoundaryGroup( 'cyl_1_ext', 'T1_EXT_I' )
125 le_cas.AddBoundaryGroup( 'cyl_1_int', 'T1_INT_O' )
126 le_cas.AddBoundaryGroup( 'cyl_1_ext', 'T1_EXT_O' )
127 le_cas.AddBoundaryGroup( 'cyl_2_int', 'T2_INT' )
128 le_cas.AddBoundaryGroup( 'cyl_2_ext', 'T2_EXT' )
137 print(". Itérations")
138 # Iteration iter_6_1 : raffinement selon les faces internes
139 iter_6_1 = le_cas.NextIteration('iter_6_1')
140 iter_6_1.SetMeshName('PIQUAGE_1')
141 iter_6_1.SetMeshFile(os.path.join(DIRCASE, "maill.01.med"))
142 iter_6_1.AssociateHypo('hypo_6_0_1')
143 erreur = iter_6_1.Compute(1, option)
144 print ("erreur = %d" % erreur)
147 # Iteration iter_6_2 : raffinement selon les faces externes
148 iter_6_2 = iter_6_1.NextIteration('iter_6_2')
149 iter_6_2.SetMeshName('PIQUAGE_2')
150 iter_6_2.SetMeshFile(os.path.join(DIRCASE, "maill.02.med"))
151 iter_6_2.AssociateHypo('hypo_6_1_2')
152 erreur = iter_6_2.Compute(1, option)
155 # Iteration iter_6_3 : second raffinement selon les faces externes
156 iter_6_3 = iter_6_2.NextIteration('iter_6_3')
157 iter_6_3.SetMeshName('PIQUAGE_3')
158 iter_6_3.SetMeshFile(os.path.join(DIRCASE, "maill.03.med"))
159 iter_6_3.AssociateHypo('hypo_6_2_3')
160 erreur = iter_6_3.Compute(1, option)
167 message += "Erreur au calcul de l'itération %d" % erreur
169 return erreur, message
171 #========================== Fin de la fonction ==================================
177 # A. Exec of HOMARD-SALOME
179 HOMARD = salome.lcc.FindOrLoadComponent('FactoryServer', 'HOMARD')
180 assert HOMARD is not None, "Impossible to load homard engine"
181 HOMARD.SetLanguageShort("fr")
183 FICMED = os.path.join(DATA_TUTORIAL, "tutorial_4.00.med")
184 FICFRMED = os.path.join(DATA_TUTORIAL, TEST_NAME+".fr.med")
186 ERREUR, MESSAGE = homard_exec("PIQUAGE", FICMED, "COURBES", FICFRMED, DEBUG)
187 except RuntimeError as eee:
189 MESSAGE = str(eee.message)
192 MESSAGE += "Pb in homard_exec"
195 # B. Test of the results
197 N_REP_TEST_FILE = N_ITER_TEST_FILE
198 DESTROY_DIR = not DEBUG
199 test_results(REP_DATA, TEST_NAME, DIRCASE, N_ITER_TEST_FILE, N_REP_TEST_FILE, DESTROY_DIR)
204 raise Exception(MESSAGE)
206 # ==================================
207 gzip_gunzip(DATA_TUTORIAL, 4, 1)
208 gzip_gunzip(DATA_TUTORIAL, 6, 1)
209 # ==================================
211 if salome.sg.hasDesktop():
212 salome.sg.updateObjBrowser(True)
213 iparameters.getSession().restoreVisualState(1)