tests/test_13.py

   1 # -*- coding: iso-8859-1 -*-
   2 # Copyright (C) 2011-2014  CEA/DEN, EDF R&D
   3 #
   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.
   8 #
   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.
  13 #
  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
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20 """
  21 Python script for HOMARD
  22 Copyright EDF-R&D 2010, 2014
  23 Test test_11 associe au tutorial 3
  24 """
  25 __revision__ = "V2.1"
  26
  27 #========================================================================
  28 Test_Name = "test_13"
  29 n_iter_test_file = 2
  30 #========================================================================
  31 import os
  32 import tempfile
  33 import sys
  34 import HOMARD
  35 import salome
  36 #
  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()
  46 os.mkdir(dircase)
  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 # ==================================
  55
  56 salome.salome_init()
  57 import iparameters
  58 ipar = iparameters.IParameters(salome.myStudy.GetCommonParameters("Interface Applicative", 1))
  59 ipar.append("AP_MODULES_LIST", "Homard")
  60 #
  61 #========================================================================
  62 #========================================================================
  63 def homard_exec(theStudy):
  64   """
  65 Python script for HOMARD
  66   """
  67   #
  68   homard.SetCurrentStudy(theStudy)
  69   #
  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__________')
  81   #
  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__________')
  94   #
  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)
 107   #
 108   # Cas
 109   # ===
 110   Case_3 = homard.CreateCase('Case_3', 'G_0', data_dir+'/tutorial_3.00.med')
 111   Case_3.SetDirName(dircase)
 112   #
 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)
 122   #
 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)
 132   #
 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)
 142   #
 143   return error
 144
 145 #========================================================================
 146
 147 homard = salome.lcc.FindOrLoadComponent('FactoryServer', 'HOMARD')
 148 assert homard is not None, "Impossible to load homard engine"
 149 homard.SetLanguageShort("fr")
 150 #
 151 # Exec of HOMARD-SALOME
 152 #
 153 try :
 154   error_main = homard_exec(salome.myStudy)
 155   if error_main :
 156     raise Exception('Pb in homard_exec at iteration %d' %error_main )
 157 except Exception, e:
 158   raise Exception('Pb in homard_exec: '+e.message)
 159 #
 160 # Test of the results
 161 #
 162 n_rep_test_file = 3
 163 test_results(Rep_Test, Test_Name, dircase, n_iter_test_file, n_rep_test_file)
 164 #
 165 # ==================================
 166 gzip_gunzip(data_dir, 3, 1)
 167 # ==================================
 168 #
 169 if salome.sg.hasDesktop():
 170   salome.sg.updateObjBrowser(1)
 171   iparameters.getSession().restoreVisualState(1)
 172