1 # -*- coding: iso-8859-1 -*-
2 # Copyright (C) 2011-2014 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 Copyright EDF-R&D 2010, 2014
23 Test test_11 associe au tutorial 3
27 #========================================================================
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 test_results
44 # Repertoire des resultats
45 dircase = tempfile.mktemp()
47 # Repertoire des donnees du tutorial
48 data_dir = os.path.join(pathHomard, "share", "doc", "salome", "gui", "HOMARD", "fr", "_downloads")
49 data_dir = os.path.normpath(data_dir)
50 sys.path.append(data_dir)
51 from tutorial_util import gzip_gunzip
52 # ==================================
53 gzip_gunzip(data_dir, 3, -1)
54 # ==================================
58 ipar = iparameters.IParameters(salome.myStudy.GetCommonParameters("Interface Applicative", 1))
59 ipar.append("AP_MODULES_LIST", "Homard")
61 #========================================================================
62 #========================================================================
63 def homard_exec(theStudy):
65 Python script for HOMARD
68 homard.SetCurrentStudy(theStudy)
70 # Hypothese "Hypo_0vers1"
71 # =======================
72 Hypo_0vers1 = homard.CreateHypothesis('Hypo_0vers1')
73 # Characterization of the field
74 Hypo_0vers1.SetField('SOLU_0__QIRE_ELEM_SIGM__________')
75 Hypo_0vers1.SetUseComp(0)
76 Hypo_0vers1.AddComp('ERREST ')
77 Hypo_0vers1.SetRefinThr(3, 1.0)
78 Hypo_0vers1.SetTypeFieldInterp(2)
79 Hypo_0vers1.AddFieldInterp('SOLU_0__DEPL____________________')
80 Hypo_0vers1.AddFieldInterp('SOLU_0__ERRE_ELEM_SIGM__________')
82 # Hypothese "Hypo_1vers2"
83 # =======================
84 Hypo_1vers2 = homard.CreateHypothesis('Hypo_1vers2')
85 # Characterization of the field
86 Hypo_1vers2.SetField('SOLU_1__QIRE_ELEM_SIGM__________')
87 Hypo_1vers2.SetUseComp(0)
88 Hypo_1vers2.AddComp('ERREST ')
89 Hypo_1vers2.SetRefinThr(3, 1.5)
90 Hypo_1vers2.SetUnRefThr(3, 6.)
91 Hypo_1vers2.SetTypeFieldInterp(2)
92 Hypo_1vers2.AddFieldInterp('SOLU_1__DEPL____________________')
93 Hypo_1vers2.AddFieldInterp('SOLU_1__QIRE_ELEM_SIGM__________')
95 # Hypothese "Hypo_1vers2_bis"
96 # ===========================
97 Hypo_1vers2_bis = homard.CreateHypothesis('Hypo_1vers2_bis')
98 # Characterization of the field
99 Hypo_1vers2_bis.SetField('SOLU_1__DEPL____________________')
100 Hypo_1vers2_bis.SetUseComp(1)
101 Hypo_1vers2_bis.AddComp('DX')
102 Hypo_1vers2_bis.AddComp('DY')
103 Hypo_1vers2_bis.AddComp('DZ')
104 Hypo_1vers2_bis.SetRefinThr(1, 0.0001)
105 Hypo_1vers2_bis.SetUnRefThr(1, 0.000001)
106 Hypo_1vers2_bis.SetTypeFieldInterp(0)
110 Case_3 = homard.CreateCase('Case_3', 'G_0', data_dir+'/tutorial_3.00.med')
111 Case_3.SetDirName(dircase)
113 # Iteration "Iter_3_1"
114 # ====================
115 Iter_3_1 = Case_3.NextIteration('Iter_3_1')
116 Iter_3_1.SetMeshName('H_1')
117 Iter_3_1.SetMeshFile(dircase+'/maill.01.med')
118 Iter_3_1.SetFieldFile(data_dir+'/tutorial_3.00.med')
119 Iter_3_1.SetTimeStepRank( 1, 1)
120 Iter_3_1.AssociateHypo('Hypo_0vers1')
121 error = Iter_3_1.Compute(1, 2)
123 # Iteration "Iter_3_2"
124 # ====================
125 Iter_3_2 = Iter_3_1.NextIteration('Iter_3_2')
126 Iter_3_2.SetMeshName('H_2')
127 Iter_3_2.SetMeshFile(dircase+'/maill.02.med')
128 Iter_3_2.SetFieldFile(data_dir+'/tutorial_3.01.med')
129 Iter_3_2.SetTimeStepRank(1, 1)
130 Iter_3_2.AssociateHypo('Hypo_1vers2')
131 error = Iter_3_2.Compute(1, 2)
133 # Iteration "Iter_3_2_bis"
134 # ========================
135 Iter_3_2_bis = Iter_3_1.NextIteration('Iter_3_2_bis')
136 Iter_3_2_bis.SetMeshName('H_2_bis')
137 Iter_3_2_bis.SetMeshFile(dircase+'/maill.02.bis.med')
138 Iter_3_2_bis.SetFieldFile(data_dir+'/tutorial_3.01.med')
139 Iter_3_2_bis.SetTimeStepRank(1, 1)
140 Iter_3_2_bis.AssociateHypo('Hypo_1vers2_bis')
141 error = Iter_3_2_bis.Compute(1, 2)
145 #========================================================================
147 homard = salome.lcc.FindOrLoadComponent('FactoryServer', 'HOMARD')
148 assert homard is not None, "Impossible to load homard engine"
149 homard.SetLanguageShort("fr")
151 # Exec of HOMARD-SALOME
154 error_main = homard_exec(salome.myStudy)
156 raise Exception('Pb in homard_exec at iteration %d' %error_main )
158 raise Exception('Pb in homard_exec: '+e.message)
160 # Test of the results
163 test_results(Rep_Test, Test_Name, dircase, n_iter_test_file, n_rep_test_file)
165 # ==================================
166 gzip_gunzip(data_dir, 3, 1)
167 # ==================================
169 if salome.sg.hasDesktop():
170 salome.sg.updateObjBrowser(1)
171 iparameters.getSession().restoreVisualState(1)