resources/ADAOSchemaCatalog.xml

   1 <?xml version='1.0' encoding='iso-8859-1' ?>
   2 <!--
   3   Copyright (C) 2010-2011 EDF R&D
   4
   5   This library is free software; you can redistribute it and/or
   6   modify it under the terms of the GNU Lesser General Public
   7   License as published by the Free Software Foundation; either
   8   version 2.1 of the License.
   9
  10   This library is distributed in the hope that it will be useful,
  11   but WITHOUT ANY WARRANTY; without even the implied warranty of
  12   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13   Lesser General Public License for more details.
  14
  15   You should have received a copy of the GNU Lesser General Public
  16   License along with this library; if not, write to the Free Software
  17   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  18
  19   See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  20
  21   Author: Andre Ribes, andre.ribes@edf.fr, EDF R&D
  22 -->
  23 <proc>
  24
  25   <objref name="computeAD" id="python:computeAD:1.0">
  26     <base>pyobj</base>
  27   </objref>
  28
  29   <!-- Types for parametric computations -->
  30   <!-- TODO On devrait pouvoir le lire depuis le KERNEL !!!-->
  31   <type name="long" kind="int"/>
  32   <struct name="SALOME_TYPES/Parameter">
  33     <member type="string" name="name"></member>
  34     <member type="string" name="value"></member>
  35   </struct>
  36   <sequence content="SALOME_TYPES/Parameter" name="SALOME_TYPES/ParameterList"></sequence>
  37   <sequence content="double" name="SALOME_TYPES/Variable"></sequence>
  38   <sequence content="SALOME_TYPES/Variable" name="SALOME_TYPES/VariableSequence"></sequence>
  39   <sequence content="SALOME_TYPES/VariableSequence" name="SALOME_TYPES/StateSequence"></sequence>
  40   <sequence content="SALOME_TYPES/StateSequence" name="SALOME_TYPES/TimeSequence"></sequence>
  41   <sequence content="string" name="SALOME_TYPES/VarList"></sequence>
  42   <struct name="SALOME_TYPES/ParametricInput">
  43     <member type="SALOME_TYPES/VarList" name="inputVarList"></member>
  44     <member type="SALOME_TYPES/VarList" name="outputVarList"></member>
  45     <member type="SALOME_TYPES/TimeSequence" name="inputValues"></member>
  46     <member type="SALOME_TYPES/ParameterList" name="specificParameters"></member>
  47   </struct>
  48   <struct name="SALOME_TYPES/ParametricOutput">
  49     <member type="SALOME_TYPES/TimeSequence" name="outputValues"></member>
  50     <member type="SALOME_TYPES/ParameterList" name="specificOutputInfos"></member>
  51     <member type="long" name="returnCode"></member>
  52     <member type="string" name="errorMessage"></member>
  53   </struct>
  54
  55
  56   <inline name="CreateAssimilationStudy">
  57     <script><code>
  58
  59 <![CDATA[
  60 import numpy, logging
  61 logging.debug("CREATE Entering in CreateAssimilationStudy")
  62 print "Entering in the assimilation study"
  63 print "Name is set to........:", Name
  64 print "Algorithm is set to...:", Algorithm
  65 print "Debug is set to.......:", Debug
  66
  67 # Create Assimilation study
  68 from daYacsIntegration.daStudy import *
  69 assim_study = daStudy(Name, Algorithm, Debug)
  70
  71 # Algorithm parameters
  72 try:
  73   AlgorithmParameters
  74 except NameError:
  75   pass
  76 else:
  77   assim_study.setAlgorithmParameters(AlgorithmParameters)
  78
  79 # Data
  80 # print "Data entered are:"
  81 # Background
  82 try:
  83   Background
  84 except NameError:
  85   pass
  86 else:
  87   logging.debug("CREATE Background is %s"%Background)
  88   logging.debug("CREATE BackgroundType is %s"%BackgroundType)
  89   assim_study.setBackgroundType(BackgroundType)
  90   assim_study.setBackground(Background)
  91
  92 # BackgroundError
  93 try:
  94   BackgroundError
  95 except NameError:
  96   pass
  97 else:
  98   logging.debug("CREATE BackgroundError is %s"%BackgroundError)
  99   logging.debug("CREATE BackgroundErrorType is %s"%BackgroundErrorType)
 100   assim_study.setBackgroundError(BackgroundError)
 101
 102 # Observation
 103 try:
 104   Observation
 105 except NameError:
 106   pass
 107 else:
 108   logging.debug("CREATE Observation is %s"%Observation)
 109   logging.debug("CREATE ObservationType is %s"%ObservationType)
 110   assim_study.setObservationType(ObservationType)
 111   assim_study.setObservation(Observation)
 112
 113 # ObservationError
 114 try:
 115   ObservationError
 116 except NameError:
 117   pass
 118 else:
 119   logging.debug("CREATE ObservationError is %s"%ObservationError)
 120   logging.debug("CREATE ObservationErrorType is %s"%ObservationErrorType)
 121   assim_study.setObservationError(ObservationError)
 122
 123 # ObservationOperator
 124 ObservationOperatorOk = 0
 125 try:
 126   ObservationOperator
 127 except NameError:
 128   pass
 129 else:
 130   logging.debug("CREATE ObservationOperator is %s"%ObservationOperator)
 131   logging.debug("CREATE ObservationOperatorType is %s"%ObservationOperatorType)
 132   assim_study.setObservationOperatorType("Matrix", ObservationOperatorType)
 133   assim_study.setObservationOperator("Matrix", ObservationOperator)
 134   ObservationOperatorOk = 1
 135
 136 if ObservationOperatorOk == 0:
 137   try:
 138     ObservationOperatorDirect
 139   except NameError:
 140     pass
 141   else:
 142     logging.debug("CREATE ObservationOperatorDirect is %s"%ObservationOperatorDirect)
 143     assim_study.setObservationOperatorType("Direct", "Function")
 144     assim_study.setObservationOperator("Direct", ObservationOperatorDirect)
 145   try:
 146     ObservationOperatorTangent
 147   except NameError:
 148     pass
 149   else:
 150     logging.debug("CREATE ObservationOperatorTangent is %s"%ObservationOperatorTangent)
 151     assim_study.setObservationOperatorType("Tangent", "Function")
 152     assim_study.setObservationOperator("Tangent", ObservationOperatorTangent)
 153   try:
 154     ObservationOperatorAdjoint
 155   except NameError:
 156     pass
 157   else:
 158     logging.debug("CREATE ObservationOperatorAdjoint is %s"%ObservationOperatorAdjoint)
 159     assim_study.setObservationOperatorType("Adjoint", "Function")
 160     assim_study.setObservationOperator("Adjoint", ObservationOperatorAdjoint)
 161
 162 # Variables
 163 for name, size in zip(InputVariablesNames, InputVariablesSizes):
 164   assim_study.setInputVariable(name, size)
 165 for name, size in zip(OutputVariablesNames, OutputVariablesSizes):
 166   assim_study.setOutputVariable(name, size)
 167
 168 if has_observers:
 169   logging.debug("CREATE Observers is %s"%observers.keys())
 170   # Adding observers to the study
 171   for observer_name in observers.keys():
 172     scheduler = ""
 173     info = ""
 174     number = str(observers[observer_name]["number"])
 175     if "scheduler" in observers[observer_name].keys():
 176       scheduler = observers[observer_name]["scheduler"]
 177     if "info" in observers[observer_name].keys():
 178       info = observers[observer_name]["info"]
 179     assim_study.addObserver(observer_name, scheduler, info, number)
 180 Study = assim_study
 181 ]]>
 182
 183 </code></script>
 184     <inport name="Name" type="string"/>
 185     <inport name="Algorithm" type="string"/>
 186     <inport name="Debug" type="bool"/>
 187     <inport name="InputVariablesNames" type="stringvec"/>
 188     <inport name="InputVariablesSizes" type="intvec"/>
 189     <inport name="OutputVariablesNames" type="stringvec"/>
 190     <inport name="OutputVariablesSizes" type="intvec"/>
 191     <inport name="has_observers" type="bool"/>
 192     <inport name="observers" type="pyobj"/>
 193     <outport name="Study" type="pyobj"/>
 194   </inline>
 195
 196   <inline name="CreateNumpyMatrixFromString">
 197     <script><code><![CDATA[
 198 import numpy, logging
 199 logging.debug("CREATE Entering in CreateNumpyMatrixFromString")
 200 matrix = numpy.matrix(matrix_in_string)
 201 type = "Matrix"
 202 logging.debug("CREATE Matrix is %s"%matrix)
 203 ]]></code></script>
 204     <inport name="matrix_in_string" type="string"/>
 205     <outport name="matrix" type="pyobj"/>
 206     <outport name="type" type="string"/>
 207   </inline>
 208
 209   <inline name="CreateNumpyMatrixFromScript">
 210     <script><code><![CDATA[
 211 import logging
 212 logging.debug("CREATE Entering in CreateNumpyMatrixFromScript")
 213 type = "Matrix"
 214
 215 # Get file path and filename
 216 import sys
 217 import os
 218 filepath = os.path.dirname(script)
 219 filename = os.path.basename(script)
 220 module_name = os.path.splitext(filename)[0]
 221 sys.path.insert(0,filepath)
 222
 223 # Import script
 224 __import__(module_name)
 225 user_script_module = sys.modules[module_name]
 226
 227 # Get Data from script
 228 ]]></code></script>
 229     <inport name="script" type="string"/>
 230     <outport name="type" type="string"/>
 231   </inline>
 232
 233   <inline name="CreateNumpyVectorFromString">
 234     <script><code><![CDATA[
 235 import numpy, logging
 236 logging.debug("CREATE Entering in CreateNumpyVectorFromString")
 237 vector = numpy.matrix(vector_in_string)
 238 type = "Vector"
 239 logging.debug("Vector is %s"%vector)
 240 ]]></code></script>
 241     <inport name="vector_in_string" type="string"/>
 242     <outport name="vector" type="pyobj"/>
 243     <outport name="type" type="string"/>
 244   </inline>
 245
 246   <inline name="CreateNumpyVectorFromScript">
 247     <script><code><![CDATA[
 248 import logging
 249 logging.debug("CREATE Entering in CreateNumpyVectorFromScript")
 250 type = "Vector"
 251
 252 # Get file path and filename
 253 import sys
 254 import os
 255 filepath = os.path.dirname(script)
 256 filename = os.path.basename(script)
 257 module_name = os.path.splitext(filename)[0]
 258 sys.path.insert(0,filepath)
 259
 260 # Import script
 261 __import__(module_name)
 262 user_script_module = sys.modules[module_name]
 263
 264 # Get Data from script
 265 ]]></code></script>
 266     <inport name="script" type="string"/>
 267     <outport name="type" type="string"/>
 268   </inline>
 269
 270   <inline name="SimpleExecuteDirectAlgorithm">
 271     <script><code><![CDATA[
 272 import logging
 273 logging.debug("EXECUTE Entering in SimpleExecuteDirectAlgorithm")
 274 from daYacsIntegration.daStudy import *
 275 ADD = Study.getAssimilationStudy()
 276 ADD.analyze()
 277 ]]></code></script>
 278     <inport name="Study" type="pyobj"/>
 279     <outport name="Study" type="pyobj"/>
 280   </inline>
 281
 282   <inline name="SimpleUserAnalysis">
 283     <script><code><![CDATA[
 284 #-*-coding:iso-8859-1-*-
 285 import logging
 286 logging.debug("TERMINATE Entering in SimpleUserAnalysis")
 287 from daYacsIntegration.daStudy import *
 288 ADD = Study.getAssimilationStudy()
 289 # User code is below
 290
 291 ]]></code></script>
 292     <inport name="Study" type="pyobj"/>
 293   </inline>
 294
 295   <inline name="FakeOptimizerLoopNode">
 296     <script><code><![CDATA[
 297 import logging
 298 logging.debug("EXECUTE Entering in FakeOptimizerLoopNode")
 299 result = None
 300 ]]></code></script>
 301     <inport name="computation" type="SALOME_TYPES/ParametricInput"/>
 302     <outport name="result" type="SALOME_TYPES/ParametricOutput"/>
 303   </inline>
 304
 305   <inline name="CreateDictFromScript">
 306     <script><code><![CDATA[
 307 import logging
 308 logging.debug("CREATE Entering in CreateDictFromScript")
 309
 310 # Get file path and filename
 311 import sys
 312 import os
 313 filepath = os.path.dirname(script)
 314 filename = os.path.basename(script)
 315 module_name = os.path.splitext(filename)[0]
 316 sys.path.insert(0,filepath)
 317
 318 # Import script
 319 __import__(module_name)
 320 user_script_module = sys.modules[module_name]
 321
 322 # Get Data from script
 323 ]]></code></script>
 324     <inport name="script" type="string"/>
 325   </inline>
 326
 327   <inline name="UserDataInitFromScript">
 328     <script><code><![CDATA[
 329 import logging
 330 logging.debug("CREATE Entering in UserDataInitFromScript")
 331
 332 # Get file path and filename
 333 import sys
 334 import os
 335 filepath = os.path.dirname(script)
 336 filename = os.path.basename(script)
 337 module_name = os.path.splitext(filename)[0]
 338 sys.path.insert(0,filepath)
 339
 340 # Import script
 341 __import__(module_name)
 342 user_script_module = sys.modules[module_name]
 343
 344 # Get Data from script
 345 ]]></code></script>
 346     <inport name="script" type="string"/>
 347     <outport name="init_data" type="pyobj"/>
 348   </inline>
 349
 350   <inline name="ReadForSwitchNode">
 351     <script><code><![CDATA[
 352 import logging
 353 logging.debug("CREATE Entering in ReadForSwitchNode")
 354 switch_value = -1
 355 for param in data["specificParameters"]:
 356   if param["name"] == "switch_value":
 357     switch_value = int(param["value"])
 358 ]]></code></script>
 359     <inport name="data" type="SALOME_TYPES/ParametricInput"/>
 360     <outport name="data" type="SALOME_TYPES/ParametricInput"/>
 361     <outport name="switch_value" type="int"/>
 362   </inline>
 363
 364   <inline name="ExtractDataNode">
 365     <script><code><![CDATA[
 366 import logging
 367 logging.debug("TERMINATE Entering in ExtractDataNode")
 368 import pickle
 369 from daCore.AssimilationStudy import AssimilationStudy
 370 var = None
 371 info = None
 372 for param in data["specificParameters"]:
 373   if param["name"] == "var":
 374     var = pickle.loads(param["value"])
 375   if param["name"] == "info":
 376     info = param["value"]
 377 ]]></code></script>
 378     <inport name="data" type="SALOME_TYPES/ParametricInput"/>
 379     <outport name="var" type="pyobj"/>
 380     <outport name="info" type="pyobj"/>
 381   </inline>
 382
 383   <inline name="ObservationNodeString">
 384     <script><code><![CDATA[
 385 #print "Entering in Observation"
 386
 387 ]]></code></script>
 388     <inport name="var" type="pyobj"/>
 389     <inport name="info" type="pyobj"/>
 390   </inline>
 391
 392   <inline name="ObservationNodeFile">
 393     <script><code><![CDATA[
 394 #print "Entering in Observation"
 395 execfile(script)
 396
 397 ]]></code></script>
 398     <inport name="var"    type="pyobj"/>
 399     <inport name="info"   type="pyobj"/>
 400     <inport name="script" type="string"/>
 401   </inline>
 402
 403   <inline name="EndObservationNode">
 404     <script><code><![CDATA[
 405 # Create a fake output object.
 406 # An observer is always successful.
 407 output = {}
 408 output["outputValues"]        = [[[[]]]]
 409 output["specificOutputInfos"] = []
 410 output["returnCode"]          = 0
 411 output["errorMessage"]        = ""
 412 ]]></code></script>
 413     <outport name="output" type="SALOME_TYPES/ParametricOutput"/>
 414   </inline>
 415
 416   <inline name="SetObserversNode">
 417     <script><code><![CDATA[
 418 #print "Setting observers"
 419 ]]></code></script>
 420     <outport name="has_observers" type="bool"/>
 421     <outport name="observers" type="pyobj"/>
 422   </inline>
 423 </proc>