1 # -*- coding: utf-8 -*-
2 # Copyright (C) 2011-2015 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 test_13 associe au tutorial 3
26 #========================================================================
30 #========================================================================
37 # ==================================
38 pathHomard = os.getenv('HOMARD_ROOT_DIR')
39 # Repertoire des donnees du test
40 Rep_Test = os.path.join(pathHomard, "share", "salome", "resources", "homard")
41 Rep_Test = os.path.normpath(Rep_Test)
42 sys.path.append(Rep_Test)
43 from test_util import remove_dir
44 from test_util import test_results
45 # Repertoire des resultats
47 dircase = os.path.join("/tmp", Test_Name)
48 if ( os.path.isdir(dircase) ) :
52 dircase = tempfile.mkdtemp()
53 # Repertoire des donnees du tutorial
54 data_dir = os.path.join(pathHomard, "share", "doc", "salome", "gui", "HOMARD", "fr", "_downloads")
55 data_dir = os.path.normpath(data_dir)
56 sys.path.append(data_dir)
57 from tutorial_util import gzip_gunzip
58 # ==================================
59 gzip_gunzip(data_dir, 3, -1)
60 # ==================================
64 ipar = iparameters.IParameters(salome.myStudy.GetCommonParameters("Interface Applicative", 1))
65 ipar.append("AP_MODULES_LIST", "Homard")
67 #========================================================================
68 #========================================================================
69 def homard_exec(theStudy):
71 Python script for HOMARD
74 homard.SetCurrentStudy(theStudy)
76 # Hypothese "Hypo_0vers1"
77 # =======================
78 Hypo_0vers1 = homard.CreateHypothesis('Hypo_0vers1')
79 # Characterization of the field
80 Hypo_0vers1.SetField('SOLU_0__QIRE_ELEM_SIGM__________')
81 Hypo_0vers1.SetUseComp(0)
82 Hypo_0vers1.AddComp('ERREST ')
83 Hypo_0vers1.SetRefinThr(3, 1.0)
84 Hypo_0vers1.SetTypeFieldInterp(2)
85 Hypo_0vers1.AddFieldInterp('SOLU_0__DEPL____________________')
86 Hypo_0vers1.AddFieldInterp('SOLU_0__ERRE_ELEM_SIGM__________')
88 # Hypothese "Hypo_1vers2"
89 # =======================
90 Hypo_1vers2 = homard.CreateHypothesis('Hypo_1vers2')
91 # Characterization of the field
92 Hypo_1vers2.SetField('SOLU_1__QIRE_ELEM_SIGM__________')
93 Hypo_1vers2.SetUseComp(0)
94 Hypo_1vers2.AddComp('ERREST ')
95 Hypo_1vers2.SetRefinThr(3, 1.5)
96 Hypo_1vers2.SetUnRefThr(3, 6.)
97 Hypo_1vers2.SetTypeFieldInterp(2)
98 Hypo_1vers2.AddFieldInterp('SOLU_1__DEPL____________________')
99 Hypo_1vers2.AddFieldInterp('SOLU_1__QIRE_ELEM_SIGM__________')
101 # Hypothese "Hypo_1vers2_bis"
102 # ===========================
103 Hypo_1vers2_bis = homard.CreateHypothesis('Hypo_1vers2_bis')
104 # Characterization of the field
105 Hypo_1vers2_bis.SetField('SOLU_1__DEPL____________________')
106 Hypo_1vers2_bis.SetUseComp(1)
107 Hypo_1vers2_bis.AddComp('DX')
108 Hypo_1vers2_bis.AddComp('DY')
109 Hypo_1vers2_bis.AddComp('DZ')
110 Hypo_1vers2_bis.SetRefinThr(1, 0.0001)
111 Hypo_1vers2_bis.SetUnRefThr(1, 0.000001)
112 Hypo_1vers2_bis.SetTypeFieldInterp(0)
116 Case_3 = homard.CreateCase('Case_3', 'G_0', data_dir+'/tutorial_3.00.med')
117 Case_3.SetDirName(dircase)
119 # Iteration "Iter_3_1"
120 # ====================
121 Iter_3_1 = Case_3.NextIteration('Iter_3_1')
122 Iter_3_1.SetMeshName('H_1')
123 Iter_3_1.SetMeshFile(dircase+'/maill.01.med')
124 Iter_3_1.SetFieldFile(data_dir+'/tutorial_3.00.med')
125 Iter_3_1.SetTimeStepRank( 1, 1)
126 Iter_3_1.AssociateHypo('Hypo_0vers1')
127 error = Iter_3_1.Compute(1, 2)
129 # Iteration "Iter_3_2"
130 # ====================
131 Iter_3_2 = Iter_3_1.NextIteration('Iter_3_2')
132 Iter_3_2.SetMeshName('H_2')
133 Iter_3_2.SetMeshFile(dircase+'/maill.02.med')
134 Iter_3_2.SetFieldFile(data_dir+'/tutorial_3.01.med')
135 Iter_3_2.SetTimeStepRank(1, 1)
136 Iter_3_2.AssociateHypo('Hypo_1vers2')
137 error = Iter_3_2.Compute(1, 2)
139 # Iteration "Iter_3_2_bis"
140 # ========================
141 Iter_3_2_bis = Iter_3_1.NextIteration('Iter_3_2_bis')
142 Iter_3_2_bis.SetMeshName('H_2_bis')
143 Iter_3_2_bis.SetMeshFile(dircase+'/maill.02.bis.med')
144 Iter_3_2_bis.SetFieldFile(data_dir+'/tutorial_3.01.med')
145 Iter_3_2_bis.SetTimeStepRank(1, 1)
146 Iter_3_2_bis.AssociateHypo('Hypo_1vers2_bis')
147 error = Iter_3_2_bis.Compute(1, 2)
151 #========================================================================
153 homard = salome.lcc.FindOrLoadComponent('FactoryServer', 'HOMARD')
154 assert homard is not None, "Impossible to load homard engine"
155 homard.SetLanguageShort("fr")
157 # Exec of HOMARD-SALOME
160 error_main = homard_exec(salome.myStudy)
162 raise Exception('Pb in homard_exec at iteration %d' %error_main )
164 raise Exception('Pb in homard_exec: '+e.message)
166 # Test of the results
169 destroy_dir = not debug
170 test_results(Rep_Test, Test_Name, dircase, n_iter_test_file, n_rep_test_file, destroy_dir)
172 # ==================================
173 gzip_gunzip(data_dir, 3, 1)
174 # ==================================
176 if salome.sg.hasDesktop():
177 salome.sg.updateObjBrowser(1)
178 iparameters.getSession().restoreVisualState(1)